-
2016-01-14 10:43:42
所有的秩1矩阵(Rank 1 matrix)都可转化为:
Am×n=um×1vTn×1
一个列向量乘以一个行向量。秩1矩阵与线性空间理论
如何描述一个空间——基与维数;
比如如下的一个线性子空间, ∀v∈R4,v1+v2+v3+v4=0 ,所有的 v 构成的子空间
S 。它的基和维数分别是多少?将 v1+v2+v3+v4=0 ,改造成方程组的形式:
Av=0⇓[1,1,1,1]⋅⎡⎣⎢⎢⎢v1v2v3v4⎤⎦⎥⎥⎥=0秩1矩阵登堂入室,此时的子空间 S 即为
A=[1,1,1,1] 的 零空间(null space)而求解null space of Am×n 的维数的一般公式是:
dimN(A)=n−rank(A)=n−r故此时 S 的维数(也即
A 的核空间的维数)为 4−1 ( rank(A)=1 秩1矩阵嘛)。我们接着给出其基的形式,一个主元(第一个位置),三个自由变量(free variable):
⎡⎣⎢⎢⎢−1100⎤⎦⎥⎥⎥,⎡⎣⎢⎢⎢−1010⎤⎦⎥⎥⎥,⎡⎣⎢⎢⎢−1001⎤⎦⎥⎥⎥更多相关内容 -
线性代数笔记15——矩阵空间和秩1矩阵
2018-10-10 23:48:01矩阵空间 矩阵空间是对向量空间的扩展,因为矩阵的本质是向量,所以与向量空间类似,也存在矩阵空间。 在向量空间中,任意两个向量的加法和数乘仍然在该空间内。类似的,所有固定大小的矩阵也组成了矩阵空间,...矩阵空间
矩阵空间是对向量空间的扩展,因为矩阵的本质是向量,所以与向量空间类似,也存在矩阵空间。
在向量空间中,任意两个向量的加法和数乘仍然在该空间内。类似的,所有固定大小的矩阵也组成了矩阵空间,在空间内的任意两个矩阵的加法和数乘也在该空间内。例如,M是所有3×3矩阵构成的空间,空间内的矩阵可以相加,也可以数乘,其结果仍然是3×3矩阵。虽然可以把它们中的两个相乘,但是没人会那么做,因为乘法与向量空间没有关系。
看起来我们需要随时注意乘法的概念,线性代数中的乘法有别于中学时代标量间的乘法。向量与一个标量相乘,被明确定义为数乘;两个相同维度的向量间存在点积和叉乘,其写法上都与九九乘法表中的乘法类似,但是没有定义两个向量的乘法。虽然矩阵间定义了乘法,但存在限制,只有满足“攘外必先安内”原则的两个矩阵才能相乘。线性代数中的概念众多,学了线性代数后,再也不能愉快地做乘法了。
“攘外必先安内”:AB两个矩阵相乘,必须满足A的列数等于B的行数:
如果AB能够相乘,必须满足n = p,看起来n和p正好夹在m和q中间,m和q在外围:
相乘的结果当然是“共同御敌,一致对外”。
以3×3矩阵构成的空间M为例,M一共有9个元素,M的维度是9。M的标准基是:
矩阵的子空间
矩阵的子空间也是对向量子空间的扩展,矩阵子空间需要满足数乘和加法仍处于同一集合内。
矩阵空间M也有一些特殊的子空间,比如将M空间内的两个对称矩阵相加或数乘,仍然得到对称矩阵;上三角矩阵子空间,将空间内的两个矩阵相加或数乘,仍然得到上三角矩阵。仍然以3×3矩阵构成的空间M为例,看看M中的几个子空间。
对称矩阵子空间
设M中的对称子空间是S,6个对称矩阵构成了S的标准基:
S维度是6,表示为:
上三角矩阵子空间
设上三角矩阵子空间是U,U的维度也是6,它的维度表示为:
很明显,上三角矩阵的基是:
对角矩阵子空间
如果一个矩阵既是对称矩阵又是上三角矩阵,则这个矩阵称为对角矩阵(diagonal matrix)。对称矩阵是一个主对角线之外的元素皆为0的矩阵,常写为diag(λ1, λ2,..., λn)。
主对角线上的元素可以为 0 或其他值,对角线上元素相等的对角矩阵称为数量矩阵;对角线上元素全为1的对角矩阵就是单位矩阵。对角矩阵的运算包括和、差运算、数乘运算、同阶对角阵的乘积运算,且这些运算的结果仍为对角阵。
对角矩阵A:
数量矩阵,对称矩阵的对角元线上的元素都相等:
单位矩阵,对称矩阵的对角元线上的元素都为1:
M3×3内的对称矩阵和上三角矩阵的交集是一个对角矩阵,所有对角矩阵也构成一个子空间其维度是3:
很明显,S∩U的标准基是:
S+U矩阵子空间
对于S∪U,它或者在对称矩阵子空间S中,或者在上三角矩阵子空间U中,或者在对角矩阵子空间S∩U中。我们对S∪U不感兴趣,主要是的原因是S∪U并不能构成子空间,可以随便举个例子:
可以看到,U1 + S1是个没什么特点的矩阵,它不属于S∪U,所以不符合子空间的加法封闭性。
现在,我们将S∪U扩大一点,变成S + U,也就是不单独地取S和U中的矩阵,而是取S中的任一矩阵和U中的任一矩阵,将二者相加:
其结果是整个M空间,即:
当然,S + U的维数也是M的维数:
矩阵子空间维数的关系
现在将上面4个矩阵子空间的维数放在一起:
可以看出:
秩1矩阵
秩1矩阵就是秩为1的矩阵,它的行空间和列空间的维度都是1:
更进一步,秩1矩阵可以表示为一列乘以一行的形式:
我们之所以对秩1矩阵感兴趣,是因为可以通过秩1矩阵搭建出任意矩阵,比如秩为4的矩阵,可以通过4个秩1矩阵搭建而成。
如果M是所有5×10矩阵的矩阵空间,那么一个由秩4矩阵组成的子集是否是一个子空间?
当然不是,因为两个矩阵之和的秩不大于两个矩阵的秩之和。设P是M中两个任意秩4矩阵之和,P的秩可能是5,不在秩4矩阵集合内。虽然P是两个秩4矩阵之和,rank(P) ≤4 + 4 = 8,但由于P仍然在M5×10内,而M5×10中的矩阵的秩不会大于5,所以rank(P)的最大值是5。同理,M中由秩1矩阵组成的子集也不是子空间。
秩1矩阵与零空间的关系
假设有一个向量x,x中的分量之和为0:
很明显x满足零空间的条件,它是某个矩阵的零空间,这个矩阵是什么呢?也就是说对于Ax = 0来说,A是什么?x的维数又是什么?
先回答最后一个,根据条件S可以确定,x的一个分量可以由另外三个分量表示:
可见x的主变量有1个,自由变量有3个,维数是3。
再看零空间所属的矩阵,可以很容易判断:
A是秩1矩阵,根据《线性代数14——行空间和左零空间》空间和维数的关系:Am×n的零空间是位于Rn下的 n – r维空间,A的零空间的维数3,x是3维的。
A的列空间和行空间都是1维的:
A的左零空间是零向量。
作者:我是8位的
出处:http://www.cnblogs.com/bigmonkey
本文以学习、研究和分享为主,如需转载,请联系本人,标明作者和出处,非商业用途!
扫描二维码关注公众号“我是8位的”
-
矩阵空间、秩1矩阵
2014-07-06 10:41:53今天要介绍一种新的向量空间,即矩阵空间,之前碰到的所有向量空间,都是n维的实数空间,现在我们将矩阵当成向量,比如说将3*3的矩阵看作向量,这相当于从原来的n维为扩展到n*n维,那么明明是矩阵为什么可以看成是...今天要介绍一种新的向量空间,即矩阵空间,之前碰到的所有向量空间,都是n维的实数空间,现在我们将矩阵当成向量,比如说将3*3的矩阵看作向量,这相当于从原来的n维为扩展到n*n维,那么明明是矩阵为什么可以看成是向量呢?因为矩阵也服从向量空间的运算,向量能相加,矩阵也能相加,向量能数乘,矩阵也可以数乘,向量可以线性组合,矩阵也可以线性组合。所以说矩阵也可以当成向量来生成空间,这个空间就是矩阵空间。比如说所有3*3的矩阵组成的空间,我们记为矩阵空间M,那么M的子空间有:所有上三角矩阵(upper triangular matrix)组成的空间U、所有的对称矩阵(symmetricmatrix)组成的空间S,所有的对角矩阵(diagonal matrix)组成的空间D(对角矩阵是上三角阵和对称阵的交集),因为两个上三角阵相加还是上三角阵,两个对称阵相加还是对称阵,两个对角阵相加还是对角阵,除此之外,我们还可以得出M空间的维数是9,9个基分别为:
S空间的维数是6,6个基分别为:
U空间的维数也是6,基分别为:
通过将U和S进行组合,我们还可以得到其他子空间,比如前面已经提过的对称阵空间S和上三角阵空间U的交集
得到的对角阵空间D,其维数是3维,但是注意并集 组成的集合不是子空间,要将
符号改成“+”,即S+U才构成子空间,如果每次任取S内任一元素加上U内任一元素,我们可以得到所有3*3矩阵,因此这个S+U构成的空间维度是9,这里我们可以发现dim(S)+dim(U)=dim(S交U)+dim(S+U),即S的维度加上U的维度等于它们交的维度加上它们和的维度,以上就是关于矩阵空间的内容,除了向量跟传统的向量有区别,其他都是一样的求法。现在我们再举个比较特殊的例子,假设有一个矩阵空间,它来自于微分方程,该微分方程为d2y/dx2+y=0,其解是y=cos( x) ,y=sin( x),这个微分方程的零空间或解空间就是这两个解的线性组合,也就是y=c1cos(x)+c2sin(x),很铭明显这个向量空间的维数是2,因为它有2个基:cos(x)和sin(x)。举这个例子的目的在于:cos(x)和sin(x)看起来并不像向量,有些像函数,但我们可以称它们为向量,因为其可以做数乘也可以作加法,所以可对它们进行线性组合,这就是线性代数中的向量,基,维数,它们不仅仅局限于我们平常所熟悉的向量形式,也可以用于矩阵甚至函数等,只要它们满足向量的特性。
秩1矩阵
关于矩阵的秩r,目前我们已经知道的是它等于矩阵主元的个数,且r<m,r<n,秩1矩阵是一类性质很好但非常简单的矩阵,因此我们在这里对其作一些简单分析。假设现有秩1矩阵
,经过消元可得知其行空间的基为A的第一行,列空间的基为A的第一列,那么A可以写成两个基相乘的形式,即
,对于所有的秩1矩阵,这个结论都是成立的,即所有的秩1矩阵都能表示为一列基乘以一行基的形式A=UVT,秩1矩阵就像搭其它矩阵的积木,任何矩阵都可以表示为若干个秩1矩阵的组合,例如如果有一个5*17的矩阵,它的秩为4,那么这个矩阵就可分解为4个秩1矩阵的组合。
-
低秩稀疏矩阵优化问题的模型与算法
2020-07-29 19:44:03低秩稀疏矩阵优化问题是一类带有组合性质的非凸非光滑优化问题. 由于零模与秩函数 的重要性和特殊性, 这类 NP-难矩阵优化问题的模型与算法研究在过去〸几年里取得了长足发展。 -
秩一矩阵的分解
2021-05-25 14:20:39a_1 \\ \vdots\\ a_n \end{bmatrix}, \quad b = \begin{bmatrix} b_1 \\ \vdots\\ b_n \end{bmatrix} a=⎣⎢⎡a1⋮an⎦⎥⎤,b=⎣⎢⎡b1⋮bn⎦⎥⎤ aTb=∑i=1naibi,abT=[a1b1⋯a1bn⋯⋯⋯anb1⋯...a = [ a 1 ⋮ a n ] , b = [ b 1 ⋮ b n ] a = \begin{bmatrix} a_1 \\ \vdots\\ a_n \end{bmatrix}, \quad b = \begin{bmatrix} b_1 \\ \vdots\\ b_n \end{bmatrix} a=⎣⎢⎡a1⋮an⎦⎥⎤,b=⎣⎢⎡b1⋮bn⎦⎥⎤
a T b = ∑ i = 1 n a i b i , a b T = [ a 1 b 1 ⋯ a 1 b n ⋯ ⋯ ⋯ a n b 1 ⋯ a n b n ] , a^T b = \sum_{i=1}^{n} a_i b_i, \quad ab^T = \begin{bmatrix} a_1b_1 & \cdots & a_1 b_n \\ \cdots & \cdots & \cdots \\ a_n b_1 & \cdots & a_n b_n \end{bmatrix}, aTb=i=1∑naibi,abT=⎣⎡a1b1⋯anb1⋯⋯⋯a1bn⋯anbn⎦⎤,
a b T ab^T abT的秩为1,其所有行或列都线性相关,只有一个非零特征值: λ ( a b T ) = a T b \lambda(ab^T) = a^Tb λ(abT)=aTb.
列变换
[ a 1 b 1 ⋯ a 1 b n ⋯ ⋯ ⋯ a n b 1 ⋯ a n b n ] [ 1 − b 2 b 1 − b 3 b 1 ⋯ − b n b 1 1 0 ⋯ 0 ⋱ ⋱ ⋮ 1 0 1 ] = [ a 1 b 1 ⋯ 0 ⋯ ⋯ ⋯ a n b 1 ⋯ 0 ] \begin{bmatrix} a_1b_1 & \cdots & a_1 b_n \\ \cdots & \cdots & \cdots \\ a_n b_1 & \cdots & a_n b_n \end{bmatrix} \begin{bmatrix} 1 & -\frac{b_2}{b_1} &-\frac{b_3}{b_1} & \cdots& -\frac{b_n}{b_1} \\ &1& 0 & \cdots & 0 \\ & & \ddots &\ddots&\vdots\\ &&&1& 0\\ &&&&1 \end{bmatrix} = \begin{bmatrix} a_1b_1 & \cdots & 0 \\ \cdots & \cdots & \cdots \\ a_n b_1 & \cdots &0 \end{bmatrix} ⎣⎡a1b1⋯anb1⋯⋯⋯a1bn⋯anbn⎦⎤⎣⎢⎢⎢⎢⎢⎡1−b1b21−b1b30⋱⋯⋯⋱1−b1bn0⋮01⎦⎥⎥⎥⎥⎥⎤=⎣⎡a1b1⋯anb1⋯⋯⋯0⋯0⎦⎤
行变换
[ 1 − a 2 a 1 1 − a 3 a 1 0 1 ⋮ ⋮ ⋱ ⋱ − a n a 1 0 ⋯ 0 1 ] [ a 1 b 1 ⋯ 0 ⋯ ⋯ ⋯ a n b 1 ⋯ 0 ] = [ a 1 b 1 0 ⋯ 0 0 0 ⋯ 0 ⋯ ⋯ ⋯ ⋯ 0 0 ⋯ 0 ] \begin{bmatrix} 1 & & & & \\ -\frac{a_2}{a_1}&1& && \\ -\frac{a_3}{a_1}& 0 & 1 && \\ \vdots &\vdots&\ddots&\ddots&\\ -\frac{a_n}{a_1} &0&\cdots&0&1 \end{bmatrix} \begin{bmatrix} a_1b_1 & \cdots & 0 \\ \cdots & \cdots & \cdots \\ a_n b_1 & \cdots &0 \end{bmatrix} = \begin{bmatrix} a_1b_1 & 0 &\cdots & 0 \\ 0 & 0 & \cdots & 0 \\ \cdots & \cdots & \cdots & \cdots \\ 0 & 0 & \cdots &0 \end{bmatrix} ⎣⎢⎢⎢⎢⎢⎡1−a1a2−a1a3⋮−a1an10⋮01⋱⋯⋱01⎦⎥⎥⎥⎥⎥⎤⎣⎡a1b1⋯anb1⋯⋯⋯0⋯0⎦⎤=⎣⎢⎢⎡a1b10⋯000⋯0⋯⋯⋯⋯00⋯0⎦⎥⎥⎤
记
P a 1 = [ 1 − a 2 a 1 1 − a 3 a 1 0 1 ⋮ ⋮ ⋱ ⋱ − a n a 1 0 ⋯ 0 1 ] , P b 1 = [ 1 − b 2 b 1 − b 3 b 1 ⋯ − b n b 1 1 0 ⋯ 0 ⋱ ⋱ ⋮ 1 0 1 ] P_{a1} = \begin{bmatrix} 1 & & & & \\ -\frac{a_2}{a_1}&1& && \\ -\frac{a_3}{a_1}& 0 & 1 && \\ \vdots &\vdots&\ddots&\ddots&\\ -\frac{a_n}{a_1} &0&\cdots&0&1 \end{bmatrix}, P_{b1} = \begin{bmatrix} 1 & -\frac{b_2}{b_1} &-\frac{b_3}{b_1} & \cdots& -\frac{b_n}{b_1} \\ &1& 0 & \cdots & 0 \\ & & \ddots &\ddots&\vdots\\ &&&1& 0\\ &&&&1 \end{bmatrix} Pa1=⎣⎢⎢⎢⎢⎢⎡1−a1a2−a1a3⋮−a1an10⋮01⋱⋯⋱01⎦⎥⎥⎥⎥⎥⎤,Pb1=⎣⎢⎢⎢⎢⎢⎡1−b1b21−b1b30⋱⋯⋯⋱1−b1bn0⋮01⎦⎥⎥⎥⎥⎥⎤特征分解与对角化
( a b T ) a = a ( b T a ) = ( a T b ) a (ab^T)a = a(b^Ta) = (a^Tb)a (abT)a=a(bTa)=(aTb)a
说明 a a a 是特征值 a T b a^Tb aTb 对应的特征向量,剩下 n − 1 n-1 n−1个特征向量只要满足和 b b b正交即可! -
19 矩阵——矩阵的相抵、相抵标准形、秩1矩阵、矩阵的满秩分解
2021-05-29 21:42:50文章目录相抵标准形相抵标准形的应用秩1矩阵矩阵的满秩分解参考 相抵标准形 命题\large\color{magenta}{\boxed{\color{brown}{命题} }}命题 对矩阵做有限次初等变换不改变矩阵的秩. 推论\large\color{magenta}{\... -
非均匀采样磁共振谱的快速低秩汉克尔矩阵分解重构方法
2021-03-08 20:43:55非均匀采样磁共振谱的快速低秩汉克尔矩阵分解重构方法 -
低秩矩阵恢复算法综述
2018-09-02 11:33:32低秩矩阵恢复算法综述,内容比较易懂,用于图像修复,推荐算法等 -
C#计算矩阵的秩实例分析
2020-09-03 10:58:04主要介绍了C#计算矩阵的秩实现方法,以实例形式较为详细的分析了C#计算矩阵秩的原理与实现技巧,具有一定参考借鉴价值,需要的朋友可以参考下 -
秩一矩阵
2018-12-01 11:17:05我们可以把一个矩阵看做一个向量 比如所有3×3的矩阵组成一个向量空间 那这个向量空间的维度dimM=9 -
秩一矩阵的优良性质
2018-02-24 10:40:10很好的从另一个角度诠释了矩阵乘法。还是很不错的。 ui→,vi→ui→,vi→\vec{u_i},\vec{v_i}都是列向量。 UVT=([u1→,u2→,...,un→])([v1→,v2→,...,vn→]T)=u1→v1→T+u2→v2→T+...,un→vn→TUVT=([u1→,u2→,... -
秩为 1 的矩阵的一些性质
2021-02-24 09:52:55从上面的分析和例题看到,对于秩为1的n阶矩阵,零是其n重或n-1重特征值,如果是n-1重,则非零特征值是矩阵的主对角线元素之和;另外还看到,秩为1的矩阵可以分解为一个非零列向量与另一个非零列向量的转置的乘积,这... -
LDPC编码器(完全或不完全秩奇偶校验矩阵):即使该矩阵不是完全秩矩阵,此程序也允许使用奇偶校验矩阵对...
2021-05-29 06:47:59LDPC 码(非满秩或满秩奇偶校验矩阵)的编码。 codeWord = ldpcEncoding(H,u): 输入变量: H:奇偶校验矩阵。 u:信息位向量。 输出变量: 码字:为信息位向量 u 生成的码字。 -
线性代数 -- 矩阵空间、秩1矩阵、小世界图
2017-08-07 19:48:42矩阵空间介绍一种新的向量空间:矩阵空间。 在这之前我们遇到的都是n维的向量空间, 对于矩阵空间先了解以下几个问题:Q:为什么矩阵也可以看成向量空间呢?A:因为矩阵可以同向量一样进行各种运算, 比如:加法, ... -
线性代数导论11——矩阵空间、秩1矩阵和小世界图
2015-04-16 19:59:00第十一课时:矩阵空间、秩1矩阵和小世界图 矩阵空间 矩阵空间,可看着是新的向量空间,比如3×3的矩阵,它们加法或数乘都停留在3×3矩阵空间,3×3矩阵有一些子空间:3×3对称矩阵的子空间(两个对称... -
P11 矩阵空间 秩1矩阵 小世界 【线性代数】
2022-01-24 21:27:17矩阵空间 秩1矩阵 小世界 -
基于截断核范数的低秩稀疏矩阵恢复新方法
2021-03-03 11:38:43基于截断核范数的低秩稀疏矩阵恢复新方法 -
用于矩阵求秩
2017-05-06 16:25:16矩阵求秩 -
关于矩阵的秩及求解Python求法
2021-08-24 01:23:09关于消元法求解线性方程组可将系数和结果转换为矩阵,并可令B为增广矩阵将A、B通过消元法求解所有的m*n的矩阵经过一系列初等变换,都可以变成如下的形式:r就是最简矩阵当中非零行的... -
C/C++求矩阵的秩
2018-10-08 19:13:07C/C++语言实现矩阵的求秩,原理是高斯消元法实现矩阵求秩。之前我在网上找的C/C++语言实现矩阵的秩的很多程序通不过大量的测试,此程序是本人经过大量的测试通过的。对于方阵和非方阵都能准确的求出其秩,如果发现有... -
【线性代数公开课MIT Linear Algebra】 第十一课 矩阵空间和秩1矩阵
2015-10-26 23:38:07本系列笔记为方便日后自己查阅而写,更多的是个人见解,也算一种学习的复习与总结,望善始善终吧~矩阵空间和之前学习的空间差不多,我们把矩阵当做向量,矩阵空间也是在空间内对一个矩阵进行加法或者scalar后仍然在... -
总结利用秩为1的矩阵相关矩阵的秩的计算问题
2016-11-28 23:12:51总结利用秩为1的矩阵相关矩阵的秩的计算问题@(线性代数)对于一个秩为1的矩阵,常常给定的是一个列向量与自己的转置之积。... -
低秩矩阵分解的概率模型研究
2021-03-12 14:42:43在过去的十年中,出现了各种低秩矩阵分解的概率模型。 低秩矩阵分解与它们相应的概率模型之间最显着的区别是后者将低秩分量视为随机变量。 本文对低秩矩阵分解的概率模型进行了调查。 首先,我们回顾了低秩矩阵分解... -
使用低秩矩阵恢复的高光谱图像恢复
2021-03-16 20:53:02本文介绍了一种基于低秩矩阵恢复(LRMR)的HSI恢复新方法,该方法可以同时去除高斯噪声,脉冲噪声,死线和条纹。 通过按字典顺序将HSI的补丁排序为二维矩阵,可以探索高光谱图像的低秩属性,这表明干净的HSI补丁... -
证明矩阵的秩=行秩=列秩
2021-11-23 22:44:27极大线性无关组:若向量组{α1,α2,…,αn}\{\alpha_1,\alpha_2,\ldots,\alpha_n\}{α1,α2,…,αn}的一个部分向量组{αi1,αi1,…,αir}\{\alpha_{i_1},\alpha_{i_1},\ldots,\alpha_{i_r}\}{αi1,αi1... -
低秩分块矩阵的核近似.pdf
2020-05-13 13:05:22低秩分块矩阵的核近似.pdf -
矩阵空间、秩1矩阵和小世界图-线性代数课时11(MIT Linear Algebra , Gilbert Strang)
2018-10-23 11:01:42这是Strang教授的第十一讲,讲解的内容是矩矩阵空间(一个新的“向量”空间)的一组基,秩1矩阵的特殊性和小世界图(small world graphs),小世界图引出图论与线性代数的关系。 矩阵空间 矩阵空间满足向量空间... -
低秩矩阵分解
2014-12-29 16:00:35低秩矩阵分解代码,inexact alm