clear,clc
f0 = 10; %czestotliwosc 10Hz
f1 = 20; %czestotliwosc 20Hz
fs = 100; %czestotliwosc probkowania 100Hz
w0 = 2*pi*f0; %pulsacja sygnalu 10Hz
w1 = 2*pi*f1; %pulsacja sygnalu 20Hz
Ts = 1/fs; %odstep czasu miedzy probkami (Ts)
t = 0:Ts:(1-Ts); %przedzial poboru probek
y0 = sin(w0*t); %generowanie funkcji sinus dla danej czestotliwosci (f0 = 10Hz)
y1 = sin(w1*t); %generowanie funkcji sinus dla danej czestotliwosci (f1 = 20Hz)
F0 = fft(y0); %wyznaczenie DFT dla y0 za pomoca Fast Fourier Transform
F1 = fft(y1); %wyznaczenie DFT dla y1 za pomoca Fast Fourier Transform
F10 = F0+F1; %wyznaczenie sumy transformat
M10 = abs(F10); %modul sumy transformat
F1_0 = fft(y0+y1); %wyznaczenie transformaty sumy
M1_0 = abs(F1_0); %wyznaczenie modulu transformaty sumy
M0 = abs(F0); %modul transformaty dla y0
M1 = abs(F1); %modul transformaty dla y1
%F0(M0<1e-6) = 0; %opcjonalne pominiecie wartosci bardzo, bardzo malych
%F1(M1<1e-6) = 0; %opcjonalne pominiecie wartosci bardzo, bardzo malych
P0 = unwrap(angle(F0)); %faza DFT dla y0
P1 = unwrap(angle(F1)); %faza DFT dla y1

figure
stem(t,M10) %stem rysuje prazki, zamiast ciaglosci jak plot
ylabel('Moduł z DFT');
xlabel('Czas t [s] - 100n/s');
grid on;
figure
stem(t,M1_0)
ylabel('Moduł z DFT');
xlabel('Czas t [s] - 100n/s');
grid on;

%kod alternatywny probkowania sinusa
%N = 100
%n = 0:1:N-1
%y1 = sin(2*pi*f/fs*n);