如何利用Python擬合函數(shù)曲線詳解
使用Python擬合函數(shù)曲線需要用到一些第三方庫：
- numpy：科學(xué)計算的基礎(chǔ)庫（例如：矩陣）
- matplotlib：繪圖庫
- scipy：科學(xué)計算庫
如果沒有安裝過這些庫，需要在命令行中輸入下列代碼進行安裝：
pip install numpy matplotlib scipy
擬合多項式
''' Author: CloudSir Date: 2021-08-01 13:40:50 LastEditTime: 2021-08-02 09:41:54 LastEditors: CloudSir Description: Python擬合多項式 https://github.com/cloudsir ''' import matplotlib.pyplot as plt import numpy as np x = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y = [2.83, 9.53, 14.52, 21.57, 38.26, 53.92, 73.15, 101.56, 129.54, 169.75, 207.59] z1 = np.polyfit(x, y, 3) #用3次多項式擬合，輸出系數(shù)從高到0 p1 = np.poly1d(z1) #使用次數(shù)合成多項式 y_pre = p1(x) plt.plot(x,y,'.') plt.plot(x,y_pre) plt.show()
函數(shù)說明
np.polyfit(x, y, n)
功能：擬合曲線
參數(shù)：
- x，y：x和y的原始數(shù)據(jù)
- n：要擬合的次數(shù)
返回值：
- 一個列表，擬合出的系數(shù)，順序為從高到底
例： n=3時，會利用
a x 3 + b x 2 + c x + d
擬合函數(shù)，并返回擬合出的系數(shù) [a, b, c, d]????
np.poly1d(li, r=False)
功能：生成多項式函數(shù)
參數(shù)：
li:
當(dāng)沒有r參數(shù)或 r=False 時，傳入一個系數(shù)列表(次數(shù)從高到低)，利用該列表生成多項式函數(shù)并返回
import numpy as np
f = np.poly1d([2, 3, 4])
“”" f ( x ) = 2 x 2 + 3 x + 4
"""print(f(2)) # 18
當(dāng)參數(shù) r=True 時，傳入一個根列表，利用該列表生成多項式函數(shù)并返回
import numpy as np
f = np.poly1d([2, 3, 4], True)
“”“ f ( x ) = ( x − 2 ) ∗ ( x − 3 ) ∗ ( x − 4 )
”""print(f(0)) # -24
返回值：
見上
擬合任意函數(shù)
'''
Author: CloudSir
Date: 2021-08-03 15:01:17
LastEditTime: 2021-08-03 15:26:05
LastEditors: CloudSir
Description: Python擬合任意函數(shù)
https://github.com/cloudsir
'''
# 引用庫函數(shù)
import numpy as np
import matplotlib.pyplot as plt
from scipy import optimize as op
plt.rcParams['font.sans-serif'] = ['SimHei'] # 用來正常顯示中文
plt.rcParams['axes.unicode_minus'] = False # 用來正常顯示負號
# 需要擬合的函數(shù)
def f_1(x, A, B, C):
return A * x**2 + B * x + C
# 需要擬合的數(shù)據(jù)組
x_group = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
y_group = [2.83, 9.53, 14.52, 21.57, 38.26, 53.92, 73.15, 101.56, 129.54, 169.75, 207.59]
# 得到返回的A，B值
A, B, C = op.curve_fit(f_1, x_group, y_group)[0]
# 數(shù)據(jù)點與原先的進行畫圖比較
plt.scatter(x_group, y_group, marker='o',label='真實值')
x = np.arange(0, 15, 0.01)
y = A * x**2 + B *x + C
plt.plot(x, y,color='red',label='擬合曲線')
plt.legend() # 顯示label
plt.show()
函數(shù)說明
op.curve_fit(f, x, y)
功能：擬合任意函數(shù)
參數(shù)：
f：要擬合的函數(shù)類型
# 構(gòu)建一個二次函數(shù)
def f(x, A, B, C):
return A * x**2 + B * x + C
op.curve_fit(f, x, y) # 進行擬合
x, y：x和y的原始數(shù)據(jù)
返回值：一個元組 (popt,pcov)
- popt是一個一維數(shù)組，表示得到的擬合方程的參數(shù)。
- pcov是一個二維數(shù)組，是在popt參數(shù)下得到的協(xié)方差。
總結(jié)
到此這篇關(guān)于如何利用Python擬合函數(shù)曲線的文章就介紹到這了,更多相關(guān)Python擬合函數(shù)曲線內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
今天小編就為大家分享一篇python 讀取文件并替換字段的實例，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2018-07-07
事實上，Python只是規(guī)定，無論是構(gòu)造方法還是實例方法，最少要包含一個參數(shù)，并沒有規(guī)定該參數(shù)的具體名稱。之所以將其命名為self，只是程序員之間約定俗成的一種習(xí)慣，遵守這個約定，可以使我們編寫的代碼具有更好的可讀性2023-01-01
本文主要介紹了python使用opencv實現(xiàn)人臉識別的相關(guān)資料。具有很好的參考價值。下面跟著小編一起來看下吧2017-04-04
使用Python編程分析火爆全網(wǎng)的魷魚游戲豆瓣影評
本文來為大家介紹如何使用Python爬取影評的操作，主要是爬取《魷魚游戲》在豆瓣上的一些影評，對數(shù)據(jù)做一些簡單的分析，用數(shù)據(jù)的角度重新審視下這部劇，有需要的朋友可以借鑒參考下2021-10-10
這篇文章主要介紹了Python實現(xiàn)簡單字典樹的方法,實例分析了Python字典樹的定義、實現(xiàn)與使用技巧,需要的朋友可以參考下2016-04-04
這篇文章主要介紹了Sanic框架流式傳輸操作,結(jié)合實例形式分析了Sanic通過流請求與響應(yīng)傳輸操作相關(guān)實現(xiàn)技巧與注意事項,需要的朋友可以參考下2018-07-07
本文主要介紹了Django配置kafka消息隊列的實現(xiàn)，文中通過示例代碼介紹的非常詳細，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-05-05
今天小編就為大家分享一篇淺談python 中類屬性共享的問題，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2019-07-07
這篇文章主要介紹了python RC4加密操作,結(jié)合實例形式分析了python實現(xiàn)RC4加密功能的具體操作步驟與相關(guān)問題解決方法,需要的朋友可以參考下2019-09-09
這篇文章主要介紹了Python代碼庫之Tuple如何append添加元素問題，具有很好的參考價值，希望對大家有所幫助。如有錯誤或未考慮完全的地方，望不吝賜教2023-01-01