【FFT】MATLAB、Pythonで株価予測 その30【IFFT⑨】

バックナンバーはこちら。
https://www.simulationroom999.com/blog/stock-predict-matlabpython-backnumber/

はじめに
登場人物
FFT出力の周波数分布をローテーション
出力結果
まとめ

はじめに

前回はMATLABによるベクトルローテーションとしてcircshift関数の挙動実験を行った。
基本は列ベクトルに対して行うもので行ベクトルだと期待する動作はしなかったが、
引数の設定次第では行ベクトルに対してもローテーション可能。
基本は列ベクトルにしておいた方が良さそうである。

今回は実際に、FFT出力の周波数分布をローテーションさせてみる。

登場人物

博識フクロウのフクさん

イラストACにて公開の「kino_k」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=iKciwKA9&area=1

エンジニア歴8年の太郎くん

イラストACにて公開の「しのみ」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=uCKphAW2&area=1

FFT出力の周波数分布をローテーション

太郎くん

circshift関数でローテーションできるのはわかったから、
早速MATLABでFFT出力をローテーションさせてみよう。

フクさん

ざっと書いたコードがこれだ。

N=1024;
L=pi;
x=linspace(-L,L,N)';
k=-N/2:N/2-1;
ft=sin(x)+sin(3*x)+sin(7*x);

% 変換前波形
subplot(2,1,1)
plot(x,ft)
title('f(t)');
grid();

% FFT後にローテーション
Fw=fft(ft);
Fw_tmp=circshift(Fw,N/2);
subplot(2,1,2);
mask = (-10<=k & k<=10);
plot(k(mask),abs(Fw_tmp(mask)));
title('F(\omega)');
grid();