%   Ardito Luca n matr 114093
%   Corso di Elettronica applicata: esercitazione n°2
%   Risposta di un passabasso al quale viene inserita un'onda quadra a frequenza 15Hz
%   di ampiezza [a,-a] e DC 50%
%
%
%   Il programma restituisce l'onda quadra approssimata in ingresso e, dopo 
%   la pressione di un tasto, l'onda in uscita al filtro passabasso.
%
%
%   Parametri:      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 homework2(arm,a,r,c,t)
    freq=15;                                                            %Viene settata la frequenza a 15Hz
    n=(1:2:arm);                                                        %Armoniche che approssimano l'onda in ingresso
    dt=(1/freq)/t;                                                      %Intorno della precisione inserita
    T=((-(1/freq)/2)-dt:dt:((1/freq)/2)+dt);                            %Definizione estremi della funzione
    oq=-(-(1/freq)/2<=T & 0>=T) + (0<=T & (1/freq)/2>=T);               %Funzione che genera un'onda quadra
    oq=oq*(a/2);                                                        %Onda quadra moltiplicata per l'ampiezza di picco
    hndl=plot(T,oq);                                                    %Visualizza l'onda quadra in ingresso
    set(hndl,'color','black');                                          %Setta il colore nero 
    'Per visualizzare l''onda in uscita premere un tasto'
    pause;                                                              %Attende la pressione del tasto invio
    hold on;                                                      
    bn=4./(n*pi);                                                       %Approssimazione dell'onda quadra
    Kv=(1/(r*c))./(j*(2*pi*freq.*n)+(1/(r*c)));                         %Funzione di trasferimento del passabasso
    modulo=abs(Kv).*bn;                                                 %Ampiezza dell'uscita
    fase=angle(Kv);                                                     %Fase dell'uscita
    for i=1:t+3
        uscita(i)=sum(modulo*(a/2).*sin(((2*pi*freq.*n)*T(i))+fase));   %Funzione d'uscita
    end
    plot(T,uscita);                                                     %Visualizza l'uscita del filtro passabasso
    hold off;
    'Per chiudere la finestra premere il tasto INVIO'
    pause;
    close all;