[½»Á÷] ¡¾ÇóÖú¡¿¸ßÈËÄܸøÖ¸µãϲ»£¿¡¾Òѽâ¾ö¡¿ ÒÑÓÐ2È˲ÎÓë

µÚÒ»´Î±à¾ØÕóÀïÓи´ÊýµÄ³ÌÐò£¡¹»·³ÈË°¡£¡ ˼·£ºÔÚÍⲿ¶¨ÒåһϵÁбäÁ¿£¬t,w,f,piΪʵÐÍ£¬fiΪ¸´ÐÍ£»¾ØÕóh0Ϊ¸´ÐÍ£¬ÒªÇóËüµÄÌØÕ÷Öµ£¬h0ÓÉÄÚ²¿×Ó³ÌÐòµ÷Óõõ½£¬ËüÓÖÓɼ¸¸öС¾ØÕó£ºc,a,b×éºÏµÃµ½¡£Ê×ÏȽ«h0¶¨ÒåΪȫ0¾ØÕó£¬È»ºóÔÙ¶ÔËüÄÚ²¿Æ¬¶Î¾ØÕó¸³Öµ£¬c,a,b°´ÕÕÒ»¶¨Ë³ÐòÅÅÁÐÔÚh0ÀïÃ棬µÃµ½h0µÄÖµºó½«Æä´«»ØÖ÷³ÌÐòÇóÌØÕ÷Öµ¡£´ó¼ÒÄÜ·ñ°ïæ¿´Ï£¬²»ÖªµÀh0¶¨ÒåΪ¸´Ðͺò»ºÃ°¡ ³ÌÐò¼ûÕâÀעÒ⣬Ç뽫ºóÃæµÄ.jpgɾ³ýºó¾ÍÊÇf90¸ñʽÁË program main use IMSL implicit none real,parameter::pi=3.1415927 complex,parameter::fi=(0.0,1.0) complex q real t,f,w    complex, allocatable::h0(:,: ) real, allocatable:: eigenvalue(: )   integer :: i,n read*,n   allocate(h0(4*n,4*n)) allocate(eigenvalue(n)) w=0.0 t=1.0 f=0.25 call hmatr(h0,n)   eigenvalue=eig(h0) do i=1,4*n    write(*,"('eigenvalue=',f10.7)" )eigenvalue(i)   end do stop contains subroutine hmatr(h0,n) integer n,i,j,n0 complex h0 complex,dimension(:,: )::c(4,4),a(4,4),b(4,4) h0=0.0 a=reshape((/0.0,0.0,0.0,t,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/),(/4,4/)) b=transpose(a) do n0=1,2*n-1,2     c=((/w,t*cmplx(exp(-fi*(-pi+2*pi/3*f*(1.5*n0+1/4)))),0,0,t*exp(fi*(-pi+2*pi/3*f*(1.5*n0+1/4))),w,t,0,0,t,w,t*exp(fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),0,0,t*exp(-fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),w/),(/4,4/))      j=2*n0-1       do j=1,4*n-3,4        h0(j:j+3,j:j+3)=c          if(j==1) then           h0(j:j+3,j+4:j+7)=b;           h0(j:j+3,4*n-3:4*n)=a;         else if(j>=5.and.j<4*n-3) then           h0(j:j+3,j+4:j+7)=b;           h0(j:j+3,j-4:j-1)=a;         else           h0(4*n-3:4*n,1:4)=b;           h0(j:j+3,j-4:j-1)=a;         end if       end do end do end subroutine hmatr end ³ö´íÐÅÏ¢£º 1£©Error: The shape matching rules of actual arguments and dummy arguments have been violated.   [H0] call hmatr(h0,n)   -----------^ 2£©E:\matlabxuexi\yyuu.f90(36) : Error: A constant or named constant is required in this context.     c=((/w,t*cmplx(exp(-fi*(-pi+2*pi/3*f*(1.5*n0+1/4)))),0,0,t*exp(fi*(-pi+2*pi/3*f*(1.5*n0+1/4))),w,t,0,0,t,w,t*exp(fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),0,0,t*exp(-fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),w/),(/4,4/)) ---------^ 3£©E:\matlabxuexi\yyuu.f90(36) : Error: An INTEGER or REAL data type is required in this context.     c=((/w,t*cmplx(exp(-fi*(-pi+2*pi/3*f*(1.5*n0+1/4)))),0,0,t*exp(fi*(-pi+2*pi/3*f*(1.5*n0+1/4))),w,t,0,0,t,w,t*exp(fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),0,0,t*exp(-fi*(-pi+2*pi/3*f*(1.5*(n0+1)+1/4))),w/),(/4,4/)) ---------^ 4£©E:\matlabxuexi\yyuu.f90(39) : Error: The syntax of this substring is invalid.   [H0]        h0(j:j+3,j:j+3)=c -------^ Error executing df.exe. yyuu.obj - 4 error(s), 0 warning(s) [ Last edited by nono2009 on 2010-8-3 at 09:59 ]

»Ø¸´´ËÂ¥