Pagina 1 di 2

Dominio della frequenza

MessaggioInviato: 25 giu 2021, 16:15
da dannywall
salve sto studiando per un esame di controlli automatici, in particolare devo analizzare l'andamento in funzione del dominio della frequenza e non so come fare dato che con matlab sono agli inizi...

grazie per l'aiuto

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 16:24
da MarcoD
Cosa cerchi in particolare?
Matlab è solo uno strumento, non è la teoria, che dovresti studiare.
I criteri per la stabilità e risposta li hanno sviluppati negli anni '30 Bode e Nyquist quando non esistevano i calcolatori numerici e le risposte in frequenza si calcolavano anche in modo grafico. O_/

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 16:33
da dannywall
come detto sono agli inizi dello studio e ho un problema di questo tipo
Immagine 2021-06-25 163048.jpg

e ho che all'esame è richiesto che i calcoli siano fatti tramite l'uso di matlab...

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 16:56
da MarcoD
Troppa fatica risolverlo.
Uno spunto: scrivi la F(s) come 10/(S x (1 + S/20)^2)
Ha un polo nell'origine e due poli per S = 20 rad/s

Secondo me, la rete correttiva, al crescere della pulsazione deve avere prima uno zero poi un polo.
Forse serve anche definire un coefficiente di guadagno per pulsazione zero.

Visto che propongono 4 risposte, fai "reverse engineering" e ricava per ogni risposta i valori dello zero e del polo.

Forse altri ti aiuteranno meglio. Auguri per l'esame

O_/

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 18:02
da g.schgor
Puoi tracciare Bode sostituendo ad s:j \cdot \omega _n (con n da 1 a 1000 e \omega _n=n/10)
Puoi così verificare le varie risposte.

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 18:20
da dannywall
g.schgor ha scritto:Puoi tracciare Bode sostituendo ad s:j \cdot \omega _n (con n da 1 a 1000 e \omega _n=n/10)
Puoi così verificare le varie risposte.


il problema è che non riesco ad assegnare la funzione, avendo scritto:

f(j*wn)=4000/(j*wn*(j*wn+20)^2)
ed avendo provato sia n come intervallo che come vettore mi da quest'errore:

Incorrect dimensions for raising a matrix to a power. Check that the matrix is square and the power is a scalar. To perform elementwise matrix powers, use '.^'.

in realtà ho anche provato a vedere la F(s) e tracciare nichols ma non so se è la stessa cosa

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 18:29
da dannywall
se fare così è valido probabilmente ho trovato la soluzione... e sarebbe la 4

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 18:44
da MarcoD
Probabilmente la soluzione 4 ( quella con zero a 3 rad/s) è quella corretta.
Metodi differenti, se corretti, devono portare allo stesso risultato.

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 18:45
da g.schgor
Ma usi le variabili 'indicizzate'?
Ecc0 come tracciare la funzione (in Mathcad)
Forum210625 - Copia.gif
Forum210625 - Copia.gif (14.87 KiB) Osservato 6842 volte

Re: Dominio della frequenza

MessaggioInviato: 25 giu 2021, 22:20
da g.schgor
Aggiungendo R (risp n.4=C2) si ottengono i diagrammi di C1 (verde), C2 (blu)
e il loro prodotto C (rosso):
Forum210625gr.gif
Forum210625gr.gif (24.13 KiB) Osservato 6796 volte

\omega =23r/s (n=230) riulta il valore di taglio
con margine di fase di 30° cdd.