Cuadros Enviar Informe 2

UNIVERSIDAD CATOLICA DE SANTA MARIA Página:1/2
FACULTAD DE CIENCIAS E INGENIERÍAS FISICAS Y FORMALES

ESCUELA PROFESIONAL
Jefe de Prácticas:
INGENIERIA MECÁNICA, MECÁNICA-ELÉCTRICA Y MECATRÓNICA Ing. Juan Carlos Cuadros
Laboratorio de Procesamiento Digital de Señales
Código: 4E08143
Práctica 2 – Tema 1: SISTEMAS DE TIEMPO DISCRETO Semestre: VIII
Grupo: 1 2
FECHA:
Apellidos y Nombres:RAMIREZ ZAMORA JOSE FELIPE BLAIR P2-Tema Nº 02 21/OCT/2020
I. OBJETIVO
I.1. Obtener la respuesta de un sistema discreto por medio de la evaluación directa de las ecuaciones de diferencias y empleando
herramientas de software.
II. MARCO TEORICO
II.1. SISTEMAS DE TIEMPO DISCRETO

 Definición: Un sistema discreto procesa una secuencia de entrada x[n] para obtener una secuencia de salida y[n]
modificando sus propiedades (salvo el sistema y[n]=x[n]).
Un sistema discreto lineal e invariante en el tiempo puede ser descrito matemáticamente por una ecuación de
diferencias de coeficientes constantes de la forma:
Esta ecuación describe la relación entre la señal de salida y la señal de entrada. Se puede ver que están
involucrados los valores de instantes actuales, así como de instantes anteriores tanto para la señal de entrada como
para la de salida. De esta expresión se puede despejar el valor actual de la señal de salida y[n]
En Matlab las ecuaciones de diferencias se representan en dos vectores: un vector contiene los coeficientes del
lado directo, bm de los términos x, y otro vector contiene los coeficientes de la parte recurrente, ak, de los términos
y. Al coeficiente a0 es costumbre darle el valor 1. (En Matlab la función filter dividirá a la salida por a0 por
lo que este coeficiente no debe ser cero)
NOTA: La función y=filter(b,a,x) implementa un filtro digital, caracterizado por los vectores de
coeficientes a y b, que filtrara los datos almacenados en x. Si x es el impulso unitario, entonces y será la
respuesta al impulso h[n]. Tenga en cuenta que la función filter devuelve tantas muestras en y como muestras
hay en x.
III. MATERIAL Y EQUIPO
III.1. Una PC con SO Windows y MATLAB con el toolbox de PDS instalado.
IV. PROCEDIMIENTO
IV.1. EDD. Cree los vectores b y a que contengan respectivamente los coeficientes de x[n] y de y[n], de la siguiente ecuación
de diferencias:
𝑦[𝑛] + 0.9𝑦[𝑛 − 2] = 0.3𝑥[𝑛] + 0.6𝑥[𝑛 − 1] + 0.3𝑥[𝑛 − 2]
IV.2. Calcule analíticamente y[n] para x[n]=δ[n]

Laboratorio de Procesamiento Digital de Señales Página: 2/2
Práctica 2 – Tema 1: SISTEMAS DE TIEMPO DISCRETO P N° 2 – T N° 01 JP: Ing. JCC
IV.3. Evalúe los 30 primeros valores muestra y haga una gráfica de su respuesta analítica.
IV.4. Basado en la EDD General, cree un script en Matlab que solucione la EDD dada en el ítem IV.1. Grafique la respuesta
obtenida con el script.
%definimos el rango de n
n=0:31
for i=1:length(n) %establesemos la condicion de conteo para el bucle
n(i)=n(i)-2
end
n;
x=(n==0).*(1)+0;%declamramos que x sera 1siempre que n sea 0
y=0;%declaramos y
for z=3:length(n)%establesemos la condicion de conteo para el bucle
y(z)=(-0.9*y(z-2)+0.3*x(z)+0.6*x(z-1)+0.3*x(z-2));
end
stem(n,y)%aplicamos stem para graficar
axis([0,30,-0.8,0.8])%acomodamos la visibilidad de la grafica

