function tiaoxing(bochang,N,m,L,handles) ^@Qi&g`lr?
global n PEBFN
global point &'7"i~pC
k=2*pi/bochang; }iSakq'
a=m*1e+006; g&4~nEp
deltax=2*a/point;%取单位长度 UNQRtR/
xx=-a:deltax:a; JQ_gM._3
L=L*1e+006; ,0Zn hS)kq
nn=length(xx); TJXraQK-=
u=ones(1,nn); nr>{ uTa
jj=1:nn; _g 4/%
for n=1:N %进行迭代 eFSC^
u7"VeTz
for b=1:nn mY`]33??v
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... {(OIu]:
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% h)_Gxe"x
?Y{^un
end Iv(Qa6(
Iz)hz9k
ma=max(uu);%归一化 D4q>R;
uu=uu./ma; bW9"0=j[{
u=uu; )M'UASB;8
axes(handles.axes1); kV ,G,wo
plot(xx,abs(u)) m*!f%}T
axes(handles.axes2); ?{: D,{+
plot(xx,angle(u)) _;:rkC fj
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 3u*hTT
End