您正在关注此问题

您将在关注的内容源中看到更新。
您可能会收到电子邮件，具体取决于您的通信预设项。

Designing a bandpass filter

5 次查看（过去 30 天）

显示更早的评论

Temisanren Okotie 2020-10-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter

⋮

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter

评论： Star Strider 2020-10-26

采纳的回答： Star Strider

How do i design a bandpass filter that is used to eliminate the low and high frequency noises of an eeg signal with a cut off frequency of 0.1- 25 Hz

采纳的回答

Star Strider 2020-10-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#answer_524283

⋮

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#answer_524283

See if the bandpass fucntion works for you.

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1086988

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1086988

Thank you.

Star Strider 2020-10-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087018

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087018

As always, my pleasure!

Filters are straightforward to design using more basic functions in MATLAB, however for a single signal and a single filter run, the bandpass and related functions are easiest. If you intend to filter more signals, request the second output of bandpass (the digital filter object) and use it with other signals and the filtfilt function to do this most efficiently, otherwise bandpass re-designs the filter each time, and that is both inefficient and slow.

Temisanren Okotie 2020-10-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087058

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087058

Thank you very much!

load Subject00_before_fp1.mat;

fs = 500;

N = length(Subject00_before_fp1);

N2 = 12000;

t = (1:N)/fs;

subplot(2,1,1);

plot(t,Subject00_before_fp1);

fc = [0.1 25]/fs/2;

y = bandpass(Subject00_before_fp1,[0.1 25],fs);

subplot(2,1,2);

plot(t,y,'k');

figure;

P = nextpow2(12000);

n = 2^P;

X = fft(y(1:N2),n);

X_mag = abs(X);

X_power = (abs(X).^2)/N;

f = (0:N-1)*fs/N;

plot(f(1:N2),20*log10(X_power(1:N2)),'k');

I am trying to plot the PSD using fft function. The value of N is 91000 which is quite large and Matlab wasn't able to process it. I used the first 12000 samples. The result i get after plotting the PSD is still quite noisy. I don't know if my code is right. I'm new to Matlab

Star Strider 2020-10-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087188

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087188

在 MATLAB Online 中打开

This worked when I tried it (increasing the lower passband frequency to 0.5 Hz):

fs = 500; 
Subject00_before_fp1 = randn(91000,1);
N = length(Subject00_before_fp1);
t = (1:N)/fs;
[y,df] = bandpass(Subject00_before_fp1,[0.5 25],fs);
figure
subplot(2,1,1)
plot(t,Subject00_before_fp1)
grid
subplot(2,1,2)
plot(t,y)
grid

This designs a very long (23093) FIR filter, according to the ‘df’ object I requested so that I could see what bandpass was designing. The original filter designed an

order elliptic filter (the sort I would design myself) however the extremely low lower passband of 0.1 Hz is a nearly impossible task for any filter, even with a sampling frequency of 500 Hz.

Temisanren Okotie 2020-10-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087333

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087333

IMG_8517.pdf
Subject00_before_fp1.mat

I took your suggestion and made the changes but the PSD plot still looks thesame, maybe that part of the code is wrong? I attached a copy of the signal and the picture of the graphs I got if that would help. Thank you very much.

Star Strider 2020-10-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087413

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087413

It appears to be the PSD of an appropriately-filtered signal.

What were you expecting?

Temisanren Okotie 2020-10-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087428

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087428

C5799098319.pdf

This is the paper I'm basing my work on, the PSD graphs look quite different.

Star Strider 2020-10-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087503

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/625888-designing-a-bandpass-filter#comment_1087503

在 MATLAB Online 中打开

They actually do not look much different.

Use:

xlim[0 25])

with your plot and they should look reasonably similar to those in the paper.

Please do not expect to reproduce the results in any specific paper exactly. There are simply too many considerations involved, specifically the filter the paper uses and its characteristics, such as transition-region width, filter type (Butterworth, elliptic, FIR, none of which I can find stated in the paper) or other details. Frequently, ‘close enough is good enough’.

Temisanren Okotie 2020-10-26