Matlab Transformada Z

1. (P4.
1) Determine la transformada z de las siguientes secuencias

usando su definicin. Indique la regin de convergencia de cada
secuencia y verifique las expresiones usando Matlab.
a.
x ( n )={3,2, 1 ,2.3 }
b.
x ( n )=( 0.8 )n u ( n2 )
c.
d.
e.
n
4
(
)
xn=
u ( 1n )
3
()
x ( n )=2|1|+
1
3
|n|
()
x ( n )=( n+1 ) ( 3 )n u (n)
2. (P4.2) Determine la transformada Z de las siguientes secuencias

usando la tabla de transformada Z y las propiedades de la
misma. Exprese X (Z) como una funcin racional en z-1. Verifique
los resultados usando Matlab. Indique la regin de convergencia
en cada caso y grafique los polos y ceros.
a.
x ( n )=2 ( n2 )+ 3u (n3)
b.
1 n (
n3
(
)
xn=
u n2 )+ ( 0.9 ) u ( n )
3
c.
d.
e.
()
x ( n )=nsin
( n3 )u ( n)+ ( 0.9) u ( n2)

n
1 n
n
cos
45 u ( n1 )
2
4
( ) (( )
x ( n )=
x ( n )=( n3 )
1
4
n2
()
cos
{( ) }
( n1 ) u(n)
2
3. (P4.3) La transformada z de
x (n)
es
X ( z )=( 1+2 z1 ) ,|z|0.
encuentre la transformada z de las siguientes secuencias e

indique su regin de convergencia.
a.
x 1 ( n )=x ( 3n ) + x ( n3 )
b.
x 2 ( n )= ( 1+ n+n2 ) x ( n )
1
x ( n2 )
2
()
c.
x 3 ( n) =
d.
x 4 ( n )=x ( n+2 )x ( n2 )
e.
x 5 ( n ) =cos
( n2 ) x( n)
4. (P4.6)
Si
las
secuencias
x 1 ( n ) , x 2 ( n ) , estnrelacionadas por x 3 ( n )=x 1 ( n )x 2 ( n ) , entonces:
x 3 ( n )=
n=
x 1(n)
n=
)(
n=
x 2(n)
a. Probar el resultado anterior sustituyendo por la definicin de la

convolucin en el lado izquierdo.
b. Probar el resultado anterior usando la propiedad de convolucin.
c. Verificar el resultado anterior usando Matlab pero eligiendo alguna de las
dos secuencias x1(n) y x2(n).
5. (P4.7) Determine los resultados de las siguientes operaciones de
polinomios usando Matlab.
a.
X 1 ( z )= ( 12 z1 +3 z24 z3 ) ( 4+3 z12 z2 + z3 )
b.
X 2 ( z ) =( z 22 z +3+2 z1+ z 2 ) ( z 3z3 )
c.
1+ z + z
X 3 ( z )=
d.
X 4 ( z ) =X 1 ( z ) X 2 ( z ) X 3 ( z )
e.
X 5 ( z ) =( z13 z3 +2 z5+5 z7z 9 ) ( z +3 z2 +2 z 3+ 4 z 4 )
f.
6. (P4.8) La funcin deconv es til cuando se dividen dos
secuencias causales. Escriba una funcin en Matlab deconv_m
para dividir dos secuencias no causales (similar a la funcin
conv). Evale su funcin con la siguiente operacin.
z + z +1+ z + z + z
3 z ++3 z
=( z 1+ 2 z12 z2) +
1
z+ 2+ z
z +2+ z1
7. (P4.11) Para los sistemas lineales e invariantes en el tiempo

descritos por la respuesta al impulso que se enlistan debajo,
determine (i) la funcin de representacin del sistema, (ii) la
representacin como ecuacin diferencial, (iii) la grfica de
polos y ceros, y (iv) la salida y(n) si la entrada es
a.
h ( n )=2
b.
h ( n )=n
1 n ( )
un
2
()
c.
1 n
u(n) .
4
()
x ( n )=
1
1 ( )
u ( n) +
un
3
4
()
h ( n )=3 ( 0.9 )n cos
( )
( n4 + 3 )u ( n )
d.
h ( n )=n [u ( n )u ( n10 ) ]
e.
h ( n )= [ 2sin ( n ) ] u( n)
8. (P4.15) Un sistema estable tiene las siguientes localizaciones de

sus polos y ceros:
1 1
1
1
z 1= j , z 2= j , p 1= + j , p 2=
j
2 2
2
2
Si tambin es conocido que la respuesta en frecuencia de la funcin
H ( e j 10 )
evaluada en =9 es igual a 0.8; esto es:
H ( e j 10 )=0.8
a. Determine la funcin del sistema H (z) e indique su regin de
convergencia.
b. Determine la representacin como ecuacin diferencial.
c. Determine la respuesta en estado estacionario Yss(n) si la entrada es
n
2
) U(n)
1
x ( n )= sin
2
d. Determine la respuesta transitoria si la entrada es :
n
2
) U(n)
1
x ( n )= sin
2
9. (P4.16) Un filtro digital es descrito por la ecuacin diferencial:
y ( n )=x ( n ) + x ( n1 ) +0.9 y ( n1 )0.81 y ( n2 )

