function tiaoxing(bochang,N,m,L,handles) {.)D)8`<d
global n FkkB#Jk4
global point 51usiOq
k=2*pi/bochang; 0 pz
X!f1~
a=m*1e+006; >FhBl\oIi
deltax=2*a/point;%取单位长度 t, YAk
?}
xx=-a:deltax:a; ,~w )@.
L=L*1e+006; T!E LH!
nn=length(xx); 40ZB;j$l
u=ones(1,nn); pr?(5{BL
jj=1:nn; Q)8t;Kx
for n=1:N %进行迭代 ( \
%y )
nwVtfsb
for b=1:nn PUArKBYM-
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... $cCB%}
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% yh!vl&8M
Fb-TCq1y#
end } 4^UVdz
iDN,}:<V
ma=max(uu);%归一化 W5{e.eI}|
uu=uu./ma; 1j(, VW
u=uu; Wn5]2D\vkT
axes(handles.axes1); ^5F/=TtE G
plot(xx,abs(u)) 548BM^^"r
axes(handles.axes2); 85]UrwlA4
plot(xx,angle(u)) mer{Jy s
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 2
{0VyLx
End
