Fortran Coder

查看: 6033|回复: 2
打印 上一主题 下一主题

[非线性] 怎么求解这个超越方程,对应的解说,我写了个word,谢谢...

[复制链接]

7

帖子

2

主题

0

精华

入门

F 币
68 元
贡献
41 点

帅哥勋章

跳转到指定楼层
楼主
发表于 2014-4-11 07:49:07 | 只看该作者 |只看大图 回帖奖励 |正序浏览 |阅读模式
本帖最后由 小范范1989 于 2014-4-11 07:52 编辑

在有限深对称势阱中,如果说固定了势阱的宽度a,势阱的深度V,那么里面的本正能量E的个数和大小就确定了。
参数为:k=sqrt(2*E*m)/hbar  . β=sqrt(2*m*(V-E))/habr .  ε=k*a/2  . η=β*a/2
第一种情况,当为偶宇称时ε*tan(ε)=η,β^2+ε^2=m*v*a*a/(2*hbar*hbar)
第二种情况,当为奇宇称时 -ε*cot(ε)=η  β^2+ε^2=m*v*a*a/(2*hbar*hbar)>=pi^2/4  才出现最低的奇宇称能级。
现在王老师让我用fortran求解上面的超越方程,做出相应的图像,画出表格,即告诉一个aV,得到E的个数和大小。
课本上用的图解法,
      奇宇称                 偶宇称

clip_image002.gif (5.63 KB, 下载次数: 305)

奇宇称

奇宇称

clip_image002.gif (6.29 KB, 下载次数: 293)

偶宇称

偶宇称
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

66

帖子

5

主题

2

精华

版主

院士级水师

F 币
481 元
贡献
273 点

管理勋章帅哥勋章爱心勋章规矩勋章

QQ
板凳
发表于 2014-6-20 11:05:14 | 只看该作者
这是量子力学导论上的问题吧,以前学过的,也是使用图解法求得,貌似这个在以前的PF论坛上还有见过(貌似哈)。你这个问题描述的太专业化了,不适宜于非专业高手理解,建议你说的在直白些,比如求解的具体方程的形式,判定的条件等,写的在明确一些,估计就会有人帮你了,最后吐槽一句,你自己编写个,添加上标注哪怕有错误在发出来,也会有更多的人帮你解决啊,你这问的就是直接叫别人帮你干活啊。当然这么多天过去了,你是不是已经解决了?发出来给大家分享下。
科研穷三代,读博毁一生

8

帖子

3

主题

0

精华

入门

F 币
30 元
贡献
18 点
沙发
发表于 2014-4-11 09:23:00 | 只看该作者
WHAT THE HELL IS IS
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-22 14:18

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表