【矩阵论】2. 矩阵分解——正规谱分解

article/2025/9/1 18:22:33

矩阵论
1. 准备知识——复数域上矩阵,Hermite变换)
1.准备知识——复数域上的内积域正交阵
1.准备知识——Hermite阵,二次型,矩阵合同,正定阵,幂0阵,幂等阵,矩阵的秩
2. 矩阵分解——SVD准备知识——奇异值
2. 矩阵分解——SVD
2. 矩阵分解——QR分解
2. 矩阵分解——正定阵分解
2. 矩阵分解——单阵谱分解
2. 矩阵分解——正规分解——正规阵
2. 矩阵分解——正规谱分解
2. 矩阵分解——高低分解
3. 矩阵函数——常见解析函数
3. 矩阵函数——谱公式,幂0与泰勒计算矩阵函数
3. 矩阵函数——矩阵函数求导
4. 矩阵运算——观察法求矩阵特征值特征向量
4. 矩阵运算——张量积
4. 矩阵运算——矩阵拉直
4.矩阵运算——广义逆——加号逆定义性质与特殊矩阵的加号逆
4. 矩阵运算——广义逆——加号逆的计算
4. 矩阵运算——广义逆——加号逆应用
4. 矩阵运算——广义逆——减号逆
5. 线性空间与线性变换——线性空间
5. 线性空间与线性变换——生成子空间
5. 线性空间与线性变换——线性映射与自然基分解,线性变换
6. 正规方程与矩阵方程求解
7. 范数理论——基本概念——向量范数与矩阵范数
7.范数理论——基本概念——矩阵范数生成向量范数&谱范不等式
7. 矩阵理论——算子范数
7.范数理论——范数估计——许尔估计&谱估计
7. 范数理论——非负/正矩阵
8. 常用矩阵总结——秩1矩阵,优阵(单位正交阵),Hermite阵
8. 常用矩阵总结——镜面阵,正定阵
8. 常用矩阵总结——单阵,正规阵,幂0阵,幂等阵,循环阵


矩阵分解可以得到简化的乘积矩阵,可以简化后续的计算与处理度

在这里插入图片描述

2.4.5 正规谱分解

a. 正规分解推导

若 A = A n × n 正规,互异根为 λ 1 , ⋯ , λ k ,则有 Q H A Q = D = ( λ 1 I 1 ⋱ λ k I k ) 其中 Q 为 U 阵, Q H = Q − 1 , I 1 , ⋯ , I k 为单位阵 ( 如 D = ( 2 ( 1 1 ) 3 ( 1 1 ) ) = ( 2 I 1 3 I 2 ) ) 可设 Q − 1 A Q = D = ( λ 1 I 1 0 ⋱ 0 λ k I k ) ( Q 为 U 阵, Q H = Q − 1 ) 写为 D = λ 1 ( I 1 0 ⋱ 0 0 ) + λ 2 ( 0 0 I 2 0 ⋱ ) + ⋯ + λ k ( 0 0 ⋱ 0 I k ) 则令 D 1 = ( I 1 0 ⋱ 0 0 ) , D 2 = ( 0 0 I 2 0 ⋱ ) , ⋯ , D k = ( 0 0 ⋱ 0 I k ) ⇒ Q − 1 A Q = D = λ 1 D 1 + λ 2 D 2 + ⋯ + λ k D k \begin{aligned} &若A=A_{n\times n} 正规,互异根为 \lambda_1,\cdots,\lambda_k,则有Q^HAQ=D=\left( \begin{matrix} \lambda_1I_1&&\\ &\ddots&\\ &&\lambda_kI_k \end{matrix} \right)\\ &其中Q为U阵,Q^H=Q^{-1},I_1,\cdots,I_k为单位阵\\ &\left(如D=\left( \begin{matrix} 2\left( \begin{matrix} 1&\\ &1 \end{matrix} \right)\\ &3\left( \begin{matrix} 1&\\ &1 \end{matrix} \right) \end{matrix} \right)=\left( \begin{matrix} 2I_1&\\ &3I_2 \end{matrix} \right)\right)\\ &可设 Q^{-1}AQ=D=\left( \begin{matrix} \lambda_1I_1&&0\\ &\ddots&\\ 0&&\lambda_kI_k \end{matrix} \right)(Q为U阵,Q^H=Q^{-1})\\ &写为D=\lambda_1\left( \begin{matrix} I_1&&0\\ &\ddots&\\ 0&&0 \end{matrix} \right)+\lambda_2\left( \begin{matrix} 0&&0\\ &I_2&\\ 0&&\ddots \end{matrix} \right)+\cdots+\lambda_k\left( \begin{matrix} 0&&0\\ &\ddots&\\ 0&&I_k \end{matrix} \right)\\ &则令D_1=\left( \begin{matrix} I_1&&0\\ &\ddots&\\ 0&&0 \end{matrix} \right),D_2=\left( \begin{matrix} 0&&0\\ &I_2&\\ 0&&\ddots \end{matrix} \right),\cdots,D_k=\left( \begin{matrix} 0&&0\\ &\ddots&\\ 0&&I_k \end{matrix} \right)\\ &\Rightarrow Q^{-1}AQ=D=\lambda_1D_1+\lambda_2D_2+\cdots+\lambda_kD_k\\ \end{aligned} A=An×n正规,互异根为λ1,,λk,则有QHAQ=D= λ1I1λkIk 其中QU阵,QH=Q1,I1,,Ik为单位阵 D= 2(11)3(11) =(2I13I2) 可设Q1AQ=D= λ1I100λkIk (QU阵,QH=Q1)写为D=λ1 I1000 +λ2 00I20 ++λk 000Ik 则令D1= I1000 ,D2= 00I20 ,,Dk= 000Ik Q1AQ=D=λ1D1+λ2D2++λkDk

