Diagrammi asintotici di bode con matlab

[enpires1]

salve a tutti! qualcuno mi sa dire come posso fare per tracciare i diagrammi asintotici di bode utilizzando matlab?
Ho provato a cercare nell'help della funzione "bode" ma non trovo nulla

[Umby2]

Girovagando per il web, ho trovato questo codice:

function [mod,fase,omega]=asbode(num,den,W_axes,M_axes,F_axes,B,R,X,Y);
% ASBODE : traccia il diagramma asintotico di Bode 
% ******************************************************************************
%                 ## SYNTAX ##
%
%  ASBODE(NUM,DEN)
%     Traccia il diagramma di Bode di una funzione di trasferimento
%
%                   NUM     b_m*s^m + b_{m-1}*s^{m-1} + ... + b_1*s + b_0
%          W(s) =  ----- = -----------------------------------------------
%                   DEN     a_n*s^n + a_{n-1}*s^{n-1} + ... + a_1*s + a_0
%
%     e scrive sullo schermo i termini che compongono la fattorizzazione di
%     Bode della W(s).
%
%     Il vettore 
%          NUM = [ b_m  b_{m-1}  ...  b_1  b_0 ]  
%     contiene i coefficienti del polinomio al numeratore mentre il vettore
%          DEN = [ a_n  a_{n-1}  ...  a_1  a_0 ]  
%     contiene i coefficienti del polinomio al denominatore della W(s).
%
%     I diagrammi asintotici dei singoli termini sono tracciati in diversi colori,
%     il diagramma asintotico complessivo e' tracciato con una linea nera spessa, 
%     il diagramma reale complessivo e' tracciato con una linea nera tratteggiata. 
%     Le scale e il vettore delle frequenze sono scelti automaticamente.
%

il link è questo:

http://www.diee.unica.it/~giua/ASD/materiale.html

[enpires1]

Si quel codice l'ho visto anche io, ma non è molto affidabile :\ ad esempio se provi a graficare con questo script il diagramma delle fasi di $(1+s)^{-1}$ utilizza un solo asintoto, mentre in realtà ce ne dovrebbero essere 3