function tiaoxing(bochang,N,m,L,handles) 4;~xRg;u&*
global n DOtz
global point /|U;_F Pmc
k=2*pi/bochang; )Mq4p'*A[
a=m*1e+006; "n3n-Y#'
deltax=2*a/point;%取单位长度 9u=]D> kb
xx=-a:deltax:a; Z]LP18m9kl
L=L*1e+006; Z5rL.a&
nn=length(xx); RoHX0
u=ones(1,nn); w\
'5lk,"
jj=1:nn; YH^U"\}i
for n=1:N %进行迭代 b/.EA'/
w@WtW8
p^
for b=1:nn -d!84_d9
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... 85"Szc-#
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% &:d`Pik6
n=rmf*,?
end um PN=0u6
^&F.T-( A
ma=max(uu);%归一化 /!&eP3^
uu=uu./ma; -;Mh|!yg
u=uu; h{yh}04P1
axes(handles.axes1); D'g@B.fXd
plot(xx,abs(u))
*W |
axes(handles.axes2); 4%v-)HGh
plot(xx,angle(u)) r8y,$Mv<)0
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 4t*VI<=<[
End