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

� Primjer2.m
clear all
close all
f1=1; % FREKVENCIJA OSNOVNE KOMPONENTE
f2=3; % FREKVENCIJA VISEG HARMONIKA
fs=15; % FREKVENCIJA ODMJERAVANJA
x=[0:1/fs:10]; % VRIJEME
y=sin(2*pi*f1*x)+sin(2*pi*f2*x) ; %SIGNAL
Pin = abs(fft(y))/length(y); % RACUNANJE SPEKTRA SIGNALA
fin =fs*(0:255)/512; % NORMALIZOVANJE FREKVENCIJE
x1=[0:1/fs:10];
y1=sin(2*pi*f1*x1)+sin(2*pi*f2*x1) ; %ODMJERENI SIGNAL
subplot(221)
stem(x1(1:2*fs),y1(1:2*fs),'.') % IZGLED SIGNALA
title('Odmjerci signala
{\itx(t)}=sin(2\pi{\itf_1t})+sin(2\pi\itf_2t)','FontSize',14)
xlabel(' vrijeme [s]')
ylabel('Amplituda signala')
text(1.55,1.55,['{\itf_s}=',num2str(fs),'Hz'],'FontSize',14)
text(1.55,1.05,['{\itf_1}=',num2str(f1),'Hz'],'FontSize',14)
text(1.55,0.5,['{\itf_2}=',num2str(f2),'Hz'],'FontSize',14)
grid on
grid on
w=fft(y1,512); % RACUNANJE SPEKTRA SIGNALA
Pyy = abs(w)/length(y1);
fyy =fs*(0:511)/512; % NORMALIZOVANJE FREKVENCIJE
subplot(222)
plot(fyy,Pyy(1:512)) % IZGLED SPEKTRA DIGITALNOG SIGNALA
title('Spektar signala nakon odmjeravanja','FontSize',14)
xlabel(' frekvencija [Hz]')
ylabel('Magnituda')
grid on
set(figure(1),'Position',[70 60 900 600],'Name','INTERPOLACIJA KAO
FILTRIRANJE')
% *****Interpolacija nultih odmjeraka*****
% Sa 'n' mijenjamo faktor konverzije,
% n=1=> konverzija navise 1:2
% n=2=> konverzija 1:3
% povecavamo frekvenciju odmjeravanja (n+1)puta
n=1;
k=length(y1);
y2(1:(n+1)*(k-1))=0;
y2(1)=y1(1);
for i=2:k
y2((n+1)*(i-1)+1)=y1(i);
end
subplot(223)
stem(x1(1:4*fs)/(n+1),y2(1:4*fs),'.')
title('Signal nakon ubacivanja nultih odmjeraka','FontSize',14)
xlabel(' vrijeme [s]')
prilog
72