MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   资料存档 (https://www.labfans.com/bbs/forumdisplay.php?f=72)
-   -   Matlab过滤问题 (https://www.labfans.com/bbs/showthread.php?t=24200)

poster 2019-12-10 20:48

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

有人知道吗?

谢谢



[B]回答:[/B]

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

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]);

[url=https://stackoverflow.com/questions/4384190]更多&回答...[/url]


所有时间均为北京时间。现在的时间是 09:24

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.