查看: 86|回复: 1|关注: 0

[已答复] 有关4个未知数,3个非线性方程组成的方程组迭代循环直到收敛

[复制链接]

新手

5 麦片

财富积分


050


1

主题

1

帖子

0

最佳答案
发表于 7 天前 | 显示全部楼层 |阅读模式
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

主题

3

帖子

0

最佳答案
发表于 昨天 19:43 | 显示全部楼层
亲,你问题解决了吗?我也遇到类似问题了:'(
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表