function tiaoxing(bochang,N,m,L,handles) wTbIS~!gF
global n '=Y ~Ir+
global point %--5bwZi
k=2*pi/bochang; / hdl
a=m*1e+006; :2gO )
'cD
deltax=2*a/point;%取单位长度 -7VV5 W
xx=-a:deltax:a; ?|YQtY
L=L*1e+006; Q</HF pE
nn=length(xx); leJ d) {
u=ones(1,nn); .JKaC>oX
jj=1:nn; BQf+1 Ly&
for n=1:N %进行迭代 @CUDD{1o
S $/3K q
for b=1:nn T )]|o+G
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... j" .6
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% hd% F7D5
XYTcG;_z
end D|#(zjl@
:e v c
ma=max(uu);%归一化 - $/{V&?t
uu=uu./ma; V7i`vo3Cc
u=uu; ?+L6o C.;
axes(handles.axes1); xF+x I6
plot(xx,abs(u)) ;H|M)z#[Z
axes(handles.axes2); `sC8ro@Fm
plot(xx,angle(u)) l6Bd<tSH
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 >;?97'M
End
