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

  • MHB
  • Thread starter mathmari
  • Start date
  • Tags
    Condition
In summary: 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?
  • #1
mathmari
Gold Member
MHB
5,049
7
Hey! :eek:

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)
 
Mathematics news on Phys.org
  • #2
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)
 
  • #3
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)
 
  • #4
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)
 

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

What is the condition number of A'A?

The condition number of A'A is a measure of how sensitive the solution of a linear system is to changes in the input data. It is calculated by taking the square root of the ratio of the largest eigenvalue to the smallest eigenvalue of A'A.

Why is the condition number of A'A important?

The condition number of A'A is important because it indicates the stability of the solution to a linear system. A lower condition number means that the solution is less sensitive to changes in the input data, while a higher condition number means that the solution is more sensitive and may be less reliable.

How is the condition number of A'A used in practice?

The condition number of A'A is often used as a measure of the quality of a matrix in numerical analysis. It can help determine the accuracy and stability of algorithms used to solve linear systems, and can also be used to identify ill-conditioned matrices that may cause problems in computations.

Can the condition number of A'A be negative?

No, the condition number of A'A cannot be negative. It is always a positive value, as it is calculated using the square root of a ratio of eigenvalues, which are always positive. A negative condition number would not make sense in the context of linear systems.

How can the condition number of A'A be improved?

The condition number of A'A can be improved by using matrix factorization techniques, such as LU decomposition or QR decomposition, to reduce the number of operations needed to solve a linear system. It can also be improved by using more accurate numerical methods and increasing the precision of the calculations.

Back
Top