close all;
clc;
%Ts = periodo di campionamento
%t = vettore dei tempi
%s = segnale discretizzato

s_real = real(s);
s_imm = imag(s);


%stampo s reale
figure;
plot(t, s_real);
grid on;
title('Parte Reale di s(kTs)');
ylabel("S Reale");
xlabel("t");

%stampo s immaginario
figure;
plot(t, s_imm);
grid on;
title('Parte Immaginaria di s(kTs)');
ylabel("S Immaginaria");
xlabel("t");


t0 = 1;
t1 = 3;
i1 = ((t0 + 10)/Ts)+1;
i2 = ((t1 + 10)/Ts)+1;
operando = (abs(s)).^2;

energia = Ts * sum(operando); %energia su tutto l'intervallo t
energia_intervallo = Ts * sum(operando(i1:i2)); 

%operando va da -10 a 10 secondi
%ho bisogno dell'intervallo tra 1 e 3 secondi