Procesamiento digital de señales

%PROGRAMA DE CONVOLUCIÓN[pic 1]

clear all;

xn = input ('Escribe los valores de la función x(n): ' ); %Entrada

xh = input ('Escribe los valores de la función h(n): ' ); %Respuesta impulso

k = input  ('Escribe el valor de k: ' ); % Definir rango desde k=0 hasta un valor k deseado

escala3 = 0;

escala2 = 0;

escala1 = 0;

figure,stem (xn , 'b');

    title ('Señal Entrada x(n)');

    ylabel ('x(n)');

    xlabel ('n');

    axis ([1 10 -20 20] );

figure,stem (xh , 'r' );

    title ('Señal  Respuesta Impulso x(h)');

    ylabel ('x(h)');

    xlabel ('h');

    axis ([1 10 -20 20])

figure,stem (y,'g' );

    title ('Señal Convolcionada y(n)');

    ylabel ('y(n)' );

    xlabel ('n');

    axis ([1 15 -50 200]);

%Estas instrucciones diseñan tres filtros digitales.

%Usados para un problema de separación de canales.

clear all;

fs = 10000 ;    %Frecuencia de Nyquist

T = 1/fs ;      %Tiempo de muestreo

fn = fs/2;     %Frecuencia de Muestro "Total"

f1n = 100/fn;   %Corte pasabajas normalizado

f2n = 500/fn ;  %Corte pasabanda izquierdo normalizado;

f3n = 1000/fn ; %Corte pasabanda derecho normalizado

f4n = 2000/fn ; %Corte pasaaltas normalizado

[B1,A1] = butter (8,f1n) ;

[B2,A2] = butter (7,[f2n,f3n]);

[B3,A3] = butter (10,f4n,'high');

% Mas grande N = mayor es la curva de respuesta de frecuencia.

[H1,wT] = freqz(B1,A1,200);

[H2,wT] = freqz(B2,A2,200);

[H3,wT] = freqz(B3,A3,200);

hertz = wT/(2*pi*T);

subplot(2,1,1),...

    plot(hertz,abs(H1),'-',hertz,abs(H2),'-',...

    hertz, abs(H3),'-'),...

    title ('Filtros para separa canales'),...

    xlabel ('Hz'), ylabel ('Magnitud'), grid

%PROGRAMA DE VENTANADO BARTLETT, HAMMING Y VON HANN

clear all;

n = -20:20;                     %Genera un indice DT

fc= 0.2; h=2*fc*sinc(2*fc*n);   %Genera la respuesta al impulso

a=[zeros(1,20) 1 zeros(1,20)];  %Arreglo del denominador de la DTFT

f= (-200:199)/400; w=2*pi*f;    %Arreglo de la frecuencia

hi= (abs (f)<=0.2 );            %Espectro ideal de filtro, FC=0.2 AZUL

wb= 1-abs (n)/20;               %Ventana Bartlett ROJA

wv = 0.5 + 0.5*cos(pi*n/20);    %Ventana de Von Hann (Hanning) NEGRO PUNTEADO

...