function tiaoxing(bochang,N,m,L,handles) p \A ^kX^5
global n PeGA+0bm
global point wi^zXcVj
k=2*pi/bochang; ?$i`K|
a=m*1e+006; uCO-f<b
deltax=2*a/point;%取单位长度 W+36"?*k3
xx=-a:deltax:a; Nd'+s>d0
L=L*1e+006; Tj7OV}:
nn=length(xx); !`"@!
u=ones(1,nn); O32p8AxEz
jj=1:nn; GEj/Z};;[b
for n=1:N %进行迭代 4scY8(1
JJ`RF
for b=1:nn &6<>hqR^
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... t0h@i`
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% 1Q"w)Ta
d.Wq@(ZoA
end ^kJ(bBY
#$7d1bx
ma=max(uu);%归一化 0NrUB
uu=uu./ma; 'X_8j` ]#
u=uu; *"Yz"PK
axes(handles.axes1); {:BAh5e|
plot(xx,abs(u)) XgL-t~_
axes(handles.axes2); Z BjyQ4h
plot(xx,angle(u)) V'?bZcRr~
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 |s[kY
End