Wiener Filter, Correlation Matrices

the_dialogue
Messages
77
Reaction score
0
I'm trying to understand the Wiener Filter, and I have a few questions.

1. How can there be such a thing as a correlation matrix of 1 vector. I read here:

R_yy = E[ y(k) * y^T(k) ]

where y(k) is a vector, and y_T(x) is the same vector transposed. I thought correlation represents the degree of correspondence between 2 variables, so how can we say there is a correlation between 2 equivalent vectors?

2. They present a matrix R_yy as the correlation matrix (mentioned above). Then they say vector "r_yy" is the first column of R_yy. If that is so, what are the other columns of R_yy?

3. The book first presents the Wiener filter formulation by saying that in minimizing the MSE criterion, we can find "h_w" (the wiener filter vector). They go on to another form of the Wiener filter by saying:

h_w = h_1 - inv(R_yy)*r_vv

where h_1 was defined as [1 0 0 0 ...]^T and r_vv, I suppose is the correlation (once again I don't know how one can have correlation between 2 equivalent vectors).



Thank you for any help.
 
Physics news on Phys.org
Any ideas?
 
the_dialogue said:
I'm trying to understand the Wiener Filter, and I have a few questions.

1. How can there be such a thing as a correlation matrix of 1 vector. I read here:

R_yy = E[ y(k) * y^T(k) ]

where y(k) is a vector, and y_T(x) is the same vector transposed. I thought correlation represents the degree of correspondence between 2 variables, so how can we say there is a correlation between 2 equivalent vectors?

ryy is the auto-correlation, witch you can think of as the correlation between the signal and shifted versions of itself. The relationship posted by you is actually ryy(0).
ryy(n)= E{y(k)*yT(k-n)}
the_dialogue said:
2. They present a matrix R_yy as the correlation matrix (mentioned above). Then they say vector "r_yy" is the first column of R_yy. If that is so, what are the other columns of R_yy?

That matrix is a symmetric toeplitz matrix. Indeed the first column is ryy, the second is [ryy(1) ryy(0) ryy(1) ryy(2) ... ryy(N-1)]T and so on the last one is ryy reversed.
the_dialogue said:
3. The book first presents the Wiener filter formulation by saying that in minimizing the MSE criterion, we can find "h_w" (the wiener filter vector). They go on to another form of the Wiener filter by saying:

h_w = h_1 - inv(R_yy)*r_vv

where h_1 was defined as [1 0 0 0 ...]^T and r_vv, I suppose is the correlation (once again I don't know how one can have correlation between 2 equivalent vectors).

I don't quite understand what you mean by your notation here.
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top