function tiaoxing(bochang,N,m,L,handles) d} ]jw4
global n LD)P.
f
global point +,ojlTVlt
k=2*pi/bochang; x ]jJ
a=m*1e+006; S9S%7pE
deltax=2*a/point;%取单位长度 ^
Nm!b
xx=-a:deltax:a; G8!* &vR/
L=L*1e+006; ,hXhcfFl
nn=length(xx); l)(
3]
u=ones(1,nn); 5Qgh \ 4
jj=1:nn; y j C@
for n=1:N %进行迭代 RL3G7 ;X
qEV>$>}
for b=1:nn 5IiZnG u
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... MFm2p?zPm
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% TbA=bkj[4
@t?uhT*Z=
end \L{V|}" X
; )J\k2
ma=max(uu);%归一化 /%w3(e
uu=uu./ma; n |f Huv
u=uu; *.F4?i2D
axes(handles.axes1); * b+ ~@o
plot(xx,abs(u)) M[7$cfp-Y~
axes(handles.axes2); `E 2HQA@
plot(xx,angle(u)) j z72~+)T
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 +LsACSB
End
