function tiaoxing(bochang,N,m,L,handles) c$R<j'7
global n ^h"F\vIpV
global point ]YwvwmZ
k=2*pi/bochang; (Ei} :6,}
a=m*1e+006; H.[t&VO
deltax=2*a/point;%取单位长度 %1%@L7wP>
xx=-a:deltax:a; M0"}>`1lJ
L=L*1e+006; ,L^ag&!4
nn=length(xx); qiyJ4^1
u=ones(1,nn); NC{8[*Kx5
jj=1:nn; {A8w~3F
for n=1:N %进行迭代 XQ$9E?|=
Mg=R**s1x%
for b=1:nn teg[l-R"7z
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... 6ragRS/'x
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% eLN[`hJ
)TyI~5>;
end qM:*!Aq0g
UpCkB}OhR1
ma=max(uu);%归一化 WOuk>
/
uu=uu./ma; Ln&CB!u
u=uu; !yf7y/qY
axes(handles.axes1); '7>Yrzq
plot(xx,abs(u)) hwd{^
axes(handles.axes2); (j884bu
plot(xx,angle(u)) 4W36VtQ@E
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 .Sjg
End