function tiaoxing(bochang,N,m,L,handles) K 6 yD64
global n <iBn-EG l>
global point H.;2o(vD
k=2*pi/bochang; p" 6ydXn%
a=m*1e+006; 'h@&rr@5
deltax=2*a/point;%取单位长度 3 Q~0b+k
xx=-a:deltax:a; 2tg 07
L=L*1e+006; :r2d%:h%2
nn=length(xx); ^ E_chx-e}
u=ones(1,nn); _f~$iY
jj=1:nn; JAM]neKiX
for n=1:N %进行迭代 *&tTiv{^
3mHP=)
for b=1:nn Vry*=X &Q
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... njaKU?6%d2
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% XSCcumde!
Z^BZH/I?
end P 'o]#Az
HRyhq ;C
ma=max(uu);%归一化 Z&4L///
uu=uu./ma; =4s x(<
u=uu; |S~$IFN4
axes(handles.axes1); 3ZN\F
plot(xx,abs(u)) d+vAm3.Dg
axes(handles.axes2); K%W;-W*'
plot(xx,angle(u)) )H`V\ H[0P
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 \= P(?!v
End
