Procesamiento de señales de audio por medio de matlab

|%% Procesamiento de una señal de audio usando MATLAB |
|%% Crear señal de audio |
|% Frecuencia fundamental|
|f0=1e3; % 1KHz |
|% Amplitud |
|a=4; % V=4|
|% Frecuencia de muestreo |
|fs=44.1e3; % Frecuencia de una señal de audio |
|% Tiempo de duración en segundos|
|T=1.5; |
|L = round(T*fs); % Número de muestras |
|% Frecuencia normalizada|
|fn=f0/fs; |
|y = a*sin(2*pi*fn*(0:L-1))+0.5*a*sin(2*pi*2*fn*(0:L-1)); |
|% Graficar la señal original|
|subplot(411) |
|plot((0:L-1)/fs,y) |
|title('SEÑAL ORIGINAL')% Título|
|xlabel('Tiempo (s)')         % Etiqueta del eje X |
|ylabel('Amplitud (V)')      % Etiqueta del eje Y |
|xlim([0 10/1000])            % Límite de la señal|
|%% Grabar y reproducir la señal de audio |
|%wavwrite(y,fs,'audio') |
|% wavplay(y,fs)|
|%% FFT de la señal |
|subplot(412) |
|% Llamado a la función que calculala FFT |
|fft_signal(y,fs);title('ESPECTRO DE LA SEÑAL ORIGINAL') |
|xlim([0 2500]) |
|%% Filtrado de la señal...