Помогите пожалуйста. Заданы параметры фильтры (частоты полосы пропускания, затухания, а так же затухания в полосах пропускания и задерживания). Написал программу которая выводит АЧХ и ФЧХ фильтра. Помогите сделать билинейное Z-преобразование на основе этой программы (необходимо получить коэффициенты для H(Z) - причем для двух звеньев)

Rs=40;
Rp=1;
Wp=[2*pi*1000*3.320 2*pi*1000*3.565];
Ws=[2*pi*1000*2.15 2*pi*1000*9];
[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,'s');
[B,A] = CHEBY1(N,0.5,Wn,'bandpass','s');
W = 0:1:40000;
H = freqs(B, A, W);
subplot(2,1,1);
plot(W, abs(H)), grid
title('АЧХ');
subplot(2,1,2);
plot(W, angle(H)), grid
title('ФЧХ');