function tiaoxing(bochang,N,m,L,handles) (yQ
5`
global n ^[TV;9I*
global point T+"f]v
k=2*pi/bochang; b|^I<7
a=m*1e+006; ~7q uTp)
deltax=2*a/point;%取单位长度 E]^5I3=O
xx=-a:deltax:a; ]kktoP|D
L=L*1e+006; #nyv+x;
nn=length(xx); #i QX6WF
u=ones(1,nn); &UAe!{E0
jj=1:nn; "YFls#4H-
for n=1:N %进行迭代 toa-Wa{
<*'cf2Q$Av
for b=1:nn (g/7yO(s
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... bggusK<
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% p5c8YfM
Y{Ap80'\6
end |oKu=/[K
"i'bTVs
ma=max(uu);%归一化 }4jC_ZAupt
uu=uu./ma; ^Uw[x\%#gD
u=uu; y93k_iq$S
axes(handles.axes1); cErI%v}v0
plot(xx,abs(u)) SVP:D3)
axes(handles.axes2); [_DPxM=V
plot(xx,angle(u)) sB}]yw
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 '|K.k6
End