function tiaoxing(bochang,N,m,L,handles) AN-qcp6=o
global n G.CkceWRn
global point d\% |!ix
k=2*pi/bochang; X?PcEAi;w
a=m*1e+006; Li[ :L
deltax=2*a/point;%取单位长度 E0 ~\ A;
xx=-a:deltax:a; Bf
{h\>q
L=L*1e+006; {G _ :#cep
nn=length(xx); J9+<9g4-t
u=ones(1,nn); ~?l>QP|o
jj=1:nn; WCf?_\cG
for n=1:N %进行迭代 [T|_J$
;
565UxG
}
for b=1:nn OjVI4@E;Xe
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... ma__LWKM,
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% v#yeiE4
9|y?jb5im
end NJ6*
7Cd
L[M`LZpJo
ma=max(uu);%归一化 9 b?Nlk8d
uu=uu./ma; JL[xrK0
u=uu; <#c/uIN
axes(handles.axes1); BO6u<cu"-
plot(xx,abs(u)) [FhFeW>
axes(handles.axes2); u">KE6um
plot(xx,angle(u)) ~!bA<q
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 )zUbMzF
End