From 975041a2a8f441475267caec588fd7d0e6867a76 Mon Sep 17 00:00:00 2001 From: YE Victor <victor.ye@imt-atlantique.net> Date: Tue, 27 May 2025 16:17:35 +0000 Subject: [PATCH] Upload New File --- src/SAR2Q1bis.m | 43 +++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 43 insertions(+) create mode 100644 src/SAR2Q1bis.m diff --git a/src/SAR2Q1bis.m b/src/SAR2Q1bis.m new file mode 100644 index 0000000..74e70ca --- /dev/null +++ b/src/SAR2Q1bis.m @@ -0,0 +1,43 @@ +% Paramètres +T = 1; % période du signal +Fe = 1000; % fréquence d'échantillonnage +Te = 1/Fe; +t = 0:Te:5*T-Te; % 5 périodes +N = length(t); + +% Signal carré +carre = square(2*pi*t/T); % amplitude ±1 + +% Signal dent de scie +scie = sawtooth(2*pi*t/T); % amplitude ±1 + +% FFT et centrage +Carre_FFT = fftshift(fft(carre)); +Scie_FFT = fftshift(fft(scie)); +f = (-N/2:N/2-1)*(Fe/N); + +% Amplitudes linéaires normalisées +Carre_amp = abs(Carre_FFT)/N; +Scie_amp = abs(Scie_FFT)/N; + +% Limite à 0–50 Hz +indices = find(f >= 0 & f <= 50); +f_zoom = f(indices); +Carre_zoom = Carre_amp(indices); +Scie_zoom = Scie_amp(indices); + +% Affichage du signal carré +figure; +plot(f_zoom, Carre_zoom); +title('Spectre du signal carré (amplitude linéaire)'); +xlabel('Fréquence (Hz)'); +ylabel('Amplitude'); +grid on; + +% Affichage du signal dent de scie +figure; +plot(f_zoom, Scie_zoom); +title('Spectre du signal dent de scie (amplitude linéaire)'); +xlabel('Fréquence (Hz)'); +ylabel('Amplitude'); +grid on; \ No newline at end of file -- GitLab