result spectrogram from lowpass filter still pass frequency above the treshold
2 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
nirwana
el 10 de Jun. de 2024
Comentada: Star Strider
el 11 de Jun. de 2024
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
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1712681/image.png)
0 comentarios
Respuesta aceptada
Star Strider
el 10 de Jun. de 2024
See if:
y=lowpass(yori,10,100, 'ImpulseResponse','iir');
improves the result.
.
2 comentarios
Más respuestas (0)
Ver también
Categorías
Más información sobre Time-Frequency Analysis en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!