很多时候,我们有从非常分散的数据列表中生成的线图,这使得图形看起来像连接点的直线,或者非常密集,这导致数据点彼此非常接近,因此图看起来很混乱。
默认情况下,matplotlib.pyplot.plot()函数通过用直线连接数据中的两个相邻点来生成曲线,因此matplotlib.pyplot.plot()函数不会为小范围的数据点生成平滑曲线。
示例:
import numpy as np
import matplotlib.pyplot as plt# Dataset
x = np.array([ 1, 2, 3, 4, 5, 6, 7, 8 ])
y = np.array([ 20, 30, 5, 12, 39, 48, 50, 3 ])# Plotting the Graph
plt.plot(x, y)
plt.title("Curve plotted using the given points")
plt.xlabel("X")
plt.ylabel("Y")
plt.show()
可以看到,由于底层数据不遵循平滑的直线,因此该图一点也不平滑。为了绘制一条平滑的曲线,我们首先将一条样条曲线拟合到曲线上,并使用该曲线来找到x值的y值,x值被一个无限小的间隙隔开。我们可以通过用一个非常小的间隙画出这些点来得到一条光滑的曲线。
我们可以使用以下方法来创建此数据集的平滑曲线:
1.使用PyPlot绘制平滑曲线:
它通过首先使用scipy.interpolate.make_interp_spline()确定样条曲线的系数来绘制平滑的样条曲线。我们使用给定的数据点来估计样条曲线的系数,然后使用这些系数来确定非常接近的x值的y值,以使曲线看起来平滑。这里我们将使用np.linspace()方法,它返回在指定时间间隔内计算的均匀间隔的样本。可选参数num是在start和stop范围内生成的样本数。默认值为50,且必须为非负数。我们希望该参数具有足够高的值以生成平滑曲线。让我们在最小值和最大值之间沿X轴沿着取500个等距样本来绘制曲线。
语法:
numpy.linspace(start, stop, num = 50, endpoint = True, retstep =
False, dtype = None, axis = 0)X_Y_Spline = scipy.interpolate.make_interp_spline(x, y)
import numpy as np
import numpy as np
from scipy.interpolate import make_interp_spline
import matplotlib.pyplot as plt # Dataset
x = np.array([1, 2, 3, 4, 5, 6, 7, 8])
y = np.array([20, 30, 5, 12, 39, 48, 50, 3])X_Y_Spline = make_interp_spline(x, y)# Returns evenly spaced numbers
# over a specified interval.
X_ = np.linspace(x.min(), x.max(), 500)
Y_ = X_Y_Spline(X_)# Plotting the Graph
plt.plot(X_, Y_)
plt.title("Plot Smooth Curve Using the scipy.interpolate.make_interp_spline() Class")
plt.xlabel("X")
plt.ylabel("Y")
plt.show()
2. 三次插值样条曲线
它使用scipy.interpolate.interp1d类生成一条三次插值曲线,然后我们使用该曲线来确定平滑曲线的密集x值的y值。这里我们也将使用np.linspace()方法,该方法返回在指定时间间隔内计算的均匀间隔的样本。让我们在最小值和最大值之间沿X轴沿着取500个等距样本来绘制曲线。根据您希望直线弯曲的程度,可以修改第三个参数(num)的值。
语法: numpy.linspace(start, stop, num=50, endpoint=True, retstep=False,
dtype=None, axis=0)cubic_interpolation_model=scipy.interpolate.interp1d(x,y,kind=”cubic”)
import numpy as np
from scipy.interpolate import interp1d
import matplotlib.pyplot as plt # Dataset
x=np.array([1, 2, 3, 4, 5, 6, 7, 8])
y=np.array([20, 30, 5, 12, 39, 48, 50, 3])cubic_interpolation_model = interp1d(x, y, kind = "cubic")# Plotting the Graph
X_=np.linspace(x.min(), x.max(), 500)
Y_=cubic_interpolation_model(X_)plt.plot(X_, Y_)
plt.title("Plot Smooth Curve Using the scipy.interpolate.interp1d Class")
plt.xlabel("X")
plt.ylabel("Y")
plt.show()