numpy矩阵的运算2(numpy最小二乘法)

网友投稿 417 2022-08-29


numpy矩阵的运算2(numpy最小二乘法)

学习资料:

​​Numpy官方英文教材​​

通过上一节的学习,我们可以了解到一部分矩阵中元素的计算和查找操作。然而在日常使用中,对应元素的索引也是非常重要的。依然,让我们先从一个脚本开始:

import numpy as npA = np.arange(2,14).reshape((3,4)) # array([[ 2, 3, 4, 5]# [ 6, 7, 8, 9]# [10,11,12,13]]) print(np.argmin(A)) # 0print(np.argmax(A)) # 11

numpy 的几种基本运算

其中的 ​​argmin()​​​ 和 ​​argmax()​​ 两个函数分别对应着求矩阵中最小元素和最大元素的索引。相应的,在矩阵的12个元素中,最小值即2,对应索引0,最大值为13,对应索引为11。

如果需要计算统计中的均值,可以利用下面的方式,将整个矩阵的均值求出来:

print(np.mean(A)) # 7.5print(np.average(A)) # 7.5

仿照着前一节中​​dot()​​​ 的使用法则,​​mean()​​函数还有另外一种写法:

print(A.mean()) # 7.5

同样的,我们可以写出求解中位数的函数:

print(A.median()) # 7.5

另外,和matlab中的​​cumsum()​​​累加函数类似,Numpy中也具有​​cumsum()​​函数,其用法如下:

print(np.cumsum(A)) # [2 5 9 14 20 27 35 44 54 65 77 90]

在​​cumsum()​​​函数中:生成的每一项矩阵元素均是从原矩阵首项累加到对应项的元素之和。比如元素9,在​​cumsum()​​生成的矩阵中序号为3,即原矩阵中2,3,4三个元素的和。

相应的有累差运算函数:

print(np.diff(A)) # [[1 1 1]# [1 1 1]# [1 1 1]]

该函数计算的便是每一行中后一项与前一项之差。故一个3行4列矩阵通过函数计算得到的矩阵便是3行3列的矩阵。

下面我们介绍一下​​nonzero()​​函数:

print(np.nonzero(A)) # (array([0,0,0,0,1,1,1,1,2,2,2,2]),array([0,1,2,3,0,1,2,3,0,1,2,3]))

这个函数将所有非零元素的行与列坐标分割开,重构成两个分别关于行和列的矩阵。

同样的,我们可以对所有元素进行仿照列表一样的排序操作,但这里的排序函数仍然仅针对每一行进行从小到大排序操作:

import numpy as npA = np.arange(14,2, -1).reshape((3,4)) # array([[14, 13, 12, 11],# [10, 9, 8, 7],# [ 6, 5, 4, 3]])print(np.sort(A)) # array([[11,12,13,14]# [ 7, 8, 9,10]# [ 3, 4, 5, 6]])

矩阵的转置有两种表示方法:

print(np.transpose(A)) print(A.T)# array([[14,10, 6]# [13, 9, 5]# [12, 8, 4]# [11, 7, 3]])# array([[14,10, 6]# [13, 9, 5]# [12, 8, 4]# [11, 7, 3]])

特别的,在Numpy中具有​​clip()​​函数,例子如下:

print(A)# array([[14,13,12,11]# [10, 9, 8, 7]# [ 6, 5, 4, 3]])print(np.clip(A,5,9)) # array([[ 9, 9, 9, 9]# [ 9, 9, 8, 7]# [ 6, 5, 5, 5]])

这个函数的格式是​​clip(Array,Array_min,Array_max)​​​,顾名思义,​​Array​​指的是将要被执行用的矩阵,而后面的最小值最大值则用于让函数判断矩阵中元素是否有比最小值小的或者比最大值大的元素,并将这些指定的元素转换为最小值或者最大值。

实际上每一个Numpy中大多数函数均具有很多变量可以操作,你可以指定行、列甚至某一范围中的元素。更多具体的使用细节请记得查阅​​Numpy官方英文教材​​。


版权声明:本文内容由网络用户投稿,版权归原作者所有,本站不拥有其著作权,亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容,请联系我们jiasou666@gmail.com 处理,核实后本网站将在24小时内删除侵权内容。

上一篇:Numpy 索引(numpy官网)
下一篇:Java实现顺序表和链表结构
相关文章

 发表评论

暂时没有评论,来抢沙发吧~