24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiangzigood

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 321.5
红花: 1
帖子: 18
在线: 11.7小时
虫号: 1891458
注册: 2012-07-13
专业: 制造系统与自动化

[求助] matlab fsolve得出的解全部为复数解

我用matlab中的FSOLVE解一个六元非线性方程组，为什么得出的解全部为复数解，可是理论上应该解为实数，本人Matlab新手，代码如下，大神们帮忙看一下吧！
function[x,x0,fval,exitflag]=fsolveyalinjieliu
d=0.045;D=0.05;k=1.4;ff=0.004;L=100;P1=1500000;T1=300;Pa=101325;M=29;R=8.314;
x0=unifrnd(0,1,1,6);
options=optimset('MaxFunEvals',20000,'MaxIter',5000);
[x,fval,exitflag]=fsolve(@f,x0,options);
function F=f(x)
F=[((k+1)/2)*log((x(2)^2*x(3))/(x(1)^2*x(4)))-(1/(x(1)^2)-1/(x(2)^2))+4*ff*L/D;
(d^2/D^2)*x(5)*sqrt((2*k*M)*((Pa/x(5))^(2/k)-(Pa/x(5))^((k+1)/k))/(R*x(6)*(k-1)))-x(2)*x(5)*sqrt((k*M)/(R*x(6)));
(d^2/D^2)*x(5)*sqrt((2*k*M)*((Pa/x(5))^(2/k)-(Pa/x(5))^((k+1)/k))/(R*x(6)*(k-1)))-x(1)*P1*sqrt((k*M)/(R*T1));
1+((k-1)*x(1)^2)/2-x(3);
1+((k-1)*x(2)^2)/2-x(4);
]
end
end

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

1楼2013-06-15 09:45:26

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

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★
jjdg: 金币+1, 感谢参与 2013-06-18 01:00:00

用1stOpt吧，其它软件得到正解的概率几乎为0：

CODE:

Parameter x(6);

Constant d=0.045,D1=0.05,k=1.4,ff=0.004,L=100,P1=1500000,T1=300,Pa=101325,M=29,R=8.314;

Function ((k+1)/2)*ln((x2^2*x3)/(x1^2*x4))-(1/(x1^2)-1/(x2^2))+4*ff*L/D1;

        (D1^2/D1^2)*x5*sqrt((2*k*M)*((Pa/x5)^(2/k)-(Pa/x5)^((k+1)/k))/(R*x6*(k-1)))-x2*x5*sqrt((k*M)/(R*x6));

        (D1^2/D1^2)*x5*sqrt((2*k*M)*((Pa/x5)^(2/k)-(Pa/x5)^((k+1)/k))/(R*x6*(k-1)))-x1*P1*sqrt((k*M)/(R*T1));

        1+((k-1)*x1^2)/2-x3;

        1+((k-1)*x2^2)/2-x4;

        x6-(x3*T1)/x4;

x1: 0.157422221916776
x2: 0.404163373392218
x3: 1.0049585658072
x4: 1.03267007875412
x5: 576359.723370304
x6: 291.949554479051