Is the Condition Number of A'A Related to its Matrix Norm?

  • Context: MHB 
  • Thread starter Thread starter mathmari
  • Start date Start date
  • Tags Tags
    Condition
Click For Summary
SUMMARY

The discussion centers on the relationship between the condition number of a matrix \( A \) and its transpose product \( A^TA \). Specifically, it establishes that the 2-norm condition number \( k_2(A^TA) \) is equal to the square of the condition number of \( A \), expressed as \( k_2(A^TA) = (k_2(A))^2 \). The participants explore the mathematical derivations involving the norms and properties of matrices, confirming that the maximum and minimum values of the norms are preserved under these transformations.

PREREQUISITES
  • Understanding of matrix norms, specifically 2-norms.
  • Familiarity with the concept of condition numbers in linear algebra.
  • Knowledge of matrix multiplication and properties of transpose matrices.
  • Basic proficiency in eigenvalues and eigenvectors, particularly in relation to symmetric matrices.
NEXT STEPS
  • Study the derivation of the condition number for various types of matrices, including singular and non-singular matrices.
  • Learn about the implications of condition numbers on numerical stability in algorithms.
  • Explore the properties of eigenvalues and eigenvectors in relation to matrix norms.
  • Investigate the relationship between matrix norms and their applications in optimization problems.
USEFUL FOR

Mathematicians, data scientists, and engineers who work with linear algebra, particularly those involved in numerical analysis and optimization. This discussion is beneficial for anyone seeking to deepen their understanding of matrix properties and their implications in computational contexts.

mathmari
Gold Member
MHB
Messages
4,984
Reaction score
7
Hey! :o

We have a matrix $A\in \mathbb{R}^{m\times n}$ which has the rank $n$. The condition number is defined as $\displaystyle{k(A)=\frac{\max_{\|x\|=1}\|Ax\|}{\min_{\|x\|=1}\|Ax\|}}$.

I want to show that $k_2(A^TA)=\left (k_2(A)\right )^2$. We have that $$k_2(A^TA)=\frac{\max_{\|x\|_2=1}\|(A^TA)x\|_2}{\min_{\|x\|_2=1}\|(A^TA)x\|_2}$$

It holds that $$\|(A^TA)x\|_2=\left ((A^TA)x, (A^TA)x\right )=\left (A^TAx\right )^T\left (A^TAx\right )=x^TA^TAA^TAx=\left (Ax\right )^T\left (AA^T\right )\left (Ax\right )$$
We also have that $$\|Ax\|_2=\left (Ax, Ax\right )=\left (Ax\right )^T\left (Ax\right )=x^TA^TAx=x^T\left (A^TA\right )x$$

How do we continue? (Wondering)
 
Physics news on Phys.org
Hey mathmari!

Shouldn't it be $\|Ax\|_2 = \sqrt{(Ax,Ax)}$ or $\|Ax\|_2^2 = (Ax,Ax)$? (Worried)

So we have:
$$\|Ax\|_2^2=x^T\left (A^TA\right )x$$

We also have:
$$\|(A^TA)x\|_2^2=x^TA^TAA^TAx=x^T(A^TA)^2x$$

Suppose $x$ is a vector of unit length for which $x^T\left (A^TA\right )x$ takes on its maximal value.
Then $x^T(A^TA)^2x$ will also take its maximal value, won't it? (Wondering)
 
Klaas van Aarsen said:
Shouldn't it be $\|Ax\|_2 = \sqrt{(Ax,Ax)}$ or $\|Ax\|_2^2 = (Ax,Ax)$? (Worried)

Ohh yes! (Blush)
Klaas van Aarsen said:
So we have:
$$\|Ax\|_2^2=x^T\left (A^TA\right )x$$

We also have:
$$\|(A^TA)x\|_2^2=x^TA^TAA^TAx=x^T(A^TA)^2x$$

Suppose $x$ is a vector of unit length for which $x^T\left (A^TA\right )x$ takes on its maximal value.
Then $x^T(A^TA)^2x$ will also take its maximal value, won't it? (Wondering)

Will $x^T(A^TA)^2x$ take also its maximal value because the term in the middle is the square of the previous one? (Wondering)
 
mathmari said:
Will $x^T(A^TA)^2x$ take also its maximal value because the term in the middle is the square of the previous one?

Well, to be fair, this is not immediately obvious. (Worried)

Let's take a slightly different angle.

According to the wiki page about matrix norms, a matrix $A$ has norm $\|A\|_2 =\sup\limits_{\|x\|=1} \|Ax\|_2 = \sqrt{\lambda_{\text{max}}(A^*A)}=\sigma_{\text{max}}(A)$, where $A^*$ is the conjugate transpose, which is just the transpose $A^T$ for a real matrix.

So $\|A^TA\|_2 = \sqrt{\lambda_{\text{max}}((A^TA)^T A^TA)}$, isn't it? (Wondering)
 

Similar threads

  • · Replies 9 ·
Replies
9
Views
3K
Replies
4
Views
2K
Replies
12
Views
4K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
10K
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
10K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 1 ·
Replies
1
Views
10K
Replies
6
Views
2K