现代数字信号处理仿真作业.
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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 算法%%%%%%%%%%%%%%%% % S=VR(:,IX(M:-1:M-K+1; % S1=S(1:M-1,:; %
S2=S(2:M,:; % fai=S1\S2; % % [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*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 % fy=exp(j*pi*sin(n3/180*pi; % tt=[(fy.^(n1'.' zeros(1,M-P]; %
Tfy=toeplitz(cc,tt; % GfTheta=1./(P^2*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 归一化功率谱 /dB -10 -15 -20 -25 -30 -100 -80 -60 -40 -20 0 20 空间角度 40 60 80 100