PWMによる正弦波生成方法(三角波比較方式) pythonによる実装例

import numpy as np
from scipy import signal
import matplotlib.pyplot as plt

f1 = 1        # 正弦波の周波数
f2 = 50      # 三角波の周波数
K = 0.015   # ローパスフィルタの時定数
v = 0.5      # 疑似正弦波の初期電圧値
v_t = []      # 電圧値格納スペース

x = np.arange(0,1,0.001)        # 時間
sin = (np.sin(2 * np.pi * f1 * x) + 1)/2        # 狙いの正弦波
tri = (signal.sawtooth(2 * np.pi * f2 * x, 0.5) + 1)/2    # 三角波
on = sin > tri        # ON/OFF信号

for i in on:                # 疑似正弦波生成
    v = K * i + (1-K)* v    # ローパスフィルタ
    v_t = np.append(v_t,v)

plt.plot(x, sin)
plt.plot(x, on)
plt.plot(x, v_t)
plt.show()