# 5 麦片

050

1主题 0最佳答案

 function [t,x,y,z] = shougao3252(ni_a_ch4,ni_a_h2o,p,uf) for t=600:0.01:1500 syms x y z eq1=((((((3*x+y-z)/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y)))^3)*(((x-y)/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y)))*(p^2))/((((ni_a_h2o-x-y+z)/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y)))*(((ni_a_ch4-x)/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y)))))-10^((-6.61395*10)+(1.95028*0.1*t)+(-2.25232*0.0001*t^2)+(1.24065*0.0000001*t^3)+(-2.63121*0.00000000001*t^4)); eq2=((((y/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y)))*(((3*x+y-z)/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y))))/((((x-y)/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y)))*(((ni_a_h2o-x-y+z)/((x-y)+(3*x+y-z)+(ni_a_ch4-x)+(ni_a_h2o-x-y+z)+y)))))-10^((1.32097*10)+(-3.915*0.01*t)+(4.63742*0.00001*t^2)+(-2.57479*0.00000001*t^3)+(5.47301*0.000000000001*t^4)); eq3=uf-z/(3*x+y); [x,y,z]=solve(eq1,eq2,eq3) if [x,y,z]==empty syms    continue else     break end end 不知道怎么设置x y z都大于零，和怎么判断是否收敛

# 5 麦片

050

1主题 0最佳答案

 亲，你问题解决了吗？我也遇到类似问题了:'(
 您需要登录后才可以回帖 登录 | 注册 本版积分规则 回帖后跳转到最后一页