python绘图-傅里叶级数

# encoding: gbk
# 傅里叶级数
import matplotlib.pyplot as plt
import numpy as np

x = np.arange(0, 6*np.pi, 0.01)

n = 1 # 频率倍数n, 取1,3,5,7,9等

y = 1/n * np.sin(n * x)

plt.plot(x, y * 4 / np.pi)

plt.title("n="+str(n))

plt.show()

import matplotlib.pyplot as plt
import numpy as np

x = np.arange(0, 6*np.pi, 0.01)

y = np.sin(x)

label = "1"

for n in range(2, 4):
    if n % 2 != 0:
        label = label+","+str(n)
        print(n)
        y = y + 1/n * np.sin(n * x)

plt.plot(x, y * 4 / np.pi)

plt.title("n="+label)

plt.ylim(-1.5, 1.5)

plt.show()

1倍频率

3倍频率

5倍频率

1+3倍数频率合成

1+3+5倍数频率合成

1，3，5，7倍数频率合成

1，3，5，7，9倍数频率合成

1，3，5.。。。9999倍频率合成