版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

signo

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 7.6
散金: 22
红花: 1
帖子: 311
在线: 374.8小时
虫号: 1183274
注册: 2011-01-04
专业: 凝聚态物性 II ：电子结构

引用回帖:

10楼: Originally posted by onesupeng at 2013-04-18 20:20:21
你编过其中某个的程序没有？...

我以前做过的是：二元二阶微分方程组（Fortran）;

另一个是有限差分，球体的折射率是半径的函数，然后求出入射光的轨迹，这个是Matlab编的。

你说的矩阵我没做过，打靶法我看有点希望。

此处的难点在于，我不知道具这个问题的边界条件。

资料已发邮箱。

赞一下

回复此楼

11楼2013-04-18 21:18:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

signo

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 7.6
散金: 22
红花: 1
帖子: 311
在线: 374.8小时
虫号: 1183274
注册: 2011-01-04
专业: 凝聚态物性 II ：电子结构

引用回帖:

4楼: Originally posted by onesupeng at 2013-04-15 03:20:12
我可以简单地说一个思路，这个用差分法，结合matlab计算，比较容易实现。首先，假定x定义在之间，且从方程看，可能还要求x1>x0>0。把你的方程记为：eps''=f(x) eps

假定将计算域划分为N，在i点差分你的方程 ...

已把详细资料和思路发给 onesupeng ，
我们将合作发表这个文章。

赞一下

回复此楼

12楼2013-04-19 03:04:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zwb256

银虫 (正式写手)

应助: 1 (幼儿园)
金币: 114.7
散金: 150
红花: 2
帖子: 329
在线: 309.5小时
虫号: 229406
注册: 2006-03-25
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

引用回帖:

10楼: Originally posted by onesupeng at 2013-04-18 20:20:21
你编过其中某个的程序没有？...

打靶法和差分法求边值问题的代码是当例子讲的。

赞一下

回复此楼

13楼2013-04-19 09:55:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

signo

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 7.6
散金: 22
红花: 1
帖子: 311
在线: 374.8小时
虫号: 1183274
注册: 2011-01-04
专业: 凝聚态物性 II ：电子结构

引用回帖:

13楼: Originally posted by zwb256 at 2013-04-19 09:55:54
打靶法和差分法求边值问题的代码是当例子讲的。...

是啊，我看例子都是讲的方势阱和谐振子势，边界条件用的是 f(0)=f(1)=0，我这个问题中是一个x的取值只能为正的，大于0，因为是从柱坐标r 通过变量代换过来的，所以边界值我还没有想出来

赞一下

回复此楼

14楼2013-04-19 15:38:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

maxiaoran

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 14
红花: 1
帖子: 2
在线: 6.3小时
虫号: 1452307
注册: 2011-10-20
专业: 传热传质学

引用回帖:

2楼: Originally posted by zwb256 at 2013-04-14 22:44:05
可以采用打靶法或者Numerov 算法或者有限差分法求解你这个本征值问题！

如何利用打靶法求解这个本征值问题呢？可以麻烦你给解释一下吗？给出思路即可，谢谢

赞一下

回复此楼

15楼2014-01-06 17:10:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 3430
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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

赞一下(1人)

回复此楼

16楼2014-01-06 19:25:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

signo

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 7.6
散金: 22
红花: 1
帖子: 311
在线: 374.8小时
虫号: 1183274
注册: 2011-01-04
专业: 凝聚态物性 II ：电子结构

引用回帖:

16楼: Originally posted by feixiaolin at 2014-01-06 19:25:22
作变量代换，令 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 ...

你这个方法我看懂了，可是怎么求那个本征值e 呢？

赞一下

回复此楼

17楼2014-01-09 21:02:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