| 6 | 1/1 | 返回列表 |
| 查看: 1101 | 回復(fù): 5 | |||
chuzhenhui金蟲(chóng) (小有名氣)
|
[交流]
【求助】用fsolve求解非線(xiàn)性一元方程 已有2人參與
|
|
我用fsolve()求解非線(xiàn)性方程,程序如下: >> 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)問(wèn)紅色部分是什么意思?我的這個(gè)解好像不是很準(zhǔn)確,怎么改進(jìn)。 要用optmiset設(shè)置參數(shù)嗎? 大俠指導(dǎo)下,感激不盡! |

鐵桿木蟲(chóng) (著名寫(xiě)手)
方丈大師
金蟲(chóng) (小有名氣)
|
謝謝,不好意思我粗心了,其實(shí)我想解的方程是有好多項(xiàng)加起來(lái)的,我圖方便拿了其中的一項(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 ] |

鐵桿木蟲(chóng) (著名寫(xiě)手)
方丈大師
|
你可以將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 ] |
鐵桿木蟲(chóng) (著名寫(xiě)手)
方丈大師
金蟲(chóng) (小有名氣)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279分求調(diào)劑 一志愿211 +3 | chaojifeixia 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽(yáng)CF 2026-03-18 | 4/200 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專(zhuān)利,省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 材料專(zhuān)碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 材料與化工專(zhuān)碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 一志愿,福州大學(xué)材料專(zhuān)碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|