function tiaoxing(bochang,N,m,L,handles) h2}am:%mC
global n ]]V^:"ne
global point #6FaIq92V
k=2*pi/bochang; z`9l<Q/
a=m*1e+006; 59qnEIi
deltax=2*a/point;%取单位长度 /2:Q6J
xx=-a:deltax:a;
,(hY%M&\
L=L*1e+006; u-/3(dKt
nn=length(xx); er\:U0fr#@
u=ones(1,nn); `]yKM0 Z
jj=1:nn; 2TCRS#z
for n=1:N %进行迭代 xucIjPi]
H5^Y->
for b=1:nn v =bv@c
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... uk{J@&F
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% KH)pJG|NY
zuj;T,R;
end }moz9a
W,'3D~g8
ma=max(uu);%归一化 /kb$p8!C".
uu=uu./ma; ~g96o81V
u=uu; +wjlAqMQ
axes(handles.axes1); 1'OD3~[R
plot(xx,abs(u)) h&'J+b
axes(handles.axes2); Dpp@*xX>
plot(xx,angle(u)) I9s$bRbT
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 9e76pP(
End