24小时热门版块排行榜    

查看: 576  |  回复: 0

qqcrw

银虫 (初入文坛)

[求助] matlab求解一非线性方程组

麻烦帮我看一下代码有没有错误?用matlab解非线性方程组除了sovle外还有什么比较好的函数?谢谢!
方程组可以看附图,由t算出第一个方程组中的KP1.KP2,然后代入第二个方程组中求解xy。m,t,p为已知值,求x,y。麻烦帮我分析下出现了什么错误。谢谢
M文件如下:
function f=HCQ(m,t,p)
lgkp1=-9864.75/t+8.3666*log10(t)-0.0020814.*t+1.8737*10^(-7).*t.^2-11.894;
lgkp2=2.183/t-0.09361*log10(t)+0.000632.*t-1.08*10^(-7).*t.^2-2.298;
kp1=10^lgkp1;kp2=10^lgkp2;
[x,y]=solve('kp1=(p/(1+m+2*x))^2*(x-y)*(3*x+y)^3/(1-x)/(m-x-y)','kp2=y*(3*x+y)/(x-y)/(m-x-y)')
f=[x,y]
运行结果如下:   
HCQ(2,800+273.15,3000000)
??? Error using ==> solve at 162
Unable to find closed form solution.

Error in ==> HCQ at 5
[x,y]=solve('kp1=(p/(1+m+2*x))^2*(x-y)*(3*x+y)^3/(1-x)/(m-x-y)','kp2=y*(3*x+y)/(x-y)/(m-x-y)','x','y')
matlab求解一非线性方程组
捕获.PNG


matlab求解一非线性方程组-1
捕获2.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 qqcrw 的主题更新
信息提示
请填处理意见