function tiaoxing(bochang,N,m,L,handles) 7'd_]e-.
global n ybcQ,e
global point ~%m-}Sxc
k=2*pi/bochang; |{<g-)
a=m*1e+006; e$fxC-sZ
deltax=2*a/point;%取单位长度 Cj,fP[p#7
xx=-a:deltax:a; U5%EQc-"P
L=L*1e+006; eLM_?9AZ!R
nn=length(xx); P*Uu)mG)G
u=ones(1,nn); Jcy
jj=1:nn; ?E=&LAI#
for n=1:N %进行迭代 ~"iCx+pr
JEGcZeq)
for b=1:nn C
3b
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... +3!um
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% JO1KkIV
U_y)p Cd
end |U>BXX P
`b+f^6SJn
ma=max(uu);%归一化 \ 2*<Pq
uu=uu./ma; 8J7xs6@
u=uu; W:ih#YW_F
axes(handles.axes1); H'P1EZtq
plot(xx,abs(u)) S2@[F\|r
axes(handles.axes2); Zm4IN3FGLv
plot(xx,angle(u)) bo4 :|Z
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 OXK?R\ E+
End