function tiaoxing(bochang,N,m,L,handles) cix36MR_
global n |xeE3,8
global point Sr"/-
k=2*pi/bochang; GS4_jvD-
a=m*1e+006; S=PJhAF
deltax=2*a/point;%取单位长度 ="3,}qR
xx=-a:deltax:a; ^yJ:+m;6K
L=L*1e+006; sOz sY7z3Z
nn=length(xx); n`.#59-Hx
u=ones(1,nn); @[v8}D
jj=1:nn; OpQ8\[X+
for n=1:N %进行迭代 %t[K36,p
{(Fe7,.S3
for b=1:nn vYdlSe=6G
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... |!}wF}iLc)
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% :1MMa6
)8'jxiGs
end OD|1c6+X
" '[hr$h3
ma=max(uu);%归一化 ~K`1
uu=uu./ma; O,irpQ
u=uu; qd8pF!u|#
axes(handles.axes1); .[(P
plot(xx,abs(u)) q% EC
axes(handles.axes2);
W_}/ O'l{
plot(xx,angle(u)) /[?Jylj
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 }N]!0Ka
End