作者：桂。

時(shí)間：2017-10-24??21:50:16

鏈接:http://www.cnblogs.com/xingshansi/p/7726082.html?

---

前言

記錄陣列信號(hào)在DOA估計(jì)中，信源個(gè)數(shù)估計(jì)的基本方法。

一、基本估計(jì)方法

參考：王永良《空間譜估計(jì)》p42:

以MDL為例：

function [num,a] = MDL_est(Rx,ssNum) %ssNum：快拍數(shù) 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(m+1:K);Tsph = mean(negv)/((prod(negv))^(1/(K-m)));a(m+1) = T*(K-m)*log(Tsph) + m*(2*K-m)*log(T)/2; % a(m+1) = T*(K-m)*log(Tsph) + m*(2*K-m)*log(T)/2; end [y,b] = min(a); num = b - 1;

　　

二、相干信號(hào)個(gè)數(shù)估計(jì)

如果采用空間平滑的思路，平滑后的矩陣直接利用一般估計(jì)思路，即可實(shí)現(xiàn)。

當(dāng)為了避免空間平滑的盲區(qū)，采用修正空間平滑（n為較大的正數(shù)，如1024）：

Ryy = n * Ry1 - (n-1)*Ry2;

此時(shí)的R利用MDL、AIC等估計(jì)錯(cuò)誤，可簡(jiǎn)單根據(jù)特征值判斷（粗糙的思路）：

其中：

D = []; for m = 0:K-1 % negv = sd(m+1:K); % Tsph = mean(negv)/((prod(negv))^(1/(K-m))); % a(m+1) = T*(K-m)*log(Tsph) + m*(2*K-m)*log(T)/2; % a(m+1) = T*(K-m)*log(Tsph) + m*(2*K-m)*log(T)/2;D(m+1) = sd(m+1)-mean(sd); end% [y,b] = min(a); % num = b - 1; num = sum(D>0);

　　

總結(jié)

以上是生活随笔為你收集整理的空间谱专题16：信号个数估计的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。