| 6 | 1/1 | 返回列表 |
| 查看: 1100 | 回復(fù): 5 | |||
chuzhenhui金蟲 (小有名氣)
|
[交流]
【求助】用fsolve求解非線性一元方程 已有2人參與
|
|
我用fsolve()求解非線性方程,程序如下: >> t=@(D)(8/pi^2)*exp(-12000000*pi^2*D)-0.092034; >> [D,f,flag,out]=fsolve(t,0) 運(yùn)行后結(jié)果如下: Optimization terminated: norm of relative change in X is less than max(options.TolX^2,eps) and sum-of-squares of function values is less than sqrt(options.TolFun). D = 1.8370e-008 f = -8.6133e-006 flag = 2 out = iterations: 22 funcCount: 36 algorithm: 'trust-region dogleg' firstorderopt: 44.0855 message: [1x167 char] 請(qǐng)問紅色部分是什么意思?我的這個(gè)解好像不是很準(zhǔn)確,怎么改進(jìn)啊? 要用optmiset設(shè)置參數(shù)嗎? 大俠指導(dǎo)下,感激不盡! |

鐵桿木蟲 (著名寫手)
方丈大師
金蟲 (小有名氣)
|
謝謝,不好意思我粗心了,其實(shí)我想解的方程是有好多項(xiàng)加起來的,我圖方便拿了其中的一項(xiàng),具體的方程是 s=0 for n=0:inf s=s+(-8/pi^2)*(1/(2*n+1)^2)*exp(-(2*n+1)^2*pi^2*12000000*D) end 要解的方程是由上面得到的 s-0.907966=0 我拿了其中幾項(xiàng)用同樣的方法運(yùn)行程序,得到上面的紅色字的提示。 t=@(D)(8/pi^2)*(exp(-12000000*pi^2*D)+(1/9)*exp(-9*12000000*pi^2*D)+(1/25)*exp(-25*12000000*pi^2*D))-0.907966; [D,f,flag,out]=fsolve(t,0) Optimization terminated: norm of relative change in X is less than max(options.TolX^2,eps) and sum-of-squares of function values is less than sqrt(options.TolFun). D = 9.2703e-011 f = -1.7886e-005 flag = 2 我想知道的是紅色字提示的是什么意思,我得到的解與精確解差的多不多,有什么提高的方法。希望得到指點(diǎn)。 [ Last edited by chuzhenhui on 2010-5-13 at 22:42 ] |

鐵桿木蟲 (著名寫手)
方丈大師
|
你可以將12000000*D作為新未知數(shù)D,結(jié)果最后除去12000000即可 >>t=@(D)(8/pi^2)*(exp(-pi^2*D)+(1/9)*exp(-9*pi^2*D)+(1/25)*exp(-25*pi^2*D))-0.907966; >>[D,f,flag,out]=fsolve(t,0) D = 0.001111597655874 f = 5.475275788313638e-011 flag = 1 out = iterations: 3 funcCount: 8 algorithm: 'trust-region dogleg' firstorderopt: 1.163031023008677e-009 message: [1x695 char] >> D/12000000 ans = 9.263313798950567e-011 [ Last edited by change0618 on 2010-5-14 at 12:47 ] |
鐵桿木蟲 (著名寫手)
方丈大師
金蟲 (小有名氣)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|---|---|---|---|---|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[基金申請(qǐng)]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|