Python應用實現(xiàn)雙指數(shù)函數(shù)及擬合代碼實例
更新時間：2020年06月19日 10:28:52 作者：kurrrr
這篇文章主要介紹了Python應用實現(xiàn)雙指數(shù)函數(shù)及擬合代碼實例,文中通過示例代碼介紹的非常詳細，對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下
雙指數(shù)函數(shù)
待擬合曲線為 y(x) = bepx + ceqx
import matplotlib.pyplot as plt x = ([0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0]) y = ([0.33, 0.26, 0.18, 0.16, 0.12, 0.09, 0.08, 0.07, 0.06, 0.06, 0.06, 0.07, 0.09, 0.1, 0.15, 0.19, 0.25, 0.36, 0.47, 0.68]) plt.scatter(x, y) plt.show()
擬合
import numpy as np import matplotlib.pyplot as plt from scipy.optimize import curve_fit def double_exp(x, b, c, p, q): x = np.array(x) return b*np.exp(p*x) + c*np.exp(q*x) x = ([0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0]) y = ([0.33, 0.26, 0.18, 0.16, 0.12, 0.09, 0.08, 0.07, 0.06, 0.06, 0.06, 0.07, 0.09, 0.1, 0.15, 0.19, 0.25, 0.36, 0.47, 0.68]) popt, pcov = curve_fit(double_exp, x, y, [1, 1, 1, 1]) print(popt) b = popt[0] c = popt[1] p = popt[2] q = popt[3] y_fit = double_exp(x, b, c, p, q) plt.scatter(x, y) plt.plot(x, y_fit, color='red', linewidth=1.0) plt.show()
numpy 庫，實現(xiàn)列表轉矩陣，得以進行數(shù)學運算。matplotlib.pyplot 庫，繪制圖像。scipy.optimize 庫，curve_fit() 函數(shù)，使用非線性最小二乘法擬合曲線。curve_fit()popt，擬合結果，在這里指b, c, p, q 的值。povc，該擬合結果對應的協(xié)方差。
擬合結果參數(shù)原函數(shù)擬合結果誤差b0.00110.00110c0.42000.42420.42%p6.39986.49881.55%q-5.1551-5.21641.19%
誤差可以滿意。
經(jīng)過測試，如果將初始參數(shù)設置為原函數(shù)參數(shù)（保留 4 位小數(shù)），擬合得到的結果并未發(fā)生變化。
經(jīng)過測試，擬合使用的三種方法，"trf"，"lm" 和 "dogbox" 對該函數(shù)擬合結果影響微乎其微。
以上就是本文的全部內容，希望對大家的學習有所幫助，也希望大家多多支持腳本之家。
相關文章
本篇文章主要介紹基于python的七種經(jīng)典排序算法（推薦），具有一定的參考價值，這里整理了詳細的代碼,有需要的小伙伴可以參考下。2016-12-12
這篇文章主要介紹了2013-11-11
今天小編就為大家分享一篇解決Python 中英文混輸格式對齊的問題，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2018-07-07
這篇文章主要介紹了python數(shù)字圖像處理之邊緣輪廓檢測示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進步，早日升職加薪2022-06-06
在Python中使用requests庫發(fā)送HTTP請求時,有時會遇到亂碼的問題,亂碼通常是由于編碼不一致或解碼錯誤導致的,這篇文章給大家介紹了Python requests亂碼的五種解決辦法,并通過代碼示例講解的非常詳細,需要的朋友可以參考下2024-04-04
迭代器工具在產(chǎn)生數(shù)據(jù)的時候將會顯得非常便捷、高效，掌握了這些基本的方法之后，通過簡單的組合就可以獲得更多迭代器工具。2016-05-05
大家好，本篇文章主要講的是Python爬蟲之urllib庫詳解，感興趣的同學趕快來看一看吧，對你有幫助的話記得收藏一下2022-02-02
這篇文章主要介紹了Django框架請求生命周期實現(xiàn)原理,文中通過示例代碼介紹的非常詳細，對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下2020-11-11](http://img.jbzj.com/images/xgimg/bcimg8.png)
jupyter notebook運行命令顯示[*](解決辦法)
這篇文章主要介紹了jupyter notebook運行命令顯示[*],文中通過示例代碼介紹的非常詳細，對大家的學習或者工作具有一定的參考學習價值，需要的朋友們下面隨著小編來一起學習學習吧2020-05-05
Pandas 數(shù)據(jù)處理,數(shù)據(jù)清洗詳解
今天小編就為大家分享一篇Pandas 數(shù)據(jù)處理,數(shù)據(jù)清洗詳解，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2018-07-07