Remerciez-le!

Remerciez @Admin pour avoir partagé cet document gratuitement, de la manière la plus simple, en partageant sur les réseaux sociaux.

Partie A % 1 - Chargement du fichier Audio [y, fs]=audioread('C:\Users\PC\Deskt

Partie A % 1 - Chargement du fichier Audio [y, fs]=audioread('C:\Users\PC\Desktop\sonbruit.wav'); % 3 - Lecture du fichier Audio sound(y, fs); pause % 4 - Affichage du signal bruité dans une figure % Transformée de fourier du signal non filtré N=256; f = linspace(0,fs/2,N); A=abs(fft(y(1:N))); % module de la transformée subplot(211) plot(f,A); xlabel('fréquences (Hz)'); ylabel('Amplitude'); title('Spectre avec bruit'); % 5 - Suppression du bruit n = 10; fm = 1000 / (fs/2); [b,a] = butter(n,fm,'low'); yFiltre = filter(b, a, y); Fourier=abs(fft(yFiltre(1:N))); subplot(212) plot(f,Fourier); xlabel('fréquences (Hz)'); ylabel('Amplitude'); title('Spectre sans bruit'); % 6 - Lecture du nouveau son ynew = audioplayer(yFiltre, fs); ynew.play Partie B % Chargement du fichier Audio [y, fs]=audioread('C:\Users\PC\Desktop\mamusic.wav'); % Lecture du fichier Audio sound(y, fs); pause % Spectre d'amplitude du signal non filtré N=256; f = linspace(0,fs/2,N); Amp=abs(fft(y(1:N))); subplot(2,1,1) plot(f,Amp); xlabel('fréquences (Hz)'); ylabel('Amplitude'); title('Spectre d''amplitude initial'); %Traitement du signal n = 10 fm = 3000 / (fs/2); [b,a] = butter(n,fm,'high'); yfiltre = filter(b, a, y); % Spectre d'amplitude du signal filtré Ampf=abs(fft(yfiltre(1:N))); subplot(2,1,2) plot(f,Ampf); xlabel('fréquences (Hz)'); ylabel('Amplitude'); title('Spectre d''amplitude filtré'); % Lecture du nouveau son y = audioplayer(yfiltre, fs); y.play uploads/Litterature/ projet-partie-a-et-b 1 .pdf