function tiaoxing(bochang,N,m,L,handles) tz&oe
global n - kVt_
global point L`Lro:E?kL
k=2*pi/bochang; ny:4L{)
a=m*1e+006; O%.c%)4Xo
deltax=2*a/point;%取单位长度 D`2c61jyc
xx=-a:deltax:a; pPNU0]/
L=L*1e+006; ivKhzU+
nn=length(xx); `$*cW1
u=ones(1,nn); wua`e <"
jj=1:nn; CnT]u U
for n=1:N %进行迭代 K(+ ~#$|-~
Ne)H*D T
for b=1:nn u"*@k^}(
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... ep- ~;?
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% 9b8ZOk'9 _
p pjS|l*`
end 0Y 8 Si^T
Vnu*+
ma=max(uu);%归一化 M7|k"iz v
uu=uu./ma; o+o'!)
u=uu; _%` )cOr
axes(handles.axes1); H]7MN Y
plot(xx,abs(u)) u!FX 0Ip
axes(handles.axes2); ~9N n8g6
plot(xx,angle(u)) f,'gQ5\ X3
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 IXaF(2>
End
