% Some filter parameters from FDATool
% Lowpass Chebyshev type I, 10th order, Fs=48000, Fpass=9600
B = [2.9022458e-005  2.9022458e-004  1.3060106e-003  3.4826950e-003 ...
     6.0947162e-003  7.3136595e-003  6.0947162e-003  3.4826950e-003 ...
     1.3060106e-003  2.9022458e-004  2.9022458e-005];
A = [1.0000000e+000 -5.3138256e+000  1.4970401e+001 -2.8187537e+001 ...
     3.8667812e+001 -4.0030401e+001  3.1547185e+001 -1.8678371e+001 ...
     7.9749478e+000 -2.2338684e+000  3.1700222e-001];

% Draw filter frequency response
fh = figure;
freqz(B,A);

j = sqrt(-1);
tv = 0:200;
x1 = real(exp(j*0.1*pi*tv));
y1 = filter(B,A,x1);
x2 = real(exp(j*0.2*pi*tv));
y2 = filter(B,A,x2);
x3 = real(exp(j*0.3*pi*tv));
y3 = filter(B,A,x3);
x4 = real(exp(j*0.4*pi*tv));
y4 = filter(B,A,x4);
