用户2769421 | 作者
腾讯云 云+社区 | 来源
Numpy是一个用python实现的科学计算的扩展程序库,包括:
- 一个强大的N维数组对象Array;
- 比较成熟的(广播)函数库;
- 用于整合C/C++和Fortran代码的工具包;
- 实用的线性代数、傅里叶变换和随机数生成函数。numpy和稀疏矩阵运算包scipy配合使用更加方便。
NumPy(Numeric Python)提供了许多高级的数值编程工具,如:矩阵数据类型、矢量处理,以及精密的运算库。专为进行严格的数字处理而产生。多为很多大型金融公司使用,以及核心的科学计算组织如:Lawrence Livermore,NASA用其处理一些本来使用C++,Fortran或Matlab等所做的任务。
本文整理了一个Numpy的小抄表,总结了Numpy的常用操作,可以收藏慢慢看。
(图片可以点开大图查看哦~)
一.安装Numpy
可以通过 Pip 或者 Anaconda安装Numpy:
$ pip install numpy
或
$ conda install numpy
二.基础
NumPy最常用的功能之一就是NumPy数组:列表和NumPy数组的最主要区别在于功能性和速度。
列表提供基本操作,但NumPy添加了FTTs、卷积、快速搜索、基本统计、线性代数、直方图等。
两者数据科学最重要的区别是能够用NumPy数组进行元素级计算。
- axis 0:通常执行
- axis 1:通常指列
![5c3d97f4f2eaf11a894c6380ef34d41e.png](https://i-blog.csdnimg.cn/blog_migrate/b7aff51eb86d08af4793afa0eee7b2d2.jpeg)
1.占位符
![6e3a2f449c5cada25f936ba5f9bc3b30.png](https://i-blog.csdnimg.cn/blog_migrate/58196f07c7103e69c51be932682f1a53.jpeg)
举例:
import numpy as np# 1 dimensionalx = np.array([1,2,3])# 2 dimensionaly = np.array([(1,2,3),(4,5,6)])x = np.arange(3)>>> array([0, 1, 2])y = np.arange(3.0)>>> array([ 0., 1., 2.])x = np.arange(3,7)>>> array([3, 4, 5, 6])y = np.arange(3,7,2)>>> array([3, 5])
2.数组属性
![fac27f579b64784d19063dd65a6672d3.png](https://i-blog.csdnimg.cn/blog_migrate/e5180937c95dff42d876e4d0a519256a.jpeg)
3.拷贝 /排序
![661a3da74273579879e068fc5b47054c.png](https://i-blog.csdnimg.cn/blog_migrate/2d96415520cad1e3df9d75ad77c2f1d7.jpeg)
举例:
import numpy as np# Sort sorts in ascending ordery = np.array([10, 9, 8, 7, 6, 5, 4, 3, 2, 1])y.sort()print(y)>>> [ 1 2 3 4 5 6 7 8 9 10]
4.数组操作例程
增加或减少元素
![088b812c6caabc0f766fba8852d4d8c2.png](https://i-blog.csdnimg.cn/blog_migrate/970e209dc579d4c4f03b37fa524b3371.jpeg)
举例:
import numpy as np# Append items to arraya = np.array([(1, 2, 3),(4, 5, 6)])b = np.append(a, [(7, 8, 9)])print(b)>>> [1 2 3 4 5 6 7 8 9]# Remove index 2 from previous arrayprint(np.delete(b, 2))>>> [1 2 4 5 6 7 8 9]
组合数组
![6e5064fc4b92082808132a6f0f6e253b.png](https://i-blog.csdnimg.cn/blog_migrate/4442e5fd5f97ecb6c70bcc51abea0cc7.jpeg)
举例:
import numpy as npa = np.array([1, 3, 5])b = np.array([2, 4, 6])# Stack two arrays row-wiseprint(np.vstack((a,b)))>>> [[1 3 5] [2 4 6]]# Stack two arrays column-wiseprint(np.hstack((a,b)))>>> [1 3 5 2 4 6]
分割数组
![16b03be4d10e90f494c90dea083ff583.png](https://i-blog.csdnimg.cn/blog_migrate/21e3bc77d7e347ff6d10661daa64d114.jpeg)
举例:
# Split array into groups of ~3a = np.array([1, 2, 3, 4, 5, 6, 7, 8])print(np.array_split(a, 3))>>> [array([1, 2, 3]), array([4, 5, 6]), array([7, 8])]
数组形状变化
- 操作
![e3aef40596123af787f257134b161f9a.png](https://i-blog.csdnimg.cn/blog_migrate/6c0f93801b619c27bcbe9e214f39d6d3.jpeg)
- 其他
![662deb3cf3a82a2d452c55edbc075b9b.png](https://i-blog.csdnimg.cn/blog_migrate/c8f1b9ad7237dd43cd02218dec4d8be6.jpeg)
举例:
# Find inverse of a given matrix>>> np.linalg.inv([[3,1],[2,4]])array([[ 0.4, -0.1], [-0.2, 0.3]])
5.数学计算
操作
![1682f9c4f7c367c1dc286dae91969f95.png](https://i-blog.csdnimg.cn/blog_migrate/c6ff27d927538ad3226b6e6aaef826ee.jpeg)
举例:
# If a 1d array is added to a 2d array (or the other way), NumPy# chooses the array with smaller dimension and adds it to the one# with bigger dimensiona = np.array([1, 2, 3])b = np.array([(1, 2, 3), (4, 5, 6)])print(np.add(a, b))>>> [[2 4 6] [5 7 9]] # Example of np.roots# Consider a polynomial function (x-1)^2 = x^2 - 2*x + 1# Whose roots are 1,1>>> np.roots([1,-2,1])array([1., 1.])# Similarly x^2 - 4 = 0 has roots as x=±2>>> np.roots([1,0,-4])array([-2., 2.])
比较
![bc08262cb448943f146248022558293d.png](https://i-blog.csdnimg.cn/blog_migrate/f8df5a051e43d362c5bda322215b510b.jpeg)
举例:
# Using comparison operators will create boolean NumPy arraysz = np.array([1, 2, 3, 4, 5, 6, 7, 8, 9, 10])c = z >> [ True True True True True False False False False False]
基本的统计
![00e1ec08fbb4c3e1d736cb583a09798a.png](https://i-blog.csdnimg.cn/blog_migrate/6665b027a9bdb0278276718cfcc3726e.jpeg)
举例:
# Statistics of an arraya = np.array([1, 1, 2, 5, 8, 10, 11, 12])# Standard deviationprint(np.std(a))>>> 4.2938910093294167# Medianprint(np.median(a))>>> 6.5
更多
![4439fd7a3a4612118a64220bf7cfb8a7.png](https://i-blog.csdnimg.cn/blog_migrate/c02c29bed6c08ba967d849861fdd57c2.jpeg)
6.切片和子集
![09999d3f4dab0ad83ae8a4bce6fd4720.png](https://i-blog.csdnimg.cn/blog_migrate/b010d36ec673e297a42356fa87c5aa53.jpeg)
举例:
b = np.array([(1, 2, 3), (4, 5, 6)])# The index *before* the comma refers to *rows*,# the index *after* the comma refers to *columns*print(b[0:1, 2])>>> [3]print(b[:len(b), 2])>>> [3 6]print(b[0, :])>>> [1 2 3]print(b[0, 2:])>>> [3]print(b[:, 0])>>> [1 4]c = np.array([(1, 2, 3), (4, 5, 6)])d = c[1:2, 0:2]print(d)>>> [[4 5]]
切片举例:
import numpy as npa1 = np.arange(0, 6)a2 = np.arange(10, 16)a3 = np.arange(20, 26)a4 = np.arange(30, 36)a5 = np.arange(40, 46)a6 = np.arange(50, 56)a = np.vstack((a1, a2, a3, a4, a5, a6))
生成矩阵和切片图示
![d83ca24862a31976fd6e8776309da2a5.png](https://i-blog.csdnimg.cn/blog_migrate/40b399f03dcc3f91c7fb40d1e91040d8.jpeg)
![d65be11b88ae6639ecc44c22b749dbb3.png](https://i-blog.csdnimg.cn/blog_migrate/031fc2ece739037d7d00d21a56dfded9.jpeg)
7.小技巧
布尔索引
# Index trick when working with two np-arraysa = np.array([1,2,3,6,1,4,1])b = np.array([5,6,7,8,3,1,2])# Only saves a at index where b == 1other_a = a[b == 1]#Saves every spot in a except at index where b != 1other_other_a = a[b != 1]
import numpy as npx = np.array([4,6,8,1,2,6,9])y = x > 5print(x[y])>>> [6 8 6 9]# Even shorterx = np.array([1, 2, 3, 4, 4, 35, 212, 5, 5, 6])print(x[x >> [1 2 3 4 4]
参考:
https://github.com/juliangaal/python-cheat-sheet
- END -
本文为转载分享&推荐阅读,若侵权请联系后台删除