1. uvod u digitalnu sliku - Laboratorija za digitalnu obradu signala

=n+1;
n_l=7; % red Lagranzovog polinoma
b_l=intfilt(r,n_l,'Lagrange');
[h_l,w_l]=freqz(b_l,1);
y_l=filter(b_l,1,y2);
figure
set(figure(3),'Position',[70 60 900 600],'Name','LAGRANZOVA
INTERPOLACIJA')
prilog
subplot(221)
plot(w_l,abs(h_l))
axis([0 3.5 0 2.2])
title('Frekvencijska karakteristika Lagranzovog filtra','FontSize',14)
xlabel(' frekvencija [rad/odmjerku]')
ylabel('Magnituda')
grid on
subplot(222)
axis off
title('KARAKTERISTIKE FILTRA','FontSize',14)
text(-0.2, 0.9,'FIR NF filtar sa linearnom fazom','FontSize',14)
text(-0.2, 0.75,' za Lagranzovu
interpolaciju','FontSize',14)
text(-0.2,0.6,['Zahtjevi: Faktor interpolacije R =
',num2str(r)],'FontSize',14)
text(-0.2,0.45,['Stepen Lagranzovog polinoma N =
',num2str(n_l)],'FontSize',14)
text(-0.2,0.3,'Napomena: Ulazni signal za ovaj filtar mora da
sadrzi','FontSize',12)
text(-0.2,0.2,'nulte odmjerke u vremenskim trenucima u kojim
','FontSize',12)
text(-0.2,0.1,'je potrebno interpolirati signal','FontSize',12)
% plot(w_l,angle(h_l))
% title('Fazna karakteristika Lagranzovog filtra','FontSize',14)
% xlabel(' frekvencija [Hz]')
% ylabel('Fazno kasnjenje filtra [rad]')
% grid on
subplot(223)
stem(x1(1:2*fs+n_l),y_l(1:2*fs+n_l),'.')
title('Signal nakon filtriranja','FontSize',14)
xlabel(' vrijeme [s]')
ylabel('Amplituda signala')
grid on
P_l=abs(fft(y_l,512))/length(y_l);
subplot(224)
plot(fyy1,P_l(1:512))
title('Spektar signala nakon filtriranja','FontSize',14)
xlabel(' frekvencija [Hz]')
ylabel('Magnituda')
grid on
74