Grupa : 06.01.2013 r.
Sprawozdanie z laboratorium nr 6
Temat : "Interpolacja i decymacja"
Nasze zadanie polegało na wygenerowaniu sygnałów sinusoidalnego i prostokątnego oraz wygenerowanie jego wykresów po interpolacji oraz decymacji. Dodatkowo mieliśmy wygenerować widmo tych sygnałów które również poddawaliśmy interpolacji i decymacji. Wyniki tych operacji zamieściliśmy poniżej :
wprowadzone komendy :
function [t,x]=sinus(A,f,fp,T,faza)
t=zeros(1,T*fp);
x=t;
for i=0:1/fp:T
t(k)=i;
x(k)=A*sin(2*pi*f*i+(faza*180/pi));
k=k+1;
end;
funkcja interp :
xi=0:0.005:0.05
yi=interp1(t,x,xi,'spline');
decymacja :
xi=rzad
yi=decimate([t,x],xi);
Wygenerowany wyjściowy sygnał sinusoidalny:
Sygnał sinusoidalny po interpolacji:
Sygnał sinusoidalny po decymacji:
Widmo sygnału sinusoidalnego wyjściowego :
Widmo sygnału sinusoidalnego po decymacji:
Widmo sygnału sinusoidalnego po interpolacji:
Funkcje dla sygnału prostokątnego:
function [t,x]=prostokat(A,f,fp,T,faza,w)
t=zeros(1,round(T*fp+1));
x=t;
proz=floor (w*fp/f);
prz=floor ((1-w)*fp/f);
k=1;
i=0;
while k <= (T*fp)
for j=1:proz
t(k)=i;
x(k)=A;
i=i+1/fp;
k=k+1;
end;
for j=1:prz
t(k)=i;
x(k)=0;
i=i+1/fp;
k=k+1;
end;
end;
Wygenerowany wyjściowy sygnał prostokątny:
Sygnał prostokątny po interpolacji:
Sygnał prostokątny po decymacji:
Widmo sygnału prostokątnego wyjściowego :
Widmo sygnału prostokątnego po decymacji:
Widmo sygnału prostokątnego interpolacji:
Wnioski:
Na tym laboratorium poznaliśmy funkcję interpolacji i decymacji. Jak można zauważyć pierwsza z tych funkcji upraszcza wykres, tzn. pobiera wartość sygnału w kilku określonych punktach zwanych węzłami. Decymacja natomiast polega na zachowaniu co M-tej próbki sygnału i odrzuceniu pozostałych co sprawia że sygnał jest nieco zniekształcony ( ma M-próbek mniej od wyjściowego). Obie funkcje są wykorzystywane do przekształcenia sygnału do postaci cyfrowej.