UNIVERSIDAD NACIONAL MAYOR

UNIVERSIDAD NACIONAL MAYOR

DE SAN MARCOS

FACULTAD DE INGENIERA ELECTRNICA

ALUMNO: TORRES CAHUANA, DANNY OMAR

CODIGO :

01114794

CURSO :

LABORATORIO DE PDS

TEMA : SOLUCION DEL EXAMEN FINAL

PROFESOR : RAFAEL BUSTAMANTE

CICLO : VIII

Examen Final de laboratorio de Procesamiento Digital de Seales (PDS)

1) Generar una secuencia

N=256; numero de muestras

T= 1/300; T= periodo de muestreo

x(k) = sen(20kT) +sen(60kT) +sen(100kT) +sen(140kT) +sen(180kT)

2) Generar el espectro de frecuencia de la pregunta 1

3) Generar mediante simunlink la seal de la pregunta 1

4) Filtrar la seal sen(100kT) de la pregunta 1

5) Generar la siguiente secuencia

Solucion1)

%pregunta 1 generar una secuencia

N=256;

T=1/300;

k=0:N-1;

x=sin(20*pi*k*T)+sin(60*pi*k*T)+sin(100*pi*k*T)+sin(140*pi*k*T)+sin(180*pi*k*T);

subplot(2,2,1);

plot(k,x);grid

subplot(2,2,2);

stem(k,x);grid;

Solucin 2)

% generar el espectro de frecuencia de la pregunta 1

X=fft(x);

magX=abs(X);

subplot(2,2,3),

stem(k(1:N/2),magX(1:N/2));grid

hertz=k*(1/(N*T));

subplot(2,2,4);

stem(hertz(1:N/2),magX(1:N/2));grid;

A continuacin se muestra los grficos correspondientes a lo problemas 1 y 2

Solucin 3)

Generar mediante simulink la seal de la pregunta No1

Donde:

Signal Generator Frecuencia =10Hz

Signal Generator1 Frecuencia =30Hz

Signal Generator2 Frecuencia =50Hz

Signal Generator3 Frecuencia =70Hz

Signal Generator4 Frecuencia= 90Hz

A la salida obtenemos el siguiente grafico:

Solucin 4)

Filtrar la seal sen (100kT) de la pregunta 1

clc

clear

N=256;

T=1/300;

k=0:N-1;

hertz=k*(1/(N*T));

x=sin(20*pi*k*T) + sin(60*pi*k*T) + sin(100*pi*k*T) + sin(140*pi*k*T) + sin(180*pi*k*T);

===================================================

figure(1)

[B,A]=butter(1,[0.332 0.335]);

%bode(B,A)

k=0:N/2;

H=freqz(B,A,k,300);

plot(k,abs(H)),grid

y=filter(B,A,x);

=====================================================

figure(2)

plot(hertz,y),grid

figure(3)

X=fft(y);

k=0:N-1;

magY=abs(X);

hertz=k*(1/(N*T));

plot(hertz(1:N/2),magY(1:N/2)),grid

Figura 1

Figura 2

Figura 3

Solucin 5)Para resolver este problema creamos las funciones bsicas escaln y rampafunction r=escalon(x,n)

r=zeros (size(x));

set1=find(x>n) ;

r(set1)=ones(size(set1));

====================================

function r=rampa(x,n)

r=zeros(size(x));

set1=find(x>=n)

r(set1)=x(set1)-n

=====================================

Ejecutamos el siguiente cdigo para obtener la grafica pedida

x=-1:0.1:10;

y=rampa(x,0)-rampa(x,1)+rampa(x,3)-rampa(x,4)-2*escalon(x,6);

stem(x,y),grid

plot(x,y);

PAGE 5