Python數學計算:Numpy和Scipy(矩陣相關)
阿新 • • 發佈:2019-01-28
Python是一種通用語言。它被解釋執行,是動態型別語言,並且非常適合互動工作和快速實現原型,然而又足夠強大用來寫大型應用。
NumPy是一個定義了數值陣列和矩陣型別和它們的基本運算的語言擴充套件。
SciPy是另一種使用NumPy來做高等數學、訊號處理、優化、統計和許多其它科學任務的語言擴充套件。
(1)建立矩陣
import numpy as np
s = np.array([[-0.707,-0.707],[-0.707,0.707]])
d = np.array([[0.02,0],[0,0]])
(2)矩陣乘積
g = np.dot(s,d)
print g
(3)矩陣SVD分解
A=np.array([[3 ,3],[4,4]])
#print A
U,sigma,VT=la.svd(A)
print U
print sigma
print VT
(4)解方程組
前提:矩陣A的行列式值不能為0
a = np.array([[2,6,12,20],[6,20,42,72], [12,42,90,156],[20,72,156,272]]) print np.linalg.det(a) b = np.array([0,0,0,0]) print b x1 = np.linalg.solve(a,b) print x1