numpy库:常用基本
转载自 http://www.cnblogs.com/smallpi/p/4550361.html侵删致歉
numpy
本文主要列出numpy模块常用方法
大部分内容来源于网络,而后经过自己的一点思考和总结,如果有侵权,请联系我
我是一名初学者,有哪些地方有错误请留言,我会及时更改的
创建矩阵(采用ndarray对象)
对于python中的numpy模块,一般用其提供的ndarray对象。
创建一个ndarray对象很简单,只要将一个list作为参数即可。
例如
![](/assets/blank.gif)
import numpy as np #引入numpy库#创建一维的narray对象 a = np.array([1,2,3,4,5])#创建二维的narray对象 a2 = np.array([[1,2,3,4,5],[6,7,8,9,10]])#创建多维对象以其类推
![](/assets/blank.gif)
获取矩阵行数列数(二维情况)
习惯了采用matlab进行数模的编程,要对矩阵进行遍历时,一般先获取矩阵的行数和列数。要获取narray对象的各维的长度,可以通过narray对象的shape属性
import numpy as np a = np.array([[1,2,3,4,5],[6,7,8,9,10]])print(a.shape) #结果返回一个tuple元组 (2L, 5L) print(a.shape[0]) #获得行数,返回 2 print(a.shape[1]) #获得列数,返回 5
矩阵的截取
按行列截取
矩阵的截取和list相同,可以通过[](方括号)来截取
![](/assets/blank.gif)
import numpy as np a = np.array([[1,2,3,4,5],[6,7,8,9,10]])print(a[0:1]) #截取第一行,返回 [[1 2 3 4 5]] print(a[1,2:5]) #截取第二行,第三、四列,返回 [8 9]print(a[1,:]) #截取第二行,返回 [ 6 7 8 9 10]
![](/assets/blank.gif)
按条件截取
按条件截取其实是在[](方括号)中传入自身的布尔语句
例如
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3,4,5],[6,7,8,9,10]]) b = a[a>6] # 截取矩阵a中大于6的元素,范围的是一维数组 print(b) # 返回 [ 7 8 9 10]# 其实布尔语句首先生成一个布尔矩阵,将布尔矩阵传入[](方括号)实现截取 print(a>6) # 返回 [[False False False False False][False True True True True]]
![](/assets/blank.gif)
按条件截取应用较多的是对矩阵中满足一定条件的元素变成特定的值。
例如将矩阵中大于6的元素变成0。
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3,4,5],[6,7,8,9,10]]) print(a) #开始矩阵为 [[ 1 2 3 4 5][ 6 7 8 9 10]]a[a>6] = 0 print(a) #大于6清零后矩阵为 [[1 2 3 4 5][6 0 0 0 0]]
![](/assets/blank.gif)
矩阵的合并
矩阵的合并可以通过numpy中的hstack方法和vstack方法实现
![](/assets/blank.gif)
import numpy as npa1 = np.array([[1,2],[3,4]]) a2 = np.array([[5,6],[7,8]])#!注意 参数传入时要以列表list或元组tuple的形式传入 print(np.hstack([a1,a2])) #横向合并,返回结果如下 [[1 2 5 6][3 4 7 8]]print(np.vstack((a1,a2))) #纵向合并,返回结果如下 [[1 2][3 4][5 6][7 8]]
![](/assets/blank.gif)
矩阵的合并也可以通过concatenatef方法。
np.concatenate( (a1,a2), axis=0 ) 等价于 np.vstack( (a1,a2) )
np.concatenate( (a1,a2), axis=1 ) 等价于 np.hstack( (a1,a2) )
通过函数创建矩阵
numpy模块中自带了一些创建ndarray对象的函数,可以很方便的创建常用的或有规律的矩阵。
arange
![](/assets/blank.gif)
import numpy as npa = np.arange(10) # 默认从0开始到10(不包括10),步长为1 print(a) # 返回 [0 1 2 3 4 5 6 7 8 9]a1 = np.arange(5,10) # 从5开始到10(不包括10),步长为1 print(a1) # 返回 [5 6 7 8 9]a2 = np.arange(5,20,2) # 从5开始到20(不包括20),步长为2 print(a2) # 返回 [ 5 7 9 11 13 15 17 19]
![](/assets/blank.gif)
linspace
linspace()和matlab的linspace很类似,用于创建指定数量等间隔的序列,实际生成一个等差数列。
import numpy as npa = np.linspace(0,10,7) # 生成首位是0,末位是10,含7个数的等差数列 print(a) # 结果 [ 0. 1.66666667 3.33333333 5. 6.66666667 8.33333333 10. ]
logspace
linspace用于生成等差数列,而logspace用于生成等比数列。
下面的例子用于生成首位是100,末位是102,含5个数的等比数列。
import numpy as npa = np.logspace(0,2,5) print(a) # 结果 [ 1. 3.16227766 10. 31.6227766 100. ]
ones、zeros、eye、empty
ones创建全1矩阵
zeros创建全0矩阵
eye创建单位矩阵
empty创建空矩阵(实际有值)
![](/assets/blank.gif)
import numpy as npa_ones = np.ones((3,4)) # 创建3*4的全1矩阵 print(a_ones) # 结果 [[ 1. 1. 1. 1.][ 1. 1. 1. 1.][ 1. 1. 1. 1.]]a_zeros = np.zeros((3,4)) # 创建3*4的全0矩阵 print(a_zeros) # 结果 [[ 0. 0. 0. 0.][ 0. 0. 0. 0.][ 0. 0. 0. 0.]]a_eye = np.eye(3) # 创建3阶单位矩阵 print(a_eye) # 结果 [ 1. 0. 0.][ 0. 1. 0.][ 0. 0. 1.]]a_empty = np.empty((3,4)) # 创建3*4的空矩阵 print(a_empty) # 结果 [[ 1.78006111e-306 -3.13259416e-294 4.71524461e-309 1.94927842e+289][ 2.10230387e-309 5.42870216e+294 6.73606381e-310 3.82265219e-297][ 6.24242356e-309 1.07034394e-296 2.12687797e+183 6.88703165e-315]]
![](/assets/blank.gif)
fromstring
fromstring()方法可以将字符串转化成ndarray对象,需要将字符串数字化时这个方法比较有用,可以获得字符串的ascii码序列。
a = "abcdef" b = np.fromstring(a,dtype=np.int8) # 因为一个字符为8为,所以指定dtype为np.int8 print(b) # 返回 [ 97 98 99 100 101 102]
fromfunction
fromfunction()方法可以根据矩阵的行号列号生成矩阵的元素。
例如创建一个矩阵,矩阵中的每个元素都为行号和列号的和。
![](/assets/blank.gif)
import numpy as npdef func(i,j): return i+ja = np.fromfunction(func,(5,6)) # 第一个参数为指定函数,第二个参数为列表list或元组tuple,说明矩阵的大小 print(a) # 返回 [[ 0. 1. 2. 3. 4. 5.][ 1. 2. 3. 4. 5. 6.][ 2. 3. 4. 5. 6. 7.][ 3. 4. 5. 6. 7. 8.][ 4. 5. 6. 7. 8. 9.]] #注意这里行号的列号都是从0开始的
![](/assets/blank.gif)
矩阵的运算
常用矩阵运算符
numpy中的ndarray对象重载了许多运算符,使用这些运算符可以完成矩阵间对应元素的运算。
运算符 | 说明 |
---|---|
+ | 矩阵对应元素相加 |
- | 矩阵对应元素相减 |
* | 矩阵对应元素相乘 |
/ | 矩阵对应元素相除,如果都是整数则取商 |
% | 矩阵对应元素相除后取余数 |
** | 矩阵每个元素都取n次方,如**2:每个元素都取平方 |
例如
![](/assets/blank.gif)
import numpy as np a1 = np.array([[4,5,6],[1,2,3]]) a2 = np.array([[6,5,4],[3,2,1]])print(a1+a2) # 相加 # 结果 [[10 10 10][ 4 4 4]]print(a1/a2) # 整数相除取商 # 结果 [[0 1 1][0 1 3]]print(a1%a2) # 相除取余数 # 结果 [[4 0 2][1 0 0]]
![](/assets/blank.gif)
常用矩阵函数
同样地,numpy中也定义了许多函数,使用这些函数可以将函数作用于矩阵中的每个元素。
表格中默认导入了numpy模块,即 import numpy as np
a为ndarray对象。
矩阵函数 | 说明 |
---|---|
np.sin(a) | 对矩阵a中每个元素取正弦,sin(x) |
np.cos(a) | 对矩阵a中每个元素取余弦,cos(x) |
np.tan(a) | 对矩阵a中每个元素取正切,tan(x) |
np.arcsin(a) | 对矩阵a中每个元素取反正弦,arcsin(x) |
np.arccos(a) | 对矩阵a中每个元素取反余弦,arccos(x) |
np.arctan(a) | 对矩阵a中每个元素取反正切,arctan(x) |
np.exp(a) | 对矩阵a中每个元素取指数函数,ex |
np.sqrt(a) | 对矩阵a中每个元素开根号√x |
例如
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]]) print(np.sin(a)) # 结果 [[ 0.84147098 0.90929743 0.14112001][-0.7568025 -0.95892427 -0.2794155 ]]print(np.arcsin(a)) # 结果 C:\Users\Administrator\Desktop\learn.py:6: RuntimeWarning: invalid value encountered in arcsinprint(np.arcsin(a)) [[ 1.57079633 nan nan][ nan nan nan]]
![](/assets/blank.gif)
当矩阵中的元素不在定义域范围内,会产生RuntimeWarning,结果为nan(not a number)。
[−1,1]
矩阵乘法(点乘)
矩阵乘法必须满足矩阵乘法的条件,即第一个矩阵的列数等于第二个矩阵的行数。
矩阵乘法的函数为 dot
例如
![](/assets/blank.gif)
import numpy as npa1 = np.array([[1,2,3],[4,5,6]]) # a1为2*3矩阵 a2 = np.array([[1,2],[3,4],[5,6]]) # a2为3*2矩阵print(a1.shape[1]==a2.shape[0]) # True, 满足矩阵乘法条件 print(a1.dot(a2)) # a1.dot(a2)相当于matlab中的a1*a2 # 而python中的a1*a2相当于matlab中的a1.*a2 # 结果 [[22 28][49 64]]
![](/assets/blank.gif)
矩阵的转置 aT
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]])print(a.transpose()) # 结果 [[1 4][2 5][3 6]]
![](/assets/blank.gif)
矩阵的转置还有更简单的方法,就是a.T
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]]) print(a.T) # 结果 [[1 4][2 5][3 6]]
![](/assets/blank.gif)
矩阵的逆 a−1
求矩阵的逆需要先导入numpy.linalg
,用linalg的inv函数来求逆。
矩阵求逆的条件是矩阵的行数和列数相同。
![](/assets/blank.gif)
import numpy as np import numpy.linalg as lga = np.array([[1,2,3],[4,5,6],[7,8,9]])print(lg.inv(a)) # 结果 [[ -4.50359963e+15 9.00719925e+15 -4.50359963e+15][ 9.00719925e+15 -1.80143985e+16 9.00719925e+15][ -4.50359963e+15 9.00719925e+15 -4.50359963e+15]]a = np.eye(3) # 3阶单位矩阵 print(lg.inv(a)) # 单位矩阵的逆为他本身 # 结果 [[ 1. 0. 0.][ 0. 1. 0.][ 0. 0. 1.]]
![](/assets/blank.gif)
矩阵信息获取(如平均值)
最大最小值
获得矩阵中元素最大最小值的函数分别是max
和min
,可以获得整个矩阵、行或列的最大最小值。
例如
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]]) print(a.max()) #获取整个矩阵的最大值 结果: 6 print(a.min()) #结果:1# 可以指定关键字参数axis来获得行最大(小)值或列最大(小)值 # axis=0 行方向最大(小)值,即获得每列的最大(小)值 # axis=1 列方向最大(小)值,即获得每行的最大(小)值 # 例如print(a.max(axis=0)) # 结果为 [4 5 6]print(a.max(axis=1)) # 结果为 [3 6]# 要想获得最大最小值元素所在的位置,可以通过argmax函数来获得 print(a.argmax(axis=1)) # 结果为 [2 2]
![](/assets/blank.gif)
平均值
获得矩阵中元素的平均值可以通过函数mean()
。同样地,可以获得整个矩阵、行或列的平均值。
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]]) print(a.mean()) #结果为: 3.5# 同样地,可以通过关键字axis参数指定沿哪个方向获取平均值 print(a.mean(axis=0)) # 结果 [ 2.5 3.5 4.5] print(a.mean(axis=1)) # 结果 [ 2. 5.]
![](/assets/blank.gif)
方差
方差的函数为var()
,方差函数var()
相当于函数mean(abs(x - x.mean())**2)
,其中x为矩阵。
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]]) print(a.var()) # 结果 2.91666666667print(a.var(axis=0)) # 结果 [ 2.25 2.25 2.25] print(a.var(axis=1)) # 结果 [ 0.66666667 0.66666667]
![](/assets/blank.gif)
标准差
标准差的函数为std()
。
std()
相当于sqrt(mean(abs(x - x.mean())**2))
,或相当于sqrt(x.var())
。
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]]) print(a.std()) # 结果 1.70782512766print(a.std(axis=0)) # 结果 [ 1.5 1.5 1.5] print(a.std(axis=1)) # 结果 [ 0.81649658 0.81649658]
![](/assets/blank.gif)
中值
中值指的是将序列按大小顺序排列后,排在中间的那个值,如果有偶数个数,则是排在中间两个数的平均值。
例如序列[5,2,6,4,2],按大小顺序排成 [2,2,4,5,6],排在中间的数是4,所以这个序列的中值是4。
又如序列[5,2,6,4,3,2],按大小顺序排成 [2,2,3,4,5,6],因为有偶数个数,排在中间两个数是3、4,所以这个序列中值是3.5。
中值的函数是median(),调用方法为numpy.median(x,[axis]),axis可指定轴方向,默认axis=None,对所有数去中值。
![](/assets/blank.gif)
import numpy as np x = np.array([[1,2,3],[4,5,6]])print(np.median(x)) # 对所有数取中值 # 结果 3.5print(np.median(x,axis=0)) # 沿第一维方向取中值 # 结果 [ 2.5 3.5 4.5]print(np.median(x,axis=1)) # 沿第二维方向取中值 # 结果 [ 2. 5.]
![](/assets/blank.gif)
求和
矩阵求和的函数是sum(),可以对行,列,或整个矩阵求和
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]])print(a.sum()) # 对整个矩阵求和 # 结果 21print(a.sum(axis=0)) # 对行方向求和 # 结果 [5 7 9]print(a.sum(axis=1)) # 对列方向求和 # 结果 [ 6 15]
![](/assets/blank.gif)
累积和
某位置累积和指的是该位置之前(包括该位置)所有元素的和。
例如序列[1,2,3,4,5],其累计和为[1,3,6,10,15],即第一个元素为1,第二个元素为1+2=3,……,第五个元素为1+2+3+4+5=15。
矩阵求累积和的函数是cumsum(),可以对行,列,或整个矩阵求累积和。
![](/assets/blank.gif)
import numpy as npa = np.array([[1,2,3],[4,5,6]])print(a.cumsum()) # 对整个矩阵求累积和 # 结果 [ 1 3 6 10 15 21]print(a.cumsum(axis=0)) # 对行方向求累积和 # 结果 [[1 2 3][5 7 9]]print(a.cumsum(axis=1)) # 对列方向求累积和 # 结果 [[ 1 3 6][ 4 9 15]]
转载于:https://www.cnblogs.com/Archimedes/articles/7353944.html
numpy库:常用基本相关推荐
- astype函数_Python库收录——数值计算numpy库常用函数(含注释)
这里对Numpy库常用的函数进行总结.详细的函数内容可以查阅.Numpy是科学计算库,是一个强大的N维数组对象ndarray,是广播功能函数.其整合C/C++.fortran代码的工具 ,更是Scip ...
- argmax函数_Python之Numpy库常用函数合集(附注释)
文末免费领取[亚马逊.阿里巴巴股票数据] 最近学习Python,才发现原来python里的各种库才是大头!于是乎找了学习资料对Numpy库常用的函数进行总结,并带了注释.在这里分享给大家,对于库的学习 ...
- numpy随机生成01矩阵_Python数据分析Numpy库常用函数详解,提到循环就该想到的库...
Python进行数据分析的核心库肯定是Pandas,该库差不多可以解决结构化数据的绝大部分处理需求.在<Python数据分析常用函数及参数详解,可以留着以备不时之需 >一文中也已经对该库的 ...
- numpy库常用函数——np.random.normal()函数
一.函数语法: np.random.normal(loc=0.0, scale=1.0, size=None) 参数解释: loc(float):此概率分布的均值(对应着整个分布的中心centre s ...
- [基础库Numpy] 技能树总结:一篇文章学会Numpy库
目录 1.基础库的下载 (1)终端窗口下载安装 (2)包管理器搜索下载 2.Numpy库常用API简介 (1)创建数组 #构造法 #特殊值法 #随机数法 #定长分割法 #重复构造法 #网格构造法 (2 ...
- Python -- Numpy库的使用
Numpy库的使用 要点:numpy是用于处理含有同种元素的多维数组运算的第三方库. 1.numpy库概述 Python标准库中提供了一个array类型,用于保存数组类型数据,然而这个类型不支 ...
- Python - Numpy库的使用(简单易懂)
目录 numpy多维数组--数组的创建 1.array函数创建数组对象 2.通过arange.linspace函数创建等差数组对象 3.通过logspace函数创建等比数列数组 函数 zeros on ...
- python数学库的使用方法_浅谈numpy库的常用基本操作方法
NumPy数组的维数称为秩(rank),一维数组的秩为1,二维数组的秩为2,以此类推.在NumPy中,每一个线性的数组称为是一个轴(axes),秩其实是描述轴的数量.比如说,二维数组相当于是一个一维数 ...
- Py之Numpy:Numpy库中常用函数的简介、应用之详细攻略
Py之Numpy:Numpy库中常用函数的简介.应用之详细攻略 目录 Numpy库中常用函数的简介.应用 1.X, Y = np.meshgrid(X, Y) 相关文章 Py之Numpy:Numpy库 ...
最新文章
- CxImage图像处理类库
- linux dns chroot配置文件,linxu chroot DNS 配置
- linux docker 限制容器可用的 cpu
- c语言个人账册报告的课题来源,C语言个人账簿管理系统报告
- 打印一个整数二进制表示中所有的奇数位和偶数位
- ASP.NET Core Blazor Webassembly 之 数据绑定
- 年度影像旗舰vivo X70系列正式发布 售价3699元起
- mysql的增加命令_如何为mysql增加自定义命令
- sammon映射 matlab实现,matlab编写的 32个降维程序
- Android基础之图片的压缩算法
- 拓端tecdat|R语言中固定与随机效应Meta分析 - 效率和置信区间覆盖
- WINRM实现端口复用后门
- [copypaste]各种渐进表示法
- win10前面板耳机无声音解决
- Java 基础类型 —— Short
- 视频教程-nodejs教程-Node.js
- 挑战华为社招:不止面试题,笔记源码统统都有,最强技术实现
- 用pandas进行数据分析实战
- 什么是web前端技术?要学什么?
- FFMPEG音视频开发: 完成摄像头、桌面本地录制与rtmp推流(windows)