| 查看: 532 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
liaoliao8545新虫 (小有名气)
|
[求助]
matlab
|
|
|
我用ezplot画函数图找零点运行后出现下面问题 Current plot held Warning: Function failed to evaluate on array inputs; vectorizing the function may speed up its evaluation and avoid the need to loop over array elements. 什么原因呢 还有我用fsolve解非线性方程组的解 x30 x40 x50 x60 为初值 运行 出现下面的情况 x30 = 0.0158 x40 = 0.6910 x50 = 1.3824e-035 x60 = 2.8181e-032 Optimization terminated: no further progress can be made. Trust-region radius less than 2*eps. Problem may be ill-conditioned or Jacobian may be inaccurate. Try using exact Jacobian or check Jacobian for errors. z = 0.0158 0.6910 -0.0000 0.0000 什么原因呢 这是非线性方程组m文件 function f=fun_fun(x) global C1 C2 C7 C8 C9 C10 d1 d2 d3 n m;%%%%这些是由主程序求得的 f=[2*C2*x(4)*sqrt(x(1))+2*x(4)+C7*sqrt(x(4))*sqrt(x(4))+C9*sqrt(x(4))-(C1*x(3)*sqrt(x(1))+x(3))*d1; 2*C1*x(3)*sqrt(x(1))+C2*x(4)*sqrt(x(3))+2*x(1)+x(3)+C7*sqrt(x(4))*sqrt(x(1))+C8*sqrt(x(1))+C10*sqrt(x(1))*sqrt(x(2))-(C1*x(3)*sqrt(x(1))+x(3))*d2; 2*x(2)+C10*sqrt(x(1))*sqrt(x(2))-(C1*x(3)*sqrt(x(1))+x(3))*d3; C1*x(3)*sqrt(x(1))+C2*x(4)*sqrt(x(1))+x(3)+x(2)+x(3)+x(4)+C7*sqrt(x(4))*sqrt(x(3))+C8*sqrt(x(3))+C9*sqrt(x(4))+C10*sqrt(x(1))*sqrt(x(2))+1/(n*(C1*x(3)*sqrt(x(1))+x(3)))-1;]; |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
liaoliao8545
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 213.5
- 沙发: 1
- 帖子: 197
- 在线: 45.8小时
- 虫号: 2119920
- 注册: 2012-11-11
- 专业: 燃烧学
4楼2013-05-12 16:12:34
2楼2013-05-12 11:03:31
liaoliao8545
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 213.5
- 沙发: 1
- 帖子: 197
- 在线: 45.8小时
- 虫号: 2119920
- 注册: 2012-11-11
- 专业: 燃烧学
|
补充一下非线性方程组如下 ：比较容易看清 2C_2 x_4 √(x_1 )+2x_4+C_7 √(x_4 ) √(x_1 )+C_9 √(x_4 )-(C_1 x_3 √(x_1 )+x_3 ) d_1=0 2C_1 x_3 √(x_1 )+C_2 x_4 √(x_1 )+2x_1+x_3+C_7 √(x_4 ) √(x_1 )+C_8 √(x_1 )+C_10 √(x_1 ) √(x_2 )-(C_1 x_3 √(x_1 )+x_3 ) d_2=0 2x_2+C_10 √(x_1 ) √(x_2 )-(C_1 x_3 √(x_1 )+x_3 ) d_3=0 C_1 x_3 √(x_1 )+C_2 x_4 √(x_1 )+x_1+x_2+x_3+x_4+C_7 √(x_4 ) √(x_1 )+C_8 √(x_1 )+C_9 √(x_4 )+C_10 √(x_1 ) √(x_2 )+1/n (C_1 x_3 √(x_1 )+x_3 )=1 |
3楼2013-05-12 16:11:00
