# -*- coding: utf-8 -*-""" Created on Sat Jul1818:13:292026@author: Thomas"""importnumpy as npimportmatplotlib.pyplot as plt# 屏幕显示和默认保存分辨率plt.rcParams.update({"figure.dpi":180,"savefig.dpi":600,"font.size":14,"axes.titlesize":17,"axes.labelsize":15,"legend.fontsize":12,"xtick.labelsize":12,"ytick.labelsize":12,"lines.linewidth":2.5,})x_m=3.0#alphas = [0.5, 1.0, 1.5, 2.0, 3.0, 5.0, 10.0]alphas=[0.5,4.0]x=np.linspace(x_m,10,3000)foralphainalphas: pdf=alpha * x_m**alpha / x**(alpha +1)plt.plot(x, pdf,linewidth=2,label=rf"$\alpha={alpha}$")plt.xlabel(r"$x$")plt.ylabel(r"$f(x)$")plt.title(rf"Pareto PDF,$x_m={x_m}$")plt.xlim(x_m,10)plt.ylim(bottom=0)plt.grid(alpha=0.3)plt.legend()plt.tight_layout()plt.savefig("pareto_pdf.png",dpi=300)plt.show()