function tiaoxing(bochang,N,m,L,handles) ~sNBklK
global n 34J*<B[Njo
global point CsjrQ-#9yn
k=2*pi/bochang; `~N jBtQ
a=m*1e+006; d@ ] N
deltax=2*a/point;%取单位长度 -\25&m! +
xx=-a:deltax:a; p &
Kfy~
L=L*1e+006; C4
-y%W"P
nn=length(xx); K C8
u=ones(1,nn); #[Rs &$vQm
jj=1:nn; s#X fu\CP
for n=1:N %进行迭代 _]L]_Bh
y-~_ W 6\
for b=1:nn w,OPM}) il
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... Ou+b ce
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% _SMi`ie#
<00=bZzX
end "1 UpoF'w
` ^;J<l
ma=max(uu);%归一化 #S[Y}-]T
uu=uu./ma; Th$xk9TK^@
u=uu; )&%Y{a#
axes(handles.axes1); E,Xl8rC
plot(xx,abs(u)) 6;(b-Dhi
axes(handles.axes2); i\x@ s>@x}
plot(xx,angle(u)) ``>z8t[ks
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 C(M ?$s `
End
