result spectrogram from lowpass filter still pass frequency above the treshold

4 ビュー (過去 30 日間)
I am trying use lowpass filter to my time series data and then i perform spectrogram to check wheater the signal still containing freq above the treshold or not.
clear all, clc, close all
yori=load("HHN__20150101T000000Z__20150102T000000Z.DAT");
y=lowpass(yori,10,100);
Nspec=256;
wspec=hamming(Nspec);
Noverlap=Nspec/2;
fs=100;
subplot 211
[S,F,T,P]=spectrogram(yori,wspec,Noverlap,Nspec,fs,'xaxis');
imagesc(T/(60*60),F,10*log10(P));
axis xy, title ("spectogram comp:HHN filter:NO")
ylabel('Frequency(Hz)')
subplot 212
[S,F,T,P]=spectrogram(y,wspec,Noverlap,Nspec,fs,'xaxis');
imagesc(T/(60*60),F,10*log10(P));
axis xy, title ("spectogram comp:HHN filter:YES (lowpass 10Hz")
ylabel('Frequency(Hz)')
Here my code, and the spectrogram shows that signal still contain freq above 10hz (that I don't want to pass). I don't understand why this happened? can someone give me suggestion to lowpass 10Hz signal? Thanks in advance

採用された回答

Star Strider
Star Strider 2024 年 6 月 10 日
See if:
y=lowpass(yori,10,100, 'ImpulseResponse','iir');
improves the result.
.
  2 件のコメント
nirwana
nirwana 2024 年 6 月 11 日
It's help sir!! Thanks a lot
Star Strider
Star Strider 2024 年 6 月 11 日
As always, my pleasure!

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeMultirate Signal Processing についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by