サンプルプログラムのソースコードです。 # -*- coding: utf-8 -*- import numpy as np import matplotlib.pyplot as plt # データのパラメータ N = 256 # サンプル数 dt = 0.01 # サンプリング間隔 f1, f2 = 10, 20 # 周波数 t = np.arange(0, N*dt, dt) # 時間軸 freq = np.linspace(0, 1.0/dt, N) # 周波数軸 # 信号を生成(周波数10の正弦波+周波数20の正弦波+ランダムノイズ) f = np.sin(2*np.pi*f1*t) + np.sin(2*np.pi*f2*t) + 0.3 * np.random.randn(N) # 高速フーリエ変換 F = np.fft.fft(f) # 振幅スペクトルを計算 Amp = n
![【NumPy】高速フーリエ変換 (FFT)で振幅スペクトルを計算 | アルゴリズム雑記](https://cdn-ak-scissors.b.st-hatena.com/image/square/c81fd1dde433189d0ac488989199b6de0c9c52cf/height=288;version=1;width=512/https%3A%2F%2Falgorithm.joho.info%2Fwp-content%2Fuploads%2F2019%2F09%2Ffft2.jpg)