kt_estimate_space_chapter3_5.asv

阵列信号处理中的KT法进行DOA估计

function kt_estimate_space_chapter3_5(y_db,M)
Y_db=zeros(24,8);
for i=1:2*(20-M)%表示所需的矩阵数目
    if i<=(20-M)
    Y_db(i,:)=y_db(1,7+i:-1:i);
% Y_0db(2,:)=y_0db(1,9:-1:2);
% Y_0db(3,:)=y_0db(1,10:-1:3);
    else
      if i>=13
            Y_db(i,:)=conj(y_db(1,26-i:33-i));
       end
    end
end
yout_db=zeros(24,1);
yout_db=[y_db(9:20).';y_db(12:-1:1)'];
%%%KT%%%%
R=Y_db'*Y_db;
[u,s,v]=svd(R)%s为特征值矩阵,对矩阵进行奇异值分解
N=2;%N用来取矩阵中有意义的特征值
u=u(:,1:N);
v=v(:,1:N);
for i=1:N
    s(i,i)=1/s(i,i);
end
s=s(1:N,1:N);
 h_guji=u*s*v'*Y_db'*yout_db;
 gm=zeros(9,1);
 gm=[1;-h_guji];
 theta=linspace(-100,100,500);
theta1=theta*pi/180;
Y_theta=zeros(1,500);
for k=1:500
    a_theta=zeros(M+1,1);
    a_theta(1,1)=1;
    for n=2:M+1
        a_theta(n,1)=exp(j*(n-1)*pi*sin(theta1(k)));
    end
    a_theta=a_theta.';
    X_theta(k)=1/abs(a_theta*gm).^2;
end
plot(theta,X_theta)