| 查看: 822 | 回复: 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个名额)地质化工冶金领域
已经有5人回复
欢迎发来filecode的Mz6后的代码验证其规律
已经有6人回复
不应该看fileCode
已经有9人回复
Filecode 又变了,巨变
已经有4人回复
分享一下我之前已中青C的计划书的filecode
已经有5人回复
FileCode能看出啥?
已经有26人回复
关于Filecode分析方法
已经有11人回复
重要来源:本周末出结果
已经有9人回复
我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人
已经有14人回复
结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。
已经有4人回复
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










回复此楼