def IndustWater(time): pop = (Gaussian(Fit.PopulationFit(time),pPara)/Gaussian(Fit.PopulationFit(years[-1]),pPara)) pcg = (Gaussian(Fit.PCGDPFit(time),gPara)/Gaussian(Fit.PCGDPFit(years[-1]),gPara)) ele = (Gaussian(Fit.ElectricityFit(time),ePara)/Gaussian(Fit.ElectricityFit(years[-1]),ePara)) spr = (Gaussian(Fit.SteelProductFit(time),sPara)/Gaussian(Fit.SteelProductFit(years[-1]),sPara)) return Data.WaterUseIndustry[years[-1]]*(pop/4+pcg/4+ele/4+spr/4)
plt.show() plt.close() plt.figure() plt.xlabel("Population /10k people") plt.ylabel("Industry Water Usage/ 100m m^3") plt.scatter([Data.Population[t] for t in years],[Data.WaterUseIndustry[t] for t in years],marker="^",s=50,label='Raw Data') plt.plot(Fit.PopulationFit(x),Gaussian(Fit.PopulationFit(x),pPara),'g-',label="Fit") plt.legend() plt.show() plt.close() plt.figure() plt.xlabel("PCGDP /CNY") plt.ylabel("Industry Water Usage/ 100m m^3") plt.scatter([Data.PCGDP[t] for t in years],[Data.WaterUseIndustry[t] for t in years],marker="^",s=50,label='Raw Data') plt.plot(Fit.PCGDPFit(x),Gaussian(Fit.PCGDPFit(x),gPara),'g-',label="Fit") plt.legend() plt.show() plt.close() plt.figure() plt.xlabel("Electricity /100M kWh") plt.ylabel("Industry Water Usage/ 100m m^3") plt.scatter([Data.Electricity[t] for t in years],[Data.WaterUseIndustry[t] for t in years],marker="^",s=50,label='Raw Data') plt.plot(Fit.ElectricityFit(x),Gaussian(Fit.ElectricityFit(x),ePara),'g-',label="Fit") plt.legend() plt.show() plt.close()