则可得出结论:

  1. 和为单位阵: D 1 + D 2 + ⋯ + D k = ( I 1 ⋱ I k ) = I ( 单位阵 ) D_1+D_2+\cdots+D_k=\left( \begin{matrix} I_1&&\\ &\ddots&\\ &&I_k \end{matrix} \right)=I(单位阵) D1+D2++Dk= I1Ik =I(单位阵)
  2. 正交: D 1 D 2 = 0 , ⋯ , D i D j = 0 ( i ≠ j ) D_1D_2=0,\cdots,D_iD_j=0(i\neq j) D1D2=0,,DiDj=0(i=j)
  3. 幂等: D 1 2 = D 1 , ⋯ , D k 2 = D k ,且 D 1 H = D 1 , ⋯ , D k H = D k D_1^2=D_1,\cdots,D_k^2=D_k,且D_1^H=D_1,\cdots,D_k^H=D_k D12=D1,,Dk2=Dk,且D1H=D1,,DkH=Dk

故可等价写为:
Q H A Q = D = λ 1 D 1 + λ 2 D 2 + ⋯ + λ k D k ⇒ A = Q D Q H = λ 1 Q D 1 Q H + λ 2 Q D 2 Q H + ⋯ + λ k Q D k Q H 可令 G 1 = Q D 1 Q H , ⋯ , G k = Q D k Q H ⇒ A = λ 1 G 1 + ⋯ + λ k G k \begin{aligned} &Q^HAQ=D=\lambda_1D_1+\lambda_2D_2+\cdots+\lambda_kD_k\\ &\Rightarrow A=QDQ^H=\lambda_1QD_1Q^H+\lambda_2QD_2Q^H+\cdots+\lambda_kQD_kQ^H\\ &可令G_1=QD_1Q^H,\cdots,G_k=QD_kQ^H\\ &\Rightarrow A=\lambda_1G_1+\cdots+\lambda_kG_k \end{aligned} QHAQ=D=λ1D1+λ2D2++λkDkA=QDQH=λ1QD1QH+λ2QD2QH++λkQDkQH可令G1=QD1QH,,Gk=QDkQHA=λ1G1++λkGk

