| 查看: 582 | 回复: 2 | ||
zidingxiang金虫 (小有名气)
|
[求助]
运行结果与设置精度不符，结果应与初值接近才对 已有1人参与
|
|
1.首先建立函数myfun,储存方程组编程如下将myfun.m保存到工作路径中: function F = myfun( x ) %this function designed for solving nonlinear equations a = x(1); b = x(2); c = x(3); d = x(4); e = x(5); f = x(6); F = [a+b*10+c*10^2+(d+e*10)/(0.6438^2-f^2)-2*1.620015/(1.620015^2-1)*(-0.87+2.6); a+b*10+c*10^2+(d+e*10)/(0.6328^2-f^2)-2*1.620507/(1.620507^2-1)*(-0.87+2.7); a+b*10+c*10^2+(d+e*10)/(0.5893^2-f^2)-2*1.622702/(1.622702^2-1)*(-0.87+2.8); a+b*10+c*10^2+(d+e*10)/(0.5461^2-f^2)-2*1.625401/(1.625401^2-1)*(-0.87+2.9); a+b*10+c*10^2+(d+e*10)/(0.4800^2-f^2)-2*1.631029/(1.631029^2-1)*(-0.88+3.3); a+b*10+c*10^2+(d+e*10)/(0.4358^2-f^2)-2*1.636369/(1.636369^2-1)*(-0.88+3.6)]; end 2.运行程序 x0=[-3e-005,-1.5e-008,-5.0e-011,3.0e-007,2.0e-010,0]; options=optimset('Display','iter','TolFun',1e-020,'TolX',1e-020,'MaxIter',1e100); format long; [x,fval]=fsolve(@myfun,x0,options) 运行结果显示如下： Norm of First-order Trust-region Iteration Func-count f(x) step optimality radius 0 7 106.837 2.5e+003 1 1 14 1.58981 0.0393664 28.2 1 2 21 0.0425606 0.0513838 36.8 1 3 28 0.0198142 0.000578732 0.414 1 4 35 0.0194798 0.000755402 0.541 1 5 42 0.0194749 8.50804e-006 0.00609 1 6 49 0.0194748 1.11053e-005 0.00795 1 7 56 0.0194748 1.25077e-007 8.95e-005 1 8 63 0.0194748 1.6325e-007 0.000117 1 9 70 0.0194748 1.84156e-009 1.32e-006 1 10 77 0.0194748 2.41651e-009 1.58e-006 1 11 84 0.0194748 2.49273e-011 1.22e-008 1 Optimizer appears to be converging to a point which is not a root. Relative function value changing by less than max(options.TolFun^2,eps) but sum-of-squares of function values is greater than or equal to sqrt(options.TolFun) Try again with a new starting guess. x = 0.000167225746880 0.001972242468804 0.019722574638040 0.006232656757062 0.062323567700574 0.000000000137224 fval = 0.060290300323624 -0.083434220551602 -0.030535103096213 0.083729176350573 -0.030583681232919 0.000533528265654 注：函数精度远远小于设定精度值；参数值应接近于设定值，但相差很远。 |
回复此楼» 猜你喜欢
重要来源：本周末出结果
已经有6人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab解非线性方程组程序纠错与正确运行
已经有6人回复
zidingxiang
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 774.8
- 散金: 69
- 红花: 1
- 帖子: 114
- 在线: 31.7小时
- 虫号: 958621
- 注册: 2010-03-01
- 性别: MM
- 专业: 玻璃材料
2楼2014-07-08 08:23:19
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2014-07-08 09:14:23
