Funzione di trasferimento con matlab

[DirtyDeeds]

Punto 1 ...il punto non l'ho capito!!!

Qui ci va il manuale! Guarda bene anche gli esempi

Punto 2 ..come spaziarle in modo logaritmico??Ho riletto il post ma.....ignorantissima io!!

Se scrivi

Codice: Seleziona tutto

f = 10:10:50;

Ottieni 5 numeri, tra 10 e 50, tali per cui (per che va da 1 a 4). I numeri sono quindi equispaziati.

Se scrivi (v. anche logspace)

Codice: Seleziona tutto

f = logspace(log10(10),log10(50),5)

Ottieni di nuovo 5 numeri compresi tra 10 e 50, con la caratteristica, però, che è il rapporto tra due numeri successivi ad essere costante: .

Ma



implica



ovvero



In questo caso è il logaritmo delle frequenze ad essere equispaziato. Immagina di voler fare un grafico con che va da a : qual è la soluzione più economica? Senza contare che tenendo conto degli andamenti asintotici dei diagrammi di Bode, spaziare linearmente le frequenze ha poco senso.

[parolina]

Ok!!!Capito!!
Dunque ricapitolando...nel mio esercizio....(perché alla fine io sto ricavando la funzione di trasferimento di quell'esercizio che ti mandai...Pll ad 1MHz con intervallo di aggancio di 400kHz) posso fare cosi':

Codice: Seleziona tutto

K=127324;
j = sqrt(-1);
f = logspace(log10(1e3),log10(10e6),100);
w = 2*pi*f;
T = K*2*pi*40.528e3./((j*w).^2 + (j*w)*2*pi*40.528e3 + K*2*pi*40.528e3);
loglog(f,abs(T))
xlabel('Hz')


per la fase...utilizzo loglog(f,angle(T)) giusto??

[DirtyDeeds]

In alternativa, se ti piace di più (e parametrizza sempre tutto),

Codice: Seleziona tutto

    K=127324;
    flpf = 40.528e3;
    j = sqrt(-1);
    f = logspace(log10(1e3),log10(10e6),100);
    s = j*2*pi*f;
    T = K*2*pi*flpf./(s.^2 + s*2*pi*flpf + K*2*pi*flpf);
    loglog(f,abs(T))
    xlabel('Hz')
    figure;
    semilogx(f,angle(T))


Nota che per la fase il grafico deve essere semilogaritmico.

[parolina]

Sei un grandeeeeeeeeeeee...è venuta proprio bellina!!!!!

[IsidoroKZ]

Ma il listato che avevo messo io non funzionava?

[DirtyDeeds]

Da quello che ho capito, IsidoroKZ, il prof. di parolina non voleva che lei usasse il comando bode. Bah!

[parolina]

Mi controllate quest'altra??
Questa è quella del filtro!!
con Klpf=1 e sempre flpf=40.528kHz
Ho fatto cosi':

Codice: Seleziona tutto

Klpf=1;
    flpf = 40.528e3;
    j = sqrt(-1);
    f = logspace(log10(1e3),log10(10e9),100);
    s = j*2*pi*f;
    T = Klpf./(1+(s./2*pi*flpf));
    loglog(f,abs(T))
    xlabel('Hz')
    figure;
    semilogx(f,angle(T))


forse ho sbaglaito qualcosina perché non viene proprio bene ahahahhaha

[DirtyDeeds]

Non ha senso fare quel diagramma di Bode fino a 10 GHz Fermati molto più in basso

[parolina]

Si si infatti era una prova che avevo fatto...quindi aggiustato quel parametro per il resto va bene??
No perché il modulo e la fase andando a compilare non vengono fuori come dovrebbero!!!

[DirtyDeeds]

Sì.

Nella riga

Codice: Seleziona tutto

T = Klpf./(1+(s./2*pi*flpf));

puoi togliere il punto dopo la s, perché 2*pi*flpf non è un vettore:

Codice: Seleziona tutto

T = Klpf./(1+s/2*pi*flpf);

Poi, tieni conto che tu gli andamenti asintotici dei diagrammi di Bode li conosci, quindi non è difficile capire se nel grafico ci sono errori.