Infinite-Impulse-Response-Systeme - Digitale Signalverarbeitung mit MATLAB - page 111

Programmbeispiel 10-1 IIR- Filter 2. Ordnung in der transponierten Direktform II
% 2nd order IIR system in transposed direct form II
% function [y,sf] = iird2t(b,a,x,si)
% x : input signal
% b : numerator coefficients b = [b0,b1,b2]
% a : denominator cefficients a = [a0,a1,a2]
% si : initial conditions si = [s1,s2]
% y : output signal
% sf : final conditions
% iirdf2t.m * mw * 04/30/2008
function [y,sf] = iird2t(b,a,x,si)
if nargin==4
SI = si; % initial conditions

else
SI = [0,0];
end
a = a/a(1); b = b/a(1); % normalize coefficients to a0 = 1

y = zeros(size(x)); % allocate memory

% Filtering

Im

A10.6-3 Bestimmen Sie die Pole und
Nullstellen des Systems H ( z )
aus (10.14)

z

z 01,2 =

Re

1

z f1,2 =

und tragen Sie die Pole und
Nullstellen in Bild 10-4 ein.

Bild 10-4 Pol-Nullstellendiagramm zu H ( z ) in (10.14)