解一個閥門的數(shù)學模型

shouce

ngsxngtd 發(fā)表于 2015-5-24 21:08
我來簡化一下：
看似是非線性方程組，實際上不過是一元二次方程求解，初中生足以。
" S) D/ o6 f1 C# B. V* F最后一個方程：S6，α ...

>> syms a b c x
>> f=sym('a*x^2+b*x+c')

f =

( W/ h& z, w- n0 ^, }8 t8 w7 s3 La*x^2 + b*x + c
3 u- V4 X& |* |$ D" D4 w* u
# U4 F3 F. H2 @4 d( ]) u" K, t* Y>> finverse(f,x)
+ k, H1 D7 c# R' k' f
ans =

-(b + (b^2 - 4*a*c + 4*a*x)^(1/2))/(2*a)
5 Z) L: |  q6 w
( G; R% o2 _) y- p5 `4 U% G' c我用matlab 推導 一元二次方程求根公式      
$ r! Z& N/ D/ }( ?) |2 _            
               2
. D! K+ ~  I* [7 P  b + sqrt(b  - 4 a c + 4 a x)
- ----------------------------
5 K* Y) F( F9 I7 c               2 a
6 d4 C" r- @) m& A0 J, u$ X
4 b( S) G+ u* H: d
* O7 A; O6 ?( F$ A5 B3 F0 b* u: W, _這個也含有X  估計  可以刪除

" J: {8 R" Y! X! t4 h* w

shouce

ngsxngtd 發(fā)表于 2015-5-24 21:46
糊涂了，既然如此，那就是迭代吧。
不管含有多少個ri,哪怕是離散的映射關(guān)系，反正有邊界條件可以控制。
4 K8 C8 G& x3 ^具 ...

不想要開平方根的S   可以用權(quán)函數(shù)  之類的表達    可以吧

明月山河

ngsxngtd 發(fā)表于 2015-5-24 21:46
( e/ P0 u0 e! V  Y0 G糊涂了，既然如此，那就是迭代吧。
0 D0 k# @* M4 j% a" N0 u" }. L不管含有多少個ri,哪怕是離散的映射關(guān)系，反正有邊界條件可以控制。
具 ...

6 ~; V  P( N, B( [# R. A流量系數(shù)只有40，但是這個數(shù)字沒有什么意義，這個閥是用來清洗物體的換能閥，流體高頻振動的，慣性力起主要作用，而且與頻率有關(guān)， 流動阻力相比之下很小，對阻抗的關(guān)注更大一些。收縮系數(shù)暫定取0.9以下。修正系數(shù)？好像用不到。
另外圖中的通道也不是主要水流流道，而是脈沖波的傳輸通道。六個格子是分隔不同相位的脈沖，防止其能量混合相消。所以格子的位置關(guān)系非常嚴格，才有前面的模型。由于脈沖波的擴散方式還沒有選定，所以Si的函數(shù)也沒有確定。所以提取了一個子問題出來單獨解決。其他的參數(shù)還在設(shè)計中。
4 S' W0 i8 c: W8 Z+ D% y8 y# o' f" U

湖北的旱鴨子

看不懂，閥門好復雜的樣子呢~

zms9439

你給的公式真復雜。
把 R6-δ 定義為R0
S6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ
8 O7 w, {# C: ~4 i# e" ~還是3個未知，再簡化下，a角把δ/2去掉，

明月山河

zms9439 發(fā)表于 2015-7-8 18:18
你給的公式真復雜。
把 R6-δ 定義為R0
S6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ

- c. h4 _* n. V; q+ k0 }謝謝關(guān)注，我已經(jīng)解決了，用級數(shù)展開到4次項。
& Z: s" }0 K/ P& s4 @" }

zms9439

明月山河 發(fā)表于 2016-1-19 14:29
謝謝關(guān)注，我已經(jīng)解決了，用級數(shù)展開到4次項。

5 ~; |& V/ J. j! Q, ^% K5 {) p呵呵，時隔半年，終于解決了。恭喜
計算的結(jié)果可以發(fā)上來看看呀？

dahai102320

好高端的樣子，真心不懂。希望以后能懂

落雪clq

牛人就是多