Python線性擬合實(shí)現(xiàn)函數(shù)與用法示例
本文實(shí)例講述了Python線性擬合實(shí)現(xiàn)函數(shù)與用法。分享給大家供大家參考，具體如下：
1. 參考別人寫的：
#-*- coding:utf-8 -*-
import math
import matplotlib.pyplot as plt
def linefit(x , y):
N = float(len(x))
sx,sy,sxx,syy,sxy=0,0,0,0,0
for i in range(0,int(N)):
sx += x[i]
sy += y[i]
sxx += x[i]*x[i]
syy += y[i]*y[i]
sxy += x[i]*y[i]
a = (sy*sx/N -sxy)/( sx*sx/N -sxx)
b = (sy - a*sx)/N
r = abs(sy*sx/N-sxy)/math.sqrt((sxx-sx*sx/N)*(syy-sy*sy/N))
return a,b,r
if __name__ == '__main__':
x=[ 1 ,2 ,3 ,4 ,5 ,6]
y=[ 2.5 ,3.51 ,4.45 ,5.52 ,6.47 ,7.51]
a,b,r=linefit(x,y)
print("X=",x)
print("Y=",y)
print("擬合結(jié)果: y = %10.5f x + %10.5f , r=%10.5f" % (a,b,r) )
plt.plot(x, y, "r:", linewidth=2)
plt.grid(True)
plt.show()
顯示圖像如下：
2. 不用擬合，直接顯示一個(gè)一元函數(shù)
#-*- coding:utf-8 -*- import numpy as np import matplotlib.pyplot as plt import math f = lambda x:5*x+4 tx = np.linspace(0,10,50) print tx plt.plot(tx, f(tx), "r-", linewidth=2) plt.grid(True) plt.show()
運(yùn)行結(jié)果：
PS：這里再為大家推薦兩款相似的在線工具供大家參考：
在線多項(xiàng)式曲線及曲線函數(shù)擬合工具：
http://tools.jb51.net/jisuanqi/create_fun
在線繪制多項(xiàng)式/函數(shù)曲線圖形工具：
http://tools.jb51.net/jisuanqi/fun_draw
更多關(guān)于Python相關(guān)內(nèi)容感興趣的讀者可查看本站專題：《Python數(shù)學(xué)運(yùn)算技巧總結(jié)》、《Python數(shù)據(jù)結(jié)構(gòu)與算法教程》、《Python函數(shù)使用技巧總結(jié)》、《Python字符串操作技巧匯總》及《Python入門與進(jìn)階經(jīng)典教程》
希望本文所述對(duì)大家Python程序設(shè)計(jì)有所幫助。
相關(guān)文章
在本篇文章里小編給大家整理的是一篇關(guān)于python寫try語句的方法相關(guān)內(nèi)容，有興趣的朋友們可以跟著學(xué)習(xí)參考下。2020-07-07
在本篇內(nèi)容里小編給大家整理了關(guān)于django教程學(xué)習(xí)流程，想自學(xué)的朋友們可以閱讀下。2020-07-07
Python 文件操作技巧(File operation) 實(shí)例代碼分析
python遍歷文件夾和文件 perl分割路徑和文件名2008-08-08
利用Django內(nèi)置的認(rèn)證視圖實(shí)現(xiàn)用戶密碼重置功能詳解
這篇文章主要給大家介紹了關(guān)于利用Django內(nèi)置的認(rèn)證視圖實(shí)現(xiàn)用戶密碼重置功能的相關(guān)資料，文中通過示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧。2017-11-11
python神經(jīng)網(wǎng)絡(luò)tfrecords文件的寫入讀取及內(nèi)容解析
這篇文章主要為大家介紹了python神經(jīng)網(wǎng)絡(luò)tfrecords文件的寫入讀取及內(nèi)容解析，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步，早日升職加薪2022-05-05
Pyqt實(shí)現(xiàn)簡(jiǎn)易計(jì)算器功能
這篇文章主要為大家詳細(xì)介紹了Pyqt實(shí)現(xiàn)簡(jiǎn)易計(jì)算器功能，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2022-05-05
這篇文章主要介紹了Python中的下劃線詳解,本文講解了單個(gè)下劃線直接做變量名、單下劃線前綴的名稱、雙下劃線前綴的名稱等內(nèi)容,需要的朋友可以參考下2015-06-06
Python實(shí)現(xiàn)輕松找出兩個(gè)列表不同之處
在日常編程中,需要比較兩個(gè)列表并找出它們之間差異是一種常見需求,在本文中,我們將深入探討Python中查找兩個(gè)列表差異值的方法,需要的小伙伴可以參考下2023-12-12
9行Python3代碼實(shí)現(xiàn)批量提取PDF文件的指定內(nèi)容
這篇文章主要為大家詳細(xì)介紹了如何通過9行Python3代碼實(shí)現(xiàn)批量提取PDF文件的指定內(nèi)容，文中的示例代碼講解詳細(xì)，感興趣的小伙伴可以嘗試一下2022-12-12
這篇文章主要為大家詳細(xì)介紹了Python yield的使用，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下，希望能夠給你帶來幫助2022-02-02