python 繪制擬合曲線并加指定點標識


import os
import numpy as np
from scipy import log
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt
import math
from sklearn.metrics import r2_score
# 字體
plt.rcParams['font.sans-serif']=['SimHei']
# 擬合函數(shù)
def func(x, a, b):
# y = a * log(x) + b
y = x/(a*x+b)
return y
# 擬合的坐標點
x0 = [2, 4, 8, 10, 24, 28, 32, 48]
y0 = [6.66,8.35,10.81,11.55,13.63,13.68,13.69,13.67]
# 擬合,可選擇不同的method
result = curve_fit(func, x0, y0,method='trf')
a, b = result[0]
# 繪制擬合曲線用
x1 = np.arange(2, 48, 0.1)
#y1 = a * log(x1) + b
y1 = x1/(a*x1+b)
x0 = np.array(x0)
y0 = np.array(y0)
# 計算r2
y2 = x0/(a*x0+b)
#y2 = a * log(x0) + b
r2 = r2_score(y0, y2)
#plt.figure(figsize=(7.5, 5))
# 坐標字體大小
plt.tick_params(labelsize=11)
# 原數(shù)據(jù)散點
plt.scatter(x0,y0,s=30,marker='o')
# 橫縱坐標起止
plt.xlim((0, 50))
plt.ylim((0, round(max(y0))+2))
# 擬合曲線
plt.plot(x1, y1, "blue")
plt.title("標題",fontsize=13)
plt.xlabel('X(h)',fontsize=12)
plt.ylabel('Y(%)',fontsize=12)
# 指定點,y=9時求x
p = round(9*b/(1-9*a),2)
#p = b/(math.log(9/a))
p = round(p, 2)
# 顯示坐標點
plt.scatter(p,9,s=20,marker='x')
# 顯示坐標點橫線、豎線
plt.vlines(p, 0, 9, colors = "c", linestyles = "dashed")
plt.hlines(9, 0, p, colors = "c", linestyles = "dashed")
# 顯示坐標點坐標值
plt.text(p, 9, (float('%.2f'% p),9),ha='left', va='top', fontsize=11)
# 顯示公式
m = round(max(y0)/10,1)
print(m)
plt.text(48, m, 'y= x/('+str(round(a,2))+'*x+'+str(round(b,2))+')', ha='right',fontsize=12)
plt.text(48, m, r'$R^2=$'+str(round(r2,3)), ha='right', va='top',fontsize=12)
# True 顯示網(wǎng)格
# linestyle 設(shè)置線顯示的類型(一共四種)
# color 設(shè)置網(wǎng)格的顏色
# linewidth 設(shè)置網(wǎng)格的寬度
plt.grid(True, linestyle = "--", color = "g", linewidth = "0.5")
plt.show()另外有需要云服務器可以了解下創(chuàng)新互聯(lián)scvps.cn,海內(nèi)外云服務器15元起步,三天無理由+7*72小時售后在線,公司持有idc許可證,提供“云服務器、裸金屬服務器、高防服務器、香港服務器、美國服務器、虛擬主機、免備案服務器”等云主機租用服務以及企業(yè)上云的綜合解決方案,具有“安全穩(wěn)定、簡單易用、服務可用性高、性價比高”等特點與優(yōu)勢,專為企業(yè)上云打造定制,能夠滿足用戶豐富、多元化的應用場景需求。
文章題目:python繪制擬合曲線并加指定點標識的實現(xiàn)-創(chuàng)新互聯(lián)
網(wǎng)頁路徑:http://www.yijiale78.com/article18/deoegp.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供全網(wǎng)營銷推廣、服務器托管、軟件開發(fā)、域名注冊、企業(yè)建站、用戶體驗
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時需注明來源: 創(chuàng)新互聯(lián)
猜你還喜歡下面的內(nèi)容