function tiaoxing(bochang,N,m,L,handles) H-pf8
global n f%EHzm/V
global point zn-=mk;W
k=2*pi/bochang; -OJ <Lf+"=
a=m*1e+006; &fB=&jc*j
deltax=2*a/point;%取单位长度 `C: 7N=9
xx=-a:deltax:a; YtvDayR>
L=L*1e+006; X:s~w#>R
nn=length(xx); OD~Q|I(j
u=ones(1,nn); Yd4X*Ua
jj=1:nn; 0+iRgnd9?
for n=1:N %进行迭代 "mBM<rEn*
d;dT4vx$[M
for b=1:nn wuXQa
wo
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... ]^"Lc~w8&
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% vNjc
-crMO57/
end o;6~pw%
y%vAEQ2j=
ma=max(uu);%归一化 |v7Je?yh
uu=uu./ma; 7u;N/@
u=uu; G AY?F
axes(handles.axes1); UY9*)pEE
plot(xx,abs(u)) >g !Z|ju
axes(handles.axes2); =aB+|E
plot(xx,angle(u)) m^XO77"
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 1P"7.{
End