# 非平稳噪声下的MUSIC.zip

• 小小放牛娃
了解作者
• matlab
开发工具
• 1.3KB
文件大小
• zip
文件格式
• 0
收藏次数
• 5 积分
下载积分
• 0
下载次数
• 2021-12-29 21:04
上传日期

• 非平稳噪声下的MUSIC
• ffmu1.m
2.2KB

%MUSIC ALGORITHM %DOA Estimation by UCA clear all; close all; clc l=0.2; % length of wave d=0.5*l; c=1500;% velocity of wave f0=7500; % center frequency N=500; % Length of Signal M=10; % number of sensor P=2; % Number of Signals f1=7540; % signal frequency f2=7550; %signal frequency fs=15000; % sample frequency t=0:1/fs:(N-1)/fs; f=[f1,f2]'; jiao1=10; jiao2=15; jiao11=pi*jiao1/180; jiao22=pi*jiao2/180; jiao=[jiao11 jiao22] SNR=[15,15];% signal and noise of ratio for k0=1:2; SN(k0)=sqrt(2*10^(SNR(k0)/10)); end for i=1:P for q=1:M A(i,q)=exp(-1j*(q-1)*2*pi*d*sin(jiao(i))/l); end end S=exp(1j*2*pi*f*t);% signal for k1=1:2 S1(k1,:)=SN(k1)*S(k1,:); end dd=[1,1,1,1,3,1,4,7,1,1]; Wn=(randn(M,N)+1j*randn(M,N)); for k=1:length(dd) Wn(k,:)=dd(k)*Wn(k,:); end x=A'*S1+Wn; %加噪声 R=x*x'; %------------------------ R1=R(P+1:2*P,1:P); R2=R(2*P+1:M,1:P); R3=R(P+1:2*P,2*P+1:M); R4=R(2*P+1:M,2*P+1:M); Q3=R4-R2*inv(R1)*R3; T=[eye(3*P,3*P) zeros(3*P,M-3*P)]; Rt=T*R*T'; R11=Rt(1:P,1:P); R22=Rt(P+1:2*P,P+1:2*P); R33=Rt(2*P+1:3*P,2*P+1:3*P); X1=Rt(1:P,P+1:2*P); X2=Rt(1:P,2*P+1:3*P); X3=Rt(P+1:2*P,1:P); X4=Rt(P+1:2*P,2*P+1:3*P); X5=Rt(2*P+1:3*P,1:P); X6=Rt(2*P+1:3*P,P+1:2*P); Q1=R11-X2*inv(X4)*X3; Q2=R22-X3*inv(X5)*X6; Q=zeros(M,M); Q(1:P,1:P)=Q1; Q(P+1:2*P,P+1:2*P)=Q2; Q(2*P+1:M,2*P+1:M)=Q3; QQ=Q^(-1/2); Rw=QQ*R*QQ; %------------------ deal with [V D]=eig(R); tezhen=diag(D) [lambda,index] = sort((diag(D))); UU=V(:,index(1:M-P)); theta=-30:0.1:40; for i = 1:length(theta); for q=1:M; atheta(1,q)=exp(-1j*(q-1)*2*pi*d*sin(theta(i)*pi/180)/l); end WW=atheta*UU*(UU)'*atheta'; Pu(i)=abs(1/WW); end Pu= 10*log10(Pu/max(Pu)); %--------------- [Vw Dw]=eig(Rw); tezhenw=diag(Dw) [lambdaw,indexw] = sort((diag(Dw))); UUw=Vw(:,indexw(1:M-P)); for i = 1:length(theta); for q=1:M; atheta(1,q)=exp(-j*(q-1)*2*pi*d*sin(theta(i)*pi/180)/l); end WW1=atheta*QQ*UUw*(QQ*UUw)'*atheta'; Pu1(i)=abs(1/WW1); end Pu1= 10*log10(Pu1/max(Pu1)); plot(theta,Pu,'k--'); hold on plot(theta,Pu1,'k-') xlabel('DOA/度') ylabel('空间谱/dB') legend('处理前','处理后')

相关推荐
• 经典MUSIC算法
经典MUSIC算法 MATLAB程序 上传给大家共享一下 莫用谢
• music算法估计参数
MUSIC算法是一种基于矩阵特征空间分解的方法。从几何角度讲，信号处理的观测空间可以分解为信号子空间和噪声子空间，显然这两个空间是正交的。信号子空间由阵列接收到的数据协方差矩阵中与信号对应的特征向量组成，...
• music算法及其改进算法
music算法....................................................................................................
• 空间平滑MUSIC算法
空间平滑MUSIC算法 有声压阵和矢量阵的
• music算法的实现
DOA经典算法系列之MUSIC MUSIC算法的深入研究改进算法众多,这里提供Music系列二:噪声子空间法-DOA
• 空间平滑MUSIC算法
空间平滑MUSIC算法 有声压阵和矢量阵的
• MUSIC算法代码
MUSIC算法利用两个互补空间之间的正交特性来估计空间信号的方位
• music算法
基于music算法的doa估计matlab仿真代码
• MUSIC算法程序
MUSIC算法 应用于DOA估计的超分辨算法 多重信号旋转子空间不变算法 波达方向估计最经典的算法
• ROOT-MUSIC算法
经典MUSIC算法 RPPT-MUSIC算法