24小时热门版块排行榜

返回列表

zgqghh

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 1041.6
帖子: 136
在线: 46.1小时
虫号: 2925926
注册: 2014-01-11
性别: GG
专业: 多相流热物理学

[求助] MATLAB中采用同fsolve求解复杂非线性方程组的问题已有1人参与

科研需要求一组比较复杂的非线性方程组，试着编写了一些程序，fun函数代码为：
function F=myfun(q)
A=2;
delta=0.75;
beta=0.35;
gama=0.5;
r=q(1);
p=q(2);
w=q(3);
F(1)=w/(1+delta*log(w))+r*(1+gama/beta)/(2*(r+2))-1;
F(2)=w/(1+delta*log(w))+((1+delta*log(A/p))*p)/A-1;
F(3)=w-2*beta*(1+r)*p/r^2;

初始值q0=[1 1.5 2]，求解r,p,w。运行程序后发现提示Solver stopped prematurely. fsolve stopped because it exceeded the function evaluation limit, options.MaxFunEvals = 300 (the default value).由于是第一次求解这么复杂的非线性方程，不知道该怎么处理了。出现这个是因为初始值不合适？或者原方程组的非线性太大？不知道有没有可以优化处理的办法，各位大神，求指导，先行拜谢！

需求解的非线性方程组.jpg

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

求用MatAB求解已经有9人回复
求助matlab高手指点：大型非线性方程组（牛顿拉佛森法）已经有7人回复
MATLAB解非线性方程组已经有9人回复
MATLAB求解五元非线性方程组，在线等已经有7人回复
求助：非线性方程组的求解（郁闷中）已经有26人回复
求助解复杂非线性方程组的好的方法已经有24人回复
matlab fsolve得出的解全部为复数解已经有15人回复
用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
matlab如何求解一个非线性微分方程组已经有8人回复
matlab求解高阶方程组已经有3人回复
MATLAB求解非线性方程组已经有5人回复
求Matlab解非线性方程代码已经有12人回复
【求助】非线性方程组的求解问题已经有6人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
使用fsolve求解非线性方程问题已经有9人回复
matlab求解方程中的参数已经有21人回复
求助matlab---fsolve解非线性方程组已经有6人回复
用matlab求解方程出问题，请帮忙看看已经有3人回复
matlab求解非线性方程组已经有16人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
请教Matlab高手解多元非线性方程组中的参数已经有6人回复

唉。。。

1楼 2015-05-17 22:26:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

1stOpt计算：

CODE:

Constant A=2, delta=0.75, beta1=0.35, gama=0.5;

Function

w/(1+delta*ln(w))+r*(1+gama/beta1)/(2*(r+2))-1;

w/(1+delta*ln(w))+((1+delta*ln(A/p))*p)/A-1;

w-2*beta1*(1+r)*p/r^2;

w: 3.18918393381689
r: -0.735048145928154
p: 9.29065081088276

赞一下