| 查看: 455 | 回复: 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;]; |
回复此楼» 猜你喜欢
焦虑
已经有13人回复
-
一志愿985，本科211，0817化学工程与技术319求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有9人回复
-
346求调剂[0856]
已经有5人回复
-
285化工学硕求调剂（081700）
已经有11人回复
-
267一志愿南京工业大学0817化工求调剂
已经有9人回复
-
274求调剂
已经有6人回复
-
求材料调剂
已经有10人回复
-
本科郑州大学物理学院，一志愿华科070200学硕，346求调剂
已经有4人回复
-
国自科面上基金字体
已经有8人回复
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
liaoliao8545
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 213.5
- 沙发: 1
- 帖子: 197
- 在线: 45.8小时
- 虫号: 2119920
- 注册: 2012-11-11
- 专业: 燃烧学
4楼2013-05-12 16:12:34
