clear all;
fp = 1000;  %czestotliwosc probkowania
fs1 = 10;  %czestotliwosc sygnalu #1
fs2 = 100;   %czestotliwosc sygnalu #2
A = 1;      %amplituda sygnalu
N = 500;   %liczba probek
phi = pi/3; %przesuniecie fazowe
t = 0:1/fp:(N-1)/fp;%wektor czasu
% otrzymujemy wektory _x_ i _x1_
x = A*sin(2*pi*fs1*t + phi);%nasz sygnal
x1= A*sin(2*pi*fs2*t);%nasz sygnal
f = 0: fp/N: fp - fp/N;

%%transformata Fouriera
X = fft(x);
X1 = fft(x1);
X2 = fft(x1+x);

figure(1)
y=x+x1; % suma sygnałów
y1=y+0.3*randn(size(y));
plot(t,y,t,y1)


%%%wykres widma amplitudowego w dB
figure(2)
plot(f,abs(20*log10(abs(X2))))
xlabel 'czestotliwosc [Hz]'
title 'widmo amplitudowe sygnalu'

%%%
faza = unwrap(angle(X));
figure(3)
plot(f,faza*180/pi)
xlabel ''
title 'wykres widma fazowego'