clc
clear all
Fp=40;
T=2;
dt=1/Fp;
N=Fp*T;
a=4;
b=33;

t=linspace(0,T,N);

s=sin(2*pi*10*t);
g=sin(2*pi*15*t);
h=s.*a+g.*b;

%transformata
S=fft(s,N);
G=fft(g,N);
H=fft(h,N);

H1=S.*a+G.*b;

%widmo
S=abs(S)/N;
G=abs(G)/N;
H=abs(H)/N;
H1=abs(H1)/N;

%skalowanie częstotliwości
fHz=linspace(0,Fp/2,N/2);

%obcięcie drugiej połowy
S(N/2+1:N)=[];
G(N/2+1:N)=[];
H(N/2+1:N)=[];
H1(N/2+1:N)=[];

%przeniesienie mocy
S(2:end)=S(2:end).*2;
G(2:end)=S(2:end).*2;


figure
subplot(2,1,1)
plot(fHz,H)
subplot(2,1,2)
plot(fHz,H1)