ITPub博客

首页 > Linux操作系统 > Linux操作系统 > Numpy初步

Numpy初步

原创 Linux操作系统 作者:Tybyq 时间:2018-11-28 16:17:53 0 删除 编辑

1,获取矩阵行列数

Import numpyasnp

#创建二维的naaray对象

a=np.array([[1,2,3,4,5],[6,7,8,9,10]])

print(a.shape)   #返回一个形状,是一个tuple

print(a.shape[0])#获得行数,试想如果是多维的呢,所以你就会明白为什么是[0]

print(a.shape[1])   #获得列数

2,矩阵的截取

importnumpyasnp

#创建二维的naaray对象

a=np.array([[1,2,3,4,5],[6,7,8,9,10]])

print(a[0:1])#这里不是很懂,看下面

print(a[1,2:4])#返回[89],返回第二行2-3个数

print(a[1,2:5])#返回[8910]证明了取第二行的2-4个数

3,按条件截取

importnumpyasnp

#创建二维的naaray对象

a=np.array([[1,2,3,4,5],[6,7,8,9,10]])

b=a[a>6]#截取矩阵a中大于6的元素,范围的是一维数组

print(b)

print(a>6)#其实布尔语句首先生成一个布尔矩阵,将布尔矩阵传入[](方括号)实现截取

4,满足一定条件的元素变成特定的值

importnumpyasnp

#创建二维的naaray对象

a=np.array([[1,2,3,4,5],[6,7,8,9,10]])

print(a)

#大于6清零后矩阵为

a[a>6]=0

print(a)

结果分别为:

[[12345]

[678910]]

[[12345]

[60000]]

5,矩阵合并

importnumpyasnp

 

a1=np.array([[1,2],[3,4]])

a2=np.array([[5,6],[7,8]])

 

print(np.hstack([a1,a2]))#Horizontal

print(a1)#所以知道为什么是1256,2478了

print(np.vstack((a1,a2)))#vertical

np.concatenate((a1,a2),axis=0)#等价于np.vstack((a1,a2))

np.conca

6,通过函数创建矩阵

numpy自带创建narray对象的函数,可以方便的创建常用的有规律的矩阵

importnumpyasnp

a=np.arange(10)#默认从0开始到10(不包括10),步长为1

print(a)#返回[0123456789]

a1=np.arange(5,10)#从5开始到10(不包括10),步长为1

print(a1)#返回[56789]

a2=np.arange(5,20,2)#从5开始到20(不包括20),步长为2

print(a2)#返回[5791113151719]

 

7,linspace

创建指定数量等间隔的序列,实际生成一个等差数列

importnumpyasnp

 

a=np.linspace(0,10,7)#生成首位是0,末位是10,含7个数的等差数列

print(a)

8,logspace

logspace用于生成等比数列。 

importnumpyasnp

 

a=np.logspace(0,10,7)#生成首位是10**0,末位是10**4,含5个数的等比数列

print(a)

9,ones,zeros,eye,empty

ones创建全1矩阵 

zeros创建全0矩阵 

eye创建单位矩阵 

empty创建空矩阵(实际有值)

importnumpyasnp

 

one=np.ones((3,4))#创建3*4的全1矩阵

print(one)

 

zero=np.zeros((3,4))#创建3*4的全0矩阵

print(zero)

 

eye=np.eye(5)#创建5阶单位矩阵

print(eye)

 

empty=np.empty((3,4))#创建3*4的空矩阵(实际有值)

print(empty)

10,fromstring ——获得字符ASCII码

fromstring()方法可以将字符串转化成ndarray对象,需要将字符串数字化时这个方法比较有用,可以获得字符串的ascii码序列,转成相应字符的阿斯卡码。

importnumpyasnp

 

a="abcdef"

b=np.fromstring(a,dtype=np.int8)#因为一个字符为8位,所以指定dtype为np.int8

print(b)#返回[979899100101102]

11,fromfunction

fromfunction()方法可以根据矩阵的行号列号生成矩阵的元素。 

例如创建一个矩阵,矩阵中的每个元素都为行号和列号的和。

importnumpyasnp

 

deffunc(i,j):

returni+j     #这里也可以是别的,比如再加一个9

 

a=np.fromfunction(func,(5,6))

#函数定义就是如此,第一个参数为指定函数,第二个参数为列表list或元组tuple,说明矩阵的大小

print(a)

12,常用矩阵函数

同样地,numpy中也定义了许多函数,使用这些函数可以将函数作用于矩阵中的每个元素。 

表格中默认导入了numpy模块,即 import numpy asnp

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

importnumpyasnp

 

