clear all
clc

%pod a) NF filtar 25 reda fg = 2000hz koriscenjem hanov prozor

%specifikacija NF filtra
M1 = 26;%impulsni odziv ogranicimo na 26 odbirka, jer je red filtra 25
Wn1 = (2000 / 11200)*2;  %kriticna ucestanost NF filtra 2*pi*2000/11200=4000/11200*pi
%odbirci prozorskih funkcija koje se koriste
hanov = hann(M1);
%projektovanje FIR filtara koriscenjem funkcije fir1
N1 = M1-1; %potrebni red filtra
disp(length(hanov));
disp(N1);
b1 = fir1 (N1, Wn1, 'low', hanov);
%racunanje spektra u 1024 tacke koriscenjem FFT
N_fft1 = 1024;
B1 = fft (b1, N_fft1);
%odredjivanje amplitudskih karakteristika
B1a = abs(B1(1:N_fft1/2));
B1b = angle(B1(1:N_fft1/2));
%crtanje amplitudskih karakteristika na linearnoj skali
n1 = 0:N_fft1/2-1;
w1 = n1*pi/(N_fft1/2-1);
subplot (3, 2, 1), plot (w1, B1a), title ('NF filtar amplitudska k-ka - Hanov prozor');
subplot (3, 2, 2), plot (w1, B1b), title ('NF filtar fazna k-ka - Hanov prozor');



%pod b) VF filtar 23 reda fg = 4000hz koriscenjem hamingov prozor

M2 = 25; %impulsni odziv ogranicimo na 24 odbirka, jer je red filtra 23
Wn2 = (4000 / 11200)*2;  %kriticna ucestanost VF filtra 2*pi*4000/11200=8000/11200*pi
%odbirci prozorskih funkcija koje se koriste
hamingov = hamming(M2);
%projektovanje FIR filtara koriscenjem funkcije fir1
N2 = M2 - 1; %potrebni red filtra
disp(length(hamingov));
disp(N2);
b2 = fir1 (N2, Wn2, 'high', hamingov);
%racunanje spektra u 1024 tacke koriscenjem FFT
N_fft2 = 1024;
B2 = fft (b2, N_fft2);
%odredjivanje amplitudskih karakteristika
B2a = abs(B2(1:N_fft2/2));
B2b = angle(B2(1:N_fft2/2));
%crtanje amplitudskih karakteristika na linearnoj skali
n2 = 0:N_fft2/2-1;
w2 = n2*pi/(N_fft2/2-1);
subplot (3, 2, 3), plot (w2, B2a), title ('VF filtar amplitudska k-ka - Hamingov prozor');
subplot (3, 2, 4), plot (w2, B2b), title ('VF filtar fazna k-ka - Hamingov prozor');

%pod c) PO filtar 33 reda fg1 = 2000hz i f2 = 4000hz koriscenjem blekmenovog prozor

M3 = 34; %impulsni odziv ogranicimo na 34 odbirka, jer je red filtra 33
Wn3 = [(2000 / 11200)*2 (4000 / 11200)*2];
%odbirci prozorskih funkcija koje se koriste
blekmenov = blackman(M3);
%projektovanje FIR filtara koriscenjem funkcije fir1
N3 = M3 - 1; %potrebni red filtra
b3 = fir1 (N3, Wn3, blekmenov);
%racunanje spektra u 1024 tacke koriscenjem FFT
N_fft3 = 1024;
B3 = fft (b3, N_fft3);
%odredjivanje amplitudskih karakteristika
B3a = abs(B3(1:N_fft3/2));
B3b = angle(B3(1:N_fft3/2));
%crtanje amplitudskih karakteristika na linearnoj skali
n3 = 0:N_fft3/2-1;
w3 = n3*pi/(N_fft3/2-1);
subplot (3, 2, 5), plot (w3, B3a), title ('PO filtar amplutudska k-ka - Blekmenov prozor');
subplot (3, 2, 6), plot (w3, B3b), title ('PO filtar fazna k-ka - Blekmenov prozor');

% n = 0 : 1023;
%
%
% %ulaz
% x = cos(2*pi*1000/11200*n) + cos(2*pi*3000/11200*n) + cos(2*pi*5000/11200*n);
%
% figure;
% subplot (3, 2, 1),stem (n, x); title ('Ulaz');
%
% C3 = fft(x, N_fft3);
% subplot (3, 2, 2),stem (n, C3); title ('fft ulaznog signala');
% C3a = abs(C3(1:n/2));
% subplot (3, 2, 3),stem (n, C3a); title ('amplitudska k-ka ulaznog signala');
% C3b = angle(C3(1:n/2));
%
% %subplot (3, 2, 2), stem (n, C3a);
% %subplot (3, 2, 3), stem (n, C3b);