Digital Signal Processing Reference
In-Depth Information
eigenvalue, then it immediately comes esp { V b } ⊂ esp { A } . Inversely, any vector
of esp { A } is an eigenvector of x Γ associated with the null eigenvalue. Hence
esp { V b } = esp { A } . This implies that the other eigenvectors of x Γ associated
with strictly positive eigenvalues, given esp { V s }, are in esp { A }. Then we have
esp { V s } ⊂ esp { A }. Now, x Γ being hermitian, we know that there is an
orthonormal basis of eigenvectors, which implies that the dimension of esp { V s } is
P. Hence esp { V s } = esp { A }.
In the presence of a white noise, the eigenvalues of:
H2
ΓΓ
=
AA
+
σ
I
xx
ss
are those of A Γ ss A H increased by σ 2 and the eigenvectors are unchanged, hence the
results of Theorem 8.1.
We indicate by signal subspace the subspace esp { V s } spanned by the
eigenvectors associated with the biggest P eigenvalues of Γ xx , and by noise subspace
the subspace esp { V b } spanned by the M - P eigenvectors associated with the
smallest eigenvalue σ 2 . Thus, it immediately results that the complex sinusoid
vectors a ( f i ) are orthogonal to the noise subspace, given:
H
()
Va
f
==
0,
i
1,
,
P
[8.16]
b
i
Thus, the signal vectors corresponding to the searched frequencies are the
vectors a ( f ) of the form given in [8.9] which are orthogonal to V b .
A number P of independent vectors of the form of [8.9] should exist in a
subspace of dimension P of M so that the resolution of the non-linear system
[8.16] gives the searched frequencies in a unique way. This is true provided that
P + 1 < M.
THEOREM 8.2. If P + l < M, any system of P vectors of the form a ( f ) given in [8.9],
for different values of f, forms a free family of
M .
DEMONSTRATION. In M , an infinity of vectors of the form a ( f ) [8.9] exist
when f scans the set of real numbers. However, any family of I < M vectors
{ a ( f 1 ) , …, a ( f I )}, where the parameters f i are different, is free because it is easy to
verify that the matrix formed by these vectors is of the Vandermonde type. Thus, in
a subspace of dimension I < M, only I vectors of this type exist. Actually, if another
one existed, noted { a ( f I+ 1 )} with f I+ 1 f i for i = 1, …, I in this subspace of dimension
I , it would be written as a linear combination of these I vectors and the system { a ( f 1 )
, …, a ( f I+ 1 )} would not be free, which is impossible if I + 1 < M. Therefore, in
Search WWH ::




Custom Search