CASE CLOSED. THANKS YOU

7 views (last 30 days)
Nabil Amirul
Nabil Amirul on 8 Jul 2021
Commented: Mathieu NOE on 12 Jul 2021
I have 1440 data of tidal elevation (24 hours recording every 1 minute), I attached my data (Book1.xlsx). My question is how can I do High Pass Filter with those data? My experience in this case is zero.

Accepted Answer

Mathieu NOE
Mathieu NOE on 8 Jul 2021
hello Nabil
the code below implements low and high pass filtering
hope it helps
T = readtable('Book1.xlsx');
time = T.Var1;
y = T.Var2;
dt = 60; % seconds (data sampled every min)
Fs = 1/dt;
samples = length(y);
t = (0:samples - 1)*dt/3600; % now in hours
%%% filter data %%%
N = 20;
y_LP = smoothdata(y, 'gaussian' , N); % low pass filtered data
y_HP = y - y_LP; % high pass filtered data
%%% plot data %%%
figure(1)
subplot(211),plot(t,y,t,y_LP);legend('Raw','low pass filtered');
title(['Data samples at Fs = ' num2str(Fs) ' Hz / filtered with smoothdata' ]);
xlabel('Time (hours) ')
ylabel('Tidal depth () ')
subplot(212),plot(t,y_HP);legend('high pass filtered');
title(['Data samples at Fs = ' num2str(Fs) ' Hz / filtered with smoothdata' ]);
xlabel('Time (hours) ')
ylabel('Tidal depth () ')
  2 Comments
Nabil Amirul
Nabil Amirul on 9 Jul 2021
Thank you very much sir, I am very appreciate it
Mathieu NOE
Mathieu NOE on 12 Jul 2021
My pleasure

Sign in to comment.

More Answers (0)

Tags

Community Treasure Hunt

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

Start Hunting!