function tiaoxing(bochang,N,m,L,handles) 9 wSl,B-
global n -}_cO|kk
global point 7wiK.99
k=2*pi/bochang; %="~\1y
a=m*1e+006; %\X P:
deltax=2*a/point;%取单位长度 QIij>!c4
xx=-a:deltax:a; S_T{L
L=L*1e+006; TV1e
bH7q
nn=length(xx); W+d9cM=
u=ones(1,nn); d7W%zg\T
jj=1:nn; fuNl4BU
for n=1:N %进行迭代 uQ:ut(
G}=`VYK
for b=1:nn #$rf-E5g-K
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... Z7/vrME6
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% ~@[(N]=q
[^?13xMb
end >vD['XN,
<CNE>@-f
ma=max(uu);%归一化 ,_.@l+BM.
uu=uu./ma; %r M-"6Q
u=uu; 9,S,NvSq
axes(handles.axes1); pG,<_N@P
plot(xx,abs(u)) HCX!P4Hj
axes(handles.axes2); 3E<aiGU
plot(xx,angle(u)) tr}KPdE
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 lQh
E]m>+
End