| 查看: 430 | 回复: 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;]; |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有9人回复
-
为什么一个当量的苯甲醛反应不完呢
已经有5人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
为什么nbs上溴 没有产物点出现呢
已经有8人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求助一下有机合成大神
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
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