% nacrtati kosinusni signal x_c(t) = cos(2*pi*60 * t) na [-0.05 0.05]
% i izvrsiti odabiranje datog signala sa raznim periodama odabiranja T
% potom nacrtati signal x_c1 = cos(2*pi*40 * t) na [-0.05 0.05]

clear all;
close all;
f = 60;  % Hz
tmin = -0.05;
tmax = 0.05;
t = linspace(tmin, tmax, 400);
x_c = cos(2*pi*f * t);
plot(t,x_c,'k'); % crtamo kosinusni signal frekvencije 60 herca
xlabel('t (seconds)');

T = 1/100; %ovde probati sa raznim periodama obabiranja 1/130, 1/150
nmin = ceil(tmin / T);
nmax = floor(tmax / T);
n = nmin:nmax;
x1 = cos(2*pi*f * n*T); % ovo su odbirci datog kosinusnog signala iznad
hold on
plot(n*T,x1,'O')% crtamo taj semplovani signal kruzicima

% nacrtajmo jos jedan kosinusni signal (od 40 herca) i zakljucicemo da
% malopredjasnji odbirci odgovaraju i ovom kosinusnom signalu
%(ne moze se na osnovu odbiraka zakljuciti o kom polaznom signalu je rec 
%jer smo nasli jos jedan signal manje frekvencije kome odgovaraju odbirci)
% Dakle, nije moguce izvrsiti rekonstukciju signala jer nije zadovoljen
% Nikvistov uslov tj. mora biti bar T=1/(2*60) da bi bila zadovoljena teorema
% o odabiranju (i tada odbirci nece pripadati ovom drugom signalu od 40
% herca)

x_c1 = cos(2*pi*40 * t); 
nmin = ceil(tmin / T);
nmax = floor(tmax / T);
n = nmin:nmax;
x1 = cos(2*pi*f * n*T);
plot(t, x_c1,'r');
legend('prvi signal', 'odbirci', 'drugi signal');
