
% Inseriamo una fdt:
G=tf([2 4 1 5],[1 3 10 2 0]);
% Tracciamento completamente automatico dei diagrammi di Bode:
% Viene aperta una finestra grafica interattiva, ove vengono
% tracciati i diagrammi e tramite la quale si
% possono rilevare dati utili, cliccando sul grafico
bode(G);
grid on;
% alternativamente, i diagrammi di modulo e fase si possono plottare come segue:
% definiamo l'asse logaritmico delle omega:
w=logspace(-2,2,150); % 150 punti, spaziati logaritmicamente tra 10^-2 e 10^2
[Mod,Fase]=bode(G,w); % restituisce moduli in unità lineari, e fasi in gradi,
                      % calcolati alle varie pulsazioni contenute nel vettore w.
Mod_dB = 20*log10(Mod); % converte i moduli in decibel
Mod_dB = squeeze(Mod_dB); % fa diventare Mod_dB un vettore colonna
Fase = squeeze(Fase);
% plottiamo il diagramma di modulo:
semilogx(w,Mod_dB);
grid on;
xlabel('\omega [rad/s]')
ylabel('Modulo [dB]')
% plottiamo il diagramma di fase:
semilogx(w,Fase);
grid on;
xlabel('\omega [rad/s]')
ylabel('Fase [gradi]')