將numpy作為np導入
將matplotlib.pyplot作為plt導入
來自scipy . interpolate import interp 1d
將scipy.stats導入為st
sim = ST . gamma（1，loc=0，scale = 0.8）#模擬
OBS = ST . gamma（2，loc=0，scale = 0.7）#觀察到的
x = NP . linspace（0，4，1000）
simpdf = sim . pdf（x）
OBS pdf = OBS . pdf（x）
PLT . plot（x，simpdf，label =‘Simulated‘）
PLT . plot（x，obspdf，‘r-‘，label =‘Observed‘）
plt.title（“觀測和模擬降水的PDF”）
PLT . legend（loc =‘best‘）
plt.show()
工廠圖（1）
sim CDF = sim . CDF（x）
OBS CDF = OBS . CDF（x）
PLT . plot（x，simcdf，label =‘Simulated‘）
PLT . plot（x，obscdf，‘r-‘，label =‘Observed‘）
plt.title（“觀測和模擬降水的CDF”）
PLT . legend（loc =‘best‘）
plt.show()
#逆CDF
inv CDF = interp 1d（OBS CDF，x）
transfer _ func = inv CDF（sim CDF）
工廠圖（2）
PLT . plot（transfer _ func，x，‘g-‘）
plt.show()