1. 程式人生 > >python之scipy庫簡單操作

python之scipy庫簡單操作

Scipy提供了真正的矩陣,以及大量基於矩陣運算的物件與函式,功能包含最優化,線性代數,積分,插值,擬合,特殊函式,快速傅立葉變換,訊號處理和影象處理,常微分方程求解等,Scipy依賴於Numpy庫,因此安裝Scipy時應先安裝Numpy庫,Scipy安裝與其他庫一樣,可通過pip install Scipy安裝,也可以自行下載原始碼,然後用pip install 路徑+檔名全稱(包括.字尾檔名)進行安裝,原始碼下載連結:https://pypi.python.org/pypi/scipy/1.0.0,選擇對應版本下載即可

from scipy.optimize import fsolve #匯入求解方程組的函式
def f(x): #定義求解的方程組
    x1 = x[0]
    x2 = x[1]
    return [2*x1-x2**2-1,x1**2-x2-2]

result = fsolve(f,[1,1]) #設定初始值[1,1]
print(result) #輸出結果

#數值積分運算
from scipy import integrate #匯入積分函式
def g(x):
    return (1-x**2)**0.5

pi_2,err = integrate.quad(g,-1,1) #積分結果和誤差
print(pi_2*2,err) #積分結果為圓周率pi的一半