final version

Main.m

@@ -45,14 +45,14 @@ frequencySpectrum(windowed_signal, Fs, 0);
 %spectrogram(windowed_signal, Fs, 1/Fs, 1000*length(signal)/Fs);
 
 % filter using filter and butter
-%f = (0:length(windowed_signal)-1)*(Fs/length(windowed_signal));     % frequency range
-%y = fft(windowed_signal, length(windowed_signal));% compute DFT of input signal
 
-power = abs(windowed_signal).^2/length(windowed_signal); 
-[val, ind] = max(windowed_signal); % find the mx value of DFT and its index
+[val, ind] = max(windowed_signal);
 
 figure;
-[b, a] = butter(6, 10/Fs)
-s = filter(b, a, 10*log10(windowed_signal/windowed_signal(ind)))
-plot(1:length(windowed_signal), s(50:300));
+[b, a] = butter(6, 10/Fs);
+s = filter(b, a, 10*log10(windowed_signal/windowed_signal(ind)));
+plot(50:300, s(50:300));
+
+
+audiowrite("sound.wav", signal, Fs);