result spectrogram from lowpass filter still pass frequency above the treshold

2 ビュー (過去 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 ExchangeTime-Frequency Analysis についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by