版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

cq0932

金蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 734.8
帖子: 44
在線: 20.2小時(shí)
蟲(chóng)號(hào): 965611
注冊(cè): 2010-03-09
性別: GG
專業(yè): 機(jī)械測(cè)試?yán)碚撆c技術(shù)

[求助] 求教fsolve如何查看每次迭代的x值？

利用fsolve可以求解非線性方程的數(shù)值解，但是只有最后的結(jié)果，沒(méi)有每步迭代的x取值。
在求解之前設(shè)置迭代顯示的選項(xiàng)
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

輸出結(jié)果只有f(x)的值，請(qǐng)問(wèn)如何設(shè)置才能看到每步迭代的x取值？

回復(fù)此樓

1樓 2012-11-21 21:03:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cq0932

金蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 734.8
帖子: 44
在線: 20.2小時(shí)
蟲(chóng)號(hào): 965611
注冊(cè): 2010-03-09
性別: GG
專業(yè): 機(jī)械測(cè)試?yán)碚撆c技術(shù)

引用回帖:

3樓: Originally posted by csgt0 at 2012-11-22 17:01:14
這樣行不，x有3個(gè)的，復(fù)數(shù)輸出如下
Norm of First-order
Iteration Func-count f(x) step optimality CG-iterations
x(1)=1.000000+0.000 ...

請(qǐng)問(wèn)這是用什么命令調(diào)出來(lái)的？

贊一下

回復(fù)此樓

4樓2012-11-24 15:41:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 7 個(gè)回答

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

專家經(jīng)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲(chóng)號(hào): 449046
注冊(cè): 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 應(yīng)助指數(shù)-1, 謝謝參與，非應(yīng)助請(qǐng)勿選擇“應(yīng)助回帖” 2012-11-24 19:45:06

不知道直接怎么設(shè)置，自己修改迭代函數(shù)吧。

贊一下

回復(fù)此樓

showmethemoney

2樓2012-11-22 16:17:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

專家經(jīng)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲(chóng)號(hào): 449046
注冊(cè): 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

★ ★
xiegangmai: 金幣+2, 謝謝參與 2012-11-24 19:45:25

這樣行不，x有3個(gè)的，復(fù)數(shù)輸出如下
                              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

贊一下

回復(fù)此樓

showmethemoney

3樓2012-11-22 17:01:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

東平Modelon

金蟲(chóng) (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 562.8
紅花: 9
帖子: 91
在線: 95小時(shí)
蟲(chóng)號(hào): 2144786
注冊(cè): 2012-11-24
專業(yè): 動(dòng)力學(xué)與控制

【答案】應(yīng)助回帖

★
xiegangmai: 金幣+1, 謝謝參與 2012-11-25 20:04:10

你在定義
function y = f(x)
的第二行加上
x = x
應(yīng)該就可以了吧。

贊一下

回復(fù)此樓

5樓2012-11-25 01:02:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 7 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風(fēng) 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 314求調(diào)劑 +7	無(wú)懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky