clc
clear all;

t=0:149;

y = wavread('D:\Eli\Mis documentos\Wav\kamel.wav',1000);
x=length(y);

% subplot (321)
% plot(y)
% title('Se馻l de audio a(n)')

z=hist(y,150);
f1 = 1000*(normpdf(0:149,75,20))

subplot (311)
plot(t,z,t,f1);
title('Histograma de a(n)')

r = .6 + sqrt(0.1) * randn(size(y));


% subplot (323)
% plot (r,'r')
% title('Se馻l de ruido r(n)')

z2=hist(r,150);


f2 = 1000*(normpdf(0:149,75,20))

subplot (312)
plot(t,z2,'r',t,f2,'b');
title('Histograma de r(n)')


ar=y+r;

% subplot(325)
% plot(ar,'g')
% title('Se馻l de audio con ruido, a(n)+ r(n)')

z3=hist(ar,150);
f3 = 1000*(normpdf(0:149,75,20))

subplot(313)
plot(t,z3,'g',t,f3,'r')
title('Histograma de a(n)+ r(n)')