应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 通信科学 » JJM算法的问题
51/1返回列表
查看: 854  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
1楼 2011-08-09 15:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingdao89

新虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
JJM算法么的研究过啊。。。
赞一下(1人) 回复此楼
生活的一半是爱情，爱情的全部是生活。
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的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:15 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
信息提示
关闭
请填处理意见
关闭