Numpy初步

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))#对行方向求累积和