function tiaoxing(bochang,N,m,L,handles) ](^$5Am
global n uItzFX*
global point 7EI5w37
k=2*pi/bochang; ^`BiA'gPPC
a=m*1e+006; q:}Q5gzZ
deltax=2*a/point;%取单位长度 eootHK
xx=-a:deltax:a; !06
!`LT
L=L*1e+006; 3e)W_P*0?
nn=length(xx); CrvL[6i
u=ones(1,nn); !+<OED=qe
jj=1:nn; eUY/H1
for n=1:N %进行迭代 %S'gDCwq
qdss(LZ
for b=1:nn (pDu
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... d*}dM"
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% ||xiKg
nBN+.RB:(
end Lo<-;;vQ
eNQQ`ll@m
ma=max(uu);%归一化 ~!t# M2Sk
uu=uu./ma; E4CyW
u=uu; )U2cS\k'7n
axes(handles.axes1); \/wbk`2
plot(xx,abs(u)) 6k4ZzQ}
axes(handles.axes2); J*!_kg)>J
plot(xx,angle(u)) &d%0[Ui`
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 ,$H[DX
End