有类似推论:

  1. G 1 + G 2 + ⋯ + G k = I G_1+G_2+\cdots+G_k=I G1+G2++Gk=I
    ∵ G 1 + G 2 + ⋯ + G k = Q ( D 1 Q − 1 + D 2 + ⋯ + D k ) Q − 1 = Q I Q − 1 = I \because G_1+G_2+\cdots+G_k=Q(D_1Q^{-1}+D_2+\cdots+D_k)Q^{-1}=QIQ^{-1}=I G1+G2++Gk=Q(D1Q1+D2++Dk)Q1=QIQ1=I

  2. G 1 G 2 = 0 , ⋯ , G i G j = 0 ( i ≠ j ) G_1G_2=0,\cdots,G_iG_j=0(i \neq j) G1G2=0,,GiGj=0(i=j)
    ∵ G 1 G 2 = ( Q D 1 Q − 1 ) ( Q D 2 Q − 1 ) = 0 \because G_1G_2=(QD_1Q^{-1})(QD_2Q^{-1})=0 G1G2=(QD1Q1)(QD2Q1)=0

  3. G 1 2 = G 1 , ⋯ , G k 2 = G k G_1^2=G_1,\cdots,G_k^2=G_k G12=G1,,Gk2=Gk,且 G 1 H = G 1 , ⋯ , G k H = G k G_1^H=G_1,\cdots,G_k^H=G_k G1H=G1,,GkH=Gk 都是Hermite阵

b. 正规阵谱分解与谱阵性质

A = A n × n A=A_{n\times n} A=An×n 正规,全体互异根为 λ 1 , ⋯ , λ k \lambda_1,\cdots,\lambda_k λ1,,λk,则有 A = λ 1 G 1 + λ 2 G 2 + ⋯ + λ k G k A=\lambda_1G_1+\lambda_2G_2+\cdots+\lambda_kG_k A=λ1G1+λ2G2++λkGk ,其中 G 1 , ⋯ , G k G_1,\cdots,G_k G1,,Gk A A A 的谱阵

性质

