function tiaoxing(bochang,N,m,L,handles) JK<[]>O
global n 3#t#NW*e
global point k*N!U[]
k=2*pi/bochang; Dch\k<Te
a=m*1e+006; uOzoE_i
deltax=2*a/point;%取单位长度 Qf^c}!I
xx=-a:deltax:a; p% mHxYP
L=L*1e+006; p=nbsS~":
nn=length(xx); K:'^f? P
u=ones(1,nn); k5aB|xo
jj=1:nn; O$6&4p*F.
for n=1:N %进行迭代 vq(#Ih2
a_P|KRl
for b=1:nn T}^3 Re`i
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... &j1-Ouy
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% *
rlVE
v%T'!(0j/
end :<zIWje
2)\->$Q(H
ma=max(uu);%归一化 nX3?7"v
uu=uu./ma; ;G8H'gM07
u=uu; VrPsy) J68
axes(handles.axes1); =AzPAN#e
plot(xx,abs(u)) _*e_?]G-
axes(handles.axes2); -V{"Lzrfug
plot(xx,angle(u)) >Vt2@Ee
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 tIDN~[1
End