24小时热门版块排行榜    

Znn3bq.jpeg
查看: 657  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

promisewind

金虫 (著名写手)

[交流] 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');
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

Annie
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingdao89

新虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖
JJM算法么的研究过啊。。。
生活的一半是爱情,爱情的全部是生活。
3楼2011-08-18 19:14:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

promisewind

金虫 (著名写手)

呜呜,没有人做相关的吗?
Annie
2楼2011-08-16 22:29:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

promisewind

金虫 (著名写手)

理论上感觉是对的,怎么仿真不出来?唉,同道中人何在?
Annie
4楼2011-08-19 11:09:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 初试分332,一志愿报考西北工业大学, +9 故人?? 2026-04-09 9/450 2026-04-09 19:52 by gong120082
[考研] 085500求调剂材料 +4 易11122 2026-04-09 4/200 2026-04-09 17:16 by 1753564080
[材料工程] 材料调剂推荐 +3 蛋糕x2 2026-04-07 3/150 2026-04-09 17:01 by Lilly_Li
[考研] 求调剂 +7 chenxrlkx 2026-04-05 9/450 2026-04-09 09:04 by wj165256
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 347材料专硕求调剂 +18 zj8215216 2026-04-06 18/900 2026-04-08 16:27 by luoyongfeng
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 285求调剂 +8 AZMK 2026-04-04 11/550 2026-04-06 13:56 by BruceLiu320
[考研] 085600,320分求调剂 +16 大馋小子 2026-04-04 17/850 2026-04-06 07:58 by MOF_Catal
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 0835学硕299求调剂 08大类可接受 +5 useryy 2026-04-03 5/250 2026-04-04 20:07 by 蓝云思雨
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 一志愿重庆大学085404,总分314分,求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 11408,284分,二战真诚求调剂 +4 12.27 2026-04-02 4/200 2026-04-03 14:14 by dxiaoxin
信息提示
请填处理意见