IV.5. Usando la función filter solucione la EDD dada en el ítem IV.1. Grafique la respuesta obtenida.
%ecuaciones diferenciales
%y[n]+0.9y[n-2]=0.3x[n]+0.6x[n-1]+0.3x[n-2]
%definimos los vectores
a=[1,0,0.9];%definimos el vector a
b=[0.3,0.6,0.3];%definimos el vector b
n=0:29;%tiempo discreto
xn=(n==0).*(1)+0;%aplicamos el vector unitario
yn=filter(b,a,xn);%respuesta al vector unitario
stem(n,yn);%graficamos
IV.6. Haga las comparaciones entre las tres respuestas obtenidas mediante distintas formas.
o Las gráficas salen iguales la de recursión se mueve un poco, pero al momento que se le aplica el axis
se le da un intervalo y con eso podemos observar que las gráficas son iguales.
IV.7. Repita los pasos IV.2 al IV.6 para los sistemas dados por las siguientes EDD:
𝜋 1
 𝑦[𝑛] + 1.8𝑐𝑜𝑠 ( ) 𝑦 + 0.81𝑦[𝑛 − 2] = 𝑥[𝑛] + 𝑥[𝑛 − 1]
16 2
 recurcion
 %definimos el rango de n
 n=0:31
 for i=1:length(n) %establecemos la condición de conteo para el bucle
 n(i)=n(i)-2
 end
 n;
 y=[zeros(1,31)];%declaramos y
 x=(n==0).*(1)+0;%declaramos que x será 1siempre que n sea 0
 for z=3:length(n)%establecemos la condición de conteo para el bucle
 y(z)=(-1.8*cos(pi/16)*y(z-1)-0.81*y(z-2)+x(z)+1/2*x(z-1));
 end
 stem(n,y)%aplicamos stem para graficar
 axis([0,30,-2,2])%acomodamos la visibilidad de la grafica
 filter
 %ecuaciones diferenciales
 %y[n]+1.8cos(pi/16)y[n-1]+0.81y[n-2]=x[n]+1/2x[n-1]
 %definimos los vectores
 a=[1,1.8*cos(pi/16),0.81];%definimos el vector a
 b=[1,1/2,0];%definimos el vector b
 n=0:29;%tiempo discreto
 xn=(n==0).*(1)+0;%aplicamos el vector unitario
 yn=filter(b,a,xn);%respuesta al vector unitario
 stem(n,yn);%graficamos
 𝑦[𝑛] − 0.25𝑦[𝑛 − 2] = (0.5)𝑛𝑢[𝑛] con 𝑦[−1] = 0, 𝑦[−2] = 0 (Entrada diferente del impulso)
 recurcion
 %definimos el rango de n
 n=0:31
 for i=1:length(n) %establesemos la condicion de conteo para el bucle
 n(i)=n(i)-2
 end
 n;
 y=[zeros(1,31)];%declaramos y
 x=0.5.^n;%declamramos que x sera siempre que n sea 0
 for z=3:length(n)%establesemos la condicion de conteo para el bucle
 y(z)=0.5^z*x(z)+0.25*y(z-2);
 end
 stem(n,y)%aplicamos stem para graficar
 axis([0,30,0,0.15])%acomodamos la visibilidad de la grafica
 filter
 %ecuaciones diferenciales
 %y[n]-0.25*y[n-2]=0.5^n*x[x]
 %definimos los vectores
 a=[1,0,-0.25];%definimos el vector a
 b=[1,0,0];%definimos el vector b
 n=0:29;%tiempo discreto
 xn=0.5.^n;%aplicamos el vector unitario
 yn=filter(b,a,xn);%respuesta al vector unitario
 stem(n,yn);%graficamos
 𝑦[𝑛] + 4𝑦[𝑛 − 1] − 8𝑦[𝑛 − 2] = cos(𝑛𝜋) 𝑢[𝑛] con 𝑦[−1] = 0, 𝑦[−2] = 1 (Entrada diferente del impulso)
 Recurcion
 clear all, close all, clc
 n=0:31;
 for i=1:length(n)
 n(i)=n(i)-2;
 end
 n;
 x=cos(pi.*n);
 y=[1 0 zeros(1,29)];
 for j=3:length(n)
 y(j)=-4*y(j-1)+8*y(j-2)+x(j);
 end
 y;
 stem(n,y,'filled')
 filter
 clear all, close all, clc
 % Ecuación de diferencias
 % y[n]+4y[n-1]-8y[n-2]=u[n]
 % a0 a1 a2 b0
 % a=[a0 a1 a2]
 % b=[b0]
 n=0:29;
 a=[1 4 -8]; % Definimos el vector a
 b=[1]; % Definimos el vector b
 u=cos(pi.*n); % Definimos la función x[n] escalón unitario
 y=filter(b,a,u); % Obtenemos la respuesta y[n] a la entrada x[n]
 stem(n,y,'filled') % Graficamos la respuesta y[n] para la entrada x[n]
