import?pandas?as?pd
import?numpy?as?np
import?matplotlib.pyplot?as?plt
import?scipy.stats?as?st

result_parameter_peak?=?pd.read_csv("result_parameter_peak.csv",?encoding="utf_8_sig")

#?設(shè)置畫布
fig?=?plt.figure(figsize=(20,?8))??#?figsize是常用的參數(shù).(寬，高)
axl?=?fig.add_subplot(1,?1,?1)

for?i?in?range(len(result_parameter_peak)):
????x?=?np.arange(0,?400,?1)
????#?繪制gamma曲線
????y661?=?st.gamma.pdf(x,?result_parameter_peak.iloc[i,?1],?scale=result_parameter_peak.iloc[i,?2])
????axl.plot(x,?y661,?'r-.',?label="α=?9.9028，β=10.4205")
????#?設(shè)置坐標(biāo)軸標(biāo)題
????axl.set_xlabel('Time')
????axl.set_ylabel('Probility')
????axl.set_title('分布')
????#?可視化
????plt.show()

import?pandas?as?pd
import?numpy?as?np
import?matplotlib.pyplot?as?plt
import?scipy.stats?as?st

result_parameter_peak?=?pd.read_csv("result_parameter_peak.csv",?encoding="utf_8_sig")

plt.figure()
for?i,?alpha,?beta?in?result_parameter_peak.itertuples():
????x?=?np.arange(0,?300,?1)
????#?繪制gamma曲線
????y661?=?st.gamma.pdf(x,?alpha,?scale=beta)
????plt.plot(x,?y661,?'-.')
????#?設(shè)置坐標(biāo)軸標(biāo)題
????plt.xlabel('Time')
????plt.ylabel('Probility')
????plt.title('分布')

#?可視化
plt.show()

import?pandas?as?pd
import?numpy?as?np
import?matplotlib.pyplot?as?plt
import?scipy.stats?as?st

result_parameter_peak?=?pd.read_csv("result_parameter_peak.csv",?encoding="utf_8_sig")

plt.figure()
for?i,?alpha,?beta?in?result_parameter_peak.itertuples():
????x?=?np.arange(0,?300,?1)
????#?繪制gamma曲線
????y661?=?st.gamma.pdf(x,?alpha,?scale=beta)
????#?plt.plot(x,?y661,?'-.')
????plt.plot(x,?y661,?'-.',?label="α:"?+?str(alpha)?+?"β:"?+?str(beta))
????#?設(shè)置坐標(biāo)軸標(biāo)題
????plt.xlabel('Time')
????plt.ylabel('Probility')
????plt.title('fenbu')
????
#?可視化
plt.legend()
plt.show()