空間譜專題16 訊號個數估計

前言

記錄陣列訊號在doa估計中，信源個數估計的基本方法。

以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-1
ne** = sd(m+1:k);
tsph = mean(ne**)/((prod(ne**))^(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;

當為了避免空間平滑的盲區，採用修正空間平滑（n為較大的正數，如1024）：

ryy = n * ry1 - (n-1)*ry2;

此時的r利用mdl、aic等估計錯誤，可簡單根據特徵值判斷（粗糙的思路）：

其中：

d = ;

for m = 0:k-1
% ne** = sd(m+1:k);
% tsph = mean(ne**)/((prod(ne**))^(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);

空間譜專題01 訊號模型

一 訊號基本模型 根據波動方程，對於球面座標 可以得出單頻解 當訊號距離接收單元距離非常遠 l 時，即可認為遠場，這時候球面波可近似看作平面波 對應的解 其中k 2pi 稱為波數向量，其大小表示單位波長的週期數。二 窄帶訊號 a 窄帶定義 關於窄帶訊號的定義不唯一，先從實際情況入手分析分析，例如一般...

空間譜專題16 間距選取分析

前言 一 問題描述 對於空間譜測向 以music演算法為例 如果間距過大 超過半波長 可能存在相位模糊問題 對應空間譜就是偽峰 干涉儀就是模糊 beamforming就是柵瓣。以一維陣列舉例，這個時候通常借助非均勻陣來降低相位模糊問題。二 原理分析 假設測角範圍為 對應導向向量為 exp j 2 p...

空間譜專題06 寬頻訊號處理思路

前言 目前分析的問題，仍然限定在布陣的環節，暫不涉及後處理及硬體實現。一 寬頻處理的一般方式前面分析的陣列訊號模型，都是建立在窄帶訊號的基礎上，對於寬頻陣列訊號，有兩個思路 1 通道化為窄帶訊號，按窄帶的思路進行處理 2 利用聚焦變換，聚焦到特定頻率，進行處理。聚焦變換的思路對相干訊號同樣有效，這樣...