a=np.array([[1,2,3],[4,5,6]])

print(np.sin(a))

 

#结果

[[0.841470980.909297430.14112001]

[-0.7568025-0.95892427-0.2794155]]

 

print(np.arcsin(a))

 

#结果

#RuntimeWarning:invalidvalueencounteredinarcsin

print(np.arcsin(a))

[[1.57079633nannan]   #nan是not a number 的意思

[nannannan]]

13,矩阵乘法(点乘)

条件:第一个矩阵的列数等于第二个矩阵的行数,函数为dot

importnumpyasnp

 

a1=np.array([[1,2,3],[4,5,6]])

a2=np.array([[1,2],[3,4],[5,6]])

ifa1.shape[1]==a2.shape[0]:#列数等于行数的话

print(a1.dot(a2))

 

14,矩阵的转置

transpose函数

importnumpyasnp

 

a=np.array([[1,2,3],[4,5,6]])

print(a.transpose())

15,矩阵的逆

求矩阵的逆需要先导入numpy.linalg,用linalg的inv函数来求逆。 

矩阵求逆的条件是矩阵的行数和列数相同。

importnumpyasnp

importnumpy.linalgaslg

 

a=np.array([[1,2,3],[4,5,6],[7,8,9]])

print(lg.inv(a))

 

#结果

[[-4.50359963e+159.00719925e+15-4.50359963e+15]

[9.00719925e+15-1.80143985e+169.00719925e+15]

[-4.50359963e+159.00719925e+15-4.50359963e+15]]

 

a=np.eye(3)#3阶单位矩阵

print(lg.inv(a))#单位矩阵的逆为他本身

 

#结果

[[1.0.0.]

[0.1.0.]

[0.0.1.]]

16,矩阵信息获取(如平均值)

获得矩阵中元素最大最小值的函数分别是max和min,可以获得整个矩阵、行或列的最大最小值。

importnumpyasnp

a=np.array([[1,3,9],[1,5,6]])

print(a.max())

print(a.min())

print(a.max(axis=0))#[456]axis=0行方向最大(小)值,即获得每列的最大(小)值

print(a.min(axis=1))#[14]axis=1列方向最大(小)值

#要想获得最大最小值元素所在的位置,可以通过argmax函数来获得

print(a.argmax(axis=1))

17,平均值mean()

获得矩阵中元素的平均值可以通过函数mean()。同样地,可以获得整个矩阵、行或列的平均值

importnumpyasnp

 

a=np.array([[1,2,3],[4,5,6]])

print(a.mean())#结果为:3.5

 

#同样地,可以通过关键字axis参数指定沿哪个方向获取平均值

print(a.mean(axis=0))#结果[2.53.54.5]

print(a.mean(axis=1))#结果[2.5.]

18,方差var()

importnumpyasnp

a=np.array([[1,2,3],[4,5,6]])

print(a.var())

print(a.var(axis=0))

print(a.var(axis=1))

19,标准差std()

importnumpyasnp

a=np.array([[1,2,3],[4,5,6]])

print(a.std())

print(a.std(axis=0))

print(a.std(axis=1))

20,中值median()

调用方法是numpy.median(x,[axis]),axis可指定轴方向,默认为axis=none,对所有数取中值

importnumpyasnp

x=np.array([[1,2,3],[4,5,6]])

print(np.median(x))#对所有数取中值

print(np.median(x,axis=0))#沿第一维方向取中值

print(np.median(x,axis=1))#沿第二维方向取中值

21,求和sum()

importnumpyasnp

a=np.array([[1,2,3],[4,5,6]])

print(a.sum())#对整个矩阵求和

print(a.sum(axis=0))#对行方向求和

print(a.sum(axis=1))#对列方向求和

22,累积和cussum()

某位置累积和指的是该位置之前(包括该位置)所有元素的和。

例如序列[1,2,3,4,5],其累计和为[1,3,6,10,15],即第一个元素为1,第二个元素为1+2=3,……,第五个元素为1+2+3+4+5=15。

矩阵求累积和的函数是cumsum(),可以对行,列,或整个矩阵求累积和。

importnumpyasnp

a=np.array([[1,2,3],[4,5,6]])

print(a.cumsum())#对整个矩阵求累积和

print(a.cumsum(axis=0))#对行方向求累积和

来自 “ ITPUB博客 ” ,链接:http://blog.itpub.net/31557424/viewspace-2221982/,如需转载,请注明出处,否则将追究法律责任。

请登录后发表评论 登录
全部评论

注册时间:2018-10-31

  • 博文量
    166
  • 访问量
    56212