function tiaoxing(bochang,N,m,L,handles) GY4yZa
global n 5Ut0I]h|z
global point U}
g%`<
k=2*pi/bochang; ~PV>3c3l=
a=m*1e+006; 5=
F-^
deltax=2*a/point;%取单位长度 CZ0 {*K:
xx=-a:deltax:a; :<jf}[w!
L=L*1e+006; g_}@/5?y
nn=length(xx); 6)?TWr'K e
u=ones(1,nn); 1%+0OmV&
jj=1:nn; KYeA=
for n=1:N %进行迭代 6HEl1FK{@
JlH|=nIaj6
for b=1:nn X~R
qv5@-
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... E O}(MXS
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% uIba{9tM"P
5 )C~L]
end G{s ,Y^
w%$n)7<*
ma=max(uu);%归一化 >U}~Hv]
uu=uu./ma; IAtZ-cM<
u=uu; aNwDMd^+
axes(handles.axes1); 8cd,SQ}y
plot(xx,abs(u)) Kp99y
axes(handles.axes2); +\E\&^ZQ
plot(xx,angle(u)) xr=f9?%R
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 A^m hPBT_
End