24小时热门版块排行榜    

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

signo

新虫 (正式写手)

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

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

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

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

微分方程.jpg
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 0703化学336分求调剂 +3 zbzihdhd 2026-03-15 4/200 2026-03-16 23:52 by zbzihdhd
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
信息提示
请填处理意见