function tiaoxing(bochang,N,m,L,handles) - Cn x!g}
global n 9 +1}8"~
global point S<*h1}V3/
k=2*pi/bochang; <Wr n/%tL
a=m*1e+006; 9}-;OJe
deltax=2*a/point;%取单位长度 ]?p 9)d=%<
xx=-a:deltax:a; u uaoBf
L=L*1e+006; |)KOy~"
nn=length(xx); ,"e n7
u=ones(1,nn); 6(pa2
jj=1:nn; 3M^s
EaUI
for n=1:N %进行迭代 }BYs.$7
idzc4jR6BT
for b=1:nn x*'2% 3C~
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... E;1QD/E$
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% >DM^/ EAG{
uKA-<nM._c
end D( \c?X"
e^=b#!}-5:
ma=max(uu);%归一化 6 uS;H]nd<
uu=uu./ma; 62.)fCQ^
u=uu; hQb3 8W[
axes(handles.axes1); ,gO(zI-1
plot(xx,abs(u)) TI5<'
U)
axes(handles.axes2); ee2k..Tq#
plot(xx,angle(u)) V$Zl]f$S
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 qtR/ K=^i
End
