clc,close all,clear,beep off;
%% Przyk³adowe generowanie przebiegu sinusoidalnego
% Czestotliwoœci próbkowania
Fs = 1600;
% Czas próbkowania
Ts = 1/Fs;
% Liczba próbek
L = 64;
% Oœ czasu - wartoœci dyskretne
t_d = (0:L-1)*Ts;
% Oœ czasu - wartoœci ci¹g³e
t_c = linspace(0,t_d(end),100);

% Parametry sygna³y - Amplituda
A = 10;
% Parametry sygna³y - Czêstotliwoœci
f = 50;
% Parametry sygna³y - Przesuniêcie fazowe
Phi = 0;

% Sygna³ - wartoœci osi Y
y_c = A*sin(2*pi*f*t_c + Phi);
y_d = A*sin(2*pi*f*t_d + Phi);

% Wykres
figure
with_standard_size
set(gcf,'DefaultLineLineWidth',1);
set(gcf,'DefaultAxesFontSize',14);
plot(t_c,y_c,'-r')
hold on;
stem(t_d,y_d,'s','LineWidth',2)
hold off
grid on
xlabel('Czas[s]')
ylabel('Amplituda')
legend('- c -','- d -')

% DFT
Y = fft(y_d);
Y = abs(Y)/L;
Y(2:end-1) = 2*Y(2:end-1);
f = Fs*(0:(L-1))/L;

% Wykres
figure
with_standard_size
set(gcf,'DefaultLineLineWidth',1);
set(gcf,'DefaultAxesFontSize',14);
stem(f(1:L/2+1),Y(1:L/2+1),'s','LineWidth',2)
grid on
xlabel('f[Hz]')
ylabel('Amplituda')
legend('- d -')