function tiaoxing(bochang,N,m,L,handles) LgxsO:mi
global n Qv
WvS9]
global point JfsvK2I
k=2*pi/bochang; $
~Ks!8'P
a=m*1e+006; 6r.#/' "
deltax=2*a/point;%取单位长度 2!_DkE
xx=-a:deltax:a; ^WkqRs
L=L*1e+006; tc0(G~.N
nn=length(xx); F};T<#
u=ones(1,nn); m-2!r*(zt
jj=1:nn; U4NH9-U'
for n=1:N %进行迭代 }nRTw2-z
/|p\l"
for b=1:nn Jsa]RA
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... YB1Jv[
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% 9~J#> C0}
xOx=Z\ c
end d=0{vsrB
Z<iK(?@O
ma=max(uu);%归一化 " ?Ux\)*
uu=uu./ma; rI OKCL?
u=uu; >a;LBQ0
axes(handles.axes1); t)5.m}
plot(xx,abs(u)) (l -l
Y
axes(handles.axes2); tlM >=s'T
plot(xx,angle(u)) wBJ|%mc3TA
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 :*ZijN*{)$
End