①和为 I : G 1 + G 2 + ⋯ + G k = I ②正交: G 1 G 2 = 0 , ⋯ , G i G j = 0 ( i ≠ j ) ③幂等: G 1 2 = G 1 , ⋯ , G k 2 = G k ( 幂等 ) , 且 G 1 H = G 1 , ⋯ , G k H = G k ④正规阵幂次: A p = λ 1 p G 1 + ⋯ + λ k G k , p = 0 , 1 , 2. ⋯ ⑤正规阵函数: f ( A ) = f ( λ 1 ) G 1 + ⋯ + f ( λ k ) G k , f ( x ) = c 0 + c 1 x 1 + ⋯ + c p x p ∵ f ( A ) = c 0 I + c 1 A + ⋯ + c k A p = c 0 ( G 1 + ⋯ + G k ) + c 1 ( λ 1 G 1 + λ 2 G 2 + ⋯ + λ k G k ) + ⋯ + c p ( λ 1 p G 1 + λ 2 p G 2 + ⋯ + λ p k G k ) = ( c 0 + c 1 λ 1 + c p λ 1 k ) G 1 + ( c 0 + c 1 λ 2 + ⋯ + c p λ 2 p k ) G 2 + ⋯ + ( c 0 + c 1 λ k + ⋯ + c p λ k p ) G k = f ( λ 1 ) G 1 + ⋯ + f ( λ k ) G k ⑥正规阵求法:设 A 正规,全体不同根为 λ 1 , λ 2 , ⋯ , λ k ,则有谱阵公式 G 1 , G 2 , ⋯ , G k G 1 = ( A − λ 2 I ) ⋯ ( A − λ k I ) ( λ 1 − λ 2 ) ⋯ ( λ 1 − λ k ) , G 2 = ( A − λ 1 I ) ( A − λ 3 I ) ⋯ ( A − λ k I ) ( λ 2 − λ 1 ) ( λ 2 − λ 3 ) ⋯ ( λ 2 − λ k ) ⋯ G k = ( A − λ 1 I ) ( A − λ 2 I ) ⋯ ( A − λ k − 1 I ) ( λ k − λ 1 ) ( λ k − λ 2 ) ⋯ ( λ k − λ k − 1 ) ⑦谱阵中列是 A 的特征向量: A G 1 = λ 1 G 1 , A G 2 = λ 2 G 2 , ⋯ , A G k = λ k G k ∵ A G 1 = ( λ 1 G 1 + ⋯ + λ k G k ) G 1 = λ 1 G 1 \begin{aligned} &①和为I:G_1+G_2+\cdots+G_k=I\\ &②正交:G_1G_2=0,\cdots,G_iG_j=0(i\neq j)\\ &③幂等: G_1^{2}=G_1,\cdots,G_k^{2}=G_k(幂等),且G_1^H=G_1,\cdots ,G_k^H=G_k\\ &④正规阵幂次:A^p=\lambda_1^pG_1+\cdots+\lambda_kG_k,p=0,1,2.\cdots\\ &⑤正规阵函数:f(A)=f(\lambda_1)G_1+\cdots+f(\lambda_k)G_k,f(x)=c_0+c_1x_1+\cdots+c_px^p\\ & \quad \because f(A)=c_0I+c_1A+\cdots+c_kA^p\\ & \quad \quad =c_0(G_1+\cdots+G_k)+c_1(\lambda_1G_1+\lambda_2G_2+\cdots+\lambda_kG_k)+\\ & \quad \quad \quad \cdots+c_p(\lambda_1^pG_1+\lambda_2^pG_2+\cdots+\lambda_p^kG_k)\\ & \quad \quad =(c_0+c_1\lambda_1+c_p\lambda_1^k)G_1+(c_0+c_1\lambda_2+\cdots+c_p\lambda_2^pk)G_2+\\ & \quad \quad \quad \cdots+(c_0+c_1\lambda_k+\cdots+c_p\lambda_k^p)G_k\\ & \quad \quad =f(\lambda_1)G_1+\cdots+f(\lambda_k)G_k\\ &⑥正规阵求法:设A正规,全体不同根为 \lambda_1,\lambda_2,\cdots,\lambda_k,则有谱阵公式G_1,G_2,\cdots,G_k\\ &\quad G_1=\frac{(A-\lambda_2I)\cdots(A-\lambda_kI)}{(\lambda_1-\lambda_2)\cdots(\lambda_1-\lambda_k)},G_2=\frac{(A-\lambda_1I)(A-\lambda_3I)\cdots(A-\lambda_kI)}{(\lambda_2-\lambda_1)(\lambda_2-\lambda_3)\cdots(\lambda_2-\lambda_k)}\\ &\quad\cdots\\ &\quad G_k=\frac{(A-\lambda_1I)(A-\lambda_2I)\cdots(A-\lambda_{k-1}I)}{(\lambda_k-\lambda_1)(\lambda_k-\lambda_2)\cdots(\lambda_k-\lambda_{k-1})}\\ &⑦谱阵中列是A的特征向量:AG_1=\lambda_1G_1,AG_2=\lambda_2G_2,\cdots,AG_k=\lambda_kG_k\\ &\quad \because AG_1=(\lambda_1G_1+\cdots+\lambda_kG_k)G_1=\lambda_1G_1 \end{aligned} 和为IG1+G2++Gk=I正交:G1G2=0,,GiGj=0(i=j)幂等:G12=G1,,Gk2=Gk(幂等),G1H=G1,,GkH=Gk正规阵幂次:Ap=λ1pG1++λkGk,p=0,1,2.正规阵函数:f(A)=f(λ1)G1++f(λk)Gk,f(x)=c0+c1x1++cpxpf(A)=c0I+c1A++ckAp=c0(G1++Gk)+c1(λ1G1+λ2G2++λkGk)++cp(λ1pG1+λ2pG2++λpkGk)=(c0+c1λ1+cpλ1k)G1+(c0+c1λ2++cpλ2pk)G2++(c0+c1λk++cpλkp)Gk=f(λ1)G1++f(λk)Gk正规阵求法:设A正规,全体不同根为λ1,λ2,,λk,则有谱阵公式G1,G2,,GkG1=(λ1λ2)(λ1λk)(Aλ2I)(AλkI),G2=(λ2λ1)(λ2λ3)(λ2λk)(Aλ1I)(Aλ3I)(AλkI)Gk=(λkλ1)(λkλ2)(λkλk1)(Aλ1I)(Aλ2I)(Aλk1I)谱阵中列是A的特征向量:AG1=λ1G1,AG2=λ2G2,,AGk=λkGkAG1=(λ1G1++λkGk)G1=λ1G1

谱阵求法

由推论⑤可知,取k个不同的函数f(x)可求出谱阵 G 1 , G 2 , ⋯ . G k G_1,G_2,\cdots.G_k G1,G2,.Gk

证明

在这里插入图片描述

A正规有2个不同根

