function tiaoxing(bochang,N,m,L,handles) @|%ICG c
global n >Eg .c
global point %MZDm&f>Kk
k=2*pi/bochang; ]},Q`n>$
a=m*1e+006; 5XO'OSdYq
deltax=2*a/point;%取单位长度 9R"N#w.U]
xx=-a:deltax:a; -Bv1}xf=6
L=L*1e+006; f?zK"
nn=length(xx); DeQ'U!?+N
u=ones(1,nn); <D;Q8
jj=1:nn; Prz+kPP
for n=1:N %进行迭代 SLOYlRGCi
iiB )/~!O
for b=1:nn )h_7 2
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... N%|Vzc
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% b(-t)5^}
LY:?OGh
end T-2p`b}hW
TxxB0
ma=max(uu);%归一化 j{Jc6U
uu=uu./ma; ieRBD6_
u=uu; ss@}Dt^
axes(handles.axes1); #s
R0*
plot(xx,abs(u)) s1 ^mk]
axes(handles.axes2); exEld
plot(xx,angle(u)) V3~a!k
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 qs1.@l("
End