溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點(diǎn)擊 登錄注冊 即表示同意《億速云用戶服務(wù)條款》

numpy庫常用基本操作方法有哪些

發(fā)布時(shí)間:2021-08-18 14:28:45 來源:億速云 閱讀:204 作者:小新 欄目:開發(fā)技術(shù)

這篇文章主要介紹了numpy庫常用基本操作方法有哪些,具有一定借鑒價(jià)值,感興趣的朋友可以參考下,希望大家閱讀完這篇文章之后大有收獲,下面讓小編帶著大家一起了解一下。

NumPy數(shù)組的維數(shù)稱為秩(rank),一維數(shù)組的秩為1,二維數(shù)組的秩為2,以此類推。在NumPy中,每一個線性的數(shù)組稱為是一個軸(axes),秩其實(shí)是描述軸的數(shù)量。比如說,二維數(shù)組相當(dāng)于是一個一維數(shù)組,而這個一維數(shù)組中每個元素又是一個一維數(shù)組。所以這個一維數(shù)組就是NumPy中的軸(axes),而軸的數(shù)量——秩,就是數(shù)組的維數(shù)。

1、創(chuàng)建矩陣

Numpy庫中的矩陣模塊為ndarray對象,有很多屬性:T,data, dtype,flags,flat,imag,real,size,

itemsize,nbytes,ndim,shape,strides,ctypes,base等等。

1.1采用ndarray對象

import numpy as np #引入numpy庫
#創(chuàng)建一維的narray對象
a = np.array([1,2,3,4,5])
#創(chuàng)建二維的narray對象
a2 = np.array([[1,2,3,4,5],[6,7,8,9,10]])
#創(chuàng)建多維對象以其類推

1.2通過函數(shù)創(chuàng)建矩陣

1.2.1 arange

import numpy as np
a = np.arange(10) # 默認(rèn)從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]

1.2.2 linspace

linspace()和matlab的linspace很類似,用于創(chuàng)建指定數(shù)量等間隔的序列,實(shí)際生成一個等差數(shù)列。

import numpy as np
a = np.linspace(0,10,5) # 生成首位是0,末位是10,含5個數(shù)的等差數(shù)列
print(a)

1.2.3 logspace


linspace用于生成等差數(shù)列,而logspace用于生成等比數(shù)列。

下面的例子用于生成首位是100,末位是102,含5個數(shù)的等比數(shù)列

import numpy as np
a = np.logspace(0,2,5)
print(a)

1.2.4 ones,zeros,eye,empty

ones創(chuàng)建全1矩陣
zeros創(chuàng)建全0矩陣
eye創(chuàng)建單位矩陣
empty創(chuàng)建空矩陣(實(shí)際有值)

import numpy as np
a_ones = np.ones((3,4)) # 創(chuàng)建3*4的全1矩陣
print(a_ones)
# 結(jié)果
[[ 1. 1. 1. 1.]
 [ 1. 1. 1. 1.]
 [ 1. 1. 1. 1.]]
a_zeros = np.zeros((3,4)) # 創(chuàng)建3*4的全0矩陣
print(a_zeros)
# 結(jié)果
[[ 0. 0. 0. 0.]
[ 0. 0. 0. 0.]
[ 0. 0. 0. 0.]]
a_eye = np.eye(3) # 創(chuàng)建3階單位矩陣
print(a_eye)
# 結(jié)果
[ 1. 0. 0.]
[ 0. 1. 0.]
[ 0. 0. 1.]]
a_empty = np.empty((3,4)) # 創(chuàng)建3*4的空矩陣 
print(a_empty)
# 結(jié)果
[[ 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]]

1.2.5 fromstring

fromstring()方法可以將字符串轉(zhuǎn)化成ndarray對象,需要將字符串?dāng)?shù)字化時(shí)這個方法比較有用,可以獲得字符串的ascii碼序列。

a = "abcdef"
b = np.fromstring(a,dtype=np.int8) # 因?yàn)橐粋€字符為8位,所以指定dtype為np.int8
print(b) # 返回 [ 97 98 99 100 101 102]

1.2.6 fromfunction

fromfunction()方法可以根據(jù)矩陣的行號列號生成矩陣的元素。

