 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3013 | 回復(fù): 18 | ||||
[求助]
如何用 數(shù)值解法 求 二階微分方程的本征值 已有1人參與
|
||||
|
如圖片中所示的二階微分方程，其中a, b 是已知參數(shù)，需要確定 本正值 E, 就是那個小寫的 ipsilon, 就是這個微分方程有本征值，以及求這個本征解。 我只聽說可以用 切比雪夫多項式 的方法解這個微分方程，可是我沒有找到相關(guān)類似方法的具體步湊，我只找到了這個多項式的具體形式，不知道這個嘗試解的收斂性是怎么操作的。 如果大家有別的 嘗試解 的方法，有參考文獻(xiàn)資料的最好了。 微分方程.jpg |
回復(fù)此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 含偏微分之絕對值項的一階偏微分方程組的數(shù)值解法
已經(jīng)有7人回復(fù)
- 2011注冊化工工程師基礎(chǔ)考試大綱 是不是和往年的有改變？？？
已經(jīng)有5人回復(fù)
- 2011年注冊化工工程師基礎(chǔ)考試大綱
已經(jīng)有25人回復(fù)
- 【求助】請教非齊次常微分方程組的解析解法
已經(jīng)有4人回復(fù)
- 注冊化工工程師
已經(jīng)有10人回復(fù)
- 【求助】求2010年注冊化工工程師基礎(chǔ)考試的考試大綱
已經(jīng)有4人回復(fù)
- 求2010年注冊化工工程師基礎(chǔ)考試的考試大綱
已經(jīng)有3人回復(fù)
- 【求助】一個二階常微分方程的解法（重獎金幣！）
已經(jīng)有16人回復(fù)
- 求張文生《科學(xué)計算中的偏微分方程有限差分法》和胡建偉，湯懷民《微分方程數(shù)值解法》
已經(jīng)有3人回復(fù)
- 【求助】偏微分方程解法
已經(jīng)有19人回復(fù)
- 【信息】環(huán)保工程師報考條件
已經(jīng)有11人回復(fù)
- 【分享】2009年注冊環(huán)保工程師基礎(chǔ)考試大綱
已經(jīng)有11人回復(fù)
2樓2013-04-14 22:44:05
onesupeng
金蟲 (職業(yè)作家)
- 數(shù)學(xué)EPI: 17
- 應(yīng)助: 256 (大學(xué)生)
- 貴賓: 1.36
- 金幣: 2336.2
- 散金: 9224
- 紅花: 92
- 帖子: 4583
- 在線: 1303.8小時
- 蟲號: 394701
- 注冊: 2007-06-07
- 專業(yè): 流體力學(xué)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
我可以簡單地說一個思路，這個用差分法，結(jié)合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 當(dāng)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)系你導(dǎo)師發(fā)表論文掛我為作者之一。 |
4樓2013-04-15 03:20:12
feixiaolin
榮譽版主 (文壇精英)
-
專家經(jīng)驗: +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3430
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時
- 蟲號: 2139575
- 注冊: 2012-11-21
- 專業(yè): 光學(xué)信息獲取與處理
- 管轄: 數(shù)學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
signo: 金幣+50, ★★★很有幫助 2014-01-09 01:20:27
signo: 金幣+50, ★★★很有幫助 2014-01-09 01:20:27
|
作變量代換，令 p(x)=y'(x)/y(x) 可得解析解。 原式變?yōu)?p'(x)+p^2(x)=x^2+a/x+b/x^2+e 再取 p(x)=p1(x)+p2(x)+p3(x)+p4(x) p1(x)+p1^2(x)=a*x^2; p2(x)+p2^2(x)=a/x; p3(x)+p3^2(x)=a/x^2; p4(x)+p4^2(x)=a; 見附件 -------------------- 最后，{lny(x)}'=p(x) y(x)=exp{Intergrate[p(x),x]} 001.JPG 002.JPG 003.JPG 004.JPG |
16樓2014-01-06 19:25:22
|
我搜了 你說的 Numerov's method， 里面有說解 Schrödinger Equation ，我這個方程就是從 薛定諤方程畫出來的。 http://en.wikipedia.org/wiki/Numerov%27s_method |
3樓2013-04-14 22:54:32
5樓2013-04-15 05:01:12
onesupeng
金蟲 (職業(yè)作家)
- 數(shù)學(xué)EPI: 17
- 應(yīng)助: 256 (大學(xué)生)
- 貴賓: 1.36
- 金幣: 2336.2
- 散金: 9224
- 紅花: 92
- 帖子: 4583
- 在線: 1303.8小時
- 蟲號: 394701
- 注冊: 2007-06-07
- 專業(yè): 流體力學(xué)
6樓2013-04-18 03:43:01
7樓2013-04-18 04:45:15
onesupeng
金蟲 (職業(yè)作家)
- 數(shù)學(xué)EPI: 17
- 應(yīng)助: 256 (大學(xué)生)
- 貴賓: 1.36
- 金幣: 2336.2
- 散金: 9224
- 紅花: 92
- 帖子: 4583
- 在線: 1303.8小時
- 蟲號: 394701
- 注冊: 2007-06-07
- 專業(yè): 流體力學(xué)
8樓2013-04-18 06:38:33
9樓2013-04-18 10:05:48
onesupeng
金蟲 (職業(yè)作家)
- 數(shù)學(xué)EPI: 17
- 應(yīng)助: 256 (大學(xué)生)
- 貴賓: 1.36
- 金幣: 2336.2
- 散金: 9224
- 紅花: 92
- 帖子: 4583
- 在線: 1303.8小時
- 蟲號: 394701
- 注冊: 2007-06-07
- 專業(yè): 流體力學(xué)
10樓2013-04-18 20:20:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 signo 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +6 | 呼呼呼。。。。 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|