clear all;
clc
% tempo discreto

n = -1:0.001:1; % intervalo das amostras
L = size(n,2); % tamanho do vetor

FS = 50000; % frequencia de amostragem

% NFFT = 2^nextpow2(L); % calculando o numero de coeficientes
NFFT = 256;
f = FS*linspace(-1,1,NFFT); % calculando a faixa de frequencias

% fases do sistema elétrico

r = 240*sin(2*pi*60*n + pi/3);
s = 240*sin(2*pi*60*n + pi);
t = 240*sin(2*pi*60*n - pi/3);
ground = r + s + t;

% Mostrando na tela

subplot(3,1,1); plot(n(1:20),r(1:20));
subplot(3,1,2); plot(n(1:20),s(1:20),'r');
subplot(3,1,3); plot(n(1:20),t(1:20),'black');

% stem(r(1:19),'DisplayName','r');figure(gcf)

figure;
subplot(3,1,1); stem(n(1:20),r(1:20));
subplot(3,1,2); stem(n(1:20),s(1:20),'r');
subplot(3,1,3); stem(n(1:20),t(1:20),'black');

% somando frequencias no matlab 

figure;
n = -1:0.00001:1; % intervalo das amostras
L = size(n, 2);
Amp = 240 - (0:20:240);
freq = 1000:1000:13000;
theta = 0:pi/6:2*pi;

FS = 100000; % frequencia de amostragem

% NFFT = 2^nextpow2(L); % calculando o numero de coeficientes
NFFT = 256;
f = FS*linspace(-1,1,NFFT); % calculando a faixa de frequencias


% somando várias frequêcias, amplitudes e fases

soma = zeros(1, L);

for i = 1:13
soma = soma +  Amp(i)*cos(2*pi*freq(i)*n + theta(i));
end


% plotando os 400 primeiros elementos
plot(soma(1:1000))

% magnitude (par)
figure;
plot(f, abs(fftshift(fft(soma,NFFT)/L)));

% fase (Ímpar)
figure;
plot(f, angle(fftshift(fft(soma,NFFT)/L)));