线性代数基础2矩阵

devtools/2024/10/22 18:34:34/

矩阵是什么
矩阵就是二维数组,下面是一个 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 class=矩阵一定是方阵。不然没有对角线!" />

矩阵运算
矩阵加减法
矩阵的加法就是矩阵的对应位置相加,减法也是一样就是对应位置相减。
在这里插入图片描述
数乘在这里插入图片描述
矩阵乘法
矩阵的乘法和一般的乘法是不太一样!
它是把第一个矩阵的每一行,和第二个矩阵的每一列拿过来做内积得到结果。
在这里插入图片描述
在这里插入图片描述
矩阵乘法运算结果
在这里插入图片描述
矩阵转置
转置的操作和向量是一样的,就是把 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 ∗称为<a class=矩阵 A 的伴随矩阵" />
代数余子式
1、代数余子式定义在这里插入图片描述
2、代数余子式计算在这里插入图片描述
计算过程如下:在这里插入图片描述
伴随矩阵性质I,E都表示单位<a class=矩阵" />
伴随矩阵与逆矩阵
在这里插入图片描述
根据上式可得:在这里插入图片描述

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. ]]
'''

http://www.ppmy.cn/devtools/5954.html

相关文章

js生成word

js生成word 下载依赖 npm install html-docx-js引入 import htmlDocx from html-docx-js/dist/html-docx;代码 //参数 html 文件名字 下载完执行回调函数 function html2word (html,fileName,callback){const converted htmlDocx.asBlob(html);const link document.cre…

袁庭新ES系列15节|Elasticsearch客户端基础操作

前言 上一章节我们介绍了搭建Elasticsearch集群相关的知识。那么又该如何来操作Elasticsearch集群呢&#xff1f;在ES官网中提供了各种语言的客户端&#xff0c;我们在项目开发过程中有多种Elasticsearch版本和连接客户端可以选择&#xff0c;那么他们有什么区别&#xff1f;这…

ccf 202104-2 邻域均值

暴力解法 import java.util.Scanner;public class NeighborAvgOld {public static void main(String[] args) {Scanner sc new Scanner(System.in);int n sc.nextInt();int L sc.nextInt();int r sc.nextInt();int t sc.nextInt();sc.nextLine();int[][] arr new int[n]…

第15届蓝桥杯题解

A题 结果&#xff1a;2429042904288 思路很简单 前20个数分别是 20 24 40 48 60 72 80 96 100 120 140 144 160 168 180 192 200 216 220 240 第2 4 6 8 12 ...n个数分别是24的 1倍 2倍 3倍 4倍 6倍 n/2倍 所以第202420242024 个数就是 24的 101210121012倍 B题 答案&am…

koa-session获取不到session踩坑记录

前端react&#xff0c;后端koa&#xff0c;然后用session作了验证码的存储。后端用koa-cors做了跨域。 场景&#xff1a;前端请求验证码&#xff0c;后端session保存验证码&#xff0c;前端发送验证码给后端&#xff0c;后端取出session校验。这里校验的时候拿不到。 查找原因…

网站结构规范对于独立站的重要性

网站结构是指通过优化网站的内部结构和布局&#xff0c;以提高其在搜索引擎结果页面&#xff08;SERP&#xff09;中的排名和可见性的过程。 涉及到网站内部链接、导航菜单、网站url路径、网站地图等细节。 网站结构的优化可以更好地帮助搜索引擎了索引网站的内容&#xff0c…

攻防演练,作为蓝方,centos的服务器 怎么检查文件是不是被修改或者被拷贝

在基于CentOS的服务器上&#xff0c;作为蓝方进行攻防演练时&#xff0c;检查文件是否被修改或拷贝可以采用以下方法&#xff1a; 1. 安装并配置 Auditd 系统 Auditd 是 Linux 下的一个用来审计系统使用情况的工具&#xff0c;可以用来监控对文件的访问、修改和拷贝操作。 安…

MySQL事务

&#x1f4df;作者主页&#xff1a;慢热的陕西人 &#x1f334;专栏链接&#xff1a;MySQL &#x1f4e3;欢迎各位大佬&#x1f44d;点赞&#x1f525;关注&#x1f693;收藏&#xff0c;&#x1f349;留言 本博客主要内容介绍了MySQL中事务的概念&#xff0c;也是MySQL中非常重…