互质阵列DOA估计，不用空间平滑

学习小混子

已于 2024-09-04 11:22:01 修改

阅读量219

点赞数 1

文章标签：机器学习算法人工智能

于 2024-09-04 11:18:46 首次发布

本文链接：https://blog.csdn.net/weixin_48627103/article/details/141888379

版权

clear all
close all
clc
M = 5;%子阵一阵元数
N = 3;%子阵二阵元数
m = 0:M-1;
n = 0:N-1;
array= union(M*n,N*m);%不用扩展扩展互质阵列
arraymax = max(array);
snap = 500;
f0 = 1e4;%载频
fs = 2*f0;%采样频率
t = [1:snap]/fs;
theta =-50:100/6:50;%多到达角
K = length(theta);
k = rand(K,1);
S = exp(1j*2*pi*(k*f0.*t));
A = exp(-1i*pi*array'*sin(theta/180*pi));
x= A*S;
snr = 30;
x = awgn (x,snr,'measured');
vir_x = vec(x*x'/snap);
%% 构造差分阵列接收单快拍数据
len = length(array);
% diff_array = zeros(1,len^2) ;
for j=1:len
 B(:,j)=array-array(j);      % 差阵列（包含重复项）
end
B = B(:);
diff_sum = unique(sort(B));
arraymax = max(diff_sum);
xv_h =zeros(2*arraymax+1,1);
for i = -arraymax:arraymax
    [a,b]=find(B == i);
    if(a)
        xv_h(i+1+arraymax,:) = mean(vir_x(a,:));
    end
end
%%取非零数
x_u = xv_h(-length(array)+arraymax+1:length(array)+arraymax+1);
L = (length(x_u)+1)/2;
R = toeplitz(x_u(L:end));
R_SS = 1/L*(R*R');
[EV,Dv] = eig(R_SS);%特征值分解
DD = diag(Dv);%将特征值变为向量形式
[DD,I] = sort(DD);%从小到大
DD = fliplr(DD');%翻转函数，从大到小
EV = fliplr(EV(:,I));
En = EV(:,K+1:end);%噪声子空间
angle=-60:0.01:60;
for m=1:length(angle)
    a=exp(1i*pi*(0:L-1)'*sin(angle(m)*pi/180));
    Pmusic(m)=1/(a'*En*En'*a);
end
Pmusic = abs(Pmusic);
Pmax = max(Pmusic);
Pmusic_db = pow2db(Pmusic/Pmax);
plot(angle,Pmusic_db,'Linewidth',1.5,'Color',[0 0 1]);
hold on
for n = 1:K
    plot([theta(n),theta(n)],ylim,'r-.','Linewidth',1);
end
grid on
xlabel('angle(degree)');
ylabel('Spatial spectrum(dB)');