numpy中的matrix與array的區別
阿新 • • 發佈:2019-01-12
Numpy matrices必須是2維的,但是 numpy arrays (ndarrays) 可以是多維的(1D,2D,3D····ND). Matrix是Array的一個小的分支,包含於Array。所以matrix 擁有array的所有特性。
在numpy中matrix的主要優勢是:相對簡單的乘法運算子號。例如,a和b是兩個matrices,那麼a*b,就是矩陣積。而不用np.dot()。如:
import numpy as np a=np.mat('4 3; 2 1') b=np.mat('1 2; 3 4') print(a) # [[4 3] # [2 1]] print(b) # [[1 2] # [3 4]] print(a*b) # [[13 20] # [ 5 8]]
matrix 和 array 都可以通過objects後面加.T
得到其轉置。但是 matrix objects 還可以在後面加 .H
f得到共軛矩陣, 加 .I
得到逆矩陣。
相反的是在numpy裡面arrays遵從逐個元素的運算,所以對於array:c 和d的c*d運算相當於matlab裡面的c.*d運算。
c=np.array([[4, 3], [2, 1]]) d=np.array([[1, 2], [3, 4]]) print(c*d) # [[4 6] # [6 4]]
而矩陣相乘,則需要numpy裡面的dot命令 :
print(np.dot(c,d)) # [[13 20] # [ 5 8]]
但是python中矩陣沒有MATLAB中的.*這個性質,對於matrix對應乘法得用np.multiply()
>>> np.multiply(a,b) matrix([[ 4, 6], [6, 4]])
當然 **
運算子的作用也不一樣 :
print(a**2) #矩陣乘法 # [[22 15] # [10 7]] print(c**2) #對應乘法 # [[16 9] # [ 4 1]]
問題就出來了,如果一個程式裡面既有matrix 又有array,會讓人腦袋大。但是如果只用array,你不僅可以實現matrix所有的功能,還減少了程式設計和閱讀的麻煩。
當然你可以通過下面的兩條命令輕鬆的實現兩者之間的轉換:np.asmatrix
和np.asarray
對我來說,numpy 中的array與numpy中的matrix,matlab中的matrix的最大的不同是,在做歸約運算時,array的維數會發生變化,但matrix總是保持為2維。例如下面求平均值的運算
>>> m = np.mat([[1,2],[2,3]]) >>> m matrix([[1, 2], [2, 3]]) >>> mm = m.mean(1) >>> mm matrix([[ 1.5], [ 2.5]]) >>> mm.shape (2, 1) >>> m - mm matrix([[-0.5, 0.5], [-0.5, 0.5]])
對array 來說
>>> a = np.array([[1,2],[2,3]]) >>> a array([[1, 2], [2, 3]]) >>> am = a.mean(1) >>> am.shape (2,) >>> am array([ 1.5, 2.5]) >>> a - am #wrong array([[-0.5, -0.5], [ 0.5, 0.5]]) >>> a - am[:, np.newaxis] #right array([[-0.5, 0.5], [-0.5, 0.5]])