线性代数基础2矩阵

矩阵是什么
矩阵就是二维数组，下面是一个 m 乘 n 的矩阵，它有 m 行，n 列，每行每列上面都有元素，每个元素都有行标i 和列标 j， a ij 。简称m × n矩阵，记作：

注意a11的索引是 A[0,0]。
这 m×n 个数称为矩阵 A 的元素，简称为元，数aij
位于矩阵 A 的第 i 行第 j 列，称为矩阵 A 的 (i,j) 元，m×n 矩阵 A 也记作 A mn 。

常见矩阵
方阵
如果 m 等于 n，那就称为方阵

对称矩阵
定义是 a ij 等于 a ji 那么就是对称矩阵，对称矩阵首先是个方阵

单位矩阵
主对角线都是 1，其它位置是 0，这称之为单位矩阵，单位矩阵写为 I，一定是方阵，等同于数字里面的 1。

对角矩阵
对角矩阵，就是主对角线非 0，其它位置是 0。

矩阵运算
矩阵加减法
矩阵的加法就是矩阵的对应位置相加，减法也是一样就是对应位置相减。

数乘
矩阵乘法
矩阵的乘法和一般的乘法是不太一样！
它是把第一个矩阵的每一行，和第二个矩阵的每一列拿过来做内积得到结果。


矩阵乘法运算结果

矩阵转置
转置的操作和向量是一样的，就是把 a ij变成 a ji ，把行和列互换一下

矩阵运算法则
矩阵加减法
满足：分配律、结合律、交换律

矩阵乘法
满足结合律
满足分配律
不满足交换律：不一定相等，甚至 AB 的尺寸和 BA 的尺寸是不同的。

矩阵乘法 + 转置

逆矩阵
逆矩阵定义
矩阵有 AB 乘法，但是没有 A/B 这么一说，只有逆矩阵。
逆矩阵怎么定义的？
假设有个矩阵 A，注意它一定是方阵（必须是方阵），乘以矩阵 B 等于单位矩阵I：
那么我们称这里的 B 为 A 的右逆矩阵，和左逆矩阵。

有个很重要的结论就是，如果这样的 B 存在的话，它的左逆和右逆一定相等，统称为 A 的逆矩阵 A −1 。则：A=B−1 B=A−1 （-1次方）
逆矩阵作用
矩阵求逆有什么用呢？它可以帮助我们解线性方程组，比如 XW=Y 。两边同时乘以 X 的逆：
就可以求解出方程的系数，它发明的目的也是干这样的事情用的。
举例说明：

# 三元一次方程
# 3x + 2y + 4z = 19
# 2x -y + 3z = 9
# x + y - z = 0
import numpy as np
X = np.array([[3,2,4],[2,-1,3],[1,1,-1]])
Y = np.array([19,9,0])

display(X,Y)

# np.linalg.inv表示矩阵求逆
# dot表示矩阵乘法
W = np.linalg.inv(X).dot(Y)
print('求解方程得x,y,z为：',W)
'''
array([[ 3,  2,  4],
       [ 2, -1,  3],
       [ 1,  1, -1]])
array([19,  9,  0])
求解方程得x,y,z为： [1. 2. 3.]
'''

从这里我们也可以看出来单位矩阵像我们乘法里面的 1。
逆矩阵相关公式：

伴随矩阵
伴随矩阵定义
设有一矩阵 A
设 A ij 是矩阵A中元素 a ij 的代数余子式，那么矩阵 A ∗
代数余子式
1、代数余子式定义
2、代数余子式计算
计算过程如下：
伴随矩阵性质
伴随矩阵与逆矩阵

根据上式可得：

import numpy as np
# 声明矩阵
A = np.array([
    [1,2,3],
    [2,2,1],
    [3,4,3]])
A_bs = [] # 伴随矩阵
n = 3 # A方阵的行、列数量
for i in range(n):
    for j in range(n):
        row = [0,1,2] # 行索引
        col = [0,1,2] # 列索引
        row.remove(i) # 去除行
        col.remove(j) # 去除列
        # 代数余子式
        A_ij = A[np.ix_(row,col)]
        A_bs.append(((-1)**(i+j)) * np.linalg.det(A_ij))
A_bs = np.array(A_bs).reshape(3,3).T
print('根据伴随矩阵求逆矩阵：\n',A_bs/np.linalg.det(A))
print('用NumPy模块求逆矩阵 ：\n',np.linalg.inv(A))
'''
根据伴随矩阵求逆矩阵：
 [[ 1.   3.  -2. ]
 [-1.5 -3.   2.5]
 [ 1.   1.  -1. ]]
用NumPy模块求逆矩阵 ：
 [[ 1.   3.  -2. ]
 [-1.5 -3.   2.5]
 [ 1.   1.  -1. ]]
'''