function tiaoxing(bochang,N,m,L,handles) u4ZOHy_O^
global n Em(Okr, 0
global point ff
6x4t
k=2*pi/bochang; SZ~lCdWad
a=m*1e+006; ~#7uNH2
deltax=2*a/point;%取单位长度 ^;wz+u4^l
xx=-a:deltax:a; ': ]Hj8t_
L=L*1e+006; huvn _
nn=length(xx); m+xub*/
u=ones(1,nn); 2IXtIE
jj=1:nn; B(Er/\-@U
for n=1:N %进行迭代 #dQFs]:F
@ ?bY,
for b=1:nn ]rji]4s
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... .z^O y_S{
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% lLLPvW[Q
g1@rY0O
end
pRA%07?W
Y I2x*t!
ma=max(uu);%归一化 RSC Q`.
uu=uu./ma; FIx|4[&>S
u=uu; \cK# /;a#
axes(handles.axes1); ~Gl5O`w(
plot(xx,abs(u)) u-*z#e_L0
axes(handles.axes2); rA}mp]
plot(xx,angle(u)) /z :1nq
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 ]d.e(yCuE
End
