| 查看: 2784 | 回复: 18 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
如何用 数值解法 求 二阶微分方程的本征值已有1人参与
|
|||||
|
如图片中所示的二阶微分方程，其中a, b 是已知参数，需要确定 本正值 E, 就是那个小写的 ipsilon, 就是这个微分方程有本征值，以及求这个本征解。 我只听说可以用 切比雪夫多项式 的方法解这个微分方程，可是我没有找到相关类似方法的具体步凑，我只找到了这个多项式的具体形式，不知道这个尝试解的收敛性是怎么操作的。 如果大家有别的 尝试解 的方法，有参考文献资料的最好了。  微分方程.jpg |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 2930
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
signo: 金币+50, ★★★很有帮助 2014-01-09 01:20:27
signo: 金币+50, ★★★很有帮助 2014-01-09 01:20:27
|
作变量代换，令 p(x)=y'(x)/y(x) 可得解析解。 原式变为 p'(x)+p^2(x)=x^2+a/x+b/x^2+e 再取 p(x)=p1(x)+p2(x)+p3(x)+p4(x) p1(x)+p1^2(x)=a*x^2; p2(x)+p2^2(x)=a/x; p3(x)+p3^2(x)=a/x^2; p4(x)+p4^2(x)=a; 见附件 -------------------- 最后，{lny(x)}'=p(x) y(x)=exp{Intergrate[p(x),x]}  001.JPG  002.JPG  003.JPG  004.JPG |
16楼2014-01-06 19:25:22
2楼2013-04-14 22:44:05
|
我搜了 你说的 Numerov's method， 里面有说解 Schrödinger Equation ，我这个方程就是从 薛定谔方程画出来的。 http://en.wikipedia.org/wiki/Numerov%27s_method |
3楼2013-04-14 22:54:32
onesupeng
金虫 (职业作家)
- 数学EPI: 17
- 应助: 256 (大学生)
- 贵宾: 1.36
- 金币: 2336.2
- 散金: 9212
- 红花: 92
- 帖子: 4583
- 在线: 1303.8小时
- 虫号: 394701
- 注册: 2007-06-07
- 专业: 流体力学
【答案】应助回帖
感谢参与，应助指数 +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]的特征值。即为你需要的情况。如果说成这样你还不会，我可以帮你编程，但是你需要联系你导师发表论文挂我为作者之一。 |
4楼2013-04-15 03:20:12