close all
clear all
clc

Tp = 5; % periodo
d = 0.4; % duty cycle
A = 2; % ampiezza onda
om0 = 2*pi/Tp; % omega zero
N = 50; % numero di armoniche considerate nel troncamento
T = 0.01; % passo di campionamento

% segnale originario x(t)
[t, x] = ondaquadra(A,Tp,d,T);

% stampa segnale
figure
plot(t,x,'r-');
xlabel('t')
ylabel('x(t)')
title(['Serie di Fourier troncata, N=', num2str(N)])
grid

%% %%%%%%%%%%%%%%%%%%%%%% serie troncata con coefficienti ak

% 1) calcolare i coefficienti della serie di Fourier 
k = -N:N;
ak = A*d*sinc(k*d);

% 2) calcolare la serie troncata 
xn = zeros(size(t));
...
    
% 3) plottare il segnale
hold on
plot(t,xn);
hold off
legend('onda quadra','serie troncata')


%% %%%%%%%%%%%%%%%%%%%%%% serie troncata con coefficienti bk

% 1) calcolare numericamente i coefficienti della serie di Fourier 
bk = zeros(size(k));
...
    
% 2) calcolare la nuova serie troncata 
xnb = zeros(size(t));
...
    
% 3) plottare il segnale
hold on
plot(t,xnb);
hold off
legend('onda quadra','serie troncata ak','serie troncata bk')

% 4) illustrare la differenza tra i coefficienti ak e bk
...