python使用插值法畫(huà)出平滑曲線(xiàn)
本文實(shí)例為大家分享了python使用插值法畫(huà)出平滑曲線(xiàn)的具體代碼,供大家參考,具體內(nèi)容如下
實(shí)現(xiàn)所需的庫(kù)
numpy、scipy、matplotlib
實(shí)現(xiàn)所需的方法
插值
- nearest:最鄰近插值法
- zero:階梯插值
- slinear:線(xiàn)性插值
- quadratic、cubic:2、3階B樣條曲線(xiàn)插值
擬合和插值的區(qū)別
簡(jiǎn)單來(lái)說(shuō),插值就是根據(jù)原有數(shù)據(jù)進(jìn)行填充,最后生成的曲線(xiàn)一定過(guò)原有點(diǎn)。
擬合是通過(guò)原有數(shù)據(jù),調(diào)整曲線(xiàn)系數(shù),使得曲線(xiàn)與已知點(diǎn)集的差別(最小二乘)最小,最后生成的曲線(xiàn)不一定經(jīng)過(guò)原有點(diǎn)。
代碼實(shí)現(xiàn)
# -*- coding: utf-8 -*- # 調(diào)用模塊 # 調(diào)用數(shù)組模塊 import numpy as np # 實(shí)現(xiàn)插值的模塊 from scipy import interpolate # 畫(huà)圖的模塊 import matplotlib.pyplot as plt # 生成隨機(jī)數(shù)的模塊 import random # random.randint(0, 10) 生成0-10范圍內(nèi)的一個(gè)整型數(shù) # y是一個(gè)數(shù)組里面有10個(gè)隨機(jī)數(shù),表示y軸的值 y = np.array([random.randint(0, 10) for _ in range(10)]) # x是一個(gè)數(shù)組,表示x軸的值 x = np.array([num for num in range(10)]) # 插值法之后的x軸值,表示從0到9間距為0.5的18個(gè)數(shù) xnew = np.arange(0, 9, 0.5) """ kind方法: nearest、zero、slinear、quadratic、cubic 實(shí)現(xiàn)函數(shù)func """ func = interpolate.interp1d(x, y, kind='cubic') # 利用xnew和func函數(shù)生成ynew,xnew的數(shù)量等于ynew數(shù)量 ynew = func(xnew) # 畫(huà)圖部分 # 原圖 plt.plot(x, y, 'ro-') # 擬合之后的平滑曲線(xiàn)圖 plt.plot(xnew, ynew) plt.show()
注意事項(xiàng)/p>
- x, y為原來(lái)的數(shù)據(jù)(少量)
- xnew為一個(gè)數(shù)組,條件:x⊆⊆xnew
- 如:x的最小值為-2.931,最大值為10.312;則xnew的左邊界要小于-2.931,右邊界要大于10.312。當(dāng)然也最好注意一下間距,最好小于x中的精度
- func為函數(shù),里面的參數(shù)x、y、kind,x,y就是原數(shù)據(jù)的x,y,kind為需要指定的方法
- ynew需要通過(guò)xnew數(shù)組和func函數(shù)來(lái)生成
- 理論上x(chóng)new數(shù)組內(nèi)的值越多,生成的曲線(xiàn)越平滑
以上就是本文的全部?jī)?nèi)容,希望對(duì)大家的學(xué)習(xí)有所幫助,也希望大家多多支持腳本之家。
相關(guān)文章
詳解Python中可散列的數(shù)據(jù)類(lèi)型
在Python中,字典(dict)是一種常用的數(shù)據(jù)類(lèi)型,其使用鍵(key)和值(value)來(lái)存儲(chǔ)和訪(fǎng)問(wèn)數(shù)據(jù),在字典中,鍵必須是可哈希(hashable)的類(lèi)型,否則會(huì)導(dǎo)致運(yùn)行時(shí)錯(cuò)誤,因此,在Python中,有一些可散列的數(shù)據(jù)類(lèi)型可以作為字典的鍵2023-06-06pycharm進(jìn)行Git關(guān)聯(lián)和取消方式
這篇文章主要介紹了pycharm進(jìn)行Git關(guān)聯(lián)和取消方式,具有很好的參考價(jià)值,希望對(duì)大家有所幫助,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教2024-06-06使用Python實(shí)現(xiàn)一鍵往Word文檔的表格中填寫(xiě)數(shù)據(jù)
在工作中,我們經(jīng)常遇到將Excel表中的部分信息填寫(xiě)到Word文檔的對(duì)應(yīng)表格中,以生成報(bào)告,方便打印,所以本文小編就給大家介紹了如何使用Python實(shí)現(xiàn)一鍵往Word文檔的表格中填寫(xiě)數(shù)據(jù),文中有詳細(xì)的代碼示例供大家參考,需要的朋友可以參考下2023-12-12Python flask框架如何顯示圖像到web頁(yè)面
這篇文章主要介紹了Python flask框架如何顯示圖像到web頁(yè)面,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-06-06