calculate_music.m

来自「适用于线列阵的常用定位算法」· M 代码 · 共 26 行

% Function of MUSIC  for linear array
% A is the direction vector  (M*1)
% RXX is the covariance matrix (M*M)

function P = Calculate_MUSIC(No_elements,No_signals,f0,D,Rxx)
    
     C=345;    

   [Eig_matrix,Eig_value]=svd(Rxx);
   
   V_matrix=Eig_matrix(:,No_signals+1:No_elements);
    %%%%% 噪声子空间MUSIC
    P = zeros(181);
    i=1;
    for theta = -90:90
        k = (0:No_elements-1)';
        A = exp(-j*2*pi*f0*k.*D*sin(theta*pi/180)/C);
        P(i)= 1/(A'*V_matrix*V_matrix'*A);
        i=i+1;
    end

%     P=abs(P)/max(abs(P));
P = 10*log10(abs(P));