python【Numpy科学计算库】连女朋友都会用的Numpy(真の能看懂~!)
文章目录
- 1 Numpy属性
- 2 Numpy的创建array
- 3 Numpy 基础运算1
- 4 Numpy 基础运算2
- 5 Numpy索引
- 6 Numpy array合并
- 7 Numpy array分割
- 8 Numpy copy & deep copy
1 Numpy属性
这次我们会介绍几种 numpy 的属性:ndim:维度
shape:行数和列数
size:元素个数
使用numpy首先要导入模块
import numpy as np #为了方便使用numpy 采用np简写
列表转化为矩阵:array = np.array([[1,2,3],[2,3,4]]) #列表转化为矩阵
print(array)
"""
array([[1, 2, 3],[2, 3, 4]])
"""
numpy 的几种属性
接着我们看看这几种属性的结果:print('number of dim:',array.ndim) # 维度
#number of dim: 2print('shape :',array.shape) # 行数和列数
#shape : (2, 3)print('size:',array.size) # 元素个数
#size: 6
2 Numpy的创建array
import numpy as npa = np.array([2,23,4]) # list 1d
print(a)a = np.array([2,23,4],dtype=np.int)
print(a.dtype)
创建特定数据 ¶
a = np.array([[2,23,4],[2,32,4]]) # 2d 矩阵 2行3列
print(a)
“”"
[[ 2 23 4]
[ 2 32 4]]
“”"
创建全零数组
a = np.zeros((3,4)) # 数据全为0,3行4列
“”"
array([[ 0., 0., 0., 0.],
[ 0., 0., 0., 0.],
[ 0., 0., 0., 0.]])
“”"
创建全一数组, 同时也能指定这些特定数据的 dtype:
a = np.ones((3,4),dtype = np.int) # 数据为1,3行4列
“”"
array([[1, 1, 1, 1],
[1, 1, 1, 1],
[1, 1, 1, 1]])
“”"
创建全空数组, 其实每个值都是接近于零的数:
a = np.empty((3,4)) # 数据为empty,3行4列
“”"
array([[ 0.00000000e+000, 4.94065646e-324, 9.88131292e-324,
1.48219694e-323],
[ 1.97626258e-323, 2.47032823e-323, 2.96439388e-323,
3.45845952e-323],
[ 3.95252517e-323, 4.44659081e-323, 4.94065646e-323,
5.43472210e-323]])
“”"
用 arange 创建连续数组:
a = np.arange(10,20,2) # 10-19 的数据,2步长
“”"
array([10, 12, 14, 16, 18])
“”"
使用 reshape 改变数据的形状
a = np.arange(12).reshape((3,4)) # 3行4列,0到11
“”"
array([[ 0, 1, 2, 3],
[ 4, 5, 6, 7],
[ 8, 9, 10, 11]])
“”"
用 linspace 创建线段型数据:
a = np.linspace(1,10,20) # 开始端1,结束端10,且分割成20个数据,生成线段
“”"
array([ 1. , 1.47368421, 1.94736842, 2.42105263,
2.89473684, 3.36842105, 3.84210526, 4.31578947,
4.78947368, 5.26315789, 5.73684211, 6.21052632,
6.68421053, 7.15789474, 7.63157895, 8.10526316,
8.57894737, 9.05263158, 9.52631579, 10. ])
“”"
同样也能进行 reshape 工作:
a = np.linspace(1,10,20).reshape((5,4)) # 更改shape
“”"
array([[ 1. , 1.47368421, 1.94736842, 2.42105263],
[ 2.89473684, 3.36842105, 3.84210526, 4.31578947],
[ 4.78947368, 5.26315789, 5.73684211, 6.21052632],
[ 6.68421053, 7.15789474, 7.63157895, 8.10526316],
[ 8.57894737, 9.05263158, 9.52631579, 10. ]])
“”"
3 Numpy 基础运算1
numpy 的几种基本运算
上述代码中的 a 和 b 是两个属性为 array 也就是矩阵的变量,而且二者都是1行4列的矩阵, 其中b矩阵中的元素分别是从0到3。 如果我们想要求两个矩阵之间的减法,你可以尝试着输入:
c=a-b # array([10, 19, 28, 37])
通过执行上述脚本,将会得到对应元素相减的结果,即[10,19,28,37]。 同理,矩阵对应元素的相加和相乘也可以用类似的方式表示:
c=a+b # array([10, 21, 32, 43])
c=a*b # array([ 0, 20, 60, 120])
有所不同的是,在Numpy中,想要求出矩阵中各个元素的乘方需要依赖双星符号 **,以二次方举例,即:
c=b**2 # array([0, 1, 4, 9])
另外,Numpy中具有很多的数学函数工具,比如三角函数等,当我们需要对矩阵中每一项元素进行函数运算时,可以很简便的调用它们(以sin函数为例):
c=10*np.sin(a)
#array([-5.44021111, 9.12945251, -9.88031624, 7.4511316 ])
除了函数应用外,在脚本中对print函数进行一些修改可以进行逻辑判断:
print(b<3)
#array([ True, True, True, False], dtype=bool)
此时由于进行逻辑判断,返回的是一个bool类型的矩阵,即对满足要求的返回True,不满足的返回False。上述程序执行后得到的结果是[True True True False]。 需要注意的是,如果想要执行是否相等的判断, 依然需要输入 == 而不是 = 来完成相应的逻辑判断。
上述运算均是建立在一维矩阵,即只有一行的矩阵上面的计算,如果我们想要对多行多维度的矩阵进行操作,需要对开始的脚本进行一些修改:
a=np.array([[1,1],[0,1]])
b=np.arange(4).reshape((2,2))
print(a)
#array([[1, 1],
#[0, 1]])
print(b)
#array([[0, 1],
#[2, 3]])
此时构造出来的矩阵a和b便是2行2列的,其中 reshape 操作是对矩阵的形状进行重构, 其重构的形状便是括号中给出的数字。 稍显不同的是,Numpy中的矩阵乘法分为两种, 其一是前文中的对应元素相乘,其二是标准的矩阵乘法运算,即对应行乘对应列得到相应元素:
c_dot = np.dot(a,b)
#array([[2, 4],
#[2, 3]])
除此之外还有另外的一种关于dot的表示方法,即:
c_dot_2 = a.dot(b)
#array([[2, 4],
#[2, 3]])
下面我们将重新定义一个脚本, 来看看关于 sum(), min(), max()的使用:
import numpy as np
a=np.random.random((2,4))
print(a)
#array([[ 0.94692159, 0.20821798, 0.35339414, 0.2805278 ],
#[ 0.04836775, 0.04023552, 0.44091941, 0.21665268]])
因为是随机生成数字, 所以你的结果可能会不一样. 在第二行中对a的操作是令a中生成一个2行4列的矩阵,且每一元素均是来自从0到1的随机数。 在这个随机生成的矩阵中,我们可以对元素进行求和以及寻找极值的操作,具体如下:
np.sum(a) # 4.4043622002745959
np.min(a) # 0.23651223533671784
np.max(a) # 0.90438450240606416
对应的便是对矩阵中所有元素进行求和,寻找最小值,寻找最大值的操作。 可以通过print()函数对相应值进行打印检验。
如果你需要对行或者列进行查找运算,就需要在上述代码中为 axis 进行赋值。 当axis的值为0的时候,将会以列作为查找单元, 当axis的值为1的时候,将会以行作为查找单元。
为了更加清晰,在刚才的例子中我们继续进行查找:
print(“a =”,a)
#a = [[ 0.23651224 0.41900661 0.84869417 0.46456022]
#[ 0.60771087 0.9043845 0.36603285 0.55746074]]
print(“sum =”,np.sum(a,axis=1))
#sum = [ 1.96877324 2.43558896]
print(“min =”,np.min(a,axis=0))
#min = [ 0.23651224 0.41900661 0.36603285 0.46456022]
print(“max =”,np.max(a,axis=1))
#max = [ 0.84869417 0.9043845 ]
4 Numpy 基础运算2
先从一个脚本开始 :import numpy as np
A = np.arange(2,14).reshape((3,4)) # array([[ 2, 3, 4, 5]
# [ 6, 7, 8, 9]
# [10,11,12,13]])print(np.argmin(A)) # 0
print(np.argmax(A)) # 11
numpy 的几种基本运算
其中的 argmin() 和 argmax() 两个函数分别对应着求矩阵中最小元素和最大元素的索引。相应的,在矩阵的12个元素中,最小值即2,对应索引0,最大值为13,对应索引为11。如果需要计算统计中的均值,可以利用下面的方式,将整个矩阵的均值求出来:print(np.mean(A)) # 7.5
print(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 np
A = 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官方英文教材。
5 Numpy索引
我们都知道,在元素列表或者数组中,我们可以用如同a[2]一样的表示方法,同样的,在Numpy中也有相对应的表示方法:import numpy as np
A = np.arange(3,15)# array([3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14])print(A[3]) # 6
让我们将矩阵转换为二维的,此时进行同样的操作:A = np.arange(3,15).reshape((3,4))
"""
array([[ 3, 4, 5, 6][ 7, 8, 9, 10][11, 12, 13, 14]])
"""print(A[2])
# [11 12 13 14]
实际上这时的A[2]对应的就是矩阵A中第三行(从0开始算第一行)的所有元素。二维索引
如果你想要表示具体的单个元素,可以仿照上述的例子:print(A[1][1]) # 8
此时对应的元素即A[1][1],在A中即横纵坐标都为1,第二行第二列的元素,即8(因为计数从0开始)。同样的还有其他的表示方法:print(A[1, 1]) # 8
在Python的 list 中,我们可以利用:对一定范围内的元素进行切片操作,在Numpy中我们依然可以给出相应的方法:print(A[1, 1:3]) # [8 9]
这一表示形式即针对第二行中第2到第4列元素进行切片输出(不包含第4列)。 此时我们适当的利用for函数进行打印:for row in A:print(row)
"""
[ 3, 4, 5, 6]
[ 7, 8, 9, 10]
[11, 12, 13, 14]
"""
此时它会逐行进行打印操作。如果想进行逐列打印,就需要稍稍变化一下:for column in A.T:print(column)
"""
[ 3, 7, 11]
[ 4, 8, 12]
[ 5, 9, 13]
[ 6, 10, 14]
"""
上述表示方法即对A进行转置,再将得到的矩阵逐行输出即可得到原矩阵的逐列输出。最后依然说一些关于迭代输出的问题:import numpy as np
A = np.arange(3,15).reshape((3,4))print(A.flatten())
# array([3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14])for item in A.flat:print(item)# 3
# 4
……
# 14
这一脚本中的flatten是一个展开性质的函数,将多维的矩阵进行展开成1行的数列。而flat是一个迭代器,本身是一个object属性。
6 Numpy array合并
np.vstack() ¶
对于一个array的合并,我们可以想到按行、按列等多种方式进行合并。首先先看一个例子:import numpy as np
A = np.array([1,1,1])
B = np.array([2,2,2])print(np.vstack((A,B))) # vertical stack
"""
[[1,1,1][2,2,2]]
"""
vertical stack本身属于一种上下合并,即对括号中的两个整体进行对应操作。此时我们对组合而成的矩阵进行属性探究:C = np.vstack((A,B))
print(A.shape,C.shape)# (3,) (2,3)
np.hstack()
利用shape函数可以让我们很容易地知道A和C的属性,从打印出的结果来看,A仅仅是一个拥有3项元素的数组(数列),而合并后得到的C是一个2行3列的矩阵。介绍完了上下合并,我们来说说左右合并:D = np.hstack((A,B)) # horizontal stackprint(D)
# [1,1,1,2,2,2]print(A.shape,D.shape)
# (3,) (6,)
通过打印出的结果可以看出:D本身来源于A,B两个数列的左右合并,而且新生成的D本身也是一个含有6项元素的序列。np.newaxis()
说完了array的合并,我们稍稍提及一下前一节中转置操作,如果面对如同前文所述的A序列, 转置操作便很有可能无法对其进行转置(因为A并不是矩阵的属性),此时就需要我们借助其他的函数操作进行转置:print(A[np.newaxis,:])
# [[1 1 1]]print(A[np.newaxis,:].shape)
# (1,3)print(A[:,np.newaxis])
"""
[[1]
[1]
[1]]
"""print(A[:,np.newaxis].shape)
# (3,1)
此时我们便将具有3个元素的array转换为了1行3列以及3行1列的矩阵了。结合着上面的知识,我们把它综合起来:import numpy as np
A = np.array([1,1,1])[:,np.newaxis]
B = np.array([2,2,2])[:,np.newaxis]C = np.vstack((A,B)) # vertical stack
D = np.hstack((A,B)) # horizontal stackprint(D)
"""
[[1 2]
[1 2]
[1 2]]
"""print(A.shape,D.shape)
# (3,1) (3,2)
np.concatenate()
当你的合并操作需要针对多个矩阵或序列时,借助concatenate函数可能会让你使用起来比前述的函数更加方便:C = np.concatenate((A,B,B,A),axis=0)print(C)
"""
array([[1],[1],[1],[2],[2],[2],[2],[2],[2],[1],[1],[1]])
"""D = np.concatenate((A,B,B,A),axis=1)print(D)
"""
array([[1, 2, 2, 1],[1, 2, 2, 1],[1, 2, 2, 1]])
"""
axis参数很好的控制了矩阵的纵向或是横向打印,相比较vstack和hstack函数显得更加方便。
7 Numpy array分割
创建数据 ¶
首先 import 模块import numpy as np
建立3行4列的ArrayA = np.arange(12).reshape((3, 4))
print(A)
"""
array([[ 0, 1, 2, 3],[ 4, 5, 6, 7],[ 8, 9, 10, 11]])
"""
纵向分割
print(np.split(A, 2, axis=1))
"""
[array([[0, 1],[4, 5],[8, 9]]), array([[ 2, 3],[ 6, 7],[10, 11]])]
"""
横向分割
print(np.split(A, 3, axis=0))# [array([[0, 1, 2, 3]]), array([[4, 5, 6, 7]]), array([[ 8, 9, 10, 11]])]
错误的分割
范例的Array只有4列,只能等量对分,因此输入以上程序代码后Python就会报错。print(np.split(A, 3, axis=1))# ValueError: array split does not result in an equal division
为了解决这种情况, 我们会有下面这种方式.不等量的分割
在机器学习时经常会需要将数据做不等量的分割,因此解决办法为np.array_split()print(np.array_split(A, 3, axis=1))
"""
[array([[0, 1],[4, 5],[8, 9]]), array([[ 2],[ 6],[10]]), array([[ 3],[ 7],[11]])]
"""
成功将Array不等量分割!其他的分割方式
在Numpy里还有np.vsplit()与横np.hsplit()方式可用。print(np.vsplit(A, 3)) #等于 print(np.split(A, 3, axis=0))# [array([[0, 1, 2, 3]]), array([[4, 5, 6, 7]]), array([[ 8, 9, 10, 11]])]print(np.hsplit(A, 2)) #等于 print(np.split(A, 2, axis=1))
"""
[array([[0, 1],[4, 5],[8, 9]]), array([[ 2, 3],[ 6, 7],[10, 11]])]
"""
8 Numpy copy & deep copy
= 的赋值方式会带有关联性 ¶
首先 import numpy 并建立变量, 给变量赋值。import numpy as npa = np.arange(4)
# array([0, 1, 2, 3])b = a
c = a
d = b
改变a的第一个值,b、c、d的第一个值也会同时改变。a[0] = 11
print(a)
# array([11, 1, 2, 3])
确认b、c、d是否与a相同。b is a # True
c is a # True
d is a # True
同样更改d的值,a、b、c也会改变。d[1:3] = [22, 33] # array([11, 22, 33, 3])
print(a) # array([11, 22, 33, 3])
print(b) # array([11, 22, 33, 3])
print(c) # array([11, 22, 33, 3])
copy() 的赋值方式没有关联性
b = a.copy() # deep copy
print(b) # array([11, 22, 33, 3])
a[3] = 44
print(a) # array([11, 22, 33, 44])
print(b) # array([11, 22, 33, 3])
此时a与b已经没有关联。
python【Numpy科学计算库】连女朋友都会用的Numpy(真の能看懂~!)相关推荐
- NumPy科学计算库-学习笔记
文章目录 NumPy科学计算库 0 准备 安装python库 jupyter 启动:jupyter notebook jupyter扩展插件 1 基本操作 1.1 数组创建 1.2 查看属性 1.3 ...
- MOOC网神经网络与深度学习TensorFlow实践2——numpy科学计算库、matplotlib
numpy科学计算库 只打乱第一维的数据: Matplotlib 散点图 折线图和柱状图 波士顿房价数据集可视化 鸢尾花数据集可视化
- python的科学计算库总结
参考地址:https://www.stxnext.com/blog/most-popular-python-scientific-libraries/?utm_campaign=Pat&utm ...
- NumPy科学计算库(1)Python数据分析
1 NumPy简介 简单来说NumPy的功能就是就是对可以一个数值集合(eg:数组)进行各种各样的操作. NumPy(Numerical Python)是Python的一种开源的数值计算扩展.提供 ...
- python【Numpy科学计算库】Numpy的超实用基础详细教程(看不懂你来打我)
文章目录 Numpy简单介绍 1.Numpy是什么 2. 多维数组 numpy数据类型设定与转换 numpy索引与切片 多维数组索引.切片 数组文件输入输出 使用numpy.arange方法 使用nu ...
- 神经网络与深度学习——TensorFlow2.0实战(笔记)(五)(NumPy科学计算库<1>python)
多维数组 形状(Shape): 是一个元组,描述数组的维度,以及各个维度的长度. 长度(Length): 某个维度中的元素个数. 数字 单门课程成绩 1*课程成绩 一维数组 多门课程成绩 n*课程成绩 ...
- Python Scipy 科学计算库
Python机器学习及分析工具:Scipy篇 原文:https://www.jianshu.com/p/6c742912047f Scipy是一个用于数学.科学.工程领域的常用软件包,可以处理插值 ...
- python的科学计算库有哪些_python科学计算:带你初探scipy库的常量模块和函数模块...
嗨,机智哥又跟大家见面了.前面几节课,我们了解了Numpy模块的一些基本功能,知道了如何运用Numpy模块中的一些函数解决一些我们生活中的问题. 当然,Numpy模块还有大部分我们没有讲到的函数,如果 ...
- python的科学计算库有哪些_Python科学计算库-Numpy
NumPy 是 Python 语言的一个扩充程序库.支持高级大量的维度数组与矩阵运算,此外也针对数组运算提供大量的数学函数库,也是学习 python 必学的一个库. 1. 读取文件 numpy.gen ...
最新文章
- 63、使用Timer类来实现定时任务
- 如何更改linux文件的拥有者及用户组(chown和chgrp)
- 【Python生成readme文件】——Markdown语法
- .Net转Java自学之路—SpringMVC框架篇八(RESTful支持)
- Android官方开发文档Training系列课程中文版:高效显示位图之加载大位图
- webx学习(二)——Webx Framework
- 怎么在云服务器上布项目,怎么在云服务器上布项目
- 【Java学习笔记之二十六】深入理解Java匿名内部类
- Jquery获得服务器控件的方法
- (二)关于ThinkPHP2.1版本操作MSSQL类的BUG--selec查询只得出1条记录
- 一、MySQL数据库优化策略
- 格雷码和二进制相互转换
- 推荐系统(一):不了解推荐系统,基本在被out的队列,岌岌可危
- 两数互素有什么性质_两个数互质是什么意思 判断互质数的方法有哪些
- 欢迎大家加入Xcode公社
- 游戏中的脚本语言原理与发展
- php和java做网站_java与php哪个做网站简单
- 100条经典C++笔试题目及答案分享
- Network problems last Friday
- Google Play 政策更新提醒与重点解读 | 2021 年第四季度
热门文章
- pybind传输list
- 人脸关键点 姿态笔记
- type torch.cuda.FloatTensor but found type torch.cuda.ByteTensor
- 车辆检测,车牌识别WPOD-NET OCR-Net
- VS2015编译Boost1.64
- 无法解析的外部符号 WinMain,该符号在函数 int __cdecl invoke_main(void) (?invoke_main@@YAHXZ) 中被引用
- java Rserve测试
- 淮海工学院软件测试技术实验三功能测试,软件工程-实验3软件测试技术实验(7页)-原创力文档...
- 向日葵win10远程linux主机,小猪为你win10系统使用向日葵远程桌面软件远程的设置方法...
- JAVA语言教学重点_《JAVA语言》教学大纲