还记得小学时背过的“九九乘法表”吗?那是刚刚学习乘法的小学生的一道难关。我们小时候用过的铁制文具盒盖子里面就印着“乘法表”,像一个上三角形矩阵。
有一个笑话:上小学二年级的明明回家告诉妈妈:“今天数学考试有一道题问3×7=?3\times 7=?3×7=?,我不会做。快下课时我不管三七二十一,写了个28上去!”这个笑话说明光会背公式是没有用的!
言归正传。大一的同学学习矩阵的乘法是也会觉得挺抽象的。所以,今天我们就来聊聊“矩阵乘法”这点事儿!
符号
AAA为nnn阶方阵.
ei=(0,⋯ ,0,1,0,⋯ ,0)Te_i=(0,\cdots,0,1,0,\cdots,0)^Tei=(0,⋯,0,1,0,⋯,0)T, 其中第iii个分量等于111.
x=(x1,x2,⋯ ,xn)Tx=(x_1,x_2,\cdots,x_n)^Tx=(x1,x2,⋯,xn)T.
小写希腊字母如α\alphaα常常表示nnn维列向量,而αT\alpha^TαT则表示nnn维行向量.
rowirow_irowi指矩阵AAA的第iii行,colicol_icoli指矩阵AAA的第iii列.
我们仅考虑实数域。
###乘法口诀
矩阵ABABAB来相乘,左右首先要分清;
乘法没有消去律,左右因子不能去;乘法没有消去律,左右因子不能去;乘法没有消去律,左右因子不能去;
矩阵乘法真有趣,且听我来说详细:矩阵乘法真有趣,且听我来说详细:矩阵乘法真有趣,且听我来说详细:
行乘以列得实数,效果等于做内积;行乘以列得实数,效果等于做内积;行乘以列得实数,效果等于做内积;
非零列乘非零行,积乃方T阵秩为1;非零列乘非零行,积乃方^T阵秩为1;非零列乘非零行,积乃方T阵秩为1;
AAA左乘以列向量,等于AAA列作组合;
AAA乘eie_iei很容易,直将iii列来提取;
行向量左乘以AAA, 等于AAA行作组合;
eie_iei转置把AAA乘,iii行取出便为积;
初等矩阵左右乘,行列变换显神奇。
注释及应用
- 行乘以列得实数,效果等于做内积;
非零列乘非零行,积乃方阵秩为1.
解释:-
αTβ\alpha^T\betaαTβ是一个实数. 且αTβ=βTα\alpha^T\beta =\beta^T\alphaαTβ=βTα.
-
αβT\alpha\beta^TαβT是一个方阵。当α≠0,β≠0\alpha\ne 0, \beta\ne0α̸=0,β̸=0时,αβT\alpha\beta^TαβT是一个秩等于1的方阵.
-
tr(αβT)=αTβtr(\alpha\beta^T)=\alpha^T\betatr(αβT)=αTβ.这里tr(αβT)tr(\alpha\beta^T)tr(αβT)表示方阵αβT\alpha\beta^TαβT的迹,就是它的主对角元的和.
-
例1 设
α=(123),βT=(212),A=αβT,\alpha=\begin{pmatrix}1\\2\\3\end{pmatrix},
\beta^T=(2\quad 1\quad 2),A=\alpha\beta^T,α=⎝⎛123⎠⎞,βT=(212),A=αβT, 计算:An.A^n.An.
解:An=αβTαβT⋯αβTαβT=α(βTα)⋯(βTα)βT=α(βTα)n−1βT=(βTα)n−1αβT=(βTα)n−1AA^n=\alpha\beta^T\alpha\beta^T\cdots\alpha\beta^T\alpha\beta^T=\alpha(\beta^T\alpha) \cdots(\beta^T\alpha)\beta^T =\alpha(\beta^T\alpha)^{n-1}\beta^T=(\beta^T\alpha)^{n-1}\alpha\beta^T=(\beta^T\alpha)^{n-1}AAn=αβTαβT⋯αβTαβT=α(βTα)⋯(βTα)βT=α(βTα)n−1βT=(βTα)n−1αβT=(βTα)n−1A
因为
βTα=1×2+2×1+3×2=10,\beta^T\alpha=1\times2+2\times1+3\times2=10,βTα=1×2+2×1+3×2=10,所以,
An=10n−1A,A^n=10^{n-1}A,An=10n−1A,
其中,
A=(212424636).A=\begin{pmatrix}2&1&2\\4&2&4\\6&3&6\end{pmatrix}.A=⎝⎛246123246⎠⎞.
- AAA左乘以列向量,等于AAA列作组合;
AAA乘eie_iei很容易,直将iii列来提取;
解释:
(1)A(x1x2⋮xn)=x1col1+x2col2+⋯+xncoln.(1)\quad A\begin{pmatrix}x_1\\x_2\\\vdots\\x_n\end{pmatrix}=x_1col_1+x_2col_2+\cdots+x_ncol_n.(1)A⎝⎜⎜⎜⎛x1x2⋮xn⎠⎟⎟⎟⎞=x1col1+x2col2+⋯+xncoln.
(2)A(0⋮1⋮0)=coli. (2)\quad A\begin{pmatrix}0\\\vdots\\1\\\vdots\\0\end{pmatrix}=col_i.(2)A⎝⎜⎜⎜⎜⎜⎜⎛0⋮1⋮0⎠⎟⎟⎟⎟⎟⎟⎞=coli.
例2 计算:
(1)(1212−14022)(111);(1)\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}1\\1\\1\end{pmatrix};(1)⎝⎛1202−12142⎠⎞⎝⎛111⎠⎞;
(2)(1212−14022)(101);(2)\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}1\\0\\1\end{pmatrix};(2)⎝⎛1202−12142⎠⎞⎝⎛101⎠⎞;
(3)(1212−14022)(001);(3)\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}0\\0\\1\end{pmatrix};(3)⎝⎛1202−12142⎠⎞⎝⎛001⎠⎞;
(4)(1212−14022)(110100111)(4)\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}1&1&0\\1&0&0\\1&1&1\end{pmatrix}(4)⎝⎛1202−12142⎠⎞⎝⎛111101001⎠⎞
解:(1)由公式(1),
(1212−14022)(111)=1⋅(120)+1⋅(2−12)+1⋅(142)=(454).\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}1\\1\\1\end{pmatrix}=1\cdot\begin{pmatrix}1\\2\\0\end{pmatrix}+1\cdot\begin{pmatrix}2\\-1\\2\end{pmatrix}+1\cdot\begin{pmatrix}1\\4\\2\end{pmatrix}=\begin{pmatrix}4\\5\\4\end{pmatrix}.⎝⎛1202−12142⎠⎞⎝⎛111⎠⎞=1⋅⎝⎛120⎠⎞+1⋅⎝⎛2−12⎠⎞+1⋅⎝⎛142⎠⎞=⎝⎛454⎠⎞.
(2) 由公式(1),
(1212−14022)(101)=1⋅(120)+0⋅(2−12)+1⋅(142)=(262).\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}1\\0\\1\end{pmatrix}=1\cdot\begin{pmatrix}1\\2\\0\end{pmatrix}+0\cdot\begin{pmatrix}2\\-1\\2\end{pmatrix}+1\cdot\begin{pmatrix}1\\4\\2\end{pmatrix}=\begin{pmatrix}2\\6\\2\end{pmatrix}.⎝⎛1202−12142⎠⎞⎝⎛101⎠⎞=1⋅⎝⎛120⎠⎞+0⋅⎝⎛2−12⎠⎞+1⋅⎝⎛142⎠⎞=⎝⎛262⎠⎞.
(3) 由公式(2),
(1212−14022)(001)=(142).\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}0\\0\\1\end{pmatrix}=\begin{pmatrix}1\\4\\2\end{pmatrix}.⎝⎛1202−12142⎠⎞⎝⎛001⎠⎞=⎝⎛142⎠⎞.
(4) 由(1)-(3)题,用左边的矩阵每次乘以右边矩阵的一列,得,
(1212−14022)(110100111)=(421564422).\begin{pmatrix}1&2&1\\2&-1&4\\0&2&2\end{pmatrix}\begin{pmatrix}1&1&0\\1&0&0\\1&1&1\end{pmatrix}=\begin{pmatrix}4&2&1\\5&6&4\\4&2&2\end{pmatrix}.⎝⎛1202−12142⎠⎞⎝⎛111101001⎠⎞=⎝⎛454262142⎠⎞.
结束的话
掌握了上面的计算矩阵乘法的方法,一次算一列,是不是感觉特棒?
(未完待续)
更多内容,欢迎用微信扫描下图中的二维码,或搜索“大哉数学之为用”,免费关注微信公众号“大哉数学之为用”进行阅读。
如果您觉得本文对您有帮助,欢迎赞赏!您的支持是作者继续下去的动力!