24小时热门版块排行榜    

查看: 3032  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

signo

新虫 (正式写手)

[求助] 如何用 数值解法 求 二阶微分方程的本征值 已有1人参与

如图片中所示的二阶微分方程,其中a, b 是已知参数,需要确定 本正值 E, 就是那个小写的 ipsilon, 就是这个微分方程有本征值,以及求这个本征解。

我只听说可以用  切比雪夫多项式 的方法解这个微分方程,可是我没有找到相关类似方法的具体步凑,我只找到了这个多项式的具体形式,不知道这个尝试解的收敛性是怎么操作的。

如果大家有别的 尝试解 的方法,有参考文献资料的最好了。

微分方程.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

signo

新虫 (正式写手)

引用回帖:
4楼: Originally posted by onesupeng at 2013-04-15 03:20:12
我可以简单地说一个思路,这个用差分法,结合matlab计算,比较容易实现。首先,假定x定义在之间,且从方程看,可能还要求x1>x0>0。把你的方程记为:eps''=f(x) eps

假定将计算域划分为N,在i点差分你的方程 ...

(1)我的导师基本处于不管我的状态,所以我想给你挂个作者问题不大,我有70%的把握导师会答应:可问题是这个论文质量不高,只能整个1分的,你还有兴趣吗?

(2)因为这是从柱坐标薛定谔方程转化出来的,所以所以 x 的取值范围是0 到正无穷;

(3)我也知道微分方程需要边界条件,我的理解是:当x趋于0时,eps 为有限值;当x趋于正无穷,eps趋于0,请问这个是不是可以当做边界条件?
5楼2013-04-15 05:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

zwb256

银虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
可以采用打靶法或者Numerov 算法或者有限差分法求解你这个本征值问题!
2楼2013-04-14 22:44:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

signo

新虫 (正式写手)

引用回帖:
2楼: Originally posted by zwb256 at 2013-04-14 22:44:05
可以采用打靶法或者Numerov 算法或者有限差分法求解你这个本征值问题!

我搜了 你说的 Numerov's method,

里面有说解 Schrödinger Equation ,我这个方程就是从 薛定谔方程画出来的。

http://en.wikipedia.org/wiki/Numerov%27s_method
3楼2013-04-14 22:54:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
我可以简单地说一个思路,这个用差分法,结合matlab计算,比较容易实现。首先,假定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]的特征值。即为你需要的情况。如果说成这样你还不会,我可以帮你编程,但是你需要联系你导师发表论文挂我为作者之一。
长期招收博士生,参见http://fsl-unsw.com
4楼2013-04-15 03:20:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 08工科 320总分 求调剂 +6 梨花珞晚风 2026-03-17 6/300 2026-03-21 03:40 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
信息提示
请填处理意见