thiết kế bộ lọc số Elliptic Matlab

Elliptic 1) Lowpass Elliptic Transfer Function

[b,a] = ellip(6,5,40,0.6);

freqz(b,a)

dataIn = randn(1000,1);

dataOut = filter(b,a,dataIn);

2) Bandstop Elliptic Filter

[b,a] = ellip(3,5,50,[0.2 0.6],'stop');

freqz(b,a)

dataIn = randn(1000,1);

dataOut = filter(b,a,dataIn);

3) Highpass Elliptic Filter

[z,p,k] = ellip(6,3,50,300/500,'high'); sos = zp2sos(z,p,k);

fvtool(sos,'Analysis','freq')

4) Bandpass Elliptic Filter

[A,B,C,D] = ellip(10,3,40,[500 560]/750);

d = designfilt('bandpassiir','FilterOrder',20,

'PassbandFrequency1',500,'PassbandFrequency2',560,

'PassbandRipple',3,

'StopbandAttenuation1',40,'StopbandAttenuation2',40,

'SampleRate',1500);

sos = ss2sos(A,B,C,D);

fvt = fvtool(sos,d,'Fs',1500);

legend(fvt,'ellip','designfilt')