企业网站开发韵茵,网站开发 经常要清理缓存,拉米拉网站建设,wordpress 菜单 导航代码线性代数 1.1 线性代数内容介绍1.1.1 线性代数介绍1.1.2 代码实现介绍 1.2 线性代数实现1.2.1 reshape运算1.2.2 转置实现1.2.3 矩阵乘法实现1.2.4 矩阵对应运算1.2.5 逆矩阵实现1.2.6 特征值与特征向量1.2.7 求行列式1.2.8 奇异值分解实现1.2.9 线性方程组求解 1.1 线性代数内… 线性代数 1.1 线性代数内容介绍1.1.1 线性代数介绍1.1.2 代码实现介绍 1.2 线性代数实现1.2.1 reshape运算1.2.2 转置实现1.2.3 矩阵乘法实现1.2.4 矩阵对应运算1.2.5 逆矩阵实现1.2.6 特征值与特征向量1.2.7 求行列式1.2.8 奇异值分解实现1.2.9 线性方程组求解 1.1 线性代数内容介绍
1.1.1 线性代数介绍
线性代数是一门被广泛运用于各工程技术领域的学科。用线性代数的相关概念和结论可以极大地简化数据挖掘中相关公式的推导和表述。线性代数将复杂的问题简单化让我们能够对问题进行高效地数学运算。
线性代数是一个数学工具它不仅提供了有助于操作数组的技术还提供了像向量和矩阵这样的数据结构用来保存数字和规则以便进行加减乘除的运算。
1.1.2 代码实现介绍
numpy是一款基于Python的数值处理模块在处理矩阵数据方面有很强大的功能与优势。因为线性代数的主要内容就是对矩阵的处理所以本章节主要的内容都是基于numpy进行展开。另外也会涉及到方程组求解所以也会用到数学科学库scipy。
1.2 线性代数实现
导入相应库
import numpy as np
import scipy as sp 1.2.1 reshape运算
在数学中并没有 reshape 运算但是在numpy运算库中是一个非常常用的运算用来改变一个张量的维度数和每个维度的大小例如一个10x10的图片在保存时直接保存为一个包含100个元素的序列在读取后就可以使用reshape将其从 1x100 变换为 10x10 。
示例如下
生成一个包含整数0~11的向量
x np.arange(12)
print(x)结果输出
[ 0 1 2 3 4 5 6 7 8 9 10 11]注意如果是直接np.arange(12)输出结果为
array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11])扩展学习 如果是想要以逗号进行隔开方式打印x则可以指定一下
new_x np.array2string(x, separator ,)
print(new_x)输出结果为
[ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11]但是需要注意的是此处打印的是str类型。
查看数组大小
x.shape结果输出
(12,)将x转换成二维矩阵其中矩阵的第一个维度为1
x x.reshape(1,12)
print(x)结果输出 [[ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]]查看数组大小
x.shape结果输出
(1, 12)将x转换3x4的矩阵
x x.reshape(3,4)
print(x)结果输出
[[ 0, 1, 2, 3],[ 4, 5, 6, 7],[ 8, 9, 10, 11]]1.2.2 转置实现
向量和矩阵的转置是交换行列顺序而三维及以上张量的转置就需要指定转换的维度。
生成3*4的矩阵
A np.arange(12).reshape(3,4)
print(A)结果输出
[[ 0, 1, 2, 3],[ 4, 5, 6, 7],[ 8, 9, 10, 11]]转置
A.T结果输出
array([[ 0, 4, 8],
[ 1, 5, 9],
[ 2, 6, 10],
[ 3, 7, 11]])1.2.3 矩阵乘法实现
矩阵乘法记两个矩阵分别为A和B两个矩阵能够相乘的条件为第一个矩阵的 列数 等于第二个矩阵的 行数 。
代码输入
A np.arange(6).reshape(3,2)
print(A)结果输出
[[0 1][2 3][4 5]]B np.arange(6).reshape(2,3)
print(B)结果输出
[[0, 1, 2],[3, 4, 5]]矩阵相乘
np.matmul(A,B)注意不要写成 A.matmul(B)也不能写成 A*B 。 结果输出
array([[ 3, 4, 5],
[ 9, 14, 19],
[15, 24, 33]])1.2.4 矩阵对应运算
元素对应运算针对形状相同矩阵的运算统称包括元素对应相乘、相加等即对两个矩阵相同位置的元素进行加减乘除等运算。
创建矩阵
A np.arange(6).reshape(3,2)矩阵相乘
print(A*A)结果输出
array([[ 0, 1],
[ 4, 9],
[16, 25]])矩阵相加
rint(A A)结果输出
array([[ 0, 2],
[ 4, 6],
[ 8, 10]])1.2.5 逆矩阵实现
只有方阵才有 逆矩阵 逆矩阵实现。
代码输入
A np.arange(4).reshape(2,2)
print(A)结果输出
array([[0, 1],
[2, 3]])求逆矩阵
np.linalg.inv(A)结果输出
array([[-1.5, 0.5],
[ 1. , 0. ]])注意此处是 np.linalg 不是scipy.linalg。可以自行了解如何求逆此处不做过多解释。 提示设A是数域上的一个n阶方阵若在相同数域上存在另一个n阶矩B使得 ABBAE 。 则我们称B是A的逆矩阵而A则被称为可逆矩阵。其中E为单位矩阵。
说明 np.linalg 和 scipy.linalg 在提供线性代数函数方面有些重叠但也存在一些区别。如果您只需要基本的线性代数功能np.linalg足以满足您的需求。如果需要更高级的线性代数功能或特定的分解方法您可能需要查看scipy.linalg中的函数。SciPy是建立在NumPy之上的科学计算库提供了更广泛的数学、科学和工程计算功能。
如果是使用scipy的linalg
代码如下
from scipy import linalg
linalg.inv(A)输出结果为
array([[-1.5, 0.5],[ 1. , 0. ]])1.2.6 特征值与特征向量
当谈论矩阵的 特征值 和 特征向量 时我们首先需要了解 线性变换 和 向量空间 的概念。 线性变换是指将一个向量空间中的向量映射到另一个向量空间中的向量的操作。在二维平面上的旋转和缩放、三维空间中的投影等都是线性变换的例子。 向量空间是指由多个向量组成的集合其中的向量可以进行加法和数乘运算。向量空间可以是二维平面、三维空间或更高维度的空间。
现在我们来解释矩阵的特征值和特征向量 特征值Eigenvalues是一个数值表示线性变换作用后的向量在同一方向上的缩放倍数。当一个向量在经过线性变换后只发生缩放而不改变方向时这个缩放倍数就是特征值。 特征向量Eigenvectors是与特征值相关联的向量。它是表示在线性变换下保持在同一方向上的向量。特征向量不会改变方向而是在线性变换后仅仅以特征值的倍数进行缩放。
具体来说对于一个n维向量空间中的线性变换我们可以表示为一个n×n的矩阵A。如果存在一个非零向量v和一个标量λ使得下式成立
A * v λ * v其中v是 特征向量 λ是 特征值 。这意味着当矩阵A作用在特征向量v上时结果只是将v进行了缩放缩放的比例由特征值λ确定。
特征值和特征向量对于矩阵的理解和分析非常重要。它们提供了关于矩阵在线性变换过程中的行为和性质的信息。通过计算矩阵的特征值和特征向量我们可以了解线性变换的缩放效果和主要方向并在许多应用中提供有用的洞察力如主成分分析、图像处理和振动分析等。
接下来进行求矩阵的 特征值 与 特征向量 并实现可视化。
导入相应库
from scipy.linalg import eig
import numpy as np
import matplotlib.pyplot as plt求特征值与特征向量
A [[1, 2],#生成一个2*2的矩阵[2, 1]]
evals, evecs eig(A) #求A的特征值evals和特征向量(evecs)
evecs evecs[:, 0], evecs[:, 1]plt.subplots() 返回一个Figure实例fig 和一个 AxesSubplot实例ax。fig代表整个图像ax代表坐标轴和画的图。 作图
fig, ax plt.subplots()输出图片如下
让坐标轴经过原点
for spine in [left, bottom]:#让在左下角的坐标轴经过原点ax.spines[spine].set_position(zero)画出网格
ax.grid(alpha0.4)设置坐标轴的范围
xmin, xmax -3, 3
ymin, ymax -3, 3
ax.set(xlim(xmin, xmax), ylim(ymin, ymax))输出结果为
[(-3.0, 3.0), (-3.0, 3.0)]画出 特征向量 。用一个箭头指向要注释的地方再写上一段话的行为叫做annotate。text是输入内容xy:箭头指向xytext文字所处的位置arrowprops通过arrowstyle表明箭头的风格或种类
for v in evecs:ax.annotate(text, xyv, xytext(0, 0),arrowpropsdict(facecolorblue,shrink0,alpha0.6,width0.5))注意 问题原因annotate()的’s’参数自Matplotlib 3.3以后已重命名为’text’不能使用s不然会报错。
画出 特征空间
x np.linspace(xmin, xmax, 3)#在指定的间隔内返回均匀间隔的数字
for v in evecs:a v[1] / v[0] #沿特征向量方向的单位向量ax.plot(x, a * x, r-, lw0.4)# 参数 lw 表示图线的粗细
plt.show() 可视化图像 蓝箭头向量为特征向量两条红色直线组成的空间为特征空间。
1.2.7 求行列式
求一个矩阵的行列式。
代码输入
E [[1, 2, 3],[4, 5, 6],[7, 8, 9]]
print(np.linalg.det(E))np.linalg.det() 是NumPy库中的一个函数用于计算矩阵的 行列式 。在这里我们将矩阵E作为参数传递给np.linalg.det()函数
结果输出
-9.51619735392994e-16扩展阅读 图说行列式几张图让你明白行列式的性质 行列式的几何意义 总结其实一个行列式的几何意义是有向线段(一阶行列式)或有向面积(二阶行列式)或有向体积(高阶行列式)。
1.2.8 奇异值分解实现
接下来利用奇异值分解Singular Value DecompositionSVD通过文章标题出现的关键词对文章进行聚类。
导入相应模块
import numpy as np
import matplotlib.pyplot as plt输入关键字
words [books,dad,stock,value,singular,estate,decomposition]设已知8个标题7个关键字。记录每个标题中每个关键字出现的次数得矩阵X。 X中每一行表示一个标题每一列表示一个关键字矩阵中的每个元素表示一个关键字中一个标题中出现的次数。
Xnp.array([[0,2,1,0,0,0,0],[2,0,0,1,0,1,0],[1,0,0,0,0,0,1],[0,0,1,0,0,0,0],[0,1,0,0,0,0,0],[0,0,0,1,1,0,1],[0,1,0,0,1,0,0],[0,0,0,0,1,1,1]])进行奇异值分解
U,s,Vhnp.linalg.svd(X)输出左奇异矩阵U
print(U,U)输出结果
U [[-1.87135757e-01 -7.93624528e-01 2.45011855e-01 -2.05404352e-01-3.88578059e-16 5.75779114e-16 -2.57394431e-01 -4.08248290e-01][-6.92896814e-01 2.88368077e-01 5.67788037e-01 2.22142537e-012.54000254e-01 -6.37019839e-16 -2.21623012e-02 2.05865892e-17][-3.53233681e-01 1.22606651e-01 3.49203461e-02 -4.51735990e-01-7.62000762e-01 1.27403968e-15 2.72513448e-01 3.80488702e-17][-2.61369658e-02 -1.33189110e-01 7.51079037e-02 -6.44727454e-015.08000508e-01 1.77635684e-15 3.68146235e-01 4.08248290e-01][-8.04993957e-02 -3.30217709e-01 8.49519758e-02 2.19661551e-01-2.54000254e-01 -4.81127681e-16 -3.12770333e-01 8.16496581e-01][-3.95029694e-01 1.56123876e-02 -5.28290830e-01 -6.82340484e-021.27000127e-01 -7.07106781e-01 -2.09360158e-01 1.55512464e-17][-2.02089013e-01 -3.80395849e-01 -2.12899198e-01 4.80790894e-018.04483689e-16 -1.60632798e-15 7.33466480e-01 1.76241226e-16][-3.95029694e-01 1.56123876e-02 -5.28290830e-01 -6.82340484e-021.27000127e-01 7.07106781e-01 -2.09360158e-01 -1.23226632e-16]]输出奇异值矩阵
print(s,s)按每个奇异值一一对应一个左奇异向量和一个右奇异向量奇异值从大到小排列输出结果
s [2.85653844 2.63792139 2.06449303 1.14829917 1. 1.0.54848559]输出右奇异矩阵Vh
print(Vh,Vh)输出结果
Vh [[-6.08788345e-01 -2.29949618e-01 -7.46612474e-02 -3.80854846e-01-3.47325416e-01 -3.80854846e-01 -4.00237243e-01][ 2.65111314e-01 -8.71088358e-01 -3.51342402e-01 1.15234846e-01-1.32365989e-01 1.15234846e-01 5.83153945e-02][ 5.66965547e-01 1.75382762e-01 1.55059743e-01 1.91316736e-02-6.14911671e-01 1.91316736e-02 -4.94872736e-01][-6.48865369e-03 2.52237176e-01 -7.40339999e-01 1.34031699e-012.99854608e-01 1.34031699e-01 -5.12239408e-01][-2.54000254e-01 -2.54000254e-01 5.08000508e-01 3.81000381e-012.54000254e-01 3.81000381e-01 -5.08000508e-01][ 0.00000000e00 -7.68640544e-16 2.33583082e-15 -7.07106781e-01-1.21802199e-15 7.07106781e-01 1.91457709e-15][ 4.16034348e-01 -1.71550021e-01 2.01922906e-01 -4.22112199e-015.73845817e-01 -4.22112199e-01 -2.66564648e-01]]规定坐标轴的范围
plt.axis([-0.8,0.2,-0.8,0.8])输出图像为
原每个关键字由 1*8 的向量表示现降维成 1*2 的向量以便进行可视化
for i in range(len(words)):plt.text(U[i,0],U[i,1],words[i])
plt.show()可视化结果 这张图是奇异值分解SVD的结果。奇异值分解是一种在线性代数中常用的矩阵分解方法它将一个矩阵分解为三个矩阵的乘积这三个矩阵分别代表了旋转、缩放和另一次旋转。在这张图中每个点代表一个单词它们的位置是通过奇异值分解得到的。这种分解方法可以帮助我们理解数据的结构和关系。
将到2维可视化后我们可以将关键词聚类如singular和decomposition距离比较近可以被划分为一组。
扩展阅读 2分钟看懂奇异值分解 什么是奇异值分解SVD–SVD如何分解时空矩阵 补充奇异值是特征值的开根。
1.2.9 线性方程组求解
求解线性方程组比较简单只需要用到一个函数 scipy.linalg.solve() 就可以了。
比如我们对胶片中矩阵章节中的部门月度跑步案例进行线性方程组求解线性方程组如下
10x_1 8x_2 12x_3 20
4x_1 4x_2 2x_3 8
2x_1 - 4x_2- 2x_3 -5代码输入
from scipy.linalg import solve
a np.array([[10, 8, 12], [4, 4, 2], [2, -4, -2]])
b np.array([20,8,-5])
x solve(a, b)
print(x)结果输出
[0.5 1.3125 0.375 ]
