登录论坛

查看完整版本 : Matlab过滤问题


poster
2019-12-10, 20:48
我有一个数据数组,它在时域中。每个数据代表幅度。这些数据以10,000Hz的频率采样。我想在两个频率f1和f2之间做一个带通滤波器。我想我需要做一个低通滤波器和一个高通滤波器。信号先进行FFT,然后经过两个滤波器,再进行反FFT。我的问题是有一种简单的方法来做低通和高通滤波器吗?我不知道如何基于两个截止频率来推导传递函数。

有人知道吗?

谢谢



回答:

基于Tristan的答案,这是一些Octave代码,可能与Matlab兼容,也可能与Matlab不兼容。黄油函数为您导出传递函数系数。 https://i.stack.imgur.com/B8RZV.png 。 https://i.stack.imgur.com/fq70n.png

hz = 8000; x = [1:1:hz*10]; t = x./hz; pi = 3.1415; % Create signal with 10 hz, 200 hz and 500 hz components raw_signal = sin(10*2*pi*t)+sin(200*2*pi*t)+sin(500*2*pi*t); % View Raw Signal Over .1 Second Window plot(t, raw_signal) title('Raw Signal with 10Hz, 200Hz, 500Hz Components') xlabel('Time (Sec)') ylabel('Amplitude') set(gca,'XLim', [5, 5.1]); % Create Band Pass Butterworth Filter [S_numer, S_denom] = butter(5, [100/hz 350/hz]); band_passed_signal = filter(S_numer, S_denom, raw_signal); % View Band Pass Filtered Signal Over .1 Second Window plot(t, band_passed_signal) title('Band Pass Filtered Signal') xlabel('Time (Sec)') ylabel('Amplitude') set(gca,'XLim', [5, 5.1]);

更多&回答... (https://stackoverflow.com/questions/4384190)