function tiaoxing(bochang,N,m,L,handles) 0Y/k /)Ul]
global n O\^D
6\ v
global point LgA > ,.
k=2*pi/bochang; dlc'=M
a=m*1e+006; D ?r% Y
deltax=2*a/point;%取单位长度 QykHB
k
xx=-a:deltax:a; ~{!,Z nO*
L=L*1e+006; n2&M?MGX
nn=length(xx); BmCBC,j<v>
u=ones(1,nn); mH!\]fmR~
jj=1:nn; ~JohcU}d
for n=1:N %进行迭代 a >Xq
I7jIA>ZZi
for b=1:nn qF4DX$$<
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... kk+:y{0V
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% *+<H4.W
H
Hv|(V3-
end l6HT}x7OiH
T$s )aM
ma=max(uu);%归一化 n16TQe"8
uu=uu./ma; i|G /x
u=uu; jx8hh}C
axes(handles.axes1); \`WAG>'l5
plot(xx,abs(u)) 1DM$FG_Z-
axes(handles.axes2); [,q^\T
plot(xx,angle(u)) q??N,
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 FSS~E [(DL
End
