24小时热门版块排行榜

返回列表

zmoyanyu

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 1283.3
散金: 40
帖子: 397
在线: 82.7小时
虫号: 3739808
注册: 2015-03-15
专业: 光学

[求助] 求高手帮忙解答这个程序问题已有1人参与

没有提示warning,也没有提示error,但是运行的时候，却说触发了一个断点，让选择继续或中断，选择继续后，根据得不到想要的数据。这到底怎么回事呢？求高手帮忙看看，不胜感激

subroutine bohanshu(n,bhs,energy)
IMPLICIT NONE
!+++++++++++++++++++++++++++++++++++++++++++++++声明
integer,parameter::num=10 !num1代表n，n表示第n个态
integer,parameter::num1=100 !num2代表x,即坐标
real, parameter::pi=3.14
real::bhs(num,num1)    !放波函数的数组
real::energy(num)    !存放能量本征值得数组
integer::n
real::x,a
real ::xx
real ::dx
integer:: j

!++++++++++++++++++++++++++++++++++++++++++++++++赋值
xx=2.0                   !空间宽度
dx=xx/num1
a=1.50                   !势阱宽度

!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++运算

do n=1,num    !++++++++++++++++++++++++++++n循环
energy(n)=pi**2*(n**2)/(2*a**2)
  do j=1,num1 !+++++++++++++++++++++++++++x循环
   x=dx*j
   bhs(n,j)=sqrt(2/a)*sin(n*pi*x/a)
end do
end do

return
end

!******************************************************求微扰矩阵Hnk以及exp(energy(n)-energy(k))*t的子程序

program bobo
IMPLICIT NONE
!+++++++++++++++++++++++++++++++++++++++++++++++声明
integer,parameter::num=10    !在该程序中选取的最高能量本征态对应的量子数
integer,parameter::num1=100    !位置的份数
integer,parameter::num2=200    !时间的份数
real ,parameter::pi=3.14
real::H                         !哈密顿量中的微扰部分
real::bhs(num,num1)          !放波函数的数组
real::WRJZnk(num)             !存放微扰矩阵的数组
real::energy(num)          !存放能量本征值的数组
complex::Wnk(num,num2)             !用来存放exp(energy(n)-energy(k))*t的数组
complex::b                      !虚数单位i
real::x                         !位置变量
real::xx                      !能到达的最大位置
real ::dx                      !表示位置的变化
integer::k                      !第k个态的量子数
integer n                      !是第n个态的量子数，同时也是能量的循环指标
integer ::j                   !位置循环指标
integer::i                      !时间循环指标
real::t                         !时间
real::tx                      !微扰作用在系统上的截止时间
real::dt                         !数值模拟的时间间隔
!++++++++++++++++++++++++++++++++++++++++++++++++赋值
tx=10.0
xx=2.0
H=1.105 !微扰
!k=4                      !代表初始状态为第4个能量本征态
dx=xx/num1                !位置的数值间隔
dt=tx/num2                !时间的数值间隔
b=(0,1)
!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++运算

open(unit=60,file='WRJZnk.txt')
open(unit=70,file='Wnk.txt')
call bohanshu(n,bhs,energy)

do  n=1,num          !#################################################################################能量循环
   do j=1,num1    !#################################################################################位置循环
   WRJZnk(n)=WRJZnk(n)+bhs(n,j)*H*bhs(k,j)*dx       !###########################################################计算出了Hnk
   end do
   write(60,*)WRJZnk(n)
end do

do n=1,num          !################################################################################能量循环
do i=1,num2       !################################################################################时间循环
t=i*dt
Wnk(n,i)=exp(b*(energy(n)-energy(k))*t)

write(70,*)Wnk(n,i)
end do
end do

end  program bobo

回复此楼