例如創(chuàng)建一個矩陣,矩陣中的每個元素都為行號和列號的和。

import numpy as np
def func(i,j): 
 return i+j
a = np.fromfunction(func,(5,6)) 
# 第一個參數(shù)為指定函數(shù),第二個參數(shù)為列表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開始的

2.矩陣的操作

numpy中的ndarray對象重載了許多運(yùn)算符,使用這些運(yùn)算符可以完成矩陣間對應(yīng)元素的運(yùn)算。

如 +,-,*,/,%,**

2.1矩陣函數(shù)操作

2.1.1 常用數(shù)學(xué)運(yùn)算(導(dǎo)入numpy模塊:import numpy as np)

矩陣函數(shù)說明
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中每個元素取指數(shù)函數(shù),ex
np.sqrt(a)對矩陣a中每個元素開根號√x

如:求矩陣每個元素的sin值

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(np.sin(array))
#結(jié)果
# [[ 0.84147098 0.90929743 0.14112001]
# [-0.7568025 -0.95892427 -0.2794155 ]]

2.1.2矩陣乘法(點(diǎn)乘)

矩陣乘法必須滿足矩陣乘法的條件,即第一個矩陣的列數(shù)等于第二個矩陣的行數(shù)。

矩陣乘法的函數(shù)為 dot

 import numpy as np
a1 = 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)相當(dāng)于matlab中的a1*a2
# 而python中的a1*a2相當(dāng)于matlab中的a1.*a2
# 結(jié)果
[[22 28]
[49 64]]

2.1.3矩陣的轉(zhuǎn)置(transpose或T)

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(array.transpose())#或np.transpose(array)
print(array.T)
# 結(jié)果
# [[1 4]
# [2 5]
# [3 6]]
#      transpose與T效果一樣
# [[1 4]
# [2 5]
# [3 6]]

2.1.4矩陣的逆

求矩陣的逆需要先導(dǎo)入numpy.linalg,用linalg的inv函數(shù)來求逆。

矩陣求逆的條件是矩陣的行數(shù)和列數(shù)相同。

import numpy as np
import numpy.linalg as lg
a = np.array([[1,2,3],[4,5,6],[7,8,9]])
print(lg.inv(a))
# 結(jié)果
[[ -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)) # 單位矩陣的逆為他本身
# 結(jié)果
[[ 1. 0. 0.]
[ 0. 1. 0.]
[ 0. 0. 1.]]

2.1.5最大值最小值

獲得矩陣中元素最大最小值的函數(shù)分別是max和min,可以獲得整個矩陣、行或列的最大最小值。

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(array.max())#結(jié)果為6
print(array.min())#結(jié)果為1
#同時(shí)還可以指定axis關(guān)鍵字,獲取行或列的最大、最小值
print(array.max(axis=0)) #x軸最大值,0,1分別代表行列

2.1.6平均值

獲得矩陣中元素的平均值可以通過函數(shù)mean()或average()。同樣地,可以獲得整個矩陣、行或列的平均值

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(array.mean())#結(jié)果為3.5
print(np.average(array))#結(jié)果為3.5
print(array.mean(axis=0))#行方向的平均值,同樣,0,1代表維度

2.1.7方差

方差的函數(shù)為var(),方差函數(shù)var()相當(dāng)于函數(shù)mean(abs(x - x.mean())**2),其中x為矩陣。

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(array.var())#結(jié)果為2.91666666667
#同樣可通過axis指定維度0,1分別代表行列,
print(array.var(axis=0))
# 結(jié)果
# [ 2.25 2.25 2.25]

2.1.8標(biāo)準(zhǔn)差

標(biāo)準(zhǔn)差的函數(shù)為std()。

std()相當(dāng)于sqrt(mean(abs(x - x.mean())**2)),或相當(dāng)于sqrt(x.var())。

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(array.std())#結(jié)果為1.70782512766
#同樣可通過axis指定維度0,1分別代表行列,
print(array.std(axis=0))
# 結(jié)果
# [ 1.5 1.5 1.5]

2.1.9中位數(shù)or中值

中值指的是將序列按大小順序排列后,排在中間的那個值,如果有偶數(shù)個數(shù),則是排在中間兩個數(shù)的平均值。

