s=tf('s');
%% Filtr dolnoprzepustowy 300Hz
f1=300;
wg_dp1=2*pi*f1;
S_dp1=s/wg_dp1
%% Filtr dolnoprzepustowy 2kHz
f2=2000;
wg_dp2=2*pi*f2;
S_dp2=s/wg_dp2
%% 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))
%% 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)))