| 查看: 3051 | 回復: 18 | |||
[求助]
如何用 數(shù)值解法 求 二階微分方程的本征值 已有1人參與
|
|||
金蟲 (職業(yè)作家)
|
我可以簡單地說一個思路,這個用差分法,結合matlab計算,比較容易實現(xiàn)。首先,假定x定義在[x0,x1]之間,且從方程看,可能還要求x1>x0>0。把你的方程記為:eps''=f(x) eps 假定將計算域劃分為N,在i點差分你的方程,大約有: eps_(i-1)-2*eps_i+eps_(i+1)=f(x_i)dx^2 eps_i -->1*eps_(i-1)-[2+f(x_i)dx^2]*eps_(i)+1*eps_(i+1)=0 當i=1時,利用邊界條件,可以把eps_(i-1)用eps_i+eps_(i+1)表示,方程可以寫為 -[2+f(x_i)dx^2+a]*eps_(i)+[1+b]*eps_(i+1)=0 的形式。同理可求i=N點的情況。 于是方程離散為: [A]*[eps]=0 其中,[A]為N*N的三對角矩陣: a1 b1 0 0 .................... b a b 0 .................... .................................... .........................b a b ............................b2 a2 于是,使用matlab求[A]的特征值。即為你需要的情況。如果說成這樣你還不會,我可以幫你編程,但是你需要聯(lián)系你導師發(fā)表論文掛我為作者之一。 |

榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
|
我搜了 你說的 Numerov's method, 里面有說解 Schrödinger Equation ,我這個方程就是從 薛定諤方程畫出來的。 http://en.wikipedia.org/wiki/Numerov%27s_method |
金蟲 (職業(yè)作家)

金蟲 (職業(yè)作家)

金蟲 (職業(yè)作家)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[碩博家園] 招收生物學/細胞生物學調劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學碩求調劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[碩博家園] 求調劑 330分 085600材料與化工 +3 | gqhhh 2026-03-22 | 3/150 |
|
|
[考研] 一志愿太原理工安全工程300分,求調劑 +5 | 0857求調劑. 2026-03-24 | 6/300 |
|
|
[考研] 070300求調劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 一志愿南京航空航天大學材料學碩求調劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 291求調劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 材料292調劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考博] 26申博 +3 | 加油沖。 2026-03-26 | 3/150 |
|
|
[考研] 274求調劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 348求調劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 085601 材料工程 313分 求調劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 286求調劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 289求調劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 一志愿 南京郵電大學 288分 材料考研 求調劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 282求調劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 292求調劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研] 280分求調劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|