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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

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
金币: 3880
散金: 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的回帖

feixiaolin

荣誉版主 (文坛精英)

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

代回原式。

回复此楼

18楼2014-01-09 21:38:32

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

wfjiangjwf

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 368.2
帖子: 25
在线: 27.2小时
虫号: 1746266
注册: 2012-04-09
性别: GG
专业: 光学和光电子学

楼主想请教一下您是如何解决的可以吗？万分感谢

发自小木虫IOS客户端

赞一下

回复此楼

19楼2022-03-28 23:25:12

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

相关版块跳转我要订阅楼主 signo 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 304求调剂 +7	c297914 2026-04-05	8/400	2026-04-05 22:13 by hemengdong
[考研] 315求调剂 +5	＆123456789 2026-04-05	5/250	2026-04-05 19:55 by nepu_uu
[考研] 材料专硕(0856) 339分求调剂 +10	哈哈哈鹅哈哈哈 2026-04-04	10/500	2026-04-05 18:51 by 蓝云思雨
[考研] 329求调剂 +17	miaodesi 2026-04-02	20/1000	2026-04-05 18:33 by 蓝云思雨
[考研] 0703化学321分求调剂 +17	三dd. 2026-03-30	18/900	2026-04-05 18:07 by 蓝云思雨
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业083200 +4	what张 2026-04-04	5/250	2026-04-05 14:07 by chw1980_0
[考研] 工科277分求调剂材料 +8	上了上了上哦 2026-04-05	9/450	2026-04-05 13:05 by wwytracy
[考研] 085410人工智能初试316分求调剂 +6	残星拂曙 2026-03-31	6/300	2026-04-05 12:15 by rainbow11
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3	5ourr 2026-04-03	3/150	2026-04-05 10:40 by JOKER0401
[考研] 301求调剂 +14	A_JiXing 2026-04-01	14/700	2026-04-03 18:31 by ls刘帅
[考研] 085501一志愿天工大，机械专硕求调剂，跨材料 +3	33上 2026-04-03	3/150	2026-04-03 14:08 by 1753564080
[考研] 一志愿安徽大学0817化学工程与技术，求调剂 +14	我不是只因 2026-04-02	15/750	2026-04-03 09:49 by 蓝云思雨
[考研] 调剂 +3	osbbx 2026-04-02	3/150	2026-04-03 07:47 by cc8418
[考研] 环境工程297分求调剂一志愿杭高院 +15	GENJIOW 2026-03-31	16/800	2026-04-02 17:56 by cyh—315
[论文投稿] chinese chemical letters英文版投稿求助 120+4	Yishengeryi 2026-03-30	6/300	2026-04-02 17:19 by Yishengeryi
[考研] 346求调剂 +5	郑诚乐 2026-04-02	5/250	2026-04-02 16:38 by SZW_UJN
[考研] 调剂 +3	好好读书。 2026-04-01	6/300	2026-04-02 15:49 by liumengping
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8	双马尾痞老板2 2026-04-02	9/450	2026-04-02 14:45 by 5896
[考研] 材料专硕322分 +11	哈哈哈吼吼吼哈 2026-04-01	11/550	2026-04-02 10:52 by lnilvy
[考研] 07生物学求调剂一志愿同济大学359分 +3	LAMC. 2026-03-30	3/150	2026-04-02 10:26 by 18828373951