Ti 5.11.2010

Za 2 tydzień koło !!! z 3 pierwszych lekcji

3. Próbkowanie

-Przetwornik cyfrowo-analogowy

-Widmo częstotliwości

-Twierdzenie Kotielnikowa-Shannona Fs > 2 Fg

-Kwantyzacja

fs=2000; % częstotliwość próbkowania

t=0:1/fs:1; % zdyskretyzowany czas

a=1; % amplituda sinusoidy

f1=128; % częstotliwość sygnału

x= a*sin(2*pi*t*f1);

% wizualizacja

plot(t(10:128),x(10:128));

fs=2000; % częstotliwość próbkowania

t=0:1/fs:1; % zdyskretyzowany czas

a=1; % amplituda sinusoidy

f1=128; % częstotliwość sygnału

x= a*sin(2*pi*t*f1);

% wizualizacja

wiz=1:128;

plot(t(wiz),x(wiz));

fs=200; % częstotliwość próbkowania

t=0:1/fs:1; % zdyskretyzowany czas

a=1; % amplituda sinusoidy

f1=10; % częstotliwość sygnału

x= a*sin(2*pi*t*f1);

% wizualizacja

wiz=1:128;

plot(t(wiz),x(wiz));

dodanie szumów

fs=1024% częstotliwość próbkowania

t=0:1/fs:1; % zdyskretyzowany czas

a=1; % amplituda sinusoidy

f1=20; % częstotliwość sygnału

x= a*sin(2*pi*t*f1)+ 0.1*randn(size(t)); % randn <-- generator liczb losowych od 0 do 1

% wizualizacja

wiz=1:128;

plot(t(wiz),x(wiz));

fs=1024% częstotliwość próbkowania

t=0:1/fs:1; % zdyskretyzowany czas

a=1; % amplituda sinusoidy

f1=20; % częstotliwość sygnału

x= a*sin(2*pi*t*f1)+ 0.1*randn(size(t));

okno=4;

xfil=fikter(ones(1,okno)/okno,1,x);

idx=1:256;

figure(1);

plot(t(idx),x(idx));

figure(2);

plot(t(idx),xfil(idx));