function tiaoxing(bochang,N,m,L,handles) )=5 ng-
global n Qj,]N@7
global point *RXbc~
H
k=2*pi/bochang; LJ*q 1
;<E
a=m*1e+006; X}tVmO?
deltax=2*a/point;%取单位长度 {7 $c8 i
xx=-a:deltax:a; "SzdDY6
L=L*1e+006; GqWB{$J;"
nn=length(xx); f J/e(t
u=ones(1,nn); Q,p}:e
jj=1:nn; '`eO\huf
for n=1:N %进行迭代 jqv- D
eln&] d;
for b=1:nn ]3VI|f$$
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... 0o+6Q8q
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% }D7I3]2>
#>%X_o-o23
end '@p['# \uI
& Rx {.9
ma=max(uu);%归一化 'f[T&o&L/
uu=uu./ma; q0$
!y!~
u=uu; LM}Ib.
axes(handles.axes1); sA'6ty
plot(xx,abs(u)) )+}]+xRWGj
axes(handles.axes2); T(e!_VY|m
plot(xx,angle(u)) c}y [[EX
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 I3 ,= 0z
End
