| 查看: 652 | 回复: 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'); |
回复此楼» 猜你喜欢
材料考研求调剂总分280
已经有3人回复
-
求调剂
已经有4人回复
-
材料调剂
已经有10人回复
-
求调剂到材料
已经有3人回复
-
化工求调剂！
已经有9人回复
-
301求调剂
已经有13人回复
-
308求调剂
已经有13人回复
-
材料工程310专硕调剂
已经有15人回复
-
304求调剂
已经有9人回复
-
081700，311，求调剂
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关遗传算法的问题求助
已经有3人回复
- FFT算法的一个问题请教
已经有18人回复
- 一个算法多做个领域的问题，算不算一稿多投
已经有17人回复
- vasp算法的一些基础问题
已经有10人回复
- 请教遗传算法三个问题
已经有10人回复
- 请教关于遗传算法编程的问题（俺是小白）
已经有4人回复
- 关于免疫算法的几个问题。（搜索无重复）
已经有8人回复
- 【求助】DLS中数据算法问题
已经有10人回复
- 【求助】DLS中数据算法问题
已经有16人回复
- 关于算法比较的问题，请大家帮忙
已经有7人回复
- 【讨论】遗传算法解决象棋问题
已经有8人回复
- 【讨论】多目标跟踪算法问题
已经有10人回复
promisewind
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 4045.7
- 散金: 2776
- 红花: 5
- 帖子: 1070
- 在线: 238.9小时
- 虫号: 1253129
- 注册: 2011-04-02
- 专业: 通信理论与系统
2楼2011-08-16 22:29:21
么的研究过啊。。。 3楼2011-08-18 19:14:04
promisewind
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 4045.7
- 散金: 2776
- 红花: 5
- 帖子: 1070
- 在线: 238.9小时
- 虫号: 1253129
- 注册: 2011-04-02
- 专业: 通信理论与系统
4楼2011-08-19 11:09:18
