1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71
| clear all;
close all;
%------------------ définition des constantes ------------------%
%Caractéristiques du signal x(t)
A=1; % amplitude (V)
fs=16000; % frequence d'echantillonage(Hz)
coef_bruit=0.1; % bruit blanc gaussien
% paramètre affichage
duree=3;
nb_echan=fs*duree;
%Génération du signal x
xt= wavread('1.wav'); % formule du signal
%------------------ Acquisition/génération des signaux ------------------%
t=linspace(0,nb_echan,length(xt)); %définition de la base de temps
%Génération du bruit blanc
bt=coef_bruit*(randn(length(xt),1));
%Génération du signal issu de la mesure
yt=xt+bt;
%------------------ Traitement des signaux ------------------%
%------------------ Visualisation des signaux ------------------%
% trace de x(t)
figure(1);
subplot(2,2,1)
plot(t,xt,'-');
title('Tracés de x(t)');
xlabel('Temps t (s)');
ylabel('Signaux mesurés (V)');
grid on;
axis('square');
% trace de b(t)
%figure(2);
subplot(2,2,2)
plot(t,bt,'-');
title('Tracés de b(t)');
xlabel('Temps t (s)');
ylabel('Signaux mesurés (V)');
grid on;
axis ('square');
% trace de y(t)
%figure(3);
subplot(2,2,3)
plot(t,yt,'-');
title('Tracés de y(t)');
xlabel('Temps t (s)');
ylabel('Signaux mesurés (V)');
legend('Signal y(t)');
grid on;
axis ('square'); |