高等代数复习:矩阵的满秩分解

2024-03-30 21:20

本文主要是介绍高等代数复习:矩阵的满秩分解,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

文章目录

  • 矩阵的满秩分解

本篇文章适合个人复习翻阅,不建议新手入门使用

矩阵的满秩分解

定义:矩阵的左逆、右逆
A A A m × n m\times n m×n 矩阵

  1. r ( A ) = n r(A)=n r(A)=n,即 A A A 列满秩,则存在秩为 n n n n × m n\times m n×m 行满秩矩阵 B B B,使得 B A = I n BA=I_n BA=In,矩阵 B B B 称为 A A A 的左逆
  2. r ( A ) = m r(A)=m r(A)=m,即 A A A 行满秩,则存在秩为 m m m n × m n\times m n×m 列满秩矩阵 C C C,使得 A C = I m AC=I_m AC=Im,矩阵 C C C 称为 A A A 的右逆

证明思路
结论1:由相抵标准型理论,存在可逆阵 P , Q P,Q P,Q,使得 P A Q = ( I n O ) ⇒ ( I n ∣ O ) P A Q = I n ⇒ ( I n ∣ O ) P A = Q − 1 ⇒ Q ( I n ∣ O ) P A = I n PAQ=\begin{pmatrix} I_n\\O\\ \end{pmatrix}\Rightarrow (I_n| O)PAQ=I_n\\\Rightarrow (I_n |O)PA=Q^{-1}\Rightarrow Q(I_n| O)PA=I_n PAQ=(InO)(InO)PAQ=In(InO)PA=Q1Q(InO)PA=In结论 2类似

定义:矩阵的满秩分解
A A A m × n m\times n m×n 矩阵,其秩为 r r r
则存在秩为 r r r m × r m\times r m×r 列满秩矩阵 B B B,和秩为 r r r r × n r\times n r×n 行满秩矩阵 C C C,使得 A = B C A=BC A=BC这称为 A A A 的满秩分解

证明思路
考虑 A A A 的相抵标准型,容易得到
A = P ( I r O O O ) Q = P ( I r O ) ( I r O ) Q A=P\begin{pmatrix} I_r&O\\ O&O\\ \end{pmatrix}Q=P\begin{pmatrix} I_r\\O\\ \end{pmatrix}\begin{pmatrix} I_r&O\\ \end{pmatrix}Q A=P(IrOOO)Q=P(IrO)(IrO)Q

命题:满秩分解的唯一性
A A A 有两个满秩分解 A = B 1 C 1 = B 2 C 2 A=B_1C_1=B_2C_2 A=B1C1=B2C2,则存在 r r r 阶非异阵 P P P,使得 B 2 = B 1 P , C 2 = P − 1 C 1 B_2=B_1P,C_2=P^{-1}C_1 B2=B1P,C2=P1C1

证明思路
考虑 B 2 B_2 B2 的左逆 P P P C 2 C_2 C2 的右逆 Q Q Q,注意到
B 2 = B 2 ( C 2 Q ) = ( B 2 C 2 ) Q = ( B 1 C 1 ) Q = B 1 ( C 1 Q ) B_2=B_2(C_2Q)=(B_2C_2)Q=(B_1C_1)Q=B_1(C_1Q) B2=B2(C2Q)=(B2C2)Q=(B1C1)Q=B1(C1Q) C 2 = ( P B 2 ) C 2 = P ( B 2 C 2 ) = P ( B 1 C 1 ) = ( P B 1 ) C 1 C_2=(PB_2)C_2=P(B_2C_2)=P(B_1C_1)=(PB_1)C_1 C2=(PB2)C2=P(B2C2)=P(B1C1)=(PB1)C1 ( P B 1 ) ( C 1 Q ) = P ( B 1 C 1 ) Q = P ( B 2 C 2 ) Q = ( P B 2 ) ( C 2 Q ) = I r (PB_1)(C_1Q)=P(B_1C_1)Q=P(B_2C_2)Q=(PB_2)(C_2Q)=I_r (PB1)(C1Q)=P(B1C1)Q=P(B2C2)Q=(PB2)(C2Q)=Ir

命题:满秩分解的几何意义
下列等价:

  1. A = B C A=BC A=BC 是满秩分解
  2. B B B r r r 个列向量是 A A A n n n 个列向量所张成的线性空间的一组基
  3. C C C r r r 个行向量是 A A A m m m 个行向量所张成的线性空间的一组基

证明思路
A = B C A=BC A=BC 写成对应线性方程组的形式即可看出

命题:幂等矩阵关于满秩分解的刻画
n n n 阶方阵 A A A,其秩为 r r r,则 A A A 是幂等矩阵当且仅当存在秩为 r r r n × r n\times r n×r 矩阵 S S S 和秩为 r r r r × n r\times n r×n 矩阵 T T T,使得 A = S T , T S = I r A=ST,TS=I_r A=ST,TS=Ir

证明思路
充分性显然,必要性:考虑 A A A 的相抵标准型
A = P ( I r O O O ) Q A=P\begin{pmatrix} I_r&O\\ O&O\\ \end{pmatrix}Q A=P(IrOOO)Q A 2 = A A^2=A A2=A 得到
( I r O O O ) = ( I r O O O ) Q P ( I r O O O ) \begin{pmatrix} I_r&O\\ O&O\\ \end{pmatrix}=\begin{pmatrix} I_r&O\\ O&O\\ \end{pmatrix}QP\begin{pmatrix} I_r&O\\ O&O\\ \end{pmatrix} (IrOOO)=(IrOOO)QP(IrOOO)则只需令
S = P ( I r O O O ) ( I r O ) , T = ( I r O ) ( I r O O O ) Q S=P\begin{pmatrix} I_r&O\\ O&O\\ \end{pmatrix}\begin{pmatrix} I_r\\O\\ \end{pmatrix},T=\begin{pmatrix} I_r&O\\ \end{pmatrix}\begin{pmatrix} I_r&O\\ O&O\\ \end{pmatrix}Q S=P(IrOOO)(IrO),T=(IrO)(IrOOO)Q

