% Author: Adrien COMBE
% Date: 24/03/2023
% Description: Signal Processing Lab 2

pkg load signal

filename = "modulator22.wav";
[y,fs] = audioread(filename);

%TEMPORAL ANALYSIS
t = [1:1:size(y)];
plot(t,y);
title ("Temporal analysis");

%PERFORMING DFT
[power, duration] = frequencySpectrum(y,fs,false);
%PERFORMING FFT
[power2, duration2] = frequencySpectrum(y,fs,true);
duration
duration2

%SPECTOGRAM
spectrogram(y,fs,5,30);
spectrogram(y,fs,5,5);