| 查看: 541 | 回复: 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能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
帮忙看看fileCode
已经有5人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
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
100