| yzktst |
2016-06-10 09:43 |
function tiaoxing(bochang,N,m,L,handles) WvArppANo
global n PC L
;Z
global point RR:%"4 M
k=2*pi/bochang; uxDM
#
a=m*1e+006; EFx>Hu/ [G
deltax=2*a/point;%取单位长度 QnP3 U
xx=-a:deltax:a; ^ Lc\{,m
L=L*1e+006; KiI+ V;o
nn=length(xx); J;^ PM:6
u=ones(1,nn); P%Vq#5
jj=1:nn; (p12=EB<
for n=1:N %进行迭代 :] U\{;q2
0,m] W )
for b=1:nn ^dk$6%0
uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*... h/|p`MP\1
sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;% "9 c=kqkX
,m]5j _< }
end 1,)
yEeHjU
T^`; wD
ma=max(uu);%归一化 * {+ {h;p
uu=uu./ma; R{bG`C8.d
u=uu; -3) jUzD
axes(handles.axes1); 1,;zX^
plot(xx,abs(u)) \$.{*f
axes(handles.axes2); 7p]Izx8][
plot(xx,angle(u)) !dGu0wE
set(handles.edit10,'string',num2str(n));%迭代次数显示到面板 u<cnz% @
End
|
|