| 查看: 4388 | 回复: 15 | |||
xiangzigood新虫 (初入文坛)
|
[求助]
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 |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教fsolve如何查看每次迭代的x值？
已经有6人回复
- 求助 MATLAB-fsolve
已经有3人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- 使用fsolve求解非线性方程问题
已经有9人回复
- 使用matlab 中的fsolve 老提示出错
已经有4人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- Matlab,解复数方程
已经有8人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 请教下：如何求一个函数在一定范围内的全部解
已经有10人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 【求助】matlab中fsolve解的精度控制【已解决】
已经有4人回复
- 一个方程组的解法
已经有8人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-06-15 17:38:26
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555220
- 散金: 27828
- 红花: 2147
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
3楼2013-06-15 17:53:25
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
4楼2013-06-16 08:25:52
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
5楼2013-06-16 08:26:27
6楼2013-06-16 11:01:05
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
7楼2013-06-17 08:26:26
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
8楼2013-06-17 08:31:02
9楼2013-06-17 10:40:40
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
|
function[x,x0,fval,exitflag]=fsolveyalinjieliu d=0.04;D=0.05;k=1.4;ff=0.004;L=1000;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)*sqrt(2*((Pa/x(5))^(2/k)-(Pa/x(5))^((k+1)/k))/(k-1))-x(2); (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); x(6)-(x(3)*T1)/x(4); ] end end 这个应该比较完整，六个方程六个未知数 |
10楼2013-06-17 11:21:58
