线性代数:矩阵

1.矩阵加法

矩阵相加就是把对应位置上的项相加 。

前提:相加的两个矩阵需要有相同的m\times n

例如:
\begin{bmatrix} 1&0\\ 2&5\\ 3&1 \end{bmatrix}+ \begin{bmatrix} 4&0.5\\ 2&5\\ 0&1 \end{bmatrix} = \begin{bmatrix} 5&0.5\\ 4&10\\ 3&2 \end{bmatrix}
矩阵相加的结果还是一个m\times n的矩阵

python实现:

import numpy as np
A = np.array([[1,0],[2,5],[3,1]])
B = np.array([[4,0.5],[2,5],[0,1]])
print(A)
print(B)
# 矩阵加法
C= A + B
print(C)

或者

import numpy as np
A = np.matrix([[1,0],[2,5],[3,1]])
B = np.matrix([[4,0.5],[2,5],[0,1]])
print(A)
print(B)
# 矩阵加法
C= A + B
print(C)

2. 标量乘法(Scalar Multiplication)

矩阵中的每个元素都与标量(一个单独的数)相乘:
3\times \begin{bmatrix} 1&0\\ 2&5\\ 3&1 \end{bmatrix} = \begin{bmatrix} 3&0\\ 6&15\\ 9&3 \end{bmatrix}

3. 向量乘法(Vector Multiplication)


m\times n
的矩阵乘以
n\times 1
的向量,结果为
m
维的向量。
y_i
的值等于A矩阵中第
i
行的元素与对应的x向量中的第
i
个元素相乘后加起来。
例如:
\begin{bmatrix} 1&3\\ 4&0\\ 2&1 \end{bmatrix} \begin{bmatrix} 1\\ 5 \end{bmatrix}

运算过程:
1\times1+3\times5=16

4\times1 + 0\times5 = 4

2\times1+1\times5=7

所以:
\begin{bmatrix} 1&3\\ 4&0\\ 2&1 \end{bmatrix} \begin{bmatrix} 1\\ 5 \end{bmatrix} = \begin{bmatrix} 16\\ 4\\ 7 \end{bmatrix}

python实现

import numpy as np
A = np.matrix([[1,3],[4,0],[2,1]])
b = np.matrix([[1],[5]])   # 或者:b = np.matrix([1,5]).T
C = A * b
print(A)
print(b)
print(C)

使用array创建的矩阵是不能用向量乘法的。可以用dot函数

import numpy as np
A = np.array([[1,3],[4,0],[2,1]])
b = np.array([[1],[5]])
# 向量乘法使用dot
c = np.dot(A,b)
print(A)
print(b)
print(c)

4. 矩阵乘法


m\times n
矩阵乘以
n\times o
的矩阵,变成
m\times o
的矩阵。矩阵C的第
i
列第
j
行的元素是B中第
i
列元素分别与A中第
j
行相乘后加起来的结果。

python实现:

import numpy as np
A = np.matrix([[1,3,2],[4,0,1]])
B = np.matrix([[1,3],[0,1],[5,2]])
C = A * B
print(A)
print(B)
print(C)

或者

import numpy as np
A = np.array([[1,3,2],[4,0,1]])
B = np.array([[1,3],[0,1],[5,2]])
# 使用dot实现矩阵乘法
C = np.dot(A,B)
print(A)
print(B)
print(C)

5. 矩阵乘法的性质

矩阵乘法:

  • 不满足交换律(commutative):A\times B \neq B\times A
  • 满足结合律(associative):A\times (B\times C)=(A\times B)\times C

矩阵的乘法中,有一种矩阵起着特殊的作用,如同数的乘法中的1,这种矩阵被称为单位矩阵。它是个方阵(一般用IE表示),从左上角到右下角的对角线(称为主对角线)上的元素均为1。除此以外全都为0。

根据单位矩阵的特点,任何矩阵与单位矩阵相乘都等于本身。单位矩阵需要满足一下条件:

  • m\times n矩阵与单位矩阵相乘时,单位矩阵需要是n\times n,如下:
    \begin{bmatrix} 1&2&3\\ 4&5&6 \end{bmatrix}\times \begin{bmatrix} 1&0&0\\ 0&1&0\\ 0&0&1 \end{bmatrix}= \begin{bmatrix} 1&2&3\\ 4&5&6 \end{bmatrix}
    python 实现:
import numpy as np
I = np.array(
    [[1,0,0],
     [0,1,0],
     [0,0,1]]
)
A=np.array(
    [[1,2,3],
    [4,5,6]]
)
C = np.dot(A,I)
print(C)
  • 当单位矩阵与m\times n矩阵相乘时,单位矩阵需要是m\times m
    \begin{bmatrix} 1&0\\ 0&1 \end{bmatrix}\times \begin{bmatrix} 1&2&3\\ 4&5&6 \end{bmatrix} = \begin{bmatrix} 1&2&3\\ 4&5&6 \end{bmatrix}
    python实现:
import numpy as np
I = np.array(
    [[1,0],
     [0,1]]
)
A=np.array(
    [[1,2,3],
    [4,5,6]]
)
C = np.dot(I,A)
print(C)

6. 逆矩阵(inverse)

如果矩阵A是一个m\times m方阵,如果有逆矩阵A^{-1},则AA^{-1}=A^{-1}A=I

I为单位矩阵

例如:
A = \begin{bmatrix} 1&2\\ 3&4 \end{bmatrix}
则A的逆矩阵
A^{-1}=\begin{bmatrix} -2&1\\ 1.5&-0.5 \end{bmatrix}
则:
A\times A^{-1} = \begin{bmatrix} 1&2\\ 3&4 \end{bmatrix}\times \begin{bmatrix} -2&1\\ 1.5&-0.5 \end{bmatrix}= \begin{bmatrix} 1&0\\ 0&1 \end{bmatrix}

如果矩阵A没有逆矩阵,则A矩阵也叫做奇异矩阵(singular)或退化矩阵(degenerate)

矩阵转置(transpose)

Am\times n阶矩阵,定义A的装置为这样一个n\times m阶矩阵B,满足B=a(j,i),即b(i,j)=a(j,i)B的第ij列元素是A的第j行第i列元素),记A^T=B(有些书记为A'=B)。直观来看,将A的所有元素绕着一条从第1行第1列元素出发的右下方45度的射线作镜面反转,即得到的转置。例:
\begin{bmatrix} a&b\\ c&d\\ e&f \end{bmatrix} ^T= \begin{bmatrix} a&c&e\\ b&d&f \end{bmatrix}

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,504评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,434评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,089评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,378评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,472评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,506评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,519评论 3 413
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,292评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,738评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,022评论 2 329
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,194评论 1 342
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,873评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,536评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,162评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,413评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,075评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,080评论 2 352

推荐阅读更多精彩内容