24小时热门版块排行榜    

查看: 643  |  回复: 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的回帖

promisewind

金虫 (著名写手)

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

bingdao89

新虫 (初入文坛)


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

promisewind

金虫 (著名写手)

理论上感觉是对的,怎么仿真不出来?唉,同道中人何在?
Annie
4楼2011-08-19 11:09:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 promisewind 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 298求调剂 +3 Vv呀! 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
信息提示
请填处理意见