martes, 13 de marzo de 2012

PRACTICA MATLAB NOTAS MUSICALES


DO
                                      

ESPECTRO DO


>> x=wavread('Do.wav');
>> fs=10000;%frecuencia de muestreo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y)) 

RE
 
ESPECTRO RE

>> x=wavread('Re.wav');
>> fs=10000;%frecuencia de muestreo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y))

MI
ESPECTRO MI

>> x=wavread('Mi.wav');
>> fs=10000;%frecuencia de muestreo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y))

FA
ESPECTRO FA
>> x=wavread('Fa.wav');
>> fs=10000;%frecuencia de muestreo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y))

SOL
 ESPECTRO SOL
>> x=wavread('Sol.wav');
>> fs=10000;%frecuencia de muestreo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y))

LA
ESPECTRO LA

>> x=wavread('La.wav');
>> fs=10000;%frecuencia de muestreo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y))

SI
 ESPECTRO SI
>> x=wavread('Si.wav');
>> fs=10000;%frecuencia de muestreo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y))
















jueves, 8 de marzo de 2012

PRACTICA DE SONIDO

>> plot(abs(y))
>> x=wavread('probando.wav');
>> fs=22050;%frecuencia de muestreo de lo que se grabo
>> sound(x,fs)
>> plot(x)
>> y=fft(x);
>> plot(abs(y))


miércoles, 22 de febrero de 2012

TRABAJO MATLAB SERIE FOURIER

ONDA COMPLETA




>> T=100;
t= 0:(3*T)-1;
p=(t/T);
s=20.*cos((2*pi).*p);
s1=15.*cos(3*(2*pi).*p);
s2=11.*cos(5*(2*pi).*p);
s3=9.*cos(7*(2*pi).*p);
s4=7.*cos(9*(2*pi).*p);
s5=1.*cos(11*(2*pi).*p);
salida=s+s1+s2+s3+s4+s5;
plot (salida)


ONDA  "S"


>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> salida=s;
>> plot (salida)

ONDA " S1 "
>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s1=15.*cos(3*(2*pi).*p);
>> salida=s1;
>> plot (salida,'r')

ONDA " S2 "



>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s2=11.*cos(5*(2*pi).*p);
>> salida=s2;
>> plot (salida, 'y')

ONDA " S3 "
>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s3=9.*cos(7*(2*pi).*p);
>> salida=s3;
>> plot (salida, 'm')

ONDA " S4 "
>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s4=7.*cos(9*(2*pi).*p);
>> salida=s4;
>> plot (salida, 'g')

ONDA" S5 "
>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s5=1.*cos(11*(2*pi).*p);
>> salida=s5;
>> plot (salida, 'k')

SUPERPOSICIÓN DE ONDAS








CON "S+S1"

>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> s1=15.*cos(3*(2*pi).*p);
>> salida=s+s1;
>> plot (salida)

CON "S+S1+S2"
>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> s1=15.*cos(3*(2*pi).*p);
>> s2=11.*cos(5*(2*pi).*p);
>> salida=s+s1+s2;
>> plot (salida)


CON "S+S1+S2+S3"

>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> salida=s+s1;
>> s2=11.*cos(5*(2*pi).*p);
>> s3=9.*cos(7*(2*pi).*p);
>> salida=s+s1+s2+s3+s4+s5;
>> salida=s+s1+s2+s3;
>> plot (salida)

CON "S+S1+S2+S3+S4"
>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> s1=15.*cos(3*(2*pi).*p);
>> s2=11.*cos(5*(2*pi).*p);
>> s3=9.*cos(7*(2*pi).*p);
>> s4=7.*cos(9*(2*pi).*p);
>> salida=s+s1+s2+s3+s4;
>> plot (salida)


VARIANDO " T "



>> T=200;
>> t= 0:(3*T)-2;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> s1=15.*cos(3*(2*pi).*p);
>> s2=11.*cos(5*(2*pi).*p);
>> s3=9.*cos(7*(2*pi).*p);
>> s4=7.*cos(9*(2*pi).*p);
>> s5=1.*cos(11*(2*pi).*p);
>> salida=s+s1+s2+s3+s4+s5;
>> plot (salida)

VARIANDO " t "
>> T=100;
>> t= 0:(3*T)-8;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> s1=15.*cos(3*(2*pi).*p);
>> s2=11.*cos(5*(2*pi).*p);
>> s3=9.*cos(7*(2*pi).*p);
>> s4=7.*cos(9*(2*pi).*p);
>> s5=1.*cos(11*(2*pi).*p);
>> salida=s+s1+s2+s3+s4+s5;
>> plot (salida)



VARIANDO " t "
>> T=100;
>> t= 0:(9*T)-8;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> s1=15.*cos(3*(2*pi).*p);
>> s2=11.*cos(5*(2*pi).*p);
>> s3=9.*cos(7*(2*pi).*p);
>> s4=7.*cos(9*(2*pi).*p);
>> s5=1.*cos(11*(2*pi).*p);
>> salida=s+s1+s2+s3+s4+s5;
>> plot (salida)


 ESPECTRO DE LA ONDA
>> T=100;
>> t= 0:(3*T)-1;
>> p=(t/T);
>> s=20.*cos((2*pi).*p);
>> s1=15.*cos(3*(2*pi).*p);
>> s2=11.*cos(5*(2*pi).*p);
>> s3=9.*cos(7*(2*pi).*p);
>> s4=7.*cos(9*(2*pi).*p);
>> s5=1.*cos(11*(2*pi).*p);
>> salida=s+s1+s2+s3+s4+s5;
>> espectro=fft(salida);
>> esp= abs(espectro);
>> plot(p,salida)
>> stem(esp);