ola

SAR.m

-[x1,fe1]=audioread("C:\Users\camil\Documents\IMT_A\semestre_6\electrical engineering\tp-audio-ee-etudiant-c24leray\src\wav\single_tone_piano1.wav");
-[x2,fe2]=audioread("C:\Users\camil\Documents\IMT_A\semestre_6\electrical engineering\tp-audio-ee-etudiant-c24leray\src\wav\single_tone_piano2.wav");
+[x1,fe1]=audioread("C:\Users\samue\OneDrive\Documents\Cours\cours IMT atlantique\A1\elec inge\traitement du signal\SAR\tp-audio-ee-etudiant-c24leray\src\wav\single_tone_piano1.wav");
+[x2,fe2]=audioread("C:\Users\samue\OneDrive\Documents\Cours\cours IMT atlantique\A1\elec inge\traitement du signal\SAR\tp-audio-ee-etudiant-c24leray\src\wav\single_tone_piano2.wav");
 soundsc(x1,fe1);
 soundsc(x2,fe2);
 
@@ -20,4 +20,31 @@ plot(f2, log(abs(X2)),'b')
 title("spectre")
 xlabel("fréquence en Hz")
 ylabel("log|X(f)|")
-legend(piano1, piano2)
+
+
+freq_r=[220;442;663;885;1108;1331;1556;1782;2009];
+
+% Durée et fréquence d'échantillonnage
+Fe = 44100;           % fréquence d'échantillonnage (standard audio)
+duree = 1;            % durée du signal en secondes
+t = 0:1/Fe:duree;     % vecteur temps
+
+% Initialiser le signal composite
+s = zeros(size(t));
+
+% Ajouter les sinusoïdes
+for k = 1:length(freq_r)
+    s = s + sin(2*pi*freq_r(k)*t);
+end
+
+% Normaliser pour éviter la saturation (valeurs entre -1 et 1)
+s = s / max(abs(s));
+
+% Jouer le son
+audioplayer(s, Fe);
+
+
+
+
+    
+