filtrospb

TRABALHO DE G2
Acadêmico : CEDINIR FABICZ

QUESTÃO 3 % a) [z,p,k] = buttap(5); % Conforme atenuaçao desejada
[bn,an] = zp2tf(z,p,k); % gera numerador e denominador do filtro normalizado
Gn = tf(bn,an); % ft do filtro normalizado f = 1:1:100; % intervalo do gráfico wc = 100; fc = wc/(2*pi); % frequencia corte do filtro PB
[b,a] =lp2lp(bn,an,fc); % desnormaliza para frequencia wc
G = tf(b,a); % plota ft na linha comando
G=freqs(b,a,f)
semilogx(f,20.*log10(abs(G))); xlabel('Frequencia em Hz - log'); ylabel('Magnitude em dB'); title('Filtro Analogico PB'); grid função transferência

Gs = 1e10 ------------------------------------------------------------ s^5 + 323.6 s^4 + 5.236e04 s^3 + 5.236e06 s^2 + 3.236e08 s + 1e10

gráfico questão 3 letra a

QUESTÃO 3 LETRA B
[z,p,k] = buttap(5); % Conforme atenuaçao desejada de
[bn,an] = zp2tf(z,p,k); % gera numerador e denominador do filtro normalizado
Gn = tf(bn,an); % ft do filtro normalizado f = 1:1:100; % intervalo do gráfico wc = 100; fc = wc/(2*pi); % frequencia corte do filtro PB
[b,a] =lp2lp(bn,an,100); % desnormaliza para frequencia wc
Gs = tf(b,a)

%discretiza filtro PB fs = 1000;
Ts = 1/fs;
[bz,az] = bilinear(b,a,fs) % plota parametros de Gz, discretiza por TUSTIN
% Gz = c2d(Gs,Ts,'tustin') % Plota ft do filtro na linha comando (potencias positivas de z!)
Gz = tf(bz,az,Ts,'variable','z^-1') % plota ft com coeficientes do filtro na linha de comando potencias negativas de z bz = [ 2.658e-07 1.329e-06 2.658e-06 2.658e-06 1.329e-06 2.658e-07 ]; % copiar os coeficientes da ft anterior aqui: numerador az = [ 1 -4.677 8.758 -8.21 3.852 -0.7236 ]; % copiar os coeficientes da ft anterior aqui: denominador f = 1:1:100
[Gz,f] = freqz(bz,az,500,fs); % 500 pontos do