Come si implementa una retroazione in MATLAB?

[anton89]

Salve a tutti,
la domanda è del tutto generale. Come si implementa un sistema retroazionato? Qualunque esso sia. Ad esempio, sto lavorando ad un sistema di trasmissione numerico e ho bisogno di implementare un semplice PLL per il recupero di timing e portante.
Trovo difficoltà nel farlo perché il PLL l'ho studiato come circuito. Conosco gli schemi a blocchi. In Matlab riesco a simulare l'uscita di ogni blocco (il filtro d'anello lo implemento come un FIR, il comparatore di fase è un moltiplicatore, ecc.). Poi, quando devo utilizzare il risultato per fare la retroazione, mi blocco perché concretamente non so come inserire questo risultato nel codice Matlab che ho scritto
Questo non vale solo per il PLL, ma per tutti i sistemi retroazionati. Questa difficoltà la sto incontrando ora nel Matlab, ma se in futuro mi trovo a programmare un DSP che implementi uno schema di ricezione numerico e anche qui dovrò programmare un PLL???
Come diavolo si fa, me lo spiegate?

P.S. Non sto utilizzando simulink né tool già implementati in Matlab.

[simo85]

Conosco gli schemi a blocchi.

Dovresti poter implementare l'espressione del sistema nel dominio di Laplace ed implementare il feedback.

Per esempio:

Sorgente FidoCadJ | Ingrandisci | Cos'è Fidocad

[anton89]

Si, l'avevo pensato anch'io.
La difficoltà è solo implementativa.

Ad esempio, avevo pensato ad una soluzione di questo tipo. Dal segnale ricevuto, estraggo la portante con un filtro passa-banda molto stretto. Dalla portante estratta, attraverso la trasformata di Hilbert, estraggo la fase istantanea.
Dopo di che progetto il PLL non sul segnale ricevuto, ma solo sulla fase istantanea.
Dovrebbe venire così: detto F(s) il filtro d'anello e K/s il VCO (l'integratore), la funzione di trasferimento è:

T(s) = F(s)*(K/s) / 1 + F(s)*(K/s)

A questo punto con la funzione "filter" filtro il segnale fase con questa funzione di trasferimento (in termini di trasformate z, credo che sia immediato ricavarsi i coefficienti del filtro) e ottengo la fase istantanea con la quale costruire la portante coerente.

Domanda:
1) Questa procedura come vi sembra?
2) Come costruisco la portante a partire dalla fase istantanea?

[dimaios]

link
link
link
link