%Ardito Luca n matr 114093
%Corso di Elettronica applicata: esercitazione n°1
%Il programma stampa il diagramma di bode di una data Kv
% parametri
% costante
% num zeri
% array zeri
% num poli
% array poli

function homework4(k,nzeri,zeri,npoli,poli)
    om=logspace(-4,4);                      %Asse x
    A=20*log10(abs(k));                     %Valore in db della costante
    z=0;                                    %Inizializzazione variabile    
    fz=0;                                   %Inizializzazione variabile
    p=0;                                    %Inizializzazione variabile
    fp=0;                                   %Inizializzazione variabile
    for i=1:nzeri                           %Per tutti gli zeri calcola modulo e fase
        z=z+20*log10(om.^2+zeri(i)^2);      %Modulo degli zeri
        fz=fz+atan(om/zeri(i));             %Fase degli zeri
    end                                     %Fine ciclo for
    for i=1:npoli                           %Per tutti i poli calcola modulo e fase
        p=p+20*log10(om.^2+poli(i)^2);      %Modulo dei poli
        fp=fp+atan(om/poli(i));             %Fase dei poli
    end                                     %Fine ciclo for
    out=A+z-p;                              %Modulo Kv
    outf=fz-fp;                             %Fase Kv
    subplot(2,1,1);                         %Divisione del grafico in 2 e utilizzo della parte superiore
    semilogx(om,out);                         %Stampa a video della funzione risultante (20*log(|kv|)
    subplot(2,1,2);                         %Divisione del grafico in 2 e utilizzo della parte inferiore
    semilogx(om,outf,'r');