ElectroYou

Sono sempre rimasto abbastanza intrigato dai problemi di teoria delle reti elettriche. Spesso sono intimamente legati a problemi non banali di natura matematica, ed è una matematica applicata molto interessante.

Nei circuiti a parametri concentrati con un numero finito di elementi, le funzioni di rete sono del tipo:

dove e sono dei polinomi nella variabile complessa . Una cosa importante da sapere è se gli zeri dei due polinomi hanno parte reale negativa oppure no, dato che queste informazioni sono di solito legate alla stabilità del sistema, oppure al fatto che la funzione è a fase minima oppure no. In campo reale le cose sono semplici da verificare utilizzando la regola di Cartesio, ma in campo complesso la situazione è molto più spinosa.

Diventa quindi importante verificare questa proprietà per un generico polinomio, se possibile senza calcolare esplicitamente gli zeri del polinomio, dato che questo potrebbe essere un problema non banale per un grado elevato del polinomio da studiare. Esiste un test algebrico chiamato prova di Hurwitz che permette di determinare se questa proprietà è verificata oppure no, senza calcolare esplicitamente gli zeri. Un polinomio è hurwitziano se tutti i suoi zeri hanno parte reale strettamente negativa, se ce n'é qualcuno immaginario è hurwitizano in senso lato.

Qualche anno fa, avevo messo a punto un piccolo script Maxima che permette di verificare se un polinomio è hurwitziano oppure no.
Ne ho già parlato sul mio sito:
http://davbucci.chez-alice.fr/index.php ... aliano#fun
ma mi sono detto che sarebbe stato carino proporlo nuovamente qui, per stimolare una discussione e perché magari potrebbe essere utile a qualcuno.



Una volta salvato lo script in un file, per esempio hurwitz.max, lo si può utilizzare all'interno di Maxima:



La funzione hurwitz(N,p) richiede il polinomio N nella variabile p ed effettua il test, che è basato su una serie di divisioni successive fra la parte pari e la parte dispari del poliniomio di partenza. Il programma restituisce una lista contenente i risultati delle divisioni.
Una cosa che mi piacerebbe fare una volta o l'altra è preparare altri script di questo tipo che applichino alcune procedure di sintesi empirica per le reti lineari passive. Che ne pensate?

Si, un interessante script DarwinNE, ma direi che una volta disponibile Maxima, forse con una fattorizzazione seguita dalla

------->>>> regola di Cartesio

si fa prima! ... normalmente poi con Maxima si riescono anche a trovare direttamente le radici

Se facciamo per esempio riferimento al polinomio usato in questo documento
, dove è riportato il procedimento manuale,
potremo risolvere con uno scanmap('factor, f(s)); e poi con regola di Cartesio ... e concludere dicendo che: ci sarà una radice doppia nell'origine, 2 coppie di radici complesse coniugate a parte reale nulla e, vista la doppia permanenza dei segni dell'ultimo fattore, ancora una coppia di radici coniugate con parte reale negativa


ed inoltre, volendo, in questo caso particolare, determinare le soluzioni con un semplice solve.

Sarebbe invece interessante, caro Davide, la tua idea di sintesi automatica con Cauer e Foster a partire dall' ammettenza /impedenza complessa, magari con un bell'articolo

RenzoDF ha scritto:Si, un interessante script DarwinNE, ma direi che una volta disponibile Maxima, forse con una fattorizzazione seguita dalla

------->>>> regola di Cartesio

si fa prima! ... normalmente poi con Maxima si riescono anche a trovare direttamente le radici

E' vero! Oggigiorno c'è da dire che ci sono strumenti talmente potenti che trovare delle radici di un polinomio di grado non banale è relativamente semplice. Credo che nei casi che si presentano nella pratica sia sempre possibile attaccarli con algoritmi numerici, ma in qualche situazione c'è da stare attenti. Prova a fare qualche prova con il polinomio seguente, che pure è "solo" di grado 8:




L'unico modo di affrontare il problema è utilizzare la funzione "allroots" che calcola le radici numericamente ed è una cosa molto più pesante da un punto di vista computazionale del semplice test di Hurwitz.
Nell'esempio qui sotto, invece, factor non funziona e allroots mi pesca addirittura delle radici reali.




C'è da dire che gli esempi di cui sopra (il primo è expand(a*b*c*d)) sono delle varianti dei polinomi di Wilkinson, che sono ben conosciuti per essere molto ostici da trattare numericamente.
Per i problemi di sintesi empirica, provo a rifletterci.

Si d'accordo ma dal punto di vista ingegnetistico "usare" coefficienti a 15 cifre significative è pura fantascienza !

RenzoDF ha scritto:Si d'accordo ma dal punto di vista ingegnetistico "usare" coefficienti a 15 cifre significative è pura fantascienza !

Sono totalmente d'accordo, però se il problema è mal condizionato, invece delle 15 cifre significative dei coefficienti del polinomio, le soluzioni ne hanno al massimo una o due (vedi difatti cosa fa allroots! ).