作业帮 > 综合 > 作业

matlab中 while循环问题

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/07/16 15:40:48
matlab中 while循环问题
t2=50;
t21=40;
while abs(t2-t21)>0.1
ta3= 30.2684
syms t21;
w11=w01(ta3);
w22=w01(t21);
mw=ma*(w22-w11);
ha1=1.005*ta1+w11*(2.051+1.86*ta1);
ha2=0.00585*(t21+273.16)^3-0.497*(t21+273.16)^2+19.87*(t21+273.16)-207.61;
L=0;
L=(ma*ha2+(m1-ma*(w22-w11))*t21*cw-ma*ha1-m1*cw*tw1);
A=solve(L,t21);
t21=A(1,end)
t2=(t21+t2)/2
end
为什么不能循环,公式没有错,能运行
其中
function [w] =w01(t)
c0=-6.0436117*10^3;c1=1.89318833*10;c2=-2.8238594*10^(-2);
c3=1.7241129*10^(-5); c4=2.858487;
E=exp(c0*(t+273.16)^(-1)+c1+c2*(t+273.16)+c3*(t+273.16)^(2)+c4*log(t+273.16));
Pa=1.013*10^5;
w=0.662*E/(Pa-E);
end
t2=50;
t21=40;
ma=3;
ta1=2;
m1=1.5;
cw=2;
tw1=2.5;
while abs(t2-t21)>0.1
ta3= 30.2684;
syms t21;
w11=w01(ta3);
w22=w01(t21);
mw=ma*(w22-w11);
ha1=1.005*ta1+w11*(2.051+1.86*ta1);
ha2=0.00585*(t21+273.16)^3-0.497*(t21+273.16)^2+19.87*(t21+273.16)-207.61;
%L=0;
L=(ma*ha2+(m1-ma*(w22-w11))*t21*cw-ma*ha1-m1*cw*tw1);
A=solve(L,t21);
t21=A(1,end)
t2=(t21+t2)/2
end
主程序按上面试试