ElectroYou - la comunità dei professionisti del mondo elettrico

da **parolina** » 29 set 2011, 17:10

Salve!!!
Dovrei riuscire ad ottenere con matlab il grafico del modulo e della fase della seguente funzione di trasferimento
$T(s)=\frac{3.24*10^{10}}{s^2+s*254,6*10^3+ 3.24*10^{10}}$
Qualcuno potrebbe aiutarmi col codice?? Purtroppo non l'ho quasi mai utilizzato.
Grazie mille

da **parolina** » 29 set 2011, 17:31

E' corretto cosi??

s = tf('s');
W = (3.24*10^-10)/(s^2+254.6*10^3*s+3.24*10^-10);

bode(W)

da **IsidoroKZ** » 29 set 2011, 17:49

Puoi anche definire separatamente numeratore e denominatore. Ma aspettiamo il tuo angelo custode

DirtyDeeds.

num=[quello che e`, solo i coefficienti];
den=[stessa roba, per te: 1 254600 3.24e-10];
sysP=tf(num,den);
w=logspace(1,5); %questo per indicare le frequenze, nell'esempio parte da 10^1 fino a 10^5
[mag,phase]=bode(sysP,w);
loglog(w,squeeze(mag)),grid
semilog(w,squeeze(phase)),grid

da **parolina** » 29 set 2011, 19:58

Grazie mille...ma col codice che mi hai passato cosa visualizzo?

da **DirtyDeeds** » 29 set 2011, 20:00

Eh, questa volta non posso essere d'aiuto: ormai da alcuni anni ho abbandonato Matlab per Scilab e la sintassi per i diagrammi di Bode è differente.

Voglio comunque dare un consiglio a

parolina: non scrivere le funzioni di trasferimento mettendo direttamente i numeri, usa dei parametri: se devi cambiare qualcosa è più semplice e sei meno soggetta ad errori. La tua funzione di trasferimento è una funzione del second'ordine che può essere scritta come

$T(s)=\frac{\omega_0^2}{s^2+2\zeta\omega_0 s+\omega_0^2}$

A questo punto, fatti uno script matlab che inizi con

Codice: Seleziona tutto: w0 = .... z = ...

E, ancora meglio, normalizza

rispetto a $\omega_0$ .

da **parolina** » 29 set 2011, 20:53

ehm..ci provo!!!!Grazie infinite!!

da **DirtyDeeds** » 29 set 2011, 21:12

Altra cosa: in matlab, se devi scrivere un numero come

$3{,}24\times10^{-10}$

non scrivere

Codice: Seleziona tutto: 3.24*10^-10

ma

Codice: Seleziona tutto: 3.24e-10

PS: +1 babà :mrgreen:

da **IsidoroKZ** » 29 set 2011, 21:16

C'era gia` nel mio esempio: una roba al cioccolato a me!

da **RenzoDF** » 2 ott 2011, 14:45

parolina ha scritto:E' corretto cosi??
s = tf('s');
W = (3.24*10^-10)/(s^2+254.6*10^3*s+3.24*10^-10);
bode(W)

Si, e' corretto, se hai il Control Toolbox.
Matlab "capisce" anche la tua scrittura, ma e' assolutamente preferibile usare la variante suggerita da

DirtyDeeds per la notazione esponenziale.

Io userei la notazione abbreviata con i soli coefficienti
>> H = tf([3.24e-10],[1 254.6e3 3.24e-10])

specificando poi il range con
>>bode(H,{1e2,1e7})

Se non disponi del Control Toolbox puoi usare dei metodi alternativi, vedi per esempio

http://www.csus.edu/indiv/t/tatror/EEE_ ... ctions.pdf

BTW ricordo che per una rappresentazione asintotica e' disponibile un Free Tool di Erik Cheever

http://www.mathworks.com/matlabcentral/ ... hange/2664
con tutorial su
http://lpsa.swarthmore.edu/Bode/BodePlotGui.html
ma anche
http://www.mathworks.com/matlabcentral/ ... root-locus

PS: per me niente baba' ... solo 1MK :mrgreen:

da **parolina** » 3 ott 2011, 14:48

Grazie!!!Come siete DOLCI!!
comunque il prof adesso mi ha detto di utilizzare la notazione VETTORIALE.Avete idea di come possa fare???
Questa che devo ricavare con Matlab è la funzione di trasferimento a ciclo chiuso del mio PLL che ha :
K=127.324*10^3

$f_{LPF}=40.528*10^3$
La funzione di trasferimento è $T(s)=\frac{K 2\pi f_{LPF}}{ s^2 + s2\pi f_{LPF}+ K 2\pi f_{LPF}}$
Mi ha detto al posto di s utilizzare $j\omega$ ...aiuto

ElectroYou - la comunità dei professionisti del mondo elettrico

Funzione di trasferimento con matlab

[1] Funzione di trasferimento con matlab

[2] Re: Funzione di trasferimento con matlab

[3] Re: Funzione di trasferimento con matlab

[4] Re: Funzione di trasferimento con matlab

[5] Re: Funzione di trasferimento con matlab

[6] Re: Funzione di trasferimento con matlab

[7] Re: Funzione di trasferimento con matlab

[8] Re: Funzione di trasferimento con matlab

[9] Re: Funzione di trasferimento con matlab

[10] Re: Funzione di trasferimento con matlab

Chi c’è in linea

Informazioni

Seguici

Pubblicità

Credits