? ? ? 在實際的應用中,我們經常會遇到一些曲線擬合得需求,盡管,實際上點對之間可能是沒有絕對意義上的函數關系的,但是為了更好地量化表示和計算,我們往往期望于從原始的數據集中盡可能地挖掘出來可能的函數關系,今天我遇上了一個問題就是如何去擬合指數函數+對數函數這種組合型函數曲線,簡單看一下具體的實現內容:
def logExpFunc():
'''
【指數+對數】混合函數擬合
'''
x=np.linspace(1,2,15)
y=[21.5,23.1,25.9,30,32.6,38,41.9,47.2,55,61,69,80,90,105,117.6]
popt,pcov=curve_fit(lambda t,a,b,c,d,e: a+b*np.log(c*t)+d*np.exp(e*t),x,y)
print('popt: ',popt)
plt.figure()
a,b,c,d,e=popt
y_pre = a+b*np.log(c*x)+d*np.exp(e*x)
plt.plot(x, y, 'ko', label="Original Data")
plt.plot(x, y_pre, 'r-', label="Fitting Curve")
plt.legend()
plt.show()
? ? ?這里我的函數原型為:??
y = 4+2*np.exp(2*x)+3*np.log(2*x)
? ? ?y為指數函數和對數函數的組合形式,其中,一共有5個參數。
? ? ? 結果如下所示:
popt: [ 5.19163529 -1.73152698 1.57951688 2.44883023 1.92488973]
? ? ? 從結果輸出的5個參數上來看,a、b、c、d和e跟我們自己預設的函數原型中的參數是不一致的,但是不影響我們擬合得到所需要的曲線,這樣的參數不一定會是唯一的,可能跟我們函數曲線擬合的時候所使用到的點數據過少有關系吧,畢竟函數本身還是比較復雜的。
?
?
更多文章、技術交流、商務合作、聯系博主
微信掃碼或搜索:z360901061

微信掃一掃加我為好友
QQ號聯系: 360901061
您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。
【本文對您有幫助就好】元
