Laboratorio De Pds

%6.1

% Vectores a y b que contienen los coeficientes de la ecuación dediferencias

a=[1 -1];

b=[1 -0.2];

% Respuesta en frecuencia

figure(1)

[H,w] = freqz(b,a,128,'whole');

plot(H,w);

xlabel('W (rad)');

ylabel('H(w)');

title('Respuesta en frecuencia del Sistema');

% Cálculo de H(z)

z=exp(i*2*pi*24/128)

Hz1=(1-(0.2/z))/(1-1/z)

% Cálculo de H(z)

z=1.2

Hz=(z.^-1)-(z.^-2)-0.48*(z.^-3)+0.288*(z.^-4)

%ejercicio 6.2

h=[0.03 0.4 0.54 0.2 -0.2 0.1 0.2]

z1 = j; z2 = 1 + j; z3 = 1 - j; z4 = 2; z5 = -2; z6 = -j;

zetas = [j, 1+j, 1-j, 2, -2, -j];

e1 = zetas(1).^(0:47);

e2 = zetas(2).^(0:47);

e3 = zetas(3).^(0:47);

e4 = zetas(4).^(0:47);

e5 = zetas(5).^(0:47);

e6 = zetas(6).^(0:47);

y1 = conv(e1,h);

y2 = conv(e2,h);

y3 = conv(e3,h);

y4 = conv(e4,h);

y5 = conv(e5,h);

y6 = conv(e6,h);

figure(1)

plot([0:length(y1)-1],real(y1))

Hz1=0.03+0.4*(z1^-1)+0.54*(z1^-2)+0.2*(z1^-3)-0.2*(z1^-4)+0.1*(z1^-5)+0.2*(z1^-6)

Hz2=0.03+0.4*(z2^-1)+0.54*(z2^-2)+0.2*(z2^-3)-0.2*(z2^-4)+0.1*(z2^-5)+0.2*(z2^-6)

Hz3=0.03+0.4*(z3^-1)+0.54*(z3^-2)+0.2*(z3^-3)-0.2*(z3^-4)+0.1*(z3^-5)+0.2*(z3^-6)

Hz4=0.03+0.4*(z4^-1)+0.54*(z4^-2)+0.2*(z4^-3)-0.2*(z4^-4)+0.1*(z4^-5)+0.2*(z4^-6)

Hz5=0.03+0.4*(z5^-1)+0.54*(z5^-2)+0.2*(z5^-3)-0.2*(z5^-4)+0.1*(z5^-5)+0.2*(z5^-6)

Hz6=0.03+0.4*(z6^-1)+0.54*(z6^-2)+0.2*(z6^-3)-0.2*(z6^-4)+0.1*(z6^-5)+0.2*(z6^-6)

%ejercicio 6.3

Yz=[1 -1.6 1 -0.208];

Xz=[1 -1 0 2];

[H,w]=freqz(Xz,Yz,128,'whole');

figure(1)

plot(H)

title('frecuencia H')

grid

ceros=roots([Yz]);

polos=roots([Xz]);

fase=(0:0.01:2*pi);

z=exp(j*fase);

axis=([-3 3 -2 2]);

figure (2)

plot(real(ceros),imag(ceros),'o',real(polos),imag(polos),'x',real(z),imag(z),'r');

xlabel('El sistema es estable al tener ceros dentro del circulo unidada ')

grid

%ejercicio 6.4

a1=-2;

Yz=[1 0];

Xz=[1 a];

ceros=roots([Yz]);

polos=roots([Xz]);

fase=(0:0.01:2*pi);

z=exp(j*fase);

figure(1)

plot(real(ceros),imag(ceros),'o',real(polos),imag(polos),'x',real(z),imag(z),'r');

grid

a2=-1

Yz=[1 0];

Xz=[1 a2];

ceros=roots([Yz]);

polos=roots([Xz]);

fase=(0:0.01:2*pi);

...