xinhaokongjianzisuanfa.txt

语音端点 检测的一种方法.叫空间信号子算法

%%%DOA经典算法系列之MUSIC MUSIC算法的深入研究改进算法众多,这里提供Music系列一:信号子空间法-DOA series of classical algorithm MUSIC MUSIC - depth study of the many improved algorithm, Music here for a series : signal-space law %%%


文件头:
%MUSIC ALGORITHM

%DOA Estimation by ULA

%信号子空间法。

clear all; 

close all;

clc

N_x=1024; % Length of Signal 

N=8; % Size of Rx Matrix 

A=[4 3 3];

M=3; % Number of Signals 

w=[pi/6 pi/10 pi/8]';%信号频率 与2不同之处

K=0:N-1;

xx=pi/3;yy=0*pi/4;zz=-pi/6; %三个信号的入射角,小于5度就难以辨别了

B=[exp(-j*pi*K'*sin(xx)) exp(-j*pi*K'*sin(yy)) exp(-j*pi*K'*sin(zz))] ; %阵列流型,信号源决定行数，阵元数决定列数

xxx=2*exp(j*w*[0:N_x-1]);%仿真信号

x=B*xxx+randn(8,N_x)+j*randn(8,N_x); %加噪声

R=x*x';

[V D]=eig(R);

[lambda,index] = sort((diag(D)));

UU=V(:,index(6:8));

theta=-90:0.1:90;%ULA估计角度变化的范围和频率选择

for i = 1:length(theta)

AA=exp(-j*pi*K*sin(theta(i)*pi/180));%方向矢量

WW=AA*(eye(size(R))-UU*UU')*AA';

Pmusic(i)=abs(1/WW);%角谱

% PCap(i) = abs(1/(AA*inv(R)*AA'));

end

Pmusic = 20*log10(Pmusic);

sita=-90:0.1:90;

plot(sita,Pmusic);

grid

xlabel('rad')