2022-09-01から1ヶ月間の記事一覧

VTF (dooley)

import math import numpy as np import matplotlib.pyplot as plt #vtf = 5.05*math.exp(-0.138 * math.pi*300/180 * u ) * ( 1 - math.exp(-0.1 * math.pi * 300/180 * u) ) u = np.linspace(0, 10, 101) #print(u) #print(u.size) vtf = np.empty(u.size)…

FFT 参考

FFT

helve-blog.com import numpy as np import matplotlib.pyplot as plt N = 1024 # サンプル数 dt = 0.001 # サンプリング周期 [s] f1, f2 = 50, 120 # 周波数 [Hz] t = np.arange(0, N*dt, dt) # 時間 [s] x = 1.5*np.sin(2*np.pi*f1*t) + np.sin(2*np.pi*f2…