| 5 | 1/1 | 返回列表 |
| 查看: 1107 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
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)下,感激不盡。 |

鐵桿木蟲 (著名寫手)
方丈大師
|
你可以將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 ] |
鐵桿木蟲 (著名寫手)
方丈大師
金蟲 (小有名氣)
|
謝謝,不好意思我粗心了,其實(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 ] |

鐵桿木蟲 (著名寫手)
方丈大師
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 工科0856求調(diào)劑 +4 | 沐析汀汀 2026-03-21 | 4/200 |
|
|
[考研] 081700 調(diào)劑 267分 +4 | 迷人的哈哈 2026-03-23 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +6 | 紙魚ly 2026-03-21 | 6/300 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +4 | niko- 2026-03-22 | 4/200 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +5 |  ̄^ ̄゜汗 2026-03-19 | 7/350 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|