clear all;

Fs = 10000;
fp = 1500;
fs = 1300;
rp = 2;
rs = 35;

[sellip,wnellip] =  ellipord( fp/(Fs/2), fs/(Fs/2), rp, rs);
[bellip,aellip]  =  ellip(sellip, rp, rs, wnellip, 'high');

figure(1)
subplot(2,1,1)
freqz(bellip,aellip)

% parametry
f = 1000;
% casova osa
N = 400;
n = 0:N-1;
t = n(:)/Fs;
% generace signalu
x = cos( (-pi/3)*f*t);

figure(2)
subplot(2,1,1)
stem(x)
subplot(2,1,2)
zplane(bellip,aellip)

%Josef Prorok