a. Usando la funcin freqz, grafique la magnitud y fase de la respuesta en
frecuencia del filtro anterior. Haga notar la magnitud y la fase en =/3
y =.
En
/3 la magnitud es 10.5215 y la fase es -58.2595.
En
la magnitud es 0 y la fase es -90.
b. Genere 200 muestras de la seal
x ( n )=sin
( n3 )+5 cos ( n)
y procsela a
travs del filtro. Compare con la porcin en estado estacionario de la

salida de x(n). Cmo son las amplitudes y las fases de las dos
sinusoides afectadas por el filtro?
La seal de salida es una senoidal de frecuencia
/3
y de magnitud
amplificada por 10. En cambio la seal cosenoidal de frecuencia
desaparece en la seal de salida. Esto significa que el filtro atenuar por

completo aquellas seales mayores o iguales a una frecuencia de
rad/s. Para comprobar que la seal de salida tiene una frecuencia de
/ 3 , se realiz un anlisis espectral mediante la transformada de
Fourier.
En la figura el valor de frecuencia para el armnico de magnitud ms

grande es de 1.05 rad/s que es aproximadamente igual a
Cdigo:
%a)Respuestaalafrecuenciade
%y(n)=x(n)+x(n1)+0.9y(n1)0.81y(n2)
clearall;close;clc;
a=[10.90.81];b=[11];
w=linspace(0,pi,500);
H1=freqz(b,a,w);
mag1=abs(H1);ang1=angle(H1);
subplot(2,1,1);plot(w/pi,mag1);
title('RespuestaenMagnitud');gridon;
ylabel('|H1|');xlabel('Frecuencia0a\pi');
subplot(2,1,2);plot(w/pi,ang1*180/pi);
/3 .
title('RespuestaenFase');gridon;
ylabel('Grados');xlabel('Frecuencia0a\pi');
%Calculandolosvaloresenpi/3ypi
w2=[pi/3pi];H2=freqz(b,a,w2);
mag2=abs(H2)
ang2=angle(H2)*180/pi
pause();
%b)Respuestadelfiltroalasealdeentrada:
%x(n)=sin(pi*n/3)+5cos(pi*n)
n=0:200;x=sin(pi*n/3)+5*cos(pi*n);
figure;subplot(2,1,1);
plot(n,x);title('Sealdeentrada');
y=filter(b,a,x);subplot(2,1,2);
plot(n,y);title('Respuestaenestadoestacionario');
pause();
%Anlisisespectraldelarespuestaen
%estadoestacionario
L=length(y);N=2^nextpow2(L);
Y=fft(y,N)/L;
f=linspace(0,pi,N/2+1);
figure;plot(f,2*abs(Y(1:N/2+1)));
title('EspectrodefrecuenciasdelaRespuestaenEstadoEstacionario');
xlabel('Frecuencia(rad/s)');xlim([f(1)f(length(f))]);
10.
Resuelva las siguientes ecuaciones diferenciales para y(n)
usando el enfoque unilateral de la transformada z.
y ( n )=0.5 y ( n1 ) + 0.25 y ( n2 )+ x ( n ) , n 0; y (1 )=1, y (2 )=2

x ( n )=( 0.8 )n u ( n )
Genere las primeras 20 muestras de y(n) usando Matlab y comprelas
con sus resultados.
1. Aplicando transformada Z:
Y ( z ) 0.5 [ y (1 ) + z1 Y ( z ) ]0.25 [ y (2 ) + z1 y (1 ) + z2 Y ( z ) ]=
Y ( z ) ( 10.5 z 0.25 z )=10.25 z +
1
10.8 z1
Y (z )=
20.55 z 10.2 z2
11.3 z1 +0.15 z2 +0.2 z3
Y (z)=
0.1298
64
65.8685
+
1
1
1+0.3090 z
10.8 z
10.8090 z 1
2. Aplicando transformacin inversa:
1
10.8 z1
y ( n )=0.1298(3090) 64 ( 0.8 ) +65.8685(0.8090)
3. Los primeros valores del 0 al 20 usando matlab son:

