Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
高级会员
注册日期: 2019-11-21
帖子: 3,006
声望力: 66 ![]() |
![]()
我有一个数据数组,它在时域中。每个数据代表幅度。这些数据以10,000Hz的频率采样。我想在两个频率f1和f2之间做一个带通滤波器。我想我需要做一个低通滤波器和一个高通滤波器。信号先进行FFT,然后经过两个滤波器,再进行反FFT。我的问题是有一种简单的方法来做低通和高通滤波器吗?我不知道如何基于两个截止频率来推导传递函数。
有人知道吗? 谢谢 回答: 基于Tristan的答案,这是一些Octave代码,可能与Matlab兼容,也可能与Matlab不兼容。黄油函数为您导出传递函数系数。 ![]() ![]() 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]); 更多&回答... |
![]() |
![]() |