24小时热门版块排行榜    

查看: 3010  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂(0805) +3 烟汐忆海 2026-03-16 8/400 2026-03-18 00:47 by 烟汐忆海
[考研] 301求调剂 +6 yy要上岸呀 2026-03-17 6/300 2026-03-17 23:58 by 星空星月
[考研] 268求调剂 +3 简单点0 2026-03-17 3/150 2026-03-17 23:15 by wangjihu
[考研] 085601专硕,总分342求调剂,地区不限 +4 share_joy 2026-03-16 4/200 2026-03-17 21:32 by hmn_wj
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 7/350 2026-03-17 20:20 by 花125533
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
信息提示
请填处理意见