function tiaoxing(bochang,N,m,L,handles) v&MU=Tcqi
global n .%0ne:5
global point <V_7|)'/A
k=2*pi/bochang; RwTzz]
M
a=m*1e+006; } IlP:
deltax=2*a/point;%取单位长度 Z# Lx_*p]Q
xx=-a:deltax:a; J%dJw}
L=L*1e+006; S9Yt 1qb
nn=length(xx); px9>:t[P
u=ones(1,nn); j:1uP^.
jj=1:nn; | D.C!/69
for n=1:N %进行迭代 n!N\zx8
/4joC9\AB
for b=1:nn 8kW /DcLE
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... CM~MoV[k7e
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% "6I[4U"@
s=EiH
end hE!7RM+Y
GF--riyfB
ma=max(uu);%归一化 iG[?
]]
uu=uu./ma; F1A1@{8bN
u=uu; ->yeJTsE9
axes(handles.axes1); r7)@M%A
plot(xx,abs(u)) ZcZ;$*
axes(handles.axes2); sS,#0Qt.
plot(xx,angle(u)) GzdgL"M[
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 {:Vf0Mhb
End