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

subplot(221)
stem(x(1:2*fs),y(1:2*fs),'.') % IZGLED SIGNALA
title('Odmjerci signala','FontSize',14)
xlabel(' vrijeme [s]')
ylabel('Amplituda signala')
text(1.52,1.6,['{\itf_s}=',num2str(fs),'Hz'],'FontSize',14)
text(1.52,1.15,['{\itf_1}=',num2str(f1),'Hz'],'FontSize',14)
text(1.52,0.5,['{\itf_2}=',num2str(f2),'Hz'],'FontSize',14)
grid on
%w=fft(y,512); % RACUNANJE SPEKTRA SIGNALA
Pyy = abs(fft(y,512))/length(y);%w.*conj(w);
fyy =fs*(0:511)/512; % NORMALIZOVANJE FREKVENCIJE
subplot(222)
plot(fyy,Pyy(1:512)) % IZGLED SPEKTRA
title('Spektar signala nakon odmjeravanja','FontSize',14)
xlabel(' frekvencija [Hz]')
ylabel('Magnituda')
grid on
p=3;
q=2;
d=gcd(p,q); % Ako p/q nije prost razlomak, skratimo ga
p=p/d;
q=q/d;
[z,b]=resample(y,p,q);
fs1=(p/q)*fs;
x1=[0:1/fs1:10];
subplot(223)
stem(x1(1:2*fs1),z(1:2*fs1),'.')
title(['Signal nakon resamplinga pri faktoru
',int2str(q),':',int2str(p)],'FontSize',14)
xlabel(' vrijeme [s]')
ylabel('Amplituda signala')
text(1.52,1.6,['{\itf_s}=',num2str(fs*p/q),'Hz'],'FontSize',14)
grid on
wz=fft(z,512); % RACUNANJE SPEKTRA SIGNALA
Pyy1=abs(fft(z,512))/length(z);
fyy1 =fs*(p/q)*(0:511+(ceil(q/p)-1)*512)/512; % NORMALIZOVANJE
FREKVENCIJE + periodicnost spektra
for i=0:ceil(q/p)-1
Pyy2(512*i+1:512*(i+1))=Pyy1;
end
a=find(fyy1