弦波疊加

sin_square_wave.py

import numpy as np
import matplotlib.pyplot as plt

freq = 10         # 頻率
stop = .25        # 取樣範圍為 0 ~ stop
sample_rate = 800 # 取樣率，單位 x 取樣幾次
n = 1000          # 疊加的波數

x = np.linspace(0, stop, int(stop * sample_rate), endpoint = False)
i = np.arange(1, 2 * n, 2).reshape((n, 1))
y = (np.sin(i * freq * 2 * np.pi * x) / i).sum(axis = 0)

plt.plot(x, y)
plt.show()