function tiaoxing(bochang,N,m,L,handles) ^&rb I,D
global n E0DE FB
global point : }IS=A
k=2*pi/bochang; gq@."wHU
a=m*1e+006; 6~/H#8Kdn
deltax=2*a/point;%取单位长度 U; q)01
xx=-a:deltax:a; X*yl% V
L=L*1e+006; #dfW 1@m
nn=length(xx); 0?h .X= G
u=ones(1,nn); V7@xr
M
jj=1:nn; T+ t-0k
for n=1:N %进行迭代 vZDQ@\HrC
&^7)yS+C
for b=1:nn 4#YklVm
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... 5k(#kyP
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% t3XMQ']
&sRJ'oc
end l&A`
yEny2q}
ma=max(uu);%归一化 ~GeYB6F
uu=uu./ma; D{l((t3 =T
u=uu; z, 7^dlT
axes(handles.axes1); %dU }GYL_
plot(xx,abs(u)) p{J_d,JH
axes(handles.axes2); ZD{srEa/a
plot(xx,angle(u)) !T{g& f
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 v8IL[g6"
End
