24小时热门版块排行榜    

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

signo

新虫 (正式写手)

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

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

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

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

微分方程.jpg
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 19 个回答

zwb256

银虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
可以采用打靶法或者Numerov 算法或者有限差分法求解你这个本征值问题!
2楼2013-04-14 22:44:05
已阅   回复此楼   关注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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 能源材料化学课题组招收硕士研究生8-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 344求调剂 +4 knight344 2026-03-16 4/200 2026-03-17 17:27 by ruiyingmiao
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
请填处理意见