推论
A A A n n n 阶幂等矩阵,则 t r ( A ) = r ( A ) tr(A)=r(A) tr(A)=r(A)

证明
t r ( A ) = t r ( S T ) = t r ( T S ) = t r ( I r ) = r = r ( A ) tr(A)=tr(ST)=tr(TS)=tr(I_r)=r=r(A) tr(A)=tr(ST)=tr(TS)=tr(Ir)=r=r(A)

参考书:《高等代数学》谢启鸿 姚慕生 吴泉水 编著

这篇关于高等代数复习:矩阵的满秩分解的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/862606

相关文章

hdu 4565 推倒公式+矩阵快速幂

题意 求下式的值: Sn=⌈ (a+b√)n⌉%m S_n = \lceil\ (a + \sqrt{b}) ^ n \rceil\% m 其中: 0<a,m<215 0< a, m < 2^{15} 0<b,n<231 0 < b, n < 2^{31} (a−1)2<b<a2 (a-1)^2< b < a^2 解析 令: An=(a+b√)n A_n = (a +

hdu 6198 dfs枚举找规律+矩阵乘法

number number number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description We define a sequence  F : ⋅   F0=0,F1=1 ; ⋅   Fn=Fn

【408数据结构】散列 (哈希)知识点集合复习考点题目

苏泽  “弃工从研”的路上很孤独,于是我记下了些许笔记相伴,希望能够帮助到大家    知识点 1. 散列查找 散列查找是一种高效的查找方法,它通过散列函数将关键字映射到数组的一个位置,从而实现快速查找。这种方法的时间复杂度平均为(

【高等代数笔记】线性空间(一到四)

3. 线性空间 令 K n : = { ( a 1 , a 2 , . . . , a n ) ∣ a i ∈ K , i = 1 , 2 , . . . , n } \textbf{K}^{n}:=\{(a_{1},a_{2},...,a_{n})|a_{i}\in\textbf{K},i=1,2,...,n\} Kn:={(a1​,a2​,...,an​)∣ai​∈K,i=1,2,...,n

线性代数|机器学习-P35距离矩阵和普鲁克问题

文章目录 1. 距离矩阵2. 正交普鲁克问题3. 实例说明 1. 距离矩阵 假设有三个点 x 1 , x 2 , x 3 x_1,x_2,x_3 x1​,x2​,x3​,三个点距离如下: ∣ ∣ x 1 − x 2 ∣ ∣ 2 = 1 , ∣ ∣ x 2 − x 3 ∣ ∣ 2 = 1 , ∣ ∣ x 1 − x 3 ∣ ∣ 2 = 6 \begin{equation} ||x

【线性代数】正定矩阵,二次型函数

本文主要介绍正定矩阵,二次型函数,及其相关的解析证明过程和各个过程的可视化几何解释(深蓝色字体)。 非常喜欢清华大学张颢老师说过的一段话:如果你不能用可视化的方式看到事情的结果,那么你就很难对这个事情有认知,认知就是直觉,解析的东西可以让你理解,但未必能让你形成直觉,因为他太反直觉了。 正定矩阵 定义 给定一个大小为 n×n 的实对称矩阵 A ,若对于任意长度为 n 的非零向量 ,有 恒成

计算机基础知识复习9.6

点对点链路:两个相邻节点通过一个链路相连,没有第三者 应用:PPP协议,常用于广域网 广播式链路:所有主机共享通信介质 应用:早期的总线以太网,无线局域网,常用于局域网 典型拓扑结构:总线型 星型(逻辑总线型) 介质访问控制  静态划分信道 信道划分介质访问控制 频分多路复用FDM 时分多路复用TDM 波分多路复用WDM 码分多路复用CDM 动态分配信道 轮询访问介质访问控

python科学计算:NumPy 线性代数与矩阵操作

1 NumPy 中的矩阵与数组 在 NumPy 中,矩阵实际上是一种特殊的二维数组,因此几乎所有数组的操作都可以应用到矩阵上。不过,矩阵运算与一般的数组运算存在一定的区别,尤其是在点积、乘法等操作中。 1.1 创建矩阵 矩阵可以通过 NumPy 的 array() 函数创建。矩阵的形状可以通过 shape 属性来访问。 import numpy as np# 创建一个 2x3 矩阵mat

特征值分解(EVD)和奇异值分解(SVD)—应用于图片压缩

特征值分解(EVD)和奇异值分解(SVD)—应用于图片压缩 目录 前言 一、特征值分解 二、应用特征值分解对图片进行压缩 三、矩阵的奇异值分解 四、应用奇异值分解对图片进行压缩 五、MATLAB仿真代码 前言         学习了特征值分解和奇异值分解相关知识,发现其可以用于图片压缩,但网上没有找到相应代码,本文在学习了之后编写出了图片压缩的代码,发现奇异值分

【UVA】10003-Cutting Sticks(动态规划、矩阵链乘)

一道动态规划题,不过似乎可以用回溯水过去,回溯的话效率很烂的。 13988658 10003 Cutting Sticks Accepted C++ 1.882 2014-08-04 09:26:49 AC代码: #include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include