24小时热门版块排行榜

返回列表

promisewind

金虫 (著名写手)

应助: 3 (幼儿园)
金币: 3511.7
散金: 2776
红花: 5
帖子: 1022
在线: 238.9小时
虫号: 1253129
注册: 2011-04-02
专业: 通信理论与系统

[交流] JJM算法的问题已有1人参与

自己编写了一个JJM算法，可是怎么出现的结果不能估计出弱信号的doa，得到的是与强信号对称的一个方向的信号，不知道那里错了，还请大家赐教。

N=256;                         %设定采样点数
M =12;                            %设定阵元数
fs=100;                            %设定采样频率
kd=1/2;                         %阵元间距与波长的比值
f=[14 6];                         %设两个信号元的频率，单位为MHz
theta = [5, 34]*pi/180;             %入射方向角向量

SNR=[-5, 40];                         % 入射信号信噪比
AMP = 10.^(SNR./20);

m=0:M-1;
n=0:N-1;
S = repmat(AMP',1,N).*exp(j*(2*pi*f.'/fs*n));  %产生两个远场正弦信号
A = exp(j*2*pi*kd*m'*sin(theta));                            %方向矩阵

X=A*S;
randn('seed',sum(100*clock));
noise  = randn(M,N) + j*randn(M,N);
X=X+noise;                                                 %接受数据矩阵

Rx=X*(X')/N;

theta_stro = [34]*pi/180;
T1=eye(M,M-1);
for l=2:M
T1(l,l-1)=-exp(-j*pi*sin(theta_stro(1)));
end

TT=T';

Rx1=TT*Rx*TT';                                                 %计算协方差矩阵
%*********************************************************************
[V1,D1]=eig(Rx1,2*TT*TT');                   %自相关矩阵的特征分解,,得到特征值矩阵和与特征值一一对应的特征向量矩阵
[Q1,index1]=sort(diag(D1));       %对特征值排序

for k=1:M-length(theta)       %将与特征值对应的特征向量也按照特征值的顺序排列
En1(:,k)=V1(:,index1(k));       %取排了序的特征向量的前M-length(theta)个来构造噪声子空间En
end

a_theta=T'*exp(-j*2*pi*kd*[0:M-1].'*sin(theta_axis*pi/180));
for k=1:MD                      %进行谱峰搜索
d1(k)=a_theta(:,k)'*a_theta(:,k)/(a_theta(:,k)'*En1*En1'*a_theta(:,k));
k=k+1;
end

plot(theta_axis,10*log10(abs(d1)),'k');

xlabel('到达角/\circ');
ylabel('归一化空间谱/dB');

回复此楼