IV.8. Considere el sistema discreto:
 Calcule las 10 primeras muestras de su respuesta impulsional

 Encuentre la relación entrada salida
 Aplique la entrada 𝑥[𝑛] = {1,1,1, … … . } y calcule las 10 primeras muestras de la salida
 Parte analítica
 𝑌[𝑛] + 0.9𝑌[𝑛 − 2] = 0.3𝑋[𝑛] + 0.6𝑋[𝑛 − 1] + 0.3𝑋[𝑛 − 2]; 𝑋[𝑛] = 𝛿[𝑛] 𝐶. 𝐼 = 0
 𝐺𝑒𝑛𝑒𝑟𝑎𝑙: 〖 𝑌〗_(0+) 0.9𝑌_0 [𝑛 − 2] = 𝑋[𝑛]; 𝑌[−1] = 𝑌[−2] = 0
 𝑌_0𝑁 [𝑛] 𝑍 2 + 0.9 = 0
 𝑍1 = 0 + 0.9487𝑖 𝑍2 = 0 − 0.9487𝑖
 𝑌_0𝑁 [𝑛] = 𝑘1 〖(0 + 0.9487𝑖)〗𝑛 + 𝑘2 〖(0 − 0.9487𝑖)〗𝑛
 𝑌_0𝐹 [𝑛] = 𝐶 ∶ 𝑌_0𝐹 [𝑛] + 0.9𝑌_0𝐹 [𝑛 − 2] = 1
 𝐶+ 0= 1
 𝐶=1
 𝑌_0 [𝑛] = 𝑘_1 (0 + 0.9487𝑖)𝑛 + 𝑘_2 (0 − 0.9487𝑖)𝑛−1
 𝑌_0 [−1] = 0 = 𝑘_1 (0 + 0.9487𝑖)^(−1) + 𝑘_2 (0 − 0.9487𝑖)^(−1) + 1
 𝑌_0 [−2] = 0 = 𝑘_1 (0 + 0.9487𝑖)^(−2) + 𝑘_2 (0 − 0.9487𝑖)^(−2) + 1
 𝑘1 = 0.4972 − 0.4246𝑖 𝑘2 = 0.4972 + 0.5241𝑖
 𝑌_0 [𝑛] = [(0.4972 − 0.4246𝑖 ) (0 + 0.9487𝑖)(𝑛) + (0.4972 + 0.5241𝑖) (0 − 0.9487𝑖)^𝑛 + 1]𝛿[𝑛]
 𝑌_𝑧𝑠 [𝑛] = 0.3𝑌_0 [𝑛] + 0.6𝑌_0 [𝑛 − 1] + 0.3𝑌_0 [𝑛 − 2]
 𝑌_𝑧𝑠 [𝑛] = [(0.14916 − 0.12738𝑖 ) (0 + 0.9487𝑖)𝑛 + (0.14916 + 0.157231𝑖) (0 − 0.9487𝑖)𝑛 + 0.3]𝛿[𝑛] +
