function tiaoxing(bochang,N,m,L,handles) xO` O$ie
global n 94xRKQ}
global point aMycvYzH
k=2*pi/bochang; (cEjC`]
a=m*1e+006; >ay%
!X@3"
deltax=2*a/point;%取单位长度 !Won<:.[0
xx=-a:deltax:a; 0 P-eC|0
L=L*1e+006; ]W>kbHImz
nn=length(xx); Ju 0
u=ones(1,nn); (}NKW
jj=1:nn; CYWL@<p,
for n=1:N %进行迭代 s.uV,E*wu
c2fbqM~
for b=1:nn bQu1L>c,Uw
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... &^!vi2$5}
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% nq"U`z@R
A5LTgGzaW
end R#i{eE*WF
W|aFEY
ma=max(uu);%归一化 n%Gk
{h5
uu=uu./ma; Y<drRK!
u=uu; ?#rejA:
axes(handles.axes1); 3BWYSJ|
plot(xx,abs(u)) PQFr4EY?i
axes(handles.axes2); .gK>O2hI
plot(xx,angle(u)) iw==q:$
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 i[d@qp!H=
End