From 507f249a7c5462dfac6c4ac0ecdf881485df5a60 Mon Sep 17 00:00:00 2001
From: Gabri6 <gabriel.lucas.64100@gmail.com>
Date: Thu, 20 Apr 2023 11:14:54 +0200
Subject: [PATCH] add frequency filter

---
 main.m | 21 ++++++++++++++++++---
 1 file changed, 18 insertions(+), 3 deletions(-)

diff --git a/main.m b/main.m
index 2860a64..c7acb7a 100644
--- a/main.m
+++ b/main.m
@@ -53,13 +53,28 @@ size(fft_signal)
 
 figure;
 plot(t, fft_signal);
-xlabel('Time (s)');
+xlabel('Freq (Hz)');
 ylabel('Sound (dB)');
 title('FFT of signal');
 
 
-%b = fir1(600, [1/fMax, 1/fMin], 'bandpass');
-%filtered_signal = filter(b, signal);
+N=size(signal)(2);
+freq=(0:N-1)*samplingFreq/N;
+%frequency from frame rate
+
+
+minfreq=30; %Hz
+maxfreq=40; %Hz
+idx_min = find(freq >= minfreq, 1);
+idx_max = find(freq <= maxfreq, 1, 'last');
+filtered_freq = freq(idx_min:idx_max);
+figure;
+plot(t, filtered_freq);
+
+
+
+%[a, b] = butter(1, [1/fMax, 1/fMin], 'bandpass');
+%filtered_signal = filter(b,1, signal);
 %figure;
 %plot(t,filtered_signal)
 %xlabel('Time (s)');