Il codice che ti ho scritto alla fine del messaggio precedente non mi convince, perché non posso prevedere il comportamento del luogo quando ci sono zeri al finito e dai valori di k sparisce Inf (magari funziona lo stesso). Per sicurezza, ti indico una variante che si basa sul fatto che, se Matlab passa direttamente da un valore finito all’infinito, probabilmente ritiene che non ci siano variazioni significative in quel tratto:
`
ni = 1; % numero di nodi intermedi tra due adiacenti
[r, k] = rlocus(FDT)
[r, k] = rlocus(FDT, [interp1(k(1:end-1), 1:1/(ni+1):length(k)-1) Inf]);
`Aumentando ni, puoi aumentare il numero di nodi intermedi.
I punti che ottieni in questo modo appartengono sicuramente al luogo delle radici, mentre, interpolando i valori di r, non puoi esserne certo e, quando si analizza la stabilità, bisogna prestare attenzione.
Ciao
P.S.: Ci puoi dire direttamente la funzione di trasferimento?