Esercizio su transitorio

[xyz]

Una possibile soluzione del problema utilizzando il Python. Questi sono i moduli Python aggiuntivi utilizzati:

https://lcapy.readthedocs.io/en/latest/index.html
https://www.sympy.org
https://numpy.org
https://matplotlib.org

Lcapy è il modulo principale che si occupa di tutta la parte circuitale, Sympy la parte di calcolo simbolico, Numpy nei calcoli numerici e Matplotlib per il disegno dei grafici.

Allego il file Python "transient.py" con tutto il codice.

transient.zip

Script in Python

(764 Byte) Scaricato 367 volte

La prima parte descrive il circuito usando solo i nomi dei componenti senza specificare i valori, in questo modo tutti i calcoli sono simbolici. Schema del circuito generato in automatico dalla sua descrizione in testo (simile a quella di SPICE):

Schematic

E' possibile generare un grafo minimo con i soli nodi e componenti:

Graph

Vengono poi calcolate le varie analisi sul circuito, la documentazione di Lcapy spiega cosa viene calcolato.

Screenshot

Una volta fatti i calcoli simbolici vengono assegnati i valori numeri ai componenti e fatta la simulazione numerica. I grafici finali contengono sia il risultato analitico, calcolato dalle espressioni simboliche e il risultato della simulazione numerica, sono praticamente coincidenti:

Plot

[micdisav]

Buongiorno,
con Matlab (R2014, v8 gira, ci vuole tb del calcolo simbolico).

Scrittura circuito in termini di trasformata di Laplace:

Sorgente FidoCadJ | Ingrandisci | Cos'è Fidocad

da [1]:
I dati sono:

R1=1 Ohm, R2=0.4 Ohm, C=10mF, L=10mH.
E' chiesto il grafico dell'andamento di VC e VL all'applicazione
di una tensione a gradino V=10V.

Sorgenti Matlab (ho inserito qualche commento):

Esercizio_su_transitorio.zip

(1.02 KiB) Scaricato 339 volte

Vari plots di tensioni e correnti:

Poca "poesia" in tutto ciò, ma...

Saluti,
Michele.

[RenzoDF]

Con Tina si fa prima

Forti questi ungheresi!

[drGremi]

Ah, fa già il calcolo della TF? Fino a che ordine arriva? Potrebbe tornarmi utile.

[Roswell1947]

Ad averlo TINA full ultima versione senza acquistarlo...

[RenzoDF]

Ah, fa già il calcolo della TF? Fino a che ordine arriva? Potrebbe tornarmi utile.

Beh, io nel test della versione 6 demo, avevo provato questa rete

2022-09-14_155248.jpg (61.13 KiB) Osservato 9538 volte

poi ,forse, si può andare anche oltre.

[RenzoDF]

Ad averlo TINA full ultima versione senza acquistarlo...

Io uso la vecchia versione 6 demo, che con qualche furbizia permette di fare un sacco di cose.

[Roswell1947]

RenzoDF
Ancora non sono riuscito a trovare un riferimento bibliografico su quel tuo metodo di trovare i poli della rete nella seguente maniera

utile ,effettivamente ,anche in questo caso.In quale università lo hai studiato?non riesco a trovare un libro in cui se ne parla.
Grazie

[RenzoDF]

Non mi ricordo dove l'ho incontrata, ma sostanzialmente non è altro che andare ad ottenere la sua impedenza o la sua ammettenza "viste" da un conveniente "taglio" che coinvolga solo due dei suoi lati; nel primo caso andando a determinare l'impedenza come somma delle impedenze delle due sottoparti (aprendo uno dei due lati), nel secondo caso andando a determinare l'ammettenza come somma delle ammettenze delle due sottoparti.
Visto che due reti, uguali quando passivate, hanno gli stessi poli, qualunque sia l'ingresso e l'uscita, il gioco è fatto.: nel primo caso ipotizziamo infatti di inserire un GIT serie su un ramo, nel secondo caso un GIC parallelo, fra due suoi punti.

PS: e scopro solo ora che te ne avevo già parlato estesamente a marzo; nel mio caso, vista l'età, è scusabile, ma per un giovane come te no di certo ... ti consiglio il Memoril ... io passerò dalle attuali 4 alle 8 al giorno.

[Roswell1947]

RenzoDF
In questo caso p.es se volessi calcolare la corrente nell'induttore metto un generatore di tensione nel ramo in cui è presense l'induttore e mi calcolo gli zeri dell'impedenza vista dai morsetti del generatore cioè R1//(1/sC)+R2+sL=0
da cui ricavo gli autovalori e poi con le condizioni iniziali determino la corrente in funzione del tempo ,corretto?