s=tf('s');
L=[1];
M5=[1 3.236 5.235 5.235 3.236 1];
H5=tf(L,M5)
%% Filtr dolnoprzepustowy 300Hz
f1=300;
wg_dp1=2*pi*f1;
s_dp1=s/wg_dp1
G_dp1=1/(s_dp1^5+3.236*s_dp1^4+5.235*s_dp1^3+5.235*s_dp1^2+3.236*s_dp1+1)
%% Filtr dolnoprzepustowy 2kHz
f2=2000;
wg_dp2=2*pi*f2;
s_dp2=s/wg_dp2
G_dp2=1/(s_dp2^5+3.236*s_dp2^4+5.235*s_dp2^3+5.235*s_dp2^2+3.236*s_dp2+1)
%% Filtr pasmowoprzepustowy
f3=1500;
f4=3000;
wg1_p=2*pi*f3;
wg2_p=2*pi*f4;
w0_p=sqrt(wg1_p*wg2_p);
beta_p=w0_p/(wg2_p-wg1_p);
s_p=beta_p*((s/w0_p)+(w0_p/s))
G_p=1/(s_p^5+3.236*s_p^4+5.235*s_p^3+5.235*s_p^2+3.236*s_p+1)
%% Filtr pasmowozaporowy
f5=1500;
f6=3000;
wg1_z=2*pi*f5;
wg2_z=2*pi*f6;
w0_z=sqrt(wg1_z*wg2_z);
beta_z=w0_z/(wg2_z-wg1_z);
s_z=1/(beta_z*((s/w0_z)+(w0_z/s)))
G_z=1/(s_z^5+3.236*s_z^4+5.235*s_z^3+5.235*s_z^2+3.236*s_z+1)