例如序列[5,2,6,4,2],按大小順序排成 [2,2,4,5,6],排在中間的數(shù)是4,所以這個序列的中值是4。

又如序列[5,2,6,4,3,2],按大小順序排成 [2,2,3,4,5,6],因?yàn)橛信紨?shù)個數(shù),排在中間兩個數(shù)是3、4,所以這個序列中值是3.5。

中值的函數(shù)是median(),調(diào)用方法為numpy.median(x,[axis]),axis可指定軸方向,默認(rèn)axis=None,對所有數(shù)去中值

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(np.median(array))#結(jié)果:3.5
#指定維度
print(np.median(array,axis=0))
#結(jié)果
# [ 2.5 3.5 4.5]

2.1.10 求和

矩陣求和的函數(shù)是sum(),可以對行,列,或整個矩陣求和

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(array.sum())#結(jié)果:21
#指定維度
print(array.sum(axis=0))
#結(jié)果
# [5 7 9]

2.1.11 累計(jì)和或累加

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

例如序列[1,2,3,4,5],其累計(jì)和為[1,3,6,10,15],即第一個元素為1,第二個元素為1+2=3,……,第五個元素為1+2+3+4+5=15。

矩陣求累積和的函數(shù)是cumsum(),可以對行,列,或整個矩陣求累積和。

import numpy as np
a = np.array([[1,2,3],[4,5,6]])
print(a.cumsum())      # 對整個矩陣求累積和
# 結(jié)果 [ 1 3 6 10 15 21]
print(a.cumsum(axis=0)) # 對行方向求累積和
# 結(jié)果
[[1 2 3]
[5 7 9]]
print(a.cumsum(axis=1)) # 對列方向求累積和
# 結(jié)果
[[ 1 3 6]
[ 4 9 15]]

2.2矩陣的截取

2.2.1 按行列截取

矩陣的截取和list相同,可以通過[](方括號)來截取

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 10]
print(a[1,:]) #截取第二行,返回 [ 6 7 8 9 10]

2.2.2按條件截取

import numpy as np
a = np.array([[1,2,3,4,5],[6,7,8,9,10]])
b = a[a>6] # 截取矩陣a中大于6的元素,范圍的是一維數(shù)組
print(b) # 返回 [ 7 8 9 10]
# 其實(shí)布爾語句首先生成一個布爾矩陣,將布爾矩陣傳入[](方括號)實(shí)現(xiàn)截取
print(a>6) 
# 返回
[[False False False False False]
[False True True True True]]

按條件截取應(yīng)用較多的是對矩陣中滿足一定條件的元素變成特定的值。

例如將矩陣中大于6的元素變成0。

import numpy as np
a = 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]]

2.2.3 clip截取

clip(矩陣,min,max)#返回值:所有小于min的值都等于min,所有大于max的值都等于max

import numpy as np
array=np.array([[1,2,3],
        [4,5,6]])
print(array.clip(2,4))
#結(jié)果
# [[2 2 3]
# [4 4 4]]

2.3 矩陣的合并

矩陣的合并可以通過numpy中的hstack方法和vstack方法實(shí)現(xiàn)

import numpy as np
a1 = np.array([[1,2],[3,4]])
a2 = np.array([[5,6],[7,8]])
#!注意 參數(shù)傳入時(shí)要以列表list或元組tuple的形式傳入
print(np.hstack([a1,a2])) 
#橫向合并,返回結(jié)果如下 
[[1 2 5 6]
[3 4 7 8]]
print(np.vstack((a1,a2)))
#縱向合并,返回結(jié)果如下
[[1 2]
[3 4]
[5 6]
[7 8]]

矩陣的合并也可以通過concatenatef方法。

np.concatenate( (a1,a2), axis=0 ) 等價(jià)于 np.vstack( (a1,a2) )

np.concatenate( (a1,a2), axis=1 ) 等價(jià)于 np.hstack( (a1,a2) )

感謝你能夠認(rèn)真閱讀完這篇文章,希望小編分享的“numpy庫常用基本操作方法有哪些”這篇文章對大家有幫助,同時(shí)也希望大家多多支持億速云,關(guān)注億速云行業(yè)資訊頻道,更多相關(guān)知識等著你來學(xué)習(xí)!

向AI問一下細(xì)節(jié)

免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。

AI