Problemi di vario tipo

[absentio]

Ciao, anzitutto grazie mille per il grandissimo lavoro svolto. Sto cercando di utilizzare il vostro programma per prendere dimestichezza con i diagrammi di Bode ma sto avendo dei problemi.
Ad esempio con una funzione di trasferimento pari a (10/s)*(1/(s+1)) sulla vostra funzione mag(sys,w0) alla pulsazione di 5 rad/s ho un guadagno di -14 dB, su Matlab di -8dB.
Tracciando il diagramma di bode con bodex il risultato è coerente con matlab.
Inoltre quando utilizzo la funzione per il calcolo per il margine di fase è come se la calcolatrice non svolgesse le funzioni abs e sign e quindi il risultato non viene calcolato.
Sapete come risolvere il problema?
Ho già provato a ripetere l'installazione anche su di un emulatore ma ottengo sempre lo stesso problema.

Grazie, Antonio

[loretoparisi]

Ciao Antonio! Grazie, molto interessante la tua verifica. In teoria sia l'ambiente emulato che sulla calcolatrice portatile reale il risultato dovrebbe essere lo stesso, in quanto l'emulatore usa la ROM della calcolatrice, per cui a meno di errori, avrai la stessa approssimazione.
Detto ciò sembrerebbe quindi una differenza tra il valore computato dal plot di BodeX e quello calcolato che direi che vedi nella schermata di output del programma.
Potresti allegarmi la schermata di entrambi (BodeX e valore computato), grazie.
Per completezza (e anche per ricordalo in quanto l'ho sviluppata molti anni fa :) aggiungo la screenshot della funzione in questione

Da The CST Reference Guide, First Edtion October 2005.

Nel frattempo che indago, ti do un suggerimento, di fare una verifica con gli esempi presenti in The CST User Guide, Fifth Edtion October 2005. Qui troverai degli esempio di risoluzione di numerosi problemi tipici del Controllo Dei Sistemi, dei Controlli Automatici e di Teoria Dei Sistemi.

Proverò a ripetere il tuo caso specifico per capirne la causa.

[absentio]

Dalla calcolatrice:
Svolgendo i calcoli a mano in pratica la calcolatrice trova il rapporto(non logaritmico) come sqrt(26)/26 invece di ottenere sqrt(26)/13, cio è dovuto al valore errati ritornato dalle funzione 6 del menù data (ritorna 1/sqrt(w^2+1))
Anche la funzione phase(SYS,w0) ritorna un valore non corretto della fase alla data pulsazione.
La funzione margin(sys) non funziona esattamente come nell'emulatore(screenshot allegati)
Dall'emulatore:
Le funzioni 6 e 7 del menù data danno risultati diversi dalle stesse funzioni eseguite nella calcolatrice.
La funzione margin invece no, riporta le stesse identiche problematiche della calcolatrice.


La funzione mag funziona correttamente e ritorna un valore esatto.
La funzione phase si blocca così e non restituisce alcun valore:

[loretoparisi]

@absentio grande! Potresti verificare anche se eseguendo alcuni degli esempi nel tutorial hai esattamente gli stessi valori riportati sulla guida? Se non è così, allora il problema può essere non la funzione in se, ma altro. E' la Titanium giusto?

[absentio]

Si si tratta di una Titanium versione hardware 3. Prima di farti impazzire ho ripetuto l'installazione sulla calcolatrice da capo e adesso ottengo gli stessi risultati dell'emulatore. Ho seguito il primo degli esercizi della tua guida e la mia calcolatrice andava di passo con gli screenshots allegati al documento. Ah,dimenticavo,utilizzo la versione di CST presente all'interno della cartella src della root di questa repository
EDIT: la parte a valori discreti non mi va perchè non trovo la funzione sampler

[absentio]

Ciao, volevo aggiungere che anche il programma nyquist() da qualche problema e si pianta con alcune fdt, ad esempio 10/[(s+1)(s+2)(s+3)]. Lasciato lavorare per diversi minuti alla fine da un Error: break

[absentio]

@loretoparisi penso il problema sia dovuto alla funzione sign(0) che non viene risolta dalla calcolatrice.

[loretoparisi]

@absentio molto strano, è possibile che sia cambiato qualcosa nel sw sulla Titanium? Devo fare una verifica più accurata.

[absentio]

@loretoparisi ti lascio quello che ho trovato nel manuale della Titanium. La mia titanium ha OS 3.10

[absentio]

@loretoparisi riesci a dargli un'occhiata a breve? Avrei la prova di controlli automatici giorno 19 :(