Take the 2-minute tour ×
MathOverflow is a question and answer site for professional mathematicians. It's 100% free, no registration required.

I'm investigating the eigenvalue ratios

$$ \frac{\lambda_1}{\sum_{j=2}^N\lambda_j} \quad\mbox{and}\quad \frac{\sum_{j=1}^N\lambda_j}{\sum_{j=2}^N\lambda_j} $$

of the NxN matrix $B=AA^T$. $\lambda_1$ denotes the largest eigenvalue. The ratios can be thought of as a measure of "rank-1-ness" of $B$.

I haven't found any mention of either ratio in literature, regardless of constraints on $B$. Has any work been done on this before, and could someone point me there?

share|improve this question
    
Not heard of this before, but what do you do to make sure the denominator isn't 0? –  Anthony Quas Jun 18 '12 at 15:17
    
@Anthony: The matrix for my application originates from a noisy sampling process, so the denominator is almost surely positive. –  Anna Jun 18 '12 at 15:37
    
@Anna: Which data are available? Is your matrix $B$ a sample covariance matrix? –  Stanislav Jun 18 '12 at 20:08
    
@Stanislav: It is indeed. The only assumption on $A$ so far is that the expected energy of each sample series is equal, i.e. the 2-norm of each column of $A$ is roughly equal. –  Anna Jun 18 '12 at 20:35
1  
@Stanislav: In my case, A is a function of a a lot of parameters, and takes ~0.1s to calculate (Monte Carlo simulation). I'm trying to optimise the ratio with regard to those parameters. But for this question, I'm just asking if some properties of these ratios have been investigated before. –  Anna Jun 19 '12 at 12:18
show 2 more comments

3 Answers

I am not sure if you are allowed to change your objective function, but a natural alternative for measuring "rank-1-ness" is $$ \frac{\lambda_1^2}{\lambda_1^2+\lambda_2^2+\dots+\lambda_n^2}. $$ This ratio is easy to compute: the denominator is the squared Frobenius norm of the matrix (i.e., sum of squares of all entries), and the numerator is the squared spectral radius, which you can easily estimate using a few iterations of the power method.

The norm $\lambda_1+\lambda_2+\dots+\lambda_n$ is called nuclear norm, and as far as I know it is not as easy to compute as the other two.

share|improve this answer
add comment

The eigenvalues you are computing are the singular values of $A.$ The largest singular value is the operator norm of the matrix, the sum of all singular values is another norm (I am not sure what it is called, it's analogous to the $L^1$ norm, while the largest s.v. is the analogue of the $L^\infty$ norm). I am not sure what you are interested in, but if you google "largest singular value", you will see a wealth of information.

share|improve this answer
add comment

Also, if you have a sample covariance matrix, you should look into the Wishart distribution, there is a lot of work on the distribution of its eigenvalues. See Muirhead: "Aspects of Multivariate Statistical Theory".

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.