![]() |
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.