function tiaoxing(bochang,N,m,L,handles) 9H/>M4RT
global n R7/S SuG6\
global point `R$i|,9)
k=2*pi/bochang; a(|6)w-
a=m*1e+006; ]:Gy]qkO
deltax=2*a/point;%取单位长度 %+H _V1F
xx=-a:deltax:a; A<ds+0
L=L*1e+006; g#=<;X2
nn=length(xx); KGYbPty}
u=ones(1,nn); =3.dgtH
jj=1:nn; GLn=*Dh#
for n=1:N %进行迭代 ' @RF
Y
{^*y
for b=1:nn E ?Mgbd3
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... bGCC?}\
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% 3'Y-~^ml|
{+WBi(=W
end -67Z!N
=I`S7oF
ma=max(uu);%归一化 |n/;x$Cb
uu=uu./ma; 8f9wUPr
u=uu; } Bf@69
axes(handles.axes1); (dq_,LI
plot(xx,abs(u)) TP
rq:"K
axes(handles.axes2); ,*J@ic7"
plot(xx,angle(u)) m ol,iM*l
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 Nvgi&iBh8
End