| 查看: 3107 | 回复: 6 | ||||
cq0932金虫 (初入文坛)
|
[求助]
求教fsolve如何查看每次迭代的x值？
|
|
利用fsolve可以求解非线性方程的数值解，但是只有最后的结果，没有每步迭代的x取值。 在求解之前设置迭代显示的选项 options=optimset('Display','iter'); 输出窗口中显示如下： Norm of First-order Trust-region Iteration Func-count f(x) step optimality radius 0 2 2.61482e-007 0.000171 1 1 3 2.61482e-007 0.00152725 0.000171 1 2 5 1.51338e-007 0.000381811 0.000109 0.000382 3 6 1.51338e-007 0.000954528 0.000109 0.000955 4 8 8.78862e-008 0.000238632 8.52e-005 0.000239 5 10 7.98377e-008 0.00059658 8.75e-005 0.000597 6 11 7.98377e-008 0.00059658 8.75e-005 0.000597 7 13 6.53061e-008 0.000149145 8.2e-005 0.000149 8 15 5.49734e-008 0.000149145 7.79e-005 0.000149 9 17 4.84382e-008 0.000149145 7.17e-005 0.000149 10 19 4.1506e-008 0.000149145 6.76e-005 0.000149 11 21 3.68301e-008 0.000149145 6.25e-005 0.000149 12 23 3.1814e-008 0.000149145 5.91e-005 0.000149 13 25 2.87067e-008 0.000149145 5.52e-005 0.000149 14 27 2.5256e-008 0.000149145 5.26e-005 0.000149 15 29 2.33751e-008 0.000149145 4.98e-005 0.000149 16 31 2.10974e-008 0.000149145 4.8e-005 0.000149 输出结果只有f(x)的值，请问如何设置才能看到每步迭代的x取值？ |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab典型案例 |
» 猜你喜欢
26博士申请
已经有3人回复
-
268求调剂
已经有3人回复
-
302求调剂
已经有5人回复
-
311求调剂
已经有10人回复
-
被我言中：新模板不强调格式了，假专家开始管格式了
已经有4人回复
-
303求调剂
已经有5人回复
-
考研求调剂
已经有4人回复
-
085601专硕，总分342求调剂，地区不限
已经有4人回复
-
085601材料工程专硕求调剂
已经有5人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助，方程组求解出数值解后，第一个解是下个循环的一个参数
已经有1人回复
- 求助，方程组求解出数值解后，第一个解是下个循环的一个参数
已经有0人回复
- 曲线的光滑连接-插值求教
已经有4人回复
2楼2012-11-22 16:17:23
【答案】应助回帖
★ ★
xiegangmai: 金币+2, 谢谢参与 2012-11-24 19:45:25
xiegangmai: 金币+2, 谢谢参与 2012-11-24 19:45:25
|
这样行不，x有3个的，复数输出如下 Norm of First-order Iteration Func-count f(x) step optimality CG-iterations x(1)=1.000000+0.000000i x(2)=1.000000+0.000000i x(3)=1.000000+0.000000i 0 4 267.268 89.2 x(1)=1.000000+0.000000i x(2)=1.000000+0.000000i x(3)=1.000000+0.000000i 1 8 267.268 7.26736 89.2 0 x(1)=0.588062+0.400105i x(2)=0.223569+0.191305i x(3)=-0.158092+-0.995244i 2 12 73.6223 1.81684 13.8 0 x(1)=1.914850+-0.478769i x(2)=-0.760707+1.025506i x(3)=-2.589293+-2.754627i 3 16 24.8735 3.63368 36.1 0 x(1)=1.473642+-0.017474i x(2)=-0.533834+1.507987i x(3)=-2.481998+-3.272860i 4 20 0.30956 0.985798 2.03 0 x(1)=1.513119+-0.047554i x(2)=-0.438356+1.425158i x(3)=-2.617916+-3.192791i 5 24 0.000498233 0.208146 0.126 0 x(1)=1.514985+-0.047818i x(2)=-0.436067+1.426855i x(3)=-2.617734+-3.195948i 6 28 1.45562e-010 0.00465524 7.53e-005 0 x(1)=1.514984+-0.047818i x(2)=-0.436070+1.426856i x(3)=-2.617733+-3.195949i 7 32 2.43937e-023 3.0415e-006 1.52e-011 0 Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. xlast = 1.514983763593827 - 0.047817849126132i -0.436069706071287 + 1.426855556443150i -2.617733027757387 - 3.195948930097455i |
3楼2012-11-22 17:01:14
cq0932
金虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 734.8
- 帖子: 44
- 在线: 20.2小时
- 虫号: 965611
- 注册: 2010-03-09
- 性别: GG
- 专业: 机械测试理论与技术
4楼2012-11-24 15:41:32
5楼2012-11-25 01:02:09
6楼2012-11-26 11:09:52
szuwusongbin
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3115.4
- 散金: 10
- 红花: 2
- 帖子: 410
- 在线: 152.8小时
- 虫号: 3130374
- 注册: 2014-04-12
- 性别: GG
- 专业: 概率论与随机分析
7楼2016-01-28 00:17:39
50