function tiaoxing(bochang,N,m,L,handles) 1n6%EC|X
global n ' qT\I8%
global point %<O~eXY
k=2*pi/bochang; ]M "U 'Z
a=m*1e+006; 5mL4Zq"
deltax=2*a/point;%取单位长度 (*vBpJyz%
xx=-a:deltax:a; :T@} CJ
L=L*1e+006; l LBzY`j
nn=length(xx); 3;uLBuZOCN
u=ones(1,nn); z wwJyy%/
jj=1:nn; # Rs5W
for n=1:N %进行迭代 M djxTr^
muK.x7zyl
for b=1:nn y=
8SD7P'
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... &Wdi
5T8
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% &oZU=CN
h^,L) E
end o7PS1qcya<
?djH!
ma=max(uu);%归一化 f}.t
uu=uu./ma; heWQPM|s
u=uu; f\RTO63|O
axes(handles.axes1); d mTZEO
plot(xx,abs(u)) ?-0, x|ul
axes(handles.axes2); 96; gzG@1!
plot(xx,angle(u)) Cd6th
F)
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 @S5HMJ2=
End