[(0.29832 − 0.25476𝑖 ) (0 + 0.9487𝑖)(𝑛+2) + (024832 + 0.31446𝑖) (0 − 0.9487𝑖)(𝑛−1) + 0.6]𝛿[𝑛 − 1] +
 [(0.14916 − 0.12738𝑖 ) (0 + 0.9487𝑖)(𝑛−2) + (0.14916 + 0.15723𝑖) (0 − 0.9487𝑖)(𝑛−2) + 0.3]𝛿[𝑛 − 2]
 𝑌_𝑧1 [𝑛]: 𝑌_𝑧1 [𝑛] + 0.9𝑌_𝑧1 [𝑛 − 2] = 0
 𝑌_𝑧1 [𝑛] = 𝑘_1 (0 + 0.9487𝑖)(𝑛) + 𝑘_2 (0 − 0.9487𝑖)(𝑛)
 𝑌_𝑧1 [−1] = 0 = 𝑘_1 (0 − 0.9487𝑖)(−1) + 𝑘_2 (0 − 0.9487𝑖)(−1) = 0
 𝑌_𝑧1 [−2] = 0 = 𝑘_1 (0 + 0.9487𝑖)−2 + 𝑘_2 (0 − 0.9487𝑖)−2 = 0
 𝑘_1 = 0.4972 − 0.4246𝑖 𝑘_2 = 0.4972 + 0.5241𝑖
 𝐹𝑖𝑛𝑎𝑙𝑚𝑒𝑛𝑡𝑒:
 𝑌[𝑛] = [(0.14916 − 0.12738𝑖 ) (0 + 0.9487𝑖)^𝑛 + (0.14916 + 0.157231𝑖) (0 − 0.9487𝑖)^𝑛 + 0.3]𝛿[𝑛] +
[(0.29832 − 0.25476𝑖 ) (0 + 0.9487𝑖)^(𝑛 − 1) + (024832 + 0.31446𝑖) (0 − 0.9487𝑖)^(𝑛 − 1) + 0.6]𝛿[𝑛 −
1] +
 [(0.14916 − 0.12738𝑖 ) (0 + 0.9487𝑖)𝑛−2 + (0.14916 + 0.15723𝑖) (0 − 0.9487𝑖)𝑛−2 + 0.3]𝛿[𝑛 − 2]
 +[(0.4972 − 0.4246𝑖 ) (0 + 0.9487𝑖)𝑛 + (0.4972 + 0.5241𝑖) (0 − 0.9487𝑖)𝑛 ]𝛿[𝑛]
VI. CONCLUSIONES y OBSERVACIONES
 Se logró solucionar los problemas propuestos con ayuda del software ya que es muy necesario al momento de realizar
este tipo de ejercicios
 Para realizar la parte analítica de los problemas es necesario ser ordenado y seguir el orden ya que es muy simple
equivocarse el arrastrar errores
 Para comenzar la mayoría de los problemas es necesario declarar los vectores y establecer el rango que tomara n para
que al momento de graficar sea mucho más simple y también nos ayudar a ver la gráfica de manera más cómoda
 Es necesario para el entendimiento correcto de la gráfica hay casos en los que la gráfica queda muy a la izquierda y no
se puede apreciar de manera óptima en estos casos se tiene que colorar la herramienta que nos brinda Matlab “axis”para
poder centrar ala grafica al intervalo que necesitamos
 Se concluye que los resultados tanto como los por recursión, filled y la parte analítica dan resultados iguales de esta
manera podemos comprobar que los métodos usados son correctos y se puede optar por cualquiera de ellos para
solucionar problemas de este tipo

Cuadros Enviar Informe 2

Cargado por

Copyright:

Formatos disponibles

Cuadros Enviar Informe 2

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Cuadros Enviar Informe 2

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD CATOLICA DE SANTA MARIA Página:1/2

FACULTAD DE CIENCIAS E INGENIERÍAS FISICAS Y FORMALES

II. MARCO TEORICO

II.1. SISTEMAS DE TIEMPO DISCRETO

III. MATERIAL Y EQUIPO

III.1. Una PC con SO Windows y MATLAB con el toolbox de PDS instalado.

𝑦[𝑛] + 0.9𝑦[𝑛 − 2] = 0.3𝑥[𝑛] + 0.6𝑥[𝑛 − 1] + 0.3𝑥[𝑛 − 2]

IV.2. Calcule analíticamente y[n] para x[n]=δ[n]

Práctica 2 – Tema 1: SISTEMAS DE TIEMPO DISCRETO P N° 2 – T N° 01 JP: Ing. JCC

obtenida con el script.

for i=1:length(n) %establesemos la condicion de conteo para el bucle

x=(n==0).*(1)+0;%declamramos que x sera 1siempre que n sea 0

for z=3:length(n)%establesemos la condicion de conteo para el bucle

stem(n,y)%aplicamos stem para graficar

axis([0,30,-0.8,0.8])%acomodamos la visibilidad de la grafica

 Calcule las 10 primeras muestras de su respuesta impulsional

VI. CONCLUSIONES y OBSERVACIONES

También podría gustarte