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

n=2; % Faktor decimacije
%******** Realizacija filtra koristenog u funkciji decimate ********
redf=8; % Red filtra
Apo=0.05; % Slabljenje u PO u dB
omega=0.8; % granicna frekvencija
[b,a]=cheby1(redf, Apo,omega/n);
[h w]=freqz(b,a);
figure
subplot(221)
plot(w,10*log10(abs(h)))
axis([0 3.5 -150 10])
title('Frekvencijska karakteristika pred filtra','FontSize',14)
xlabel(' frekvencija [rad/odmjerku]')
ylabel('Slabljenje filtra [dB]')
text(0.2, -50,'Digitalni Chebishev NF filtar I tipa','FontSize',14)
text(0.2,-75,['Red filtra: N = ',num2str(redf)],'FontSize',14)
text(0.2,-100,['{\bf{\omega}_c}= ',num2str(omega/n),'{\pi}','
rad/odmjerku'],'FontSize',14)
text(0.2,-125,['A_p_o = ',num2str(Apo),'dB'],'FontSize',14)
grid on
c=filtfilt(b,a,y);
subplot(222)
stem(x(1:2*fs),c(1:2*fs),'.')
title('Signal nakon filtriranja','FontSize',14)
xlabel(' vrijeme [s]')
ylabel('Amplituda signala')
grid on
z=decimate(y,n);
x1=10*[0.5:1:length(z)-0.5]/length(z);
subplot(223)
stem(x1(1:fs),z(1:fs),'.')
xlabel(' vrijeme [s]')
ylabel('Amplituda signala')
title(['Signal nakon decimacije pri n = ',int2str(n)],'FontSize',14)
text(1.52,1.25,['{\itf_s}=',num2str(fs/n),'Hz'],'FontSize',14)
grid on
w1=fft(z,512); % RACUNANJE SPEKTRA SIGNALA
L1=length(z);
Pyy1 = abs(fft(z))/L1;
Pyy1=[Pyy1,Pyy1];
subplot(224)
plot([1:length(Pyy1)]*fs/(length(Pyy1)),Pyy1)
title('Spektar signala nakon decimacije','FontSize',14)
xlabel(' frekvencija [Hz]')
ylabel('Magnituda')
grid on
set(figure(2),'Position',[70 60 900 600],'Name','DECIMACIJA')
prilog
71