vpasolve function incorrect solution

2020 7 17
1 个回答
回答已采纳
更新时间：2020 7 18
1 次查看（30 天）

0 个投票

Hello
I am trying to solve this set of equations, however am getting incorrect solution:
syms Q1 Q2 f1 f2 Re1 Re2
eq1= Re1==(31.77E6)*Q1;
eq2= Re2==(15.88E6)*Q2;
eq3= 1/(f1^(1/2))==-1.8*log10(6.9/(Re1)+((1.125E-3)/3.7)^(1.11));
eq4= 1/(f2^(1/2))==-1.8*log10(6.9/(Re2)+((5.625E-4)/3.7)^(1.11));
eq5= 0.572==8*Q1*(1+sqrt(32*f1/f2))+(29.05E6)*f1*(Q1^2);
eq6= Q2==sqrt(32*f1/f2)*Q1;
S=vpasolve(eq1,eq2,eq3,eq4,eq5,eq6);
Matlab answer for f1 for example is 415671.45, however the correct answer should be approx 0.022. What am i doing wrong?

6 个评论
显示 4更早的评论 隐藏 4更早的评论

madhan ravi
madhan ravi 2020-7-17
编辑：madhan ravi 2020-7-17
You can specify the range into the solver.
Samer Ali
Samer Ali 2020-7-17
How can i specify the range for just one variable, lets say i know the solution of f1 to be between 0 and 1.
Walter Roberson
Walter Roberson 2020-7-17
There seem to be at least 3 valid solutions. The solution that MATLAB produces is one of them. Another has f1 near 0.02653046284, and a third has f1 near 384805.6353. The f2 values are fairly similar (but not exactly the same) for all three of those.
Walter Roberson
Walter Roberson 2020-7-17
S = vpasolve([eq1,eq2,eq3,eq4,eq5,eq6],[f1 f2 Q1 Q2 Re1 Re2], [0 1; [-inf(5,1), inf(5,1)]] );
Samer Ali
Samer Ali 2020-7-17
编辑：Samer Ali 2020-7-17
Hello Walter
executing your command:
S = vpasolve([eq1,eq2,eq3,eq4,eq5,eq6],[f1 f2 Q1 Q2 Re1 Re2], [0 1; [-inf(5,1), inf(5,1)]] );
Matlab gives me: Empty sym: 0-by-1?
The solution associated with f1 = 0.02653 should be the correct one but how i would get Matlab to search for this solution?
Alex Sha
Alex Sha 2020-7-18
编辑：Alex Sha 2020-7-18
1:
q1: 0.000823993387369594
q2: 0.00526536192506034
f1: 0.0265304628607812
f2: 0.0207914989607889
re1: 26178.269916732
re2: 83613.9473699583
feval:
-3.63797880709171E-12
1.16415321826935E-10
-1.77635683940025E-15
7.99360577730113E-15
2.22044604925031E-16
1.11022302462516E-16
2:
q1: 0.000859390439731193
q2: 4.34840412474342E-7
f1: 0.0263399277193314
f2: 3292196.85150343
re1: 27302.83427026
re2: 6.90526575009254
feval:
-3.63797880709171E-12
-1.06581410364015E-14
3.5527136788005E-15
4.41487124636097E-16
0
-1.70167860312205E-17

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

Walter Roberson
Walter Roberson 2020-7-17

0 个投票

S1 = solve([eq1,eq2,eq5,eq6],[Q1,Q2,Re1,Re2]);
E2 = subs([eq3,eq4], S1);
sol1 = vpasolve(E2(1,:), [0 1; 0 1]);
The other branch, by solving E2(2,:), eventually gives up.

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Samer Ali
Samer Ali 2020-7-17
Thank you Walter for the help. I solved it using Iterative Gauss Seidel method, which i guess was more reliable.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by