empty sym error, how to solve

2020 8 20
1 个回答
回答已采纳
更新时间：2020 8 22
9 次查看（30 天）

0 个投票

syms P1 P2 V1 V2 z1 z2 g p debi Wpump Wturb Eloss Wpumpshaft npump nturb alfa1 alfa2
g=9.807
p=1000
V1=0
z1=0
alfa1=1;
P2=0
V2=1
z2=0
alfa2=1;
debi=70 %kg/sn
Wpumpshaft=20100
npump=0.683074626865671
Wpump=Wpumpshaft*npump
Wturb=0;
nturb=1
Eloss=0;
ar=((P1/p+V1^2/(2)*alfa1+z1*g)*debi+Wpump==(P2/p+V2^2/(2)*alfa2+z2*g)*debi+Wturb*nturb+Eloss);
vpa(solve(ar,P1),8)
when i try to solve this it give error. it is ridiculuous. Where is the error.

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

 采纳的回答

Alan Stevens
Alan Stevens 2020-8-22

0 个投票

You could use fzero:
P10 = -100;
P1 = fzero(@fn, P10);
function F = fn(P1)
g=9.807;
p=1000;
V1=0;
z1=0;
alfa1=1;
P2=0;
V2=1;
z2=0;
alfa2=1;
debi=70; %kg/sn
Wpumpshaft=20100;
npump=0.683074626865671;
Wpump=Wpumpshaft*npump;
Wturb=0;
nturb=1;
Eloss=0;
F = ((P1/p+V1^2/(2)*alfa1+z1*g)*debi+Wpump ...
- (P2/p+V2^2/(2)*alfa2+z2*g)*debi+Wturb*nturb+Eloss);
end

更多回答（0 个）

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by