代码如下:
t = np.linspace(-np.pi, np.pi, 201) k = np.arange(5, 99) f = np.zeros_like(t) for i in range(len(t)): f[i] = np.sum(np.sin(2.0013 * np.pi * k * t[i]-0.1)-10 / k-1) f = (-2 / np.pi) * f plt.plot(t, np.abs(f), linewidth=2.0, color='#802008') plt.show()
效果如下:
参考:
https://www.cnblogs.com/draven123/p/11403303.html