質問: 理想的な低域通過フィルタにかけた信号のグラフ(横軸:時間、縦軸:電圧)について trace_1は「trace_1(:,1):時間、trace_1(:,2):電圧」になっています。 trace_1を理想的なlpf( 周波数が10mhz以上なら利得が0 )にかけた時の「trace_1のグラフ(横軸:時間、縦軸:電圧)」を出力したいです。 butter(https://jp.mathworks.com/help/signal/ref/butter.html)やlowpass(https://jp.mathworks.com/help/signal/ref/lowpass.html)を使うと理想的なlpfが出来ないので次のような方法で理想的なlpfを作成し、lpfに通した後の信号を逆フーリエ変換することを検討しています。 手順1 trace1をフーリエ変換する。 手順2 周波数が10mhz以上なら利得0にする。(理想的なlpfの作成) 手順3 lpfにかけた信号を逆フーリエ変換する。(横軸:時間、縦軸:電圧) どのようにプログラムをかけば、理想的なlpfにかけたあとの信号を逆フーリエ変換することが出来ますか?ご回答お待ちしております。
回答: 理想的なlpfは、物理的に実現できないので、どのような意味があるかはわかりませんが、実現できたと想定した場合の特性は計算はできると思います。 x = trace_1(:,2)として 手順1 xのフーリエ変換x = fft(x)を求めます 手順2 周波数信号xの内、10mz ~ サンプリング周波数/2 に相当する要素に0をいれ、xfilterとします。 手順3 この信号の逆フーリエ変換をy = ifft(xfilter)として求めます。 フーリエ変換を参考にされてください。 注意する点は手順2で、フーリエ変換xには正の周波数と負の周波数の部分があるので、負の(10mz ~ サンプリング周波数/2)の部分にも忘れずに0を入れないと、結果の信号yが実数になりません。