Numpy科学计算从放弃到入门

什么是NumPy

NumPy是Python科学计算的基础包,不仅是python中使用最多的第三方库,还是SciPy、Pandas等数据科学的基础库。所提供的结构比Python自身的更高级、更高效。
可以说,NumPy所提供的数据结果是Python数据分析的基础。
标准的Python用列表list保存数组的值,list的元素在系统中内存是分散存储的,而Numpy的数组结构是存储在一个均匀连续的内存块中,遍历所有的元素更高效,节省了计算资源。
主要提供了以下功能:

  • 快速高效的多维数组对象的ndarray
  • 用于对数组执行元素级计算以及直接对数组执行数学运算的函数
  • 用于读写硬盘上基于数组的数据集的工具
  • 线性代数运算、傅立叶变换,以及随机数组生成
  • 用于将C、C++、Fortran代码集成到Python的工具
  • 作为算法之间传递数据的容器,对于数值型数据,Numpy数组在存储和处理数据时要比内置的Python数据结构高效得多

numpy提供了一个ndarray对象,多维数组的含义,维度称为rank(秩),每个线性的数组称为一个axes(轴)。数组是有属性的,可以通过函数shape获得数组的大小,通过dtype获得元素的属性

相关帮助文档

一、如何创建数组

  • array
    将输入数据(列表、元组、数组或其他序列类型)转换为ndarray。要么推断出dtype,要么显示指定dtype

  • asarray
    将输入转换为ndarray,如果输入本身就是一个ndarray就不再复制
    array和asarray都可以将结构数据转化为ndarray,但是主要区别就是当数据源是ndarray时,array仍然会copy出一个副本,占用新的内存,但asarray不会

  • arange
    类似内置的range,返回一个ndarray

  • linespace
    生成等差数列,可以指定包含名都不包含终点值

  • zeros、zeros_like
    根据指定的形状和dtype创建一个全是0的数组。
    zeros_like以另一个数组为参数,并根据其形状和dtype创建一个全0数组

  • ones、ones_like
    根据指定的形状和dtype创建一个全是1的数组。
    ones_like以另一个数组为参数,并根据其形状和dtype创建一个全1数组

  • empty、empty_like
    类似ones和ones_like,只是产生的是全0数组

  • eye、identity
    创建一个正方的N*N单位矩阵,对角线为1,其余为0

二、如何访问数组

  • 下标索引
    arr2d=np.array([[1,2,3],[4,5,6],[7,8,9]])
    arr2d[0][2]和arr2d[0,2]是等价的

  • 切片索引
    可以一次传入多个切片,
    arr2d[:2, 1:]

  • 布尔型索引

将data中所有负值都设置为0

data[data < 0] = 0
  • 选出所有名字为“cqh”的所有行
data[names == 'cqh']
  • 整数索引
    利用整数数组进行索引
arr = np.empty((8,4))
for i in range(8):
    arr[i] = i
传入指定顺序
arr[[4,3,0,6]]
* 负数索引 
arr[[-3,-5,-7]]
  • 方形索引
    使用np.ix_方法,可以将两个一维整数数组转换为一个用于选取方形区域的索引
arr = np.arange(32).reshape((8,4))

获取行[1,5,7,2],列[0,3,1,2],为

arr[np.ix_([1,5,7,2],[0,3,1,2])]

相当于执行了

arr[[1,5,7,2]][:,[0,3,1,2]]

错误方法:

arr[[1,5,7,2],[0,3,1,2]]

三、如何做算数运算

x1 = np.arange(1,11,2)
x2 = np.linspace(1,9,5)
np.add(x1, x2)
np.subtract(x1, x2)
np.multiply(x1, x2)
np.divide(x1, x2)
np.power(x1, x2)
np.remainder(x1, x2)

四、如何使用统计函数

默认axis=0代表沿着列统计,axis=1代表沿着行统计

  • amax:矩阵中的最大值
  • amin:矩阵中的最小值
  • ptp:统计最大值与最小值的差
  • percentile:统计数组的百分位数
  • median:统计中位数
  • mean:平均数
  • average:加权平均值(需要传入各个元素的权重)
  • std:标准差
  • var:方差

五、如何转置数组和反转

  • reshape
arr = np.arange(15).reshape((3,5))

  • T
arr.T

六、如何去重、计算差集并集等集合操作

  • unique(x):计算x中的唯一元素,并返回有序结果
  • intersect1d(x,y):计算x和y中的公共元素,并返回有序结果
  • union(x,y):计算x和y的并集,并返回有序结果
  • in1d(x,y):得到一个表示“x的元素是否包含在y”的布尔型数组
  • setdiff1d(x,y) 集合的差,元素在x中且不在y中
  • setxor1d(x,y) 集合的对称差,存在于一个数组中但不同时存在两个数组中的元素

七、如何生成随机数

  • rand:产生[0,1]均匀分布的样本值

  • uniform:产生在[0,1]中均匀分页的样本值

  • randn:产生正态分布(平均值0,标准差1)的样本值,类似MATLAB接口

  • normal:生成标准正态分布的样本数组

  • randint:在给定的上下限范围内随机选取整数

  • seed:确定随机数生成器的种子

  • permutation:返回一个序列的随机排列

  • shuffle:对一个序列就地随机排列

  • bionmial:产生二项分布的样本值

  • beta:生产Beta分布的样本值

  • chisquare:产生卡方分布的样本值

  • gamma:产生Gamma分布的样本值

八、如何输入输出文件

以二进制格式保存到磁盘

  • np.save:保存到磁盘
  • np.savez:将多个数组保存到一个压缩文件中,加载时得到一个字典对象
  • np.load:从磁盘加载

存取文本文件

  • loadtxt:指定各种分隔符
  • savetxt:将数组写到以某种分隔符隔开的文本文件
posted @ 2020-01-21 17:49  雪山飞猪  阅读(770)  评论(0编辑  收藏  举报