diff --git a/src/wav/question_2_1.m b/src/wav/question_2_1.m new file mode 100644 index 0000000000000000000000000000000000000000..def5523d5a5a1b5f23b747246d4d9b1be10c2a99 --- /dev/null +++ b/src/wav/question_2_1.m @@ -0,0 +1,50 @@ +% Paramètres +A = 1; % Amplitude +T = 1; % Période +N = 1000; % Nombre d'échantillons +Fs = N / T; % Fréquence d'échantillonnage +t = linspace(0, T, N); % Temps sur une période +f = (-N/2:N/2-1) * Fs/N; % Axe des fréquences + +% ----- Signal carré centré ----- +square_signal = A * square(2*pi*t/T); % par défaut, valeur entre -1 et 1 + +% ----- Signal dent de scie centré ----- +sawtooth_signal = A * sawtooth(2*pi*t/T); % par défaut entre -1 et 1 + +% ----- FFT ----- +fft_square = fftshift(fft(square_signal)); +fft_saw = fftshift(fft(sawtooth_signal)); + +% Magnitude normalisée +fft_square_mag = abs(fft_square)/N; +fft_saw_mag = abs(fft_saw)/N; + +% ----- Tracés ----- +figure; + +subplot(2,2,1); +plot(t, square_signal); +title('Signal carré (centré)'); +xlabel('Temps (s)'); +ylabel('Amplitude'); + +subplot(2,2,2); +plot(f, fft_square_mag); +xlim([-50 50]); +title('Spectre du signal carré'); +xlabel('Fréquence (Hz)'); +ylabel('|X(f)|'); + +subplot(2,2,3); +plot(t, sawtooth_signal); +title('Signal dent de scie (centré)'); +xlabel('Temps (s)'); +ylabel('Amplitude'); + +subplot(2,2,4); +plot(f, fft_saw_mag); +xlim([-50 50]); +title('Spectre du signal dent de scie'); +xlabel('Fréquence (Hz)'); +ylabel('|X(f)|');