新手求救:帮我看看这个程序的错误在那

litao814

新手求救:帮我看看这个程序的错误在那,谢谢<BR>function[N,G]=solution(ALAR,ALAS,FLS1);<BR>G=cell(1,20),A=cell(1,20),B=cell(1,20);<BR>X=exp(0.1.*ALAR.*log(10))-1;<BR>Y=exp(0.1.*ALAS.*log(10))-1;<BR>AC=sqrt(X./Y);<BR>AA=AC+sqrt(AC.*AC-1.0);<BR>AB=FLS1+sqrt(FLS1.*FLS1-1.0);<BR>BC=log(AA)./log(AB)+nan;<BR>if fix(BC)==BC;<BR>    N=BC;<BR>else;<BR>    N=BC+1;<BR>end;<BR>M=fix(N/2).*2;<BR>if  M==N;<BR>    N=N+1;<BR>else;<BR>    N=N<BR>end;<BR>G(0)=1.0;<BR>G(N+1)=1.0;<BR>C=tanh(ALAR./17.37);<BR>BETA=log(1.0./c);<BR>GAMA=sinh(BETA./2.0./N);<BR>for  J=1:N;<BR>    A(J)=sin((2.0.*J-1.0).*PI./(2.0.*N));<BR>    B(J)=GAMA.*GAMA+sin(J.*PI./N).* sin(J.*PI./N); <BR>end<BR>for  J=2:N;<BR>    G(J)=4.0.*A(J-1).*A(J)./(B(J-1).*G(J-1))<BR>end;

happy

<P>把ALAR,ALAS,FLS1参数给出来，好调试<BR>一句句看没工夫</P>