G 1 = A − λ 2 I λ 1 − λ 2 , G 2 = I − G 1 = A − λ 1 I λ 2 − λ 1 \begin{aligned} G_1=\frac{A-\lambda_2I}{\lambda_1-\lambda_2},G_2=I-G_1=\frac{A-\lambda_1I}{\lambda_2-\lambda_1} \end{aligned} G1=λ1λ2Aλ2I,G2=IG1=λ2λ1Aλ1I

A正规且有3个不同根

G 1 = ( A − λ 2 I ) ( A − λ 3 I ) ( λ 1 − λ 2 ) ( λ 1 − λ 3 ) , G 2 = ( A − λ 1 I ) ( A − λ 3 I ) ( λ 2 − λ 1 ) ( λ 2 − λ 3 ) , G 3 = I − G 1 − G 2 = ( A − λ 1 I ) ( A − λ 2 I ) ( λ 3 − λ 1 ) ( λ 3 − λ 2 ) \begin{aligned} &G_1=\frac{(A-\lambda_2I)(A-\lambda_3I)}{(\lambda_1-\lambda_2)(\lambda_1-\lambda_3)},G_2=\frac{(A-\lambda_1I)(A-\lambda_3I)}{(\lambda_2-\lambda_1)(\lambda_2-\lambda_3)},\\ &G_3=I-G_1-G_2=\frac{(A-\lambda_1I)(A-\lambda_2I)}{(\lambda_3-\lambda_1)(\lambda_3-\lambda_2)} \end{aligned} G1=(λ1λ2)(λ1λ3)(Aλ2I)(Aλ3I),G2=(λ2λ1)(λ2λ3)(Aλ1I)(Aλ3I),G3=IG1G2=(λ3λ1)(λ3λ2)(Aλ1I)(Aλ2I)

正规分解性质的应用例
1

A = ( 1 0 − 2 0 0 0 − 2 0 4 ) ( 实对称正规阵 ) ,求 A 与 f ( A ) 的谱分解,与 A 100 \begin{aligned} A=\left( \begin{matrix} 1&0&-2\\ 0&0&0\\ -2&0&4 \end{matrix} \right)(实对称正规阵),求A与f(A)的谱分解,与A^{100} \end{aligned} A= 102000204 (实对称正规阵),求Af(A)的谱分解,与A100

可知 A 为秩 1 矩阵, r ( A ) = 1 , 特征根为 λ ( A ) = 5 , 0 , 0 , 不同根为 λ 1 = 5 , λ 2 = 0 G 1 = A − λ 2 I λ 1 − λ 2 = 1 5 ( 1 0 − 2 0 0 0 − 2 0 4 ) , G 2 = A − λ 1 I λ 2 − λ 1 = I − G 1 = 1 5 ( 4 0 2 0 5 0 2 0 1 ) 得谱分解: A = λ 1 G 1 + λ 2 G 2 , 且 f ( A ) = f ( λ 1 ) G 1 + f ( λ 2 ) G 2 即 f ( A ) = f ( 5 ) G 1 + f ( 0 ) G 2 令 f ( x ) = x 100 ⇒ A 100 = f ( 5 ) G 1 + f ( 0 ) G 2 = 5 100 G 1 = 5 100 1 5 A = 5 99 A 由于 G 中列向量都是 A 的特征向量,所以 λ 1 = 5 的特征向量 α 1 = ( 1 0 − 2 ) λ 2 = 0 的特征向量 α 2 = ( 2 0 1 ) , α 3 = ( 0 1 0 ) \begin{aligned} &可知A为秩1矩阵,r(A)=1,特征根为\lambda(A)={5,0,0},不同根为\lambda_1=5,\lambda_2=0\\ &G_1=\frac{A-\lambda_2I}{\lambda_1-\lambda_2}=\frac{1}{5}\left( \begin{matrix} 1&0&-2\\ 0&0&0\\ -2&0&4 \end{matrix} \right),G_2=\frac{A-\lambda_1I}{\lambda_2-\lambda_1}=I-G_1=\frac{1}{5}\left( \begin{matrix} 4&0&2\\ 0&5&0\\ 2&0&1 \end{matrix} \right)\\ &得谱分解:A=\lambda_1G_1+\lambda_2G_2,且f(A)=f(\lambda_1)G_1+f(\lambda_2)G_2\\ &即f(A)=f(5)G_1+f(0)G_2\\ &令f(x)=x^{100}\Rightarrow A^{100}=f(5)G_1+f(0)G_2=5^{100}G_1=5^{100}\frac{1}{5}A=5^{99}A\\ &由于G中列向量都是A的特征向量,所以\lambda_1=5的特征向量\alpha_1=\left( \begin{matrix} 1\\0\\-2 \end{matrix} \right)\\ &\lambda_2=0的特征向量\alpha_2=\left( \begin{matrix} 2\\0\\1 \end{matrix} \right),\alpha_3=\left( \begin{matrix} 0\\1\\0 \end{matrix} \right) \end{aligned} 可知A为秩1矩阵,r(A)=1,特征根为λ(A)=5,0,0,不同根为λ1=5,λ2=0G1=λ1λ2Aλ2I=51 102000204 ,G2=λ2λ1Aλ1I=IG1=51 402050201 得谱分解:A=λ1G1+λ2G2,f(A)=f(λ1)G1+f(λ2)G2f(A)=f(5)G1+f(0)G2f(x)=x100A100=f(5)G1+f(0)G2=5100G1=510051A=599A由于G中列向量都是A的特征向量,所以λ1=5的特征向量α1= 102 λ2=0的特征向量α2= 201 ,α3= 010

