function tiaoxing(bochang,N,m,L,handles) 9Y~A2C
global n >Efv?8$E\
global point }9jy)gF*e
k=2*pi/bochang; Ig]Gg/1G
a=m*1e+006; 'Exj|Y&
deltax=2*a/point;%取单位长度 Q[S""P.Z|
xx=-a:deltax:a; 1S <V,9(
L=L*1e+006; T0v;8Ee
nn=length(xx); JhIgqW2
u=ones(1,nn); $TWt[
jj=1:nn; F.K7w
for n=1:N %进行迭代 1)vdM(y3j
GYZzWN}U
for b=1:nn ,qyH B2v
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... q*,];j/>k
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% yX?& K}JI
UV$v:>K#
end o`7Bvh2
?$v#;n?@I
ma=max(uu);%归一化 Ar N *9
uu=uu./ma; K)N7Y=C3
u=uu; | x/,
axes(handles.axes1); 9PIm/10pP^
plot(xx,abs(u)) A2!7a}*1(
axes(handles.axes2); @u#Tx%
plot(xx,angle(u)) y#o ,Vg*V
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 j!xt&t4D
End