Python實(shí)現(xiàn)的擬合二元一次函數(shù)功能示例【基于scipy模塊】
本文實(shí)例講述了Python實(shí)現(xiàn)的擬合二元一次函數(shù)功能。分享給大家供大家參考，具體如下：
背景：
使用scipy擬合一元二次函數(shù)。
參考：
HYRY Studio-《用Python做科學(xué)計(jì)算》
代碼：
# -*- coding:utf-8 -*-
#! python3
import numpy as np
from scipy.optimize import leastsq
import pylab as pl
def func(x,p):
"""
數(shù)組擬合函數(shù)
"""
A,k,theta = p
return A*(x-k)**2+theta
def residuals(p,y,x):
"""
殘差
"""
return y-func(x,p)
x = np.linspace(0,2,100)
A,k,theta = 10.,1,2. #真實(shí)數(shù)據(jù)參數(shù)
y0 = func(x,[A,k,theta]) #真實(shí)數(shù)據(jù)
y1 = y0 + 2 * np.random.randn(len(x)) #加入噪聲序列
p0 = [7.,0.2,1.]
plsq = leastsq(residuals,p0,args = (y1,x))
print("真實(shí)參數(shù):",[A,k,theta])
print("擬合參數(shù):",plsq[0]) #試驗(yàn)數(shù)據(jù)擬合后的參數(shù)
pl.plot(x,y0,label = "real")
pl.plot(x,y1,label = "real+noise")
pl.plot(x,func(x,plsq[0]),label = "fitting")
pl.legend()
pl.show()
結(jié)果：
（貌似這里的求解方法用了智能算法，因?yàn)槊看蔚慕Y(jié)果都有細(xì)小差異。具體資料沒見到，以后有精力再找）
真實(shí)參數(shù): [10.0, 1, 2.0]
擬合參數(shù): [ 10.83391995 0.98950039 1.63356065]
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)典教程》
希望本文所述對大家Python程序設(shè)計(jì)有所幫助。
相關(guān)文章
python量化之搭建Transformer模型用于股票價(jià)格預(yù)測
這篇文章主要介紹了python量化之搭建Transformer模型用于股票價(jià)格預(yù)測，文章圍繞主題展開基于python搭建Transformer，需要的小伙伴可以參考一下2022-05-05
python進(jìn)程管理工具supervisor的安裝與使用教程
supervisor是用python寫的一個(gè)進(jìn)程管理工具，用來啟動(dòng)，重啟，關(guān)閉進(jìn)程。下面這篇文章主要給大家介紹了關(guān)于python實(shí)現(xiàn)的進(jìn)程管理工具supervisor的安裝與使用的相關(guān)資料，需要的朋友可以參考借鑒，下面來一起看看吧。2017-09-09
用python實(shí)現(xiàn)監(jiān)控視頻人數(shù)統(tǒng)計(jì)
今天教各位小伙伴學(xué)習(xí)怎么用python實(shí)現(xiàn)監(jiān)控視頻人數(shù)統(tǒng)計(jì),文中有非常詳細(xì)的代碼示例,對正在學(xué)習(xí)python的小伙伴有很大的幫助,需要的朋友可以參考下2021-05-05
這篇文章主要為大家介紹了Langchain集成管理prompt功能示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步，早日升職加薪2023-03-03
Python使用方法重載實(shí)現(xiàn)訪問者模式
這篇文章主要為大家詳細(xì)介紹了Python如何使用方法重載實(shí)現(xiàn)訪問者模式,文中的示例代碼講解詳細(xì),感興趣的小伙伴可以跟隨小編一起學(xué)習(xí)一下2024-12-12
Python實(shí)現(xiàn)計(jì)算字符串中出現(xiàn)次數(shù)最多的字符示例
這篇文章主要介紹了Python實(shí)現(xiàn)計(jì)算字符串中出現(xiàn)次數(shù)最多的字符,涉及Python針對字符串的遍歷、統(tǒng)計(jì)等相關(guān)操作技巧,需要的朋友可以參考下2019-01-01
OpenCV+Python幾何變換的實(shí)現(xiàn)示例
這篇文章主要介紹了OpenCV+Python幾何變換的實(shí)現(xiàn)示例，圖像的幾何變換是指將一幅圖像映射到另一幅圖像內(nèi)。有縮放、翻轉(zhuǎn)、仿射變換、透視、重映射等操作。感興趣的可以了解一下2021-03-03
這篇文章主要介紹了python中管道用法,實(shí)例分析了管道的創(chuàng)建、輸入、輸出及關(guān)閉等技巧,需要的朋友可以參考下2015-06-06
解決pycharm每次新建項(xiàng)目都要重新安裝一些第三方庫的問題
今天小編就為大家分享一篇解決pycharm每次新建項(xiàng)目都要重新安裝一些第三方庫的問題，具有很好的參考價(jià)值，希望對大家有所幫助。一起跟隨小編過來看看吧2019-01-01
利用pyinstaller或virtualenv將python程序打包詳解
這篇文章主要給大家介紹了利用pyinstaller將python程序打包的相關(guān)資料，文中介紹的非常詳細(xì)，相信對大家具有一定的參考價(jià)值，需要的朋友們下面來一起看看吧。2017-03-03