1、 xlabel('空間角度' % % % %RootMUSIC 算法% % syms z % pz=z.(0:M-1' % pz1=(z(-1.(0:M-1; % fz=z(M-1*pz1*G*G'*pz; % a=sym2poly(fz; % r=roots(a; % r1=abs(r; % for i=1:2*K % Y,I(i=min(abs(r1-1; % r1(I(i=inf; % end % for i=1:2*K % theta_esti1(i=asin(-angle(r(I(i/pi*180/pi; % end % % % % %ESPRIT 算

2、法% % S=VR(:,IX(M:-1:M-K+1; % S1=S(1:M-1,:; % S2=S(2:M,:; % fai=S1S2; % % U_fai,V_fai=eig(fai; % for i=1:K % theta_esti2(i=asin(-angle(V_fai(i,i/pi*180/pi; % end % % % %MVDR 算法% % P=; % for n=-pi/2:pi/180:pi/2 % a=exp(-j*0:M-1'*pi*sin(n; % P=P,1/(a'*inv(R*a; % end % P=abs(P/max(abs(P; % P=10*

3、log10(P; % figure(2; % plot(t,P; % title('MVDR 算法' % ylabel('歸一化功率譜/dB' % xlabel('空間角度' % % % %F-SAPES 算法% % P=6; % L=M+1-P; % Rf=zeros(L,L; % for i=1:P % Rf=Rf+X(i:i+L-1,:*X(i:i+L-1,:'/N; % end % Rf=Rf/P; % n1=0:P-1; % n2=0:L-1; % cc=1 zeros(1,L-1; % for n3=-90:.5:90 % f

4、y=exp(j*pi*sin(n3/180*pi; % tt=(fy.(n1'.' zeros(1,M-P; % Tfy=toeplitz(cc,tt; % GfTheta=1./(P2*Tfy*R*Tfy' % Qf=Rf-Rf-GfTheta; % aTheta=fy.(-n2' % Wof=(inv(Qf*aTheta./(aTheta'*inv(Qf*aTheta; % sigma2Theta(n3+90/.5+1=Wof'*GfTheta*Wof; % end % sigma2Theta=abs(sigma2Theta/max(abs(sigma2Theta; % sigma2Theta=10*log10(sigma2Theta; % figure(3; % t1=-90:.5:90; % plot(t1,sigma2Theta; % title('F-SAPES 算法' % ylabel('歸一化功率譜/dB' % xlabel('空間角度' MVDR算 法 0 -5 -10 歸 一 化 功 率 譜 /dB -15 -20 -25 -30 -35 -100 -80 -60 -40 -20 0 20 空間角度 40 60 80 100 F-SAPES 算 法 0 -5 歸