python多項(xiàng)式擬合之np.polyfit 和 np.polyld詳解
python數(shù)據(jù)擬合主要可采用numpy庫(kù),庫(kù)的安裝可直接用pip install numpy
等。
1. 原始數(shù)據(jù):假如要擬合的數(shù)據(jù)yyy來(lái)自sin函數(shù),np.sin
import numpy as np import matplotlib.pyplot as plt xxx = np.arange(0, 1000) # x值,此時(shí)表示弧度 yyy = np.sin(xxx*np.pi/180) #函數(shù)值,轉(zhuǎn)化成度
2. 測(cè)試不同階的多項(xiàng)式,例如7階多項(xiàng)式擬合,使用np.polyfit擬合,np.polyld得到多項(xiàng)式系數(shù)
z1 = np.polyfit(xxx, yyy, 7) # 用7次多項(xiàng)式擬合,可改變多項(xiàng)式階數(shù); p1 = np.poly1d(z1) #得到多項(xiàng)式系數(shù),按照階數(shù)從高到低排列 print(p1) #顯示多項(xiàng)式
3. 求對(duì)應(yīng)xxx的各項(xiàng)擬合函數(shù)值
yvals=p1(xxx) # 可直接使用yvals=np.polyval(z1,xxx)
4. 繪圖如下
plt.plot(xxx, yyy, '*',label='original values') plt.plot(xxx, yvals, 'r',label='polyfit values') plt.xlabel('x axis') plt.ylabel('y axis') plt.legend(loc=4) # 指定legend在圖中的位置,類(lèi)似象限的位置 plt.title('polyfitting') plt.show()
5. np.polyfit函數(shù):采用的是最小二次擬合,numpy.polyfit(x, y, deg, rcond=None, full=False, w=None, cov=False),前三個(gè)參數(shù)是必須的
官方文檔:https://docs.scipy.org/doc/numpy-1.13.0/reference/generated/numpy.polyfit.html
6. np.polyld函數(shù):得到多項(xiàng)式系數(shù),主要有三個(gè)參數(shù)
A one-dimensional polynomial class. A convenience class, used to encapsulate "natural" operations on polynomials so that said operations may take on their customary form in code (see Examples). Parameters ---------- c_or_r : array_like The polynomial's coefficients, in decreasing powers, or if the value of the second parameter is True, the polynomial's roots (values where the polynomial evaluates to 0). For example, ``poly1d([1, 2, 3])`` returns an object that represents :math:`x^2 + 2x + 3`, whereas ``poly1d([1, 2, 3], True)`` returns one that represents :math:`(x-1)(x-2)(x-3) = x^3 - 6x^2 + 11x -6`. r : bool, optional If True, `c_or_r` specifies the polynomial's roots; the default is False. variable : str, optional Changes the variable used when printing `p` from `x` to `variable` (see Examples).
參數(shù)1表示:在沒(méi)有參數(shù)2(也就是參數(shù)2默認(rèn)False時(shí)),參數(shù)1是一個(gè)數(shù)組形式,且表示從高到低的多項(xiàng)式系數(shù)項(xiàng),例如參數(shù)1為[4,5,6]表示:
參數(shù)2表示:為T(mén)rue時(shí),表示將參數(shù)1中的參數(shù)作為根來(lái)形成多項(xiàng)式,即參數(shù)1為[4,5,6]時(shí)表示:(x-4)(x-5)(x-6)=0,也就是:
參數(shù)3表示:換參數(shù)標(biāo)識(shí),用慣了x,可以用 t,s之類(lèi)的
用法:
1. 直接進(jìn)行運(yùn)算,例如多項(xiàng)式的平方,分別得到
xx=np.poly1d([1,2,3]) print(xx) yy=xx**2 #求平方,或者用 xx * xx print(yy)
2. 求值:
yy(1) = 36
3. 求根:即等式為0時(shí)的未知數(shù)值
yy.r
4. 得到系數(shù)形成數(shù)組:
yy.c 為:array([ 1, 4, 10, 12, 9])
5. 返回最高次冪數(shù):
yy.order = 4
6. 返回系數(shù):
yy[0] —— 表示冪為0的系數(shù)
yy[1] —— 表示冪為1的系數(shù)
總結(jié)
以上所述是小編給大家介紹的python多項(xiàng)式擬合之np.polyfit 和 np.polyld詳解,希望對(duì)大家有所幫助,也非常感謝大家對(duì)腳本之家網(wǎng)站的支持!
- python 對(duì)任意數(shù)據(jù)和曲線進(jìn)行擬合并求出函數(shù)表達(dá)式的三種解決方案
- python之?dāng)M合的實(shí)現(xiàn)
- python實(shí)現(xiàn)最小二乘法線性擬合
- python如何實(shí)現(xiàn)數(shù)據(jù)的線性擬合
- 在python中利用最小二乘擬合二次拋物線函數(shù)的方法
- 對(duì)python指數(shù)、冪數(shù)擬合curve_fit詳解
- Python 做曲線擬合和求積分的方法
- Python應(yīng)用實(shí)現(xiàn)雙指數(shù)函數(shù)及擬合代碼實(shí)例
相關(guān)文章
Pycharm 實(shí)現(xiàn)下一個(gè)文件引用另外一個(gè)文件的方法
今天小編就為大家分享一篇Pycharm 實(shí)現(xiàn)下一個(gè)文件引用另外一個(gè)文件的方法,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-01-01python實(shí)現(xiàn)柵欄加解密 支持密鑰加密
這篇文章主要為大家詳細(xì)介紹了python實(shí)現(xiàn)柵欄加解密,支持密鑰加密,具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2019-03-03Python按條件篩選、剔除表格數(shù)據(jù)并繪制剔除前后的直方圖(示例代碼)
本文介紹基于Python語(yǔ)言,讀取Excel表格文件數(shù)據(jù),以其中某一列數(shù)據(jù)的值為標(biāo)準(zhǔn),對(duì)于這一列數(shù)據(jù)處于指定范圍的所有行,再用其他幾列數(shù)據(jù)的數(shù)值,加以數(shù)據(jù)篩選與剔除,感興趣的朋友跟隨小編一起看看吧2024-07-07用Python實(shí)現(xiàn)一個(gè)簡(jiǎn)單的線程池
這篇文章主要介紹了用Python實(shí)現(xiàn)一個(gè)簡(jiǎn)單的線程池,通過(guò)這個(gè)小程序可以幫助更好地理解Python中線程的運(yùn)行機(jī)制,需要的朋友可以參考下2015-04-04淺談Python中range與Numpy中arange的比較
這篇文章主要介紹了淺談Python中range與Numpy中arange的比較,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-03-03python numpy--數(shù)組的組合和分割實(shí)例
這篇文章主要介紹了python numpy--數(shù)組的組合和分割實(shí)例,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-02-02基于python實(shí)現(xiàn)藍(lán)牙通信代碼實(shí)例
這篇文章主要介紹了基于python實(shí)現(xiàn)藍(lán)牙通信代碼實(shí)例,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-11-11在Python web中實(shí)現(xiàn)驗(yàn)證碼圖片代碼分享
這篇文章主要介紹了在Python web中實(shí)現(xiàn)驗(yàn)證碼圖片代碼分享,具有一定參考價(jià)值,需要的朋友可以了解下。2017-11-11