VRBA GEJ

% Parametri filtera
f_sampling = 3000; % Hz, frekvencija odmeravanja
f_cutoff = 500;    % Hz, frekvencija odsecanja
Rp = 5;            % dB, maksimalno dozvoljeno slabljenje u propusnom opsegu
Rs = 30;           % dB, minimalno dozvoljeno slabljenje u nepropusnom opsegu

% Pretvaranje frekvencija u normalizovane frekvencije
Wp = (2 * f_cutoff) / f_sampling;

% Projektovanje eliptičkog filtera
[n, Wn] = ellipord(Wp, 1.5 * Wp, Rp, Rs);
[b, a] = ellip(n, Rp, Rs, Wn, 'high');

% Prikazivanje karakteristika filtera
freqz(b, a, 1024, f_sampling);

% Filtriranje ulaznog signala
t = 0:1/f_sampling:1; % vreme od 0 do 1 sekunde
x = sin(2*pi*300*t) + sin(2*pi*1000*t); % primer ulaznog signala

y = filter(b, a, x);

% Prikazivanje rezultata
figure;
subplot(2,1,1);
plot(t, x);
title('Ulazni signal');
subplot(2,1,2);
plot(t, y);
title('Filtrirani signal');