[1.9983 2.0475 2.1621 2.1039 2.0013 1.8536
1.3481
1.1873
1.0378
0.9014
0.7787
0.4892 0.4160 0.3528 0.2984 0.2517 0.2120]
1.6887
0.6696
1.5171
0.5734
1.6915
0.5744
1.5196
0.4901
4. La solucin numrica es:

[2.0000 2.0500 2.1650 2.1070 2.0044
1.3505 1.1894
1.0397
0.9031
0.7802
0.4168 0.3535 0.2989 0.2522 0.2124]
1.8566
0.6708
Al comparar ambas salidas, el error absoluto fue de 0.0031.

Cdigo:
%Solucinanalticaynumricadeunaecuacinendiferencias
clearall;close;clc;
n=0:20;
%Solucinanaltica
y1=0.1298*(0.3090).^n64*(0.8).^n+65.8685*(0.8090).^n
%Solucinnumrica
b=1;a=[10.50.25];y0=[12];
x=(0.8).^n;
%Condicionesinicialesequivalentes
V=filtic(b,a,y0);
%Secuenciadesalida
y2=filter(b,a,x,V)
%Clculodelerrorabsoluto
abs(max(y2y1))

Matlab Transformada Z

Cargado por

Copyright:

Formatos disponibles

Matlab Transformada Z

Cargado por

Información del documento

Descripción 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

Matlab Transformada Z

Cargado por

Copyright:

Formatos disponibles

1. (P4.

1) Determine la transformada z de las siguientes secuencias

x ( n )=( n+1 ) ( 3 )n u (n)

2. (P4.2) Determine la transformada Z de las siguientes secuencias

( n3 )u ( n)+ ( 0.9) u ( n2)

X ( z )=( 1+2 z1 ) ,|z|0.

encuentre la transformada z de las siguientes secuencias e

x 1 ( n ) , x 2 ( n ) , estnrelacionadas por x 3 ( n )=x 1 ( n )x 2 ( n ) , entonces:

a. Probar el resultado anterior sustituyendo por la definicin de la

X 1 ( z )= ( 12 z1 +3 z24 z3 ) ( 4+3 z12 z2 + z3 )

X 2 ( z ) =( z 22 z +3+2 z1+ z 2 ) ( z 3z3 )

X 5 ( z ) =( z13 z3 +2 z5+5 z7z 9 ) ( z +3 z2 +2 z 3+ 4 z 4 )

7. (P4.11) Para los sistemas lineales e invariantes en el tiempo

h ( n )=3 ( 0.9 )n cos

8. (P4.15) Un sistema estable tiene las siguientes localizaciones de

evaluada en =9 es igual a 0.8; esto es:

d. Determine la respuesta transitoria si la entrada es :

y ( n )=x ( n ) + x ( n1 ) +0.9 y ( n1 )0.81 y ( n2 )

/3 la magnitud es 10.5215 y la fase es -58.2595.

la magnitud es 0 y la fase es -90.

b. Genere 200 muestras de la seal

travs del filtro. Compare con la porcin en estado estacionario de la

La seal de salida es una senoidal de frecuencia

amplificada por 10. En cambio la seal cosenoidal de frecuencia

desaparece en la seal de salida. Esto significa que el filtro atenuar por

rad/s. Para comprobar que la seal de salida tiene una frecuencia de

/ 3 , se realiz un anlisis espectral mediante la transformada de

En la figura el valor de frecuencia para el armnico de magnitud ms

y ( n )=0.5 y ( n1 ) + 0.25 y ( n2 )+ x ( n ) , n 0; y (1 )=1, y (2 )=2

Y ( z ) ( 10.5 z 0.25 z )=10.25 z +

2. Aplicando transformacin inversa:

y ( n )=0.1298(3090) 64 ( 0.8 ) +65.8685(0.8090)

3. Los primeros valores del 0 al 20 usando matlab son:

4. La solucin numrica es:

Al comparar ambas salidas, el error absoluto fue de 0.0031.

También podría gustarte