Calcolo sfasamento con Matlab

[Filippo]

Buongiorno a tutti
Volevo chiedervi se conoscete un buon metodo per il calcolo dello sfasamento( gradi o rad) tra 2 sinusoidi aventi la stessa frequenza ma ampiezze diverse.
Ho già provato con arcoseni ma non riesco ad avere in uscita una lettura corretta. Possibile che non esista un blocco già predisposto in Matlab??

[g.schgor]

Che c'entra la differenza di ampiezza?
In che forma hai le frequenze?

[Russell]

Possibile che non esista un blocco già predisposto in Matlab??

ci saranno tecniche migliori, ma se i segnali sono puliti allora
trova il max di xcorr(s1, s2), il suo indice è relativo allo sfasamento

[Filippo]

Che c'entra la differenza di ampiezza?
In che forma hai le frequenze?

Con la stessa ampiezza potrei confrontale..
In che senso forma delle frequenze? La frequenza è di 50 Hz..

[g.schgor]

Intendo: come sono espresse le frequenze?

[Filippo]

Intendo: come sono espresse le frequenze?

Sto lavoranto in simulink. Praticamente devo misurare lo sfasamento tra tensione e corrente, dove la tensione è quella di rete, che la setto io con un generatore, mentre la corrente la misuro sul carico qualsiasi.
Quindi in uscita ho 2 sinusoidi sfasate.

[Russell]

ribadisco

Codice: Seleziona tutto

[M, indice]=max( xcorr(s1, s2) );


indice, a parte un valore costante, è lo sfasamento (in campioni)


la costante di cui sopra è relativa alla lunghezza dei segnali, il modo piu' semplice/stupido per trovarla è

Codice: Seleziona tutto

[M, costante]=max( xcorr(s1, s1) );


non è difficile... no?

[Filippo]

ok.. Grazie mille delle risposte, proverò a usare la correlazione..

[dimaios]

Con l'ausilio della trigonometria e di un filtro passa basso si risolve il problema in modo elegante.




Moltiplicando le due sinusoidi si ha :




Il primo termine è una costante che dipende dallo sfasamento delle due sinusoidi mentre il secondo può essere filtrato con un filtro LP adeguato.

A questo punto basta ricavare .

Nel caso si abbia molto rumore sovrapposto e si desideri stimare più velocemente la differenza di fase bisogna impiegare la teoria degli stimatori.