Ejercicios Transformada Z MT228
Transformada z
Ejercicios
Ricardo Rodriguez Bustinza
robust@[Link]
Motivación
Prof. Ricardo Rodríguez Bustinza 1
Ejercicios Transformada Z MT228
o La identidad de Euler es una expresión que se encuentra en el corazón de la teoría de
números compleja. La expresión plantea una pregunta tentadora: ¿qué pasa con un
número positivo cuando lo elevamos a una potencia imaginaria?
o La expresión posee el número “e” de Euler, la
base de logaritmos naturales que se recluta
extensamente en el cálculo. Es un número
trascendental cuyo valor es 2.71828... Posee “i”,
la unidad compleja o imaginaria, que es la raíz
cuadrada de -1 o la solución de la función
“x²+1=0”. Es de suma importancia en la
ingeniería y ha proporcionado grandes
conocimientos en mecánica cuántica.
e = cos θ + i sin θ e = (−1) + 0
e = cos π + i sin π e +1 =0
Ejercicio 1
o Considere una señal sinusoidal x(t) = 3cos(2πt + π/4), −∞ < t < ∞. Determine un
período de muestreo apropiado T de acuerdo con el muestreo de Nyquist se desea
obtener la señal de tiempo discreto x[kT] correspondiente al mayor período de
muestreo permitido.
Solución:
o Para muestrear x(t) para que no se pierda información, la tasa de muestreo de Nyquist
es una condición indica por el período de muestreo.
o Una señal de banda limitada x(t), es decir, su espectro de paso bajo X(ω) es tal que
donde ω es la frecuencia máxima en x(t), se puede muestrear uniformemente y sin
alias de frecuencia usando una frecuencia de muestreo llamado condición de frecuencia
de muestreo de Nyquist.
2π
ω = ≥ 2ω
T
Prof. Ricardo Rodríguez Bustinza 2
Ejercicios Transformada Z MT228
f=1;
T=1/f;
n=2;
w=2*pi*f;
t=0:0.001:n*T;
x=3*cos(w*t+pi/4);
% Ts<= pi/wmax
wmax=w;
Ts=pi/wmax;
fprintf('Ts=%1.2f\n',Ts)
tk=0:Ts:n*T;
xk=3*cos(wmax*tk+pi/4);
2π π
ω = ≥ 2ω T ≤
T ω
Ejercicio 2
o Un sistema en tiempo–discreto de primer orden tiene la ecuación en diferencias:
a y k +a y k−1 = b x k +b x k−1
o Asumiendo las condiciones iniciales diferentes de cero, encontramos la transformada-z
definida según:
b +b z b x −1 − a y(−1)
Y z = X z +
a +a z a +a z
o la función de transferencia del sistema discreto presenta un cero en z = −b /b , un
polo en z = −a /a , condiciones iniciales x(−1) y y(−1).
o Asuma los coeficientes a = 1, a = −0.9, b = 0, b = 0.1 y las condiciones iniciales
x(−1) = 1 y y(−1) = 0. Programe la solución para encontrar la respuesta del sistema
debido a una entrada escalón unitario. Use herramientas de MATLAB.
Prof. Ricardo Rodríguez Bustinza 3
Ejercicios Transformada Z MT228
% -----------------------------------------
% Fcn Recursivo
% a0y(k) + a1y(k-1) = b0u(k) + b1u(k-1)
% ----------------------------------------
% coeficientes % Fcn FILTER
T=0.01; tau=0.1; % ----------------------------------------
a0=1; a1=-(1-T/tau); zi=b1*xm1-a1*ym1;
b0=0; b1=T/tau; yf=filter(n,d,x,zi);
a=a1; subplot(412)
b=[b0 b1]; stem(k,yf,'r','filled')
% condiciones iniciales title('Fcn Filter')
xm1=1; ym1=0; xlabel('k')
% Indice de tiempo % -----------------------------------------
k=0:10; % Superposicion
% Entrada % -----------------------------------------
x=ones(1,length(k)); % Y(z) = Ya(z) + Yb(z) --> Y(z) = Gz1*X(z)
yr=recursivo(a,b,k,x,xm1,ym1); + Gz2*H(z)
%
y1=lsim(Gz1,x,k);
Gz2=tf([zi 0],d,-1);
xi=[1 zeros(1,length(k)-1)];
y2=lsim(Gz2,xi,k);
y=y1+y2;
Prof. Ricardo Rodríguez Bustinza 4
Ejercicios Transformada Z MT228
Función recursivo.m
function y = recursivo(a,b,n,x,x0,y0)
% y = recursivo(a,b,n,x,x0,y0)
% resuelve para y[n] desde:
% y[n] + a1*y[n-1] + a2*y[n-2]... + an*y[n-N]
% = b0*x[n] + b1*x[n-1] + ... + bm*x[n-M]
%
% a, b, n, x, x0 and y0 son vectores
% a = [a1 a2 ... aN]
% b = [b0 b1 ... bM]
% n contiene los valores d tiempo para que la solución sea computada
% y0 contiene las condiciones iniciales para y, en orden,
% i.e., y0 = [y[n0-N], y[n0-N+1], ...,y[n0-1]]
% donde n0 representa el primer elemento de n
% x0 contiene la condicion inicila en, en orden
% i.e., x0 = [x[n0-M],...,x[n0-1]]
% la salida, y, tiene length(n)
%
N = length(a);
M = length(b)-1;
if length(y0) ~= N,
error('Longitudes de a y y0 deben ser iguales')
end
if length(x0) ~= M,
error('Longitudes de x0 debe ser igual a la longitud de b-
1')
end
y = [y0 zeros(1,length(n))];
x = [x0 x];
a1 = a(length(a):-1:1); % reversa de elementos en a
b1 = b(length(b):-1:1);
for i=N+1:N+length(n),
% ecuacion en diferencias general
y(i) = -a1*y(i-N:i-1)' + b1*x(i-N:i-N+M)';
end
y = y(N+1:N+length(n));
Prof. Ricardo Rodríguez Bustinza 5
Ejercicios Transformada Z MT228
Ejercicio 3
o Similar al cálculo de la transformada de Laplace, el cálculo de la transformación Z se
puede hacer usando las herramientas de matemática simbólica de MATLAB.
(a) h k = 0.9u k syms k w0
% (a)
(b) h k = u k − u k − 10 h1=0.9.^k;
(c) h k = cos(ω k)u k H1 = ztrans(h1);
% H1 = z/(z-9/10);
% (b)
h2=heaviside(k)-heaviside(k-10);
H2= ztrans(h2);
H2=simplify(H2);
% 1/(z - 1) - (1/(z - 1) + 1/2)/z^10 + 1/2
% (c)
h3=cos(w0*k)*heaviside(k);
h3=ztrans(h3);
% (z*(z-cos(w0)))/(z^2-2*cos(w0)*z+1);
Ejercicio 4
o Halle la transformada z de x(k) use la identidad de Euler. x k = cos ω k + θ u(k)
o Aplicando la propiedad de linealidad y usando la transformación Z anterior cuando =
y su conjugado ∗ = obtenemos:
e( ) e ( )
x k = + u(k)
2 2
1 e e
X z = +
2 1−e z 1−e z
cos θ − cos ω − θ z
X z =
1 − 2 cos ω z + z
Prof. Ricardo Rodríguez Bustinza 6
Ejercicios Transformada Z MT228
o Analizar la transformada z en MATLAB cuando θ = 0 y es variable de 0.1 a 1.5.
Prof. Ricardo Rodríguez Bustinza 7
