{"id":772,"date":"2012-04-19T21:57:14","date_gmt":"2012-04-19T13:57:14","guid":{"rendered":"http:\/\/www.hsyyf.me\/?p=772"},"modified":"2012-04-24T20:23:39","modified_gmt":"2012-04-24T12:23:39","slug":"%e5%9b%9b%e9%98%b6%e9%be%99%e6%a0%bc%e5%ba%93%e5%a1%94%e6%b3%95%e8%a7%a3%e5%90%ab%e6%97%b6%e8%96%9b%e5%ae%9a%e8%b0%94%e6%96%b9%e7%a8%8b","status":"publish","type":"post","link":"https:\/\/www.hsyyf.me\/?p=772","title":{"rendered":"\u56db\u9636\u9f99\u683c\u5e93\u5854\u6cd5\u89e3\u542b\u65f6\u859b\u5b9a\u8c14\u65b9\u7a0b"},"content":{"rendered":"

\u56db\u9636\u9f99\u683c\u5e93\u5854\u6cd5\u89e3\u542b\u65f6\u859b\u5b9a\u8c14\u65b9\u7a0b\uff0c\u4e00\u7ef4\u7ebf\u6027\u8c10\u632f\u5b50\u9644\u52a0\u542b\u65f6\u9a71\u52a8\u529b\u3002<\/p>\n

external f
\nparameter(n=100001,m=11)
\nDimension Y(m),D(m),z(m,n),b(m),c(m,n)
\ncomplex(8) y,d,z,b,h
\nreal(8) t,a
\nt=0.0d0<\/p>\n

do j=1,m
\nY(j)=(0.0,0.0)
\nenddo
\nY(1)=(1.0d0,0.0d0)<\/p>\n

h=2500.d0\/(n-1)
\ncall grkt1(t,y,m,h,n,z,f,d,b)
\nwrite(*,*)
\ncall Psi(z)
\nend<\/p>\n

subroutine Psi(z)
\nparameter(n=100001,m=11,pi=3.1415926d0)
\ncomplex(8) z(m,n),p(20,1000)
\nreal(8) a(m,n)
\ninteger t<\/p>\n

open(200,File='test.txt')
\nopen(201,File='test1.txt')
\nopen(202,File='test2.txt')
\nopen(203,File='test3.txt')<\/p>\n

do k=1,20,1
\ndo j=1,1000
\nx=-5.0d0+j*0.01d0
\nP(k,j)=(0.d0,0.d0)
\ndo i=1,m
\np(k,j)=p(k,j)+z(i,k*5000)\/sqrt((sqrt(pi)*(2**i)*JC(i)))*exp(-x*x\/2.0d0)*H(i,x*1.0)*exp((0.d0,-1.d0)*(2.d0*i+1.d0)*k*5000*(2500.d0\/(n-1)))
\n!p(k,j)=p(k,j)+z(i,k)\/sqrt((sqrt(pi)*(2**i)*JC(i)))*exp(-x*x\/2.0d0)*H(i,x*1.0)*exp((0.d0,-1.d0)*(2.d0*i+1.d0)*k*(2500.d0\/(n-1)))
\nenddo
\nenddo
\nenddo<\/p>\n

do j=1,1000
\nwrite(200,*)(-5.d0+j*0.01d0),(abs(p(i,j))**2,i=1,2)
\nwrite(201,*)(abs(p(i,j))**2,i=3,5)
\nwrite(202,*)(abs(p(i,j))**2,i=6,8)
\nwrite(203,*)(abs(p(i,j))**2,i=9,11)
\nenddo<\/p>\n

end subroutine<\/p>\n

subroutine F(t,y,m,d)<\/p>\n

Dimension Y(m),D(m)
\ncomplex(8) y,d
\nreal(8) t
\nD(1)=-(0.d0,1.d0)*g(t)*(Y(2)*sqrt(1.d0\/2.d0)*exp(-2.d0*(0,1)*t))
\ndo j=2,m-1
\nD(j)=-(0,1)*g(t)*(Y(j-1)*sqrt((j-1)\/2.d0)*exp(2.d0*(0,1)*t)+Y(j+1)*sqrt(j\/2.)*exp(-2.d0*(0,1)*t))
\nenddo
\nD(m)=-(0,1)*g(t)*(Y(m-1)*sqrt((m-1)\/2.d0)*exp(-2.d0*(0,1)*t))
\nreturn
\nend<\/p>\n

function g(t)
\nreal(8) g,t
\npi=3.1415926535d0
\nif(t<=2000.d0)then
\ng=0.5*(1+sin(t\/2000.*pi-pi\/2.))
\n!\u00a0\u00a0 \u00a0g=t\/2000.d0
\nelse
\ng=1.00d0
\nendif
\nend<\/p>\n

real recursive function H(i,x) result(Hi)
\nselect case(i)
\ncase(1)
\nHi=1.0
\ncase(2)
\nHi=2.0*x
\ncase(3:)
\nHi=2*x*h(i-1,x)-2*(i-1)*H(i-2,x)
\nend select
\nend<\/p>\n

integer Recursive function JC(i) result(J)
\nif(i==0)then
\nJ=1
\nelse
\nJ=i*JC(i-1)
\nendif
\nend<\/p>\n

SUBROUTINE GRKT1(T,Y,M,H,N,Z,F,D,B)
\nDIMENSION Y(M),D(M),Z(M,N),A(4),B(M)
\ncomplex(8) Y,D,Z,A,B,H,X,TT
\nreal(8) t
\nA(1)=H\/2.0d0
\nA(2)=A(1)
\nA(3)=H
\nA(4)=H
\nDO I=1,M
\nZ(I,1)=Y(I)
\nenddo
\nX=T
\nDO J=2,N
\nCALL F(T,Y,M,D)
\nDO I=1,M
\nB(I)=Y(I)
\nenddo
\nDO K=1,3
\nDO I=1,M
\nY(I)=Z(I,J-1)+A(K)*D(I)
\nB(I)=B(I)+A(K+1)*D(I)\/3.0
\nenddo
\nTT=T+A(K)
\nCALL F(TT,Y,M,D)
\nenddo
\nDO I=1,M
\nY(I)=B(I)+H*D(I)\/6.0
\nenddo
\nDO I=1,M
\nZ(I,J)=Y(I)
\nenddo
\nT=T+H
\nenddo
\nT=X
\nRETURN
\nEND<\/code><\/p>\n","protected":false},"excerpt":{"rendered":"

\u56db\u9636\u9f99\u683c\u5e93\u5854\u6cd5\u89e3\u542b\u65f6\u859b\u5b9a\u8c14\u65b9\u7a0b\uff0c\u4e00\u7ef4\u7ebf\u6027\u8c10\u632f\u5b50\u9644\u52a0\u542b\u65f6\u9a71\u52a8\u529b\u3002 external f parameter(n= […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"ngg_post_thumbnail":0,"footnotes":""},"categories":[16],"tags":[108,103],"class_list":["post-772","post","type-post","status-publish","format-standard","hentry","category-16","tag-fortran","tag-103"],"_links":{"self":[{"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=\/wp\/v2\/posts\/772","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=772"}],"version-history":[{"count":0,"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=\/wp\/v2\/posts\/772\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=772"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=772"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.hsyyf.me\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=772"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}