ryamada https://blog.hatena.ne.jp/ryamada/ ryamadaのコンピュータ・数学メモ https://ryamada.hatenadiary.jp/ python フーリエ変換 高速フーリエ変換 from scipy.fftpack import fft, ifft # 周期性データを作る n = 1000 # 点の数 t = np.linspace(-1,1,n) * 5*np.pi # 時刻点 k = 50 # 重ね合わせる周期関数の数 # a cos(bt+c)の係数a,b,c a = np.random.randn(k) b = np.random.randn(len(a)) c = np.random.randn(len(a)) # ゼロベクトル x = np.zeros_like(t) for i in arange(len(a)) : x = x + a… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada.hatenadiary.jp%2Fentry%2F20150101%2F1419899282" title="フーリエ変換〜数値解析諸法に慣れる〜scipy - ryamadaのコンピュータ・数学メモ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/r/ryamada/20141230/20141230055959.png Hatena Blog https://hatena.blog 2015-01-01 09:28:02 rich https://ryamada.hatenadiary.jp/entry/20150101/1419899282 1.0 100%