高埗东莞网站建设,南阳网(网站).,重庆网站建设论坛,有什么做礼品的卖家网站作者#xff1a;桂。 时间#xff1a;2017-10-24 21:50:16 链接:http://www.cnblogs.com/xingshansi/p/7726082.html 前言 记录阵列信号在DOA估计中#xff0c;信源个数估计的基本方法。 一、基本估计方法 参考#xff1a;王永良《空间谱估计》p42: 以MDL为例#xff1a… 作者桂。 时间2017-10-24 21:50:16 链接:http://www.cnblogs.com/xingshansi/p/7726082.html 前言 记录阵列信号在DOA估计中信源个数估计的基本方法。 一、基本估计方法 参考王永良《空间谱估计》p42: 以MDL为例 function [num,a] MDL_est(Rx,ssNum)
%ssNum快拍数
T ssNum;
arrNum size(Rx,1);
K arrNum;
[u,s,v] svd(Rx);
sd diag(s);
a zeros(1,K);
for m 0:K-1negv sd(m1:K);Tsph mean(negv)/((prod(negv))^(1/(K-m)));a(m1) T*(K-m)*log(Tsph) m*(2*K-m)*log(T)/2;
% a(m1) T*(K-m)*log(Tsph) m*(2*K-m)*log(T)/2;
end
[y,b] min(a);
num b - 1;二、相干信号个数估计 如果采用空间平滑的思路平滑后的矩阵直接利用一般估计思路即可实现。 当为了避免空间平滑的盲区采用修正空间平滑n为较大的正数如1024 Ryy n * Ry1 - (n-1)*Ry2; 此时的R利用MDL、AIC等估计错误可简单根据特征值判断粗糙的思路 其中 D [];
for m 0:K-1
% negv sd(m1:K);
% Tsph mean(negv)/((prod(negv))^(1/(K-m)));
% a(m1) T*(K-m)*log(Tsph) m*(2*K-m)*log(T)/2;
% a(m1) T*(K-m)*log(Tsph) m*(2*K-m)*log(T)/2;D(m1) sd(m1)-mean(sd);
end% [y,b] min(a);
% num b - 1;
num sum(D0);
