| 查看: 823 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】fsovle问题【已完成】
|
|||
|
改变初值，结果会变的很厉害，这是为什么？请高人指导~ ----------------------------------------------------------- function xFsovle clear all; clc; global x Mn Mw format short g x0=[6.3323e-007 60.515 0.18628 0.007572 0.0014898 0.42654 243.81 2.6725e-020 2.6725e-020 2.6725e-015]; [x,FVAL,EXITFLAG] =fsolve(@NonlinEqs,x0) Mn=42*(x(6)+x(9))/(x(5)+x(8)) Mw=42*(x(7)+x(10))/(x(6)+x(9)) PDI=Mw/Mn %------------------------------------- function f =NonlinEqs(x) f(1)=108.85*x(1)*x(2)-(4.8*x(3).^0.5+0.12*x(4)+0.001)*x(5)+(0-0.038)/(90*3600); f(2)=108.85*x(1)*x(2)+108.85*x(2)*x(5)-(4.8*x(3).^0.5+0.012*x(2)+0.12*x(4)+0.001)*x(6)+0.012*x(2)*x(5)+(0-28.6)/(90*3600); f(3)=108.85*x(1)*x(2)+2*108.85*x(2)*x(6)+108.85*x(2)*x(5)-(4.8*x(3).^0.5+0.012*x(2)+0.12*x(4)+0.001)*x(7)+0.012*x(2)*x(5)+(0-42901)/(90*3600); f(4)=(4.8*x(3).^0.5+0.012*x(2)+0.12*x(4)+0.001)*x(8)+(0-153.32)/(90*3600); f(5)=(4.8*x(3).^0.5+0.012*x(2)+0.12*x(4)+0.001)*x(9)+(0-114967.8)/(90*3600); f(6)=(4.8*x(3).^0.5+0.012*x(2)+0.12*x(4)+0.001)*x(10)+(0-172300000)/(90*3600); f(7)=-108.85*x(1)*x(2)-(108.85*x(2)+0.012*x(2))*x(5)+(817500000-796900000-4334380)/(90*42*3600); f(8)=-4.8*x(5)*x(3).^0.5+(4000000-3998000)/(2*3600*90); f(9)=0.0013*0.001146+(4.8*x(3).^0.5+0.12*x(4))*x(5)-108.85*x(1)*x(2); f(10)=-0.12*x(4)*x(5)+(10000-9950)/(90*3600*114); [ Last edited by nono2009 on 2009-9-24 at 07:33 ] |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有28人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有23人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
5楼2009-06-25 09:43:27
2楼2009-06-23 13:50:56
3楼2009-06-23 14:31:36
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
4楼2009-06-25 00:05:05
