function tiaoxing(bochang,N,m,L,handles) 7vO3+lT/Y;
global n 9jW/"
global point MG6taOO!
k=2*pi/bochang; hDBVL"
a=m*1e+006; P(AcDG6K
deltax=2*a/point;%取单位长度 whFaL}2C
xx=-a:deltax:a; 0}v_usP
L=L*1e+006; _voU^-
nn=length(xx); f/+UD-@%m
u=ones(1,nn); zv/owK
jj=1:nn; o^HzE;L}
for n=1:N %进行迭代 R8ZI}C1
Tk1U
for b=1:nn mz)Z
=`hy
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... $}_a`~u
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% fP 4
X|]&K
end gLp7<gx6
X'[93
C|K
ma=max(uu);%归一化 &IYSoA"Nz
uu=uu./ma; o}QtKf)W
u=uu; w
K)/m`{g
axes(handles.axes1); oMdqg4HUF
plot(xx,abs(u)) QxUsdF?p
axes(handles.axes2); * F[;D7sZ~
plot(xx,angle(u)) i_<