2

在这里插入图片描述

2.分块法

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-QSnf9MdK-1671548454484)(2.矩阵分解/image-20221031231951042.png)]

c. 新平方根公式

A = λ 1 G 1 + λ 2 G 2 + ⋯ + λ n G n 条件: A 正规,且 A = λ 1 G 1 + ⋯ + λ k G k \begin{aligned} &\sqrt{A}=\sqrt{\lambda_1}G_1+\sqrt{\lambda_2}G_2+\cdots+\sqrt{\lambda_n}G_n\\ &条件:A正规,且A=\lambda_1G_1+\cdots+\lambda_kG_k \end{aligned} A =λ1 G1+λ2 G2++λn Gn条件:A正规,且A=λ1G1++λkGk

证明:

( A ) 2 = ( λ 1 G 1 + λ 2 G 2 + ⋯ + λ k G k ) 2 = λ 1 G 1 2 + λ 2 G 2 2 + ⋯ + λ k G k 2 = λ 1 G 1 + λ 2 G 2 + ⋯ + λ k G k \begin{aligned} (\sqrt{A})^2&=(\sqrt{\lambda_1}G_1+\sqrt{\lambda_2}G_2+\cdots+\sqrt{\lambda_k}G_k)^2\\ &=\lambda_1G_1^2+\lambda_2G_2^2+\cdots+\lambda_kG_k^2\\ &=\lambda_1G_1+\lambda_2G_2+\cdots+\lambda_kG_k \end{aligned} (A )2=(λ1 G1+λ2 G2++λk Gk)2=λ1G12+λ2G22++λkGk2=λ1G1+λ2G2++λkGk

eg

A = ( 5 4 4 5 ) , 求 A A=\left( \begin{matrix} 5&4\\ 4&5 \end{matrix} \right),求\sqrt{A} A=(5445),A


由于 A 为行和相等矩阵, ∴ λ ( A ) = { 9 , t r ( A ) − 9 } = { 9 , 1 } , 令 λ 1 = 9 , λ 2 = 1 ⇒ G 1 = A − λ 2 I λ 1 − λ 2 = 1 2 ( 1 1 1 1 ) , G 2 = A − λ 1 I λ 2 − λ 1 = 1 2 ( 1 − 1 − 1 1 ) ⇒ A = λ 1 G 1 + λ 2 G 2 = 3 G 1 + G 2 = ( 2 1 1 2 ) \begin{aligned} &由于A为行和相等矩阵,\therefore \lambda(A)=\{9,tr(A)-9\}=\{9,1\},令\lambda_1=9,\lambda_2=1\\ &\Rightarrow G_1=\frac{A-\lambda_2I}{\lambda_1-\lambda_2}=\frac{1}{2}\left( \begin{matrix} 1&1\\1&1 \end{matrix} \right),G_2=\frac{A-\lambda_1I}{\lambda_2-\lambda_1}=\frac{1}{2}\left( \begin{matrix} 1&-1\\-1&1 \end{matrix} \right)\\ &\Rightarrow \sqrt{A}=\sqrt{\lambda_1}G_1+\sqrt{\lambda_2}G_2=3G_1+G_2=\left( \begin{matrix} 2&1\\1&2 \end{matrix} \right) \end{aligned} 由于A为行和相等矩阵,λ(A)={9,tr(A)9}={9,1},λ1=9,λ2=1G1=λ1λ2Aλ2I=21(1111),G2=λ2λ1Aλ1I=21(1111)A =λ1 G1+λ2 G2=3G1+G2=(2112)


