Finding Contract Transformation of Block Matrix by Elementary Transformation Method

1. Use the "elementary transformation method" to find the diagonal matrix contracted by the symmetric matrix

In order to make this article complete, here is how to use the "elementary transformation method" to find the diagonal matrix contracted by the symmetric matrix.

When learning linear algebra, we know that for any n × nn \times nn×Real symmetric matrixAA of nA , you can use the "elementary transformation method" to find the invertible matrixCCC and Diagonal MatrixDDD , such thatAAA andDDD contract, ieCTAC = DC^TAC=DCTAC=D , whereCTC^TCT stands forCCTranspose of C.

Specific method: make 2 n × n 2n \times n2 n×n- matrix

[ AI ] → perform the same elementary column transformation on 2 n × n matrix, perform elementary row transformation on A [ DC ] \left[ \begin{matrix} A\\ I \end{matrix} \right] \xrightarrow[for 2n \times n Matrix performs the same elementary column transformation]{Perform elementary row transformation on A} \left[ \begin{matrix} D\\ C \end{matrix} \right][AI]Perform elementary row transformation on A Apply the same elementary column transformation to a 2 n × n matrix[DC]

C T A C = D C^TAC=D CTAC=D

Note that an elementary row transformation needs to be applied first, and then the same elementary column transformation is applied to the entire column at once.

Example 1

Known real symmetric matrix
A = [ 1 1 1 1 2 3 1 3 5 ] A= \left[ \begin{matrix} 1 & 1 & 1\\ 1 & 2 & 3\\ 1 & 3 & 5\\ \ end{matrix} \right]A= 111123135
Finding the Reversible Matrix CC by Elementary Transformation MethodC and diagonal matrixDDD , such thatAAA andDDD contract.
[ A I ] = [ 1 1 1 1 2 3 1 3 5 1 0 0 0 1 0 0 0 1 ] → c 2 − c 1 r 2 − r 1 [ 1 0 1 0 1 2 1 2 5 1 − 1 0 0 1 0 0 0 1 ] → c 3 − c 1 r 3 − r 1 [ 1 0 0 0 1 2 0 2 4 1 − 1 − 1 0 1 0 0 0 1 ] → c 3 − c 2 r 3 − r 2 [ 1 0 0 0 1 0 0 0 0 1 − 1 1 0 1 − 2 0 0 1 ] \left[ \begin{matrix} A\\ I \end{matrix} \right]= \left[ \begin{matrix} 1 & 1 & 1\\ 1 & 2 & 3\\ 1 & 3 & 5\\ 1 & 0 & 0\\ 0 & 1 & 0\\ 0 & 0 & 1\\ \end{matrix} \right] \xrightarrow[c_2-c_1]{r_2-r_1} \left[ \begin{matrix} 1 & 0 & 1\\ 0 & 1 & 2\\ 1 & 2 & 5\\ 1 & -1 & 0\\ 0 & 1 & 0\\ 0 & 0 & 1\\ \end{matrix} \right] \xrightarrow[c_3-c_1]{r_3-r_1} \left[ \begin{matrix} 1 & 0 & 0\\ 0 & 1 & 2\\ 0 & 2 & 4\\ 1 & -1 & -1\\ 0 & 1 & 0\\ 0 & 0 & 1\\ \end{matrix} \right] \xrightarrow[c_3-c_2]{r_3-r_2} \left[ \begin{matrix} 1 & 0 & 0\\ 0 & 1 & 0\\ 0 & 0 & 0\\ 1 & -1 & 1\\ 0 & 1 & -2\\ 0 & 0 & 1\\ \end{matrix} \right] [AI]= 111100123010135001 r2r1 c2c1 101100012110125001 r3r1 c3c1 100100012110024101 r3r2 c3c2 100100010110000121

The required invertible matrix CCC and diagonal matrixDDD also
C = [ 1 − 1 1 0 1 − 2 0 0 1 ] , D = [ 1 0 0 0 1 0 0 0 0 ] C= \left[ \begin{matrix} 1 & -1 & 1\\ 0 & 1 & -2\\ 0 & 0 & 1\\ \end{matrix} \right], D= \left[\begin{matrix} 1 & 0 & 0\\ 0 & 1 & 0\\ 0& 0&0\\\end{matrix}\right]C= 100110121 ,D= 100010000

C T A C = D C^TAC=D CTAC=D

2. Use the "elementary transformation method" to find the block diagonal matrix contracted by the block symmetric matrix

Next, we extend the real symmetric matrix to the block symmetric matrix, and the elements in the matrix are no longer limited to real numbers, but can be complex numbers. By analogy to the above elementary transformation method without proof, we give the elementary transformation method for finding the block diagonal matrix contracted by the block symmetric matrix.

