%   Ardito Luca n matr 114093
%   Corso di Elettronica applicata: esercitazione n°1
%   Risposta di un passabasso al quale viene inserita un'onda quadra di
%   ampiezza [a,-a] e DC 50%
%
%
%   Il programma restituisce alla pressione di un tasto l'onda quadra
%   approssimata in ingresso e l'onda in uscita al filtro.
%
%
%   Parametri:      freq=frequenza fondamentale
%                   arm=numero di armoniche considerate               
%                   amp=ampiezza picco-picco del segnale in ingresso
%                   r=valore resistenza
%                   c=capacità del condensatore
%                   t=numero di intervalli temporali in cui viene diviso il
%                   tempo di osservazione
%



function uscita=homework1(freq,arm,amp,r,c,t)

    n=[1:2:arm];                        % Indice (dispari) delle armoniche considerate
    bn=4./(n.*pi);                      % Coefficienti serie di fourier
    om=2.*freq.*n.*pi;                  % Frequenze delle singole armoniche
    T=1/freq;                           % Periodo onda quadra in ingresso
    tempo=[0:T/t:T];                    % Istanti di osservazione
    kv=(1/(r*c))./(j.*om+(1/(r*c)));    % Funzione di trasferimento del filtro
    modulo=abs(kv);                     % Modulo kv
    fase=angle(kv);                     % Fase kv
    for indice=1:t+1                    % Ricostruzione del segnale in ingresso e di quello in uscita
        singr(indice)=sum((amp/2).*bn.*(amp/2)*sin(om*tempo(indice)));
        uscita(indice)=sum((amp/2).*bn.*modulo.*sin(om*tempo(indice)+fase));
    end
    plot(singr);                        % Visualizzazione segnale in ingresso
    grid on;                            % Attivazione griglia
    'Per visualizzare l''onda in uscita premere un tasto'
    pause;                              % Attesa pressione di un tasto per stampare la funzione di uscita
    plot(uscita);                       % Visualizzazione segnale in uscita
    grid on;