24小时热门版块排行榜    

查看: 3006  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285化工学硕求调剂(081700) +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 302求调剂 +8 负心者当诛 2026-03-11 8/400 2026-03-17 09:05 by 哦哦123
[考研] 085601求调剂 +3 Du.11 2026-03-16 3/150 2026-03-16 20:42 by 无际的草原
[文学芳草园] 伙伴们,祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 309求调剂 +5 花与叶@ 2026-03-10 5/250 2026-03-16 14:13 by 哦哦123
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
请填处理意见