For an n × nn \times nn×Block Symmetric MatrixAA of nA , you can use the "elementary transformation method" to find the invertible matrixCCC and Block Diagonal MatrixDDD , such thatAAA andDDD contract, namelyCHAC = DC^HAC=DCHAC=D , whereCHC^HCH meansCCThe conjugate transpose of C.

Specific method: make 2 n × n 2n \times n2 n×Matrix of n
[ AI ] → perform the same elementary column transformation on 2 n × n matrix [ DC ] \left[ \begin{matrix} A\\ I \end{matrix} \right] \xrightarrow [Perform the same elementary column transformation on 2n\times n matrix]{Perform elementary row transformation on A} \left[ \begin{matrix} D\\ C \end{matrix} \right][AI]Perform elementary row transformation on A Apply the same elementary column transformation to a 2 n × n matrix[DC]
C H A C = D C^HAC=D CHAC=D

In addition, it is also necessary to understand the elementary transformation rules of the block matrix, as follows:

  • Elementary Row Transformation Rules for Block Matrix

    • Put a block row to the left PPP times (PPP is the matrix) added to another block row;
    • Swap the position of two block lines
    • Left -multiply a block of rows by an invertible matrix
  • Elementary Column Transformation Rules of Block Matrix

    • Put a block column on the right PPP times (PPP is matrix) added to another block column;
    • Swap the positions of two block columns
    • Right- multiply a block of columns by an invertible matrix
  • It should be noted that there is a difference between "left" and "right" in row transformation and column transformation.

  • And (emphasis added), if you put a block row on the left PPP times (PPP is a matrix) is added to another block row, and the corresponding elementary column transformation is to convert theright PHP^HPH (note the conjugate transpose here) is multiplied to another block column.

It will be easier to understand directly by looking at the example below.

Example 2

Environmental value
B = [ A α α H 1 λ 0 ] B= \left[ \begin{matrix} A & \alpha\\ \alpha^H & \frac{1}{\lambda_0}\end {matrix}\right]B=[AaHal01]
Among themAAA isn × nn \times nn×The Hermite matrix of n (ieAH = AA^H=AAH=A ), andAAA is positive definite,α \alphaα isnnn- dimensional unit column vector. Finding the Reversible Matrix CCby Elementary Transformation MethodC and diagonal matrixDDD , such thatBBB andDDD contract.

[ B I ] = [ A α α H 1 λ 0 I n O O 1 ] → c 2 − c 1 ( α H A − 1 ) H r 2 − α H A − 1 r 1 [ A O O 1 λ 0 − α H A − 1 α I n − ( A H ) − 1 α O 1 ] \left[ \begin{matrix} B\\ I\\ \end{matrix} \right]= \left[ \begin{matrix} A & \alpha\\ \alpha^H & \frac{1}{\lambda_0}\\ I_n & O\\ O & 1 \end{matrix} \right] \xrightarrow[c_2-c_1(\alpha^HA^{-1})^H]{r_2-\alpha^HA^{-1}r_1} \left[ \begin{matrix} A & O\\ O & \frac{1}{\lambda_0}-\alpha^HA^{-1}\alpha\\ I_n & -(A^H)^{-1}\alpha\\ O & 1 \end{matrix} \right] [BI]= AaHInOal01O1 r2aHA1r1 c2c1( aHA1)H AOInOOl01aHA1 a(AH)1 a1

The required invertible matrix CCC and diagonal matrixDDD are

C = [ I n − ( A H ) − 1 α O 1 ] , D = [ A O O 1 λ 0 − α H A − 1 α ] C= \left[ \begin{matrix} I_n & -(A^H)^{-1}\alpha\\ O & 1 \end{matrix} \right], D= \left[ \begin{matrix} A & O\\ O & \frac{1}{\lambda_0}-\alpha^HA^{-1}\alpha\\ \end{matrix} \right] C=[InO(AH)1 a1],D=[AOOl01aHA1 a]
C H A C = D C^HAC=D CHAC=D

It should be noted that because BBB is a block matrix, so the identity matrixIII also score blocks, and the block method should be the same as thatof BBB agrees.

3. Positive definiteness of application-judgment matrix

Look at an engineering matrix (may also come from advanced algebra) topic:

Find a block diagonal matrix, which is consistent with the original block symmetric matrix. It is often used to analyze the positive definiteness of a block matrix, such as the following question. insert image description here
Answer:
insert image description here

insert image description here
Use EE in the answer to the above questionE represents the unit matrix, and IIis used in the previous articleI need to pay attention.

The first step in this question is to find the matrix of the conjugate contract with the original matrix, which is to use the method of this article.
In addition, if the two matrices are conjugated, their positive definiteness will remain consistent, so a new matrix that is conjugated to the original matrix can be obtained through elementary transformations, and the positive definiteness of this new matrix is ​​easy to prove, then the original The positive definiteness of the matrix is ​​then proved.

Guess you like

Origin blog.csdn.net/weixin_44286126/article/details/131440632