% czas
t = 0:0.0001:10;

% sygnal
s1 = 2*sin(2*pi*2*t);
%s2 = 3*cos(2*pi*7*t);
%s3 = 4*sin(2*pi*2*t);

% suma sygnalow
signal = s1;

figure(1);
subplot(2,2,[1,2]);
plot(t,signal);
ylabel('Amplituda');
xlabel('Czas');
suptitle('Sygnał');


centerOfMass = zeros(2);
i= 1;

max_f_zwiniecia = 3;
for fzwiniecia=0.1:0.001:max_f_zwiniecia

    funczwinieta = exp(-2*pi*j*t*fzwiniecia).*signal;

    subplot(2,2,3);
    plot(real(funczwinieta),imag(funczwinieta));
    ylabel('Im');
    xlabel('Re');
    title('Sygnał zawinięty');

    avg = mean(funczwinieta);
    centerOfMass(i,:) = [fzwiniecia,abs(avg)];
    subplot(2,2,4);
    plot(centerOfMass(:,1),centerOfMass(:,2))
    ylabel('Amplituda');
    xlabel('Częstotliwość');
    title('prawie FFT');
    axis([0 max_f_zwiniecia  -inf inf])

    pause(0.001);
    i=i+1;
end