Python線性擬合實(shí)現(xiàn)函數(shù)與用法示例
本文實(shí)例講述了Python線性擬合實(shí)現(xiàn)函數(shù)與用法。分享給大家供大家參考，具體如下：
1. 參考別人寫(xiě)的：
#-*- 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)文章
這篇文章主要介紹了python禁用鍵鼠與提權(quán)代碼實(shí)例,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-08-08
Prometheus開(kāi)發(fā)中間件Exporter過(guò)程詳解
這篇文章主要介紹了Prometheus開(kāi)發(fā)中間件Exporter過(guò)程詳解,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-11-11
這篇文章主要介紹了python打包多類型文件的方法，本文給大家介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值，需要的朋友可以參考下2020-09-09
python基礎(chǔ)學(xué)習(xí)之生成器與文件系統(tǒng)知識(shí)總結(jié)
本文是參考《python數(shù)據(jù)分析》的附錄對(duì)生成器和文件系統(tǒng)結(jié)合案例的一個(gè)簡(jiǎn)單回顧,文中對(duì)python生成器與文件系統(tǒng)作了非常詳細(xì)的介紹,對(duì)正在學(xué)習(xí)python的小伙伴們有很好地幫助,需要的朋友可以參考下2021-05-05
keras自定義回調(diào)函數(shù)查看訓(xùn)練的loss和accuracy方式
這篇文章主要介紹了keras自定義回調(diào)函數(shù)查看訓(xùn)練的loss和accuracy方式，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-05-05
解決pycharm 遠(yuǎn)程調(diào)試 上傳 helpers 卡住的問(wèn)題
今天小編就為大家分享一篇解決pycharm 遠(yuǎn)程調(diào)試 上傳 helpers 卡住的問(wèn)題，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-06-06
pygame實(shí)現(xiàn)俄羅斯方塊游戲（AI篇1）
這篇文章主要為大家詳細(xì)介紹了pygame實(shí)現(xiàn)俄羅斯方塊游戲AI的第1篇，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2019-10-10
這篇文章主要介紹了Python 使用office365郵箱的示例，幫助大家利用python進(jìn)行高效辦公，感興趣的朋友可以了解下2020-10-10
這篇文章主要介紹了Python @property裝飾器的用法，在Python中，可以通過(guò)@property裝飾器將一個(gè)方法轉(zhuǎn)換為屬性，從而實(shí)現(xiàn)用于計(jì)算的屬性，下面文章圍繞主題展開(kāi)更多相關(guān)詳情，感興趣的小伙伴可以參考一下2023-01-01
Python數(shù)字圖像處理之霍夫線變換實(shí)現(xiàn)詳解
這篇文章主要介紹了Python數(shù)字圖像處理之霍夫線變換實(shí)現(xiàn)詳解，具有一定借鑒價(jià)值,需要的朋友可以參考下2018-01-01