7四个基本子空间

1.1四种空间的定义AA

m×nm×n的矩阵。

零空间N(A):n 维向量,是Ax=0Ax=0 的解,所以N(A)在RnRn里。

列空间C(A): 列向量是m维的,所以C(A)在RmRm里。

行空间C(AT)C(AT):A的行的所有线性组合,即A转置的列的线性组合(因为我们不习惯处理行向量C(AT)C(AT)RnRn里。

A的左零空间NATN(AT):A转置的零空间N(AT)N(AT)RmRm里面。

1.2四种空间的维数四种空间的维数是多少?

列空间C(A)AA的主列就是列空间的一组基,CA=RankA=r维数就是秩的大小

行空间CATC(AT) dimCAT=RankA=rdim(C(AT))=Rank(A)=r,有一个重要的性质:行空间和列空间维数相同,都等于秩的大小

零空间N(A):一组基就是一组特殊解,rr 是主变量的个数,nr 是自由变量的个数,零空间的维数等于nrn−r,即dim(N(A))=nrdim(N(A))=n−r

左零空间N(AT)N(AT):矩阵ATATmm列,而其秩为rr,因此AT自由列数目为mrm−r,所以dimNAT=mrdim(N(AT))=m−r

行空间和零空间在RnRn里,他们的维数加起来等于nn,列空间和左零空间在RmRm里,他们的维数加起来等于mm 。如下图所示:

这里写图片描述

2.正交(正交)向量和正交子空间(orthogonal subcpaces )2个正交向量,它们的点积为零。如有2个正交向量

vvw,那么 vTw=0vTw=0. 知道了什么是正交向量,下面就是正交子空间的定义:

Subspace S is orthogonal to subspace T means: every vector in S is orthogonal to every vector in T

下图中有2个子空间,即2个平面V和W,这个2个子空间就不是正交的,虽然2个平面是垂直的,但是你会看到它们之间有个相交的线,而这条直线中的向量即属于V又属于W,它们平行却不垂直,因此V和W不是正确的。因此一个人和你说有1个向量在2个正交的子空间中,那么这个向量一定是零向量,因为只有它才垂直它本身,非零向量是不可能做到这一点的。

这里写图片描述

前面的图已经展示了四个子空间间的正交关系,这里我只解释一下为什么零空间和行空间是正交的? 实际上只要根据 Ax=0 就能得出这样的结论了。由于 Ax=0,那么矩阵的每个行向量乘以 x 都会等于0,如下图所示。有一个问题,之前说过我们既可以用行来表示线性变化也可以用列来表示线性变换,这里怎么只与行空间正交呢?其实我们在做矩阵化简的时候AA与其行最简形式的列空间是不相同的,这样研究的时候就不能用列空间来看啦。

这里写图片描述

C(AT)为啥叫左零空间呢?

AT×y=0AT×y=0,将等式左右两边都转置,得:yT×A=0TyT×A=0T,如下,所以叫左零空间。

这里写图片描述

但我们一般还是习惯用ATy=0ATy=0,因为希望yy 是列向量。

所以,求矩阵的左零空间,就试着寻找一个产生零行向量的行组合,求矩阵的零空间,就试着寻找一个产生零列向量的列组合

那么如何求解这个左零空间呢?
考虑我们前面学过的Gauss-Jordan消元法,我们在A矩阵后面加一个单位阵:

这里写图片描述

这样单位矩阵Im,mIm,m记录了AA的行变换过程。也就是说我们必然存在一个矩阵可以通过作用于A来得到它的行最简形式,那么这个矩阵是谁呢?其实就是Em,mEm,m,因为Em,mIm,m=Em,mEm,mIm,m=Em,m。所以来看下面这个变换的例子:

这里写图片描述

观察RR中第三行为零行,相应的E中的第三行也为左零空间的一组基[1,0,1][−1,0,1]。即

这里写图片描述

### 正交矩阵的定义 正交矩阵是一种特殊的方阵,其主要特性在于它的列向量(或行向量)形成一组标准正交基。具体来说,如果 \( Q \) 是一个 \( n \times n \) 的实矩阵,则 \( Q \) 被称为正交矩阵,当且仅当下述条件成立: \[ Q^TQ = QQ^T = I_n \] 其中,\( Q^T \) 表示矩阵 \( Q \) 的转置,而 \( I_n \) 则代表 \( n \times n \) 的单位矩阵[^1]。 这意味着对于任意两个不同的列向量 \( q_i, q_j \in \mathbb{R}^n \),有 \( q_i^Tq_j = 0 \) (即两向量相互垂直),而对于同一个列向量 \( q_k \),则有 \( q_k^Tq_k = 1 \) (即该向量是一个单位向量)。因此,正交矩阵可以看作是由一组互相垂直且长度为一的标准正交基组成的矩阵。 ### 学性质 #### 性质 1:行列式的绝对值为 1 正交矩阵的一个重要性质是其行列式的模恒等于 1。换句话说,若 \( Q \) 是一个正交矩阵,则必有: \[ |\det(Q)| = 1 \] 这表明正交矩阵要么保持空间定向不变(对应于行列式为 +1 的情况),要么反转空间定向(对应于行列式为 -1 的情况)。前者被称为旋转矩阵,后者可能涉及镜像反射操作[^2]。 #### 性质 2:逆矩阵为其转置 由于满足关系 \( Q^TQ = I_n \),所以正交矩阵的逆非常容易计算,它实际上就是自己的转置形式: \[ Q^{-1} = Q^T \] 这一简单的关系使得在实际应用中处理正交矩阵变得极为便利[^4]。 #### 性质 3:特征值具有特殊约束 正交矩阵的所有特征值均位于复平面上的单位圆周上,也就是说这些特征值的模皆为 1。更进一步地,在实域内讨论时,除了 ±1 外不会有其他任何实值作为特征根;而在复范围内考虑的话,则可能出现成对出现的形式如 e^(iθ)[^2]。 #### 性质 4:构成特定类型的群结构 全体 n×n 实正交矩阵组成了所谓的“正交群”,记做 O(n,R) 或者简称 O(n) 。在这个群里包含了所有能够通过刚体运动实现的空间变换所对应的矩阵集合 —— 这些变化既包括单纯的转动也允许存在翻折成分的存在。特别值得注意的是子集 SO(n), 它由那些决定因子正好取正值 (+1) 的成员组成,代表着纯粹意义上的三欧几里得几何里的绕轴自旋动作[^2]。 ```python import numpy as np def check_orthogonality(matrix): """ Check whether the given square matrix is orthogonal. Parameters: matrix (numpy.ndarray): Input matrix to be checked. Returns: bool: True if input matrix is an orthogonal one; False otherwise. """ identity_matrix = np.eye(len(matrix)) product_with_transpose = np.dot(matrix.T, matrix) return np.allclose(identity_matrix, product_with_transpose) # Example usage of function with a sample orthogonal and non-orthogonal matrices sample_ortho_mat = np.array([[np.cos(np.pi/4), -np.sin(np.pi/4)], [np.sin(np.pi/4), np.cos(np.pi/4)]]) print(f"Is this matrix orthogonal? {check_orthogonality(sample_ortho_mat)}") non_ortho_sample = np.random.rand(2, 2) print(f"And what about another random matrix? {check_orthogonality(non_ortho_sample)}") ```
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值