%% Drugie zakłócenie
Wpc22=[2475 2490]/(0.5*fs)
Wsc22=[2470 2495]/(0.5*fs)
Rpc22=3;
Rsc22=45;
[nc22,Wnc22]=cheb2ord(Wpc22,Wsc22,Rpc22,Rsc22)
[bc22,ac22]=cheby2(nc22,Rsc22,Wnc22,'stop')
GC22=tf(bc22,ac22,1/fs)
figure(4)
options = bodeoptions;
options.FreqUnits = 'Hz';
options.FreqScale = 'linear'
bode(GC22, options)
xlim([0 4700])