24小时热门版块排行榜    

查看: 3102  |  回复: 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典型案例

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东平Modelon

金虫 (小有名气)

【答案】应助回帖


xiegangmai: 金币+1, 谢谢参与 2012-11-25 20:04:10
你在定义
function y = f(x)
的第二行加上
x  = x
应该就可以了吧。
5楼2012-11-25 01:02:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图


感谢参与,应助指数 +1
xiegangmai: 金币+1, 应助指数-1, 谢谢参与,非应助请勿选择“应助回帖” 2012-11-24 19:45:06
不知道直接怎么设置,自己修改迭代函数吧。
showmethemoney
2楼2012-11-22 16:17:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★
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
showmethemoney
3楼2012-11-22 17:01:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cq0932

金虫 (初入文坛)

引用回帖:
3楼: Originally posted by csgt0 at 2012-11-22 17:01:14
这样行不,x有3个的,复数输出如下
                                Norm of      First-order
Iteration  Func-count         f(x)          step          optimality   CG-iterations
x(1)=1.000000+0.000 ...

请问这是用什么命令调出来的?
4楼2012-11-24 15:41:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见