http://chatgpt.dhexx.cn/article/0WvmSDOb.shtml

相关文章

【线性代数】矩阵的特征值分解(对角化、谱分解)

目录 1 前言2 矩阵的特征值分解2.1 从定义的角度理解2.2 从变换的角度理解(来自参考文献[3]) 3 对角矩阵(补充)3.1 对角矩阵的定义3.2 对角矩阵线性变换的几何意义 4 矩阵对角化5 相似矩阵与特征值6 参考文献 1 前言 矩阵的特征值分解又可以称作矩阵的对…

矩阵的谱分解 (详细推导步骤~~~特征值分解特征向量

所谓矩阵的分解,就是将一个矩阵写成结构比较简单的或性质比较熟悉的另一些矩阵的乘积。矩阵的分解方法有很多种,包括三角分解、QR(正交三角)分解、最大秩分解、奇异值分解和谱分解,所有这些分解在数值代数和最优化问题…

c语言-gotoxy实现先全部输出再做全部输入操作

需要用到的头文件&#xff1a; #include<windows.h> #include <iostream>代码&#xff1a; gotoxy(a,b)光标控制函数 a为行&#xff0c;b为列&#xff0c;坐标原点在左上角向右是行正方向&#xff0c;向下为列正方向 中文符号汉字在列方向为2个空间&#xff0c;英…

C语言 时钟模拟(gotoxy函数的运用)

时钟模拟&#xff0c;运用gotoxy()函数和Sleep()函数。 效果&#xff1a; #include <stdio.h> #include <windows.h> #include <time.h> #define XHOUR 40 //打印小时的起始x坐标&#xff0c;即a&#xff0c;g交点横坐标 #define YHOUR 27 #define HOUR 1 …

一些神奇的小函数(一)——gotoxy篇

一、gotoxy 1.作用&#xff08;1&#xff09;控制输出的位置① 样例 2.实现&#xff08;1&#xff09;c版 1.作用 &#xff08;1&#xff09;控制输出的位置 ① 样例 : 在代码中写上一句gotoxy(1,1)&#xff0c;然后cout<<“噢&#xff01;这个函数真有用&#xff01;…

Matlab sim函数的用法

一、引言 最近用Simulink做仿真的时候&#xff0c;需要从m文件里运行Simulink模型&#xff0c;而且需要在m文件中传递一些参数到Simulink模型&#xff0c;也需要将Simulink模型的运行结果发送回m文件&#xff0c;所以要用到sim函数。 查看了sim函数的help文档&#xff0c;并百…

CUDA10.0官方文档的翻译与学习之硬件实现

背景 在之前的文章中&#xff0c;我翻译了cuda10.0官方文档的前三章&#xff0c;这次就来翻译第四章——硬件实现 英伟达GPU架构通过由多线程的流式多处理器(SM)组成的可变数组编译&#xff0c;当一个主机CPU上的cuda程序调用了一个核网格&#xff0c;网格内的线程块将会被枚…

近距离看GPU计算(3)

在先前文章《近距离看GPU计算&#xff08;2&#xff09;》中&#xff0c;我们谈到现代GPU发展出SIMT(Single Instruction Multiple Thread)的执行结构&#xff0c;硬件线程池的线程们有相对独立的运行上下文&#xff0c;以Warp为单位分发到一组处理单元按SIMD的模式运行。这些W…

CPU线程与超线程技术

文章目录 一、CPU线程与OS线程1. CPU中的thread2. OS中的thread 二、HT/SMT技术1. 定义2. 原理3. 带来的问题 三、SIMT与SIMD1. SIMT2. SIMD3. 对比 一、CPU线程与OS线程 1. CPU中的thread CPU中的线程来自同步多线程&#xff08;SMT&#xff0c;Simultaneous Multi-threadin…

GPU 软硬件基本概念

术语: single instruction, multiple thread (SIMT)&#xff1a; a single instruction is executed on several function units in parallel GPU的硬件结构中与CUDA相关的几个概念&#xff1a;thread block grid warp sp sm streaming processor(sp): 最基本的处理单元&…

GPU异构计算基础知识

CUDA Toolkit Documentation (nvidia.com) host CPU和内存 (host memory)Device GPU和显存 (device memory) SIMT模型与SIMD模型的区别 SIMD(Single Instruction Multi Data&#xff0c;单指令多数据)模型要求同一个向量中的所有元素要在统一的同步组中一起执行&#xff0c;…

GPU硬件架构以及运行机制笔记

本文是对向往大神的文章的一个笔记。 想阅读原文章移步博客园搜索向往 原文章比较长&#xff0c;这是一个精简和自己的一些理解 这篇文章要带着下面的问题去看 1、GPU是如何与CPU协调工作的&#xff1f; 2、GPU也有缓存机制吗&#xff1f;有几层&#xff1f;它们的速度差异多…

4. CUDA编程手册中文版---硬件实现

第四章 硬件实现 更多精彩内容&#xff0c;请扫描下方二维码或者访问https://developer.nvidia.com/zh-cn/developer-program 来加入NVIDIA开发者计划 NVIDIA GPU 架构围绕可扩展的多线程流式多处理器 (SM: Streaming Multiprocessors) 阵列构建。当主机 CPU 上的 CUDA 程序调…

SMIT介绍

System Management Interface Tool(系统管理界面工具) 软件安装与维护&#xff08;Sofeware Installation and Maintenance&#xff09; 软件许可管理&#xff08;Sofeware License Management&#xff09; 版本管理&#xff08;Manage Editions&#xff09; 设备管理&#…

GPU硬件结构和编程模型(源于nvidia的CUDA文档)

GPU的硬件结构 GPU通过一个可扩展的多线程流式多处理器(SMs)构建。一个multiprocessor可以在同一时间处理上百个线程。为了管理这些线程&#xff0c;使用一个特殊的结构SIMT。利用单线程中指令级的并行&#xff0c;以及同步硬件多线程实现的广泛线程级并行性。 SIMT Architec…

实例分割文献阅读笔记(一)SimT

阅读 SimT: Handling Open-set Noise for Domain Adaptive Semantic Segmentation 原作者知乎文章链接&#xff1a;知乎文章 GitHub链接&#xff1a;开源数据 SimT (CVPR22)&#xff1a;为了解决域自适应&#xff08;包含UDA和SFDA&#xff09;任务中目标域数据伪标签中存在…

第三章 SIMT 内核:指令和寄存器数据流

在本章和下一章中&#xff0c;我们将研究现代 GPU 的架构和微架构。 我们将对 GPU 架构的讨论分为两部分&#xff1a;(1) 在本章中研究实现计算部分的 SIMT 内核&#xff0c;然后 (2) 在下一章中研究内存系统。 在其传统的图形渲染角色中&#xff0c;GPU 访问数据集&#xff0…

从GPU编程到SIMT核心

本文转自&#xff1a;从GPU编程到SIMT核心 - 知乎 (zhihu.com) 1、前言&本文重点 在 GPGPU 显得愈发重要的今天&#xff0c;仅凭 nVidia, AMD 提供的编程接口来了解 GPU 未免显得太单薄了些。时至今日&#xff0c; GPU 内部如何执行一条指令的对程序员来说依然是透明的、…

并行计算范式-SIMD vs SIMT vs SMT: What’s the Difference Between Parallel Processing Models?

Modern processor architectures utilize various execution models. Out of these, two are most popular: SIMD (Single Instruction Multiple Data) and SIMT (Single Instruction Multiple Threads). There’s also SMT (Simultaneous Multithreading), but that’s someth…

SIMD<SIMT<SMT: NVIDIA GPU的并行机制

原文出处&#xff1a; SIMD < SIMT < SMT: parallelism in NVIDIA GPUs 目录 1、概述 1.1、SIMD 2、SIMD vs SIMT 2.1 单指令、多套寄存器组 2.2 单指令、多个数据访问单元 2.3 单指令、多种运算逻辑路径 3、SIMD vs SIMT 3.1 GPU通过多thread来实现高thro…