Equazione alle differenze
# Equazione alle differenze

Le equazioni alle differenze rappresentano la formulazione discreta della controparte continua, costituita dalle equazioni differenziali ordinarie (ODE), qualora si sia effettuata una discretizzazione del dominio di definizione della funzione incognita che costituisce la soluzione all'equazione data.

## Introduzione al problema

Consideriamo delle equazioni differenziali ordinarie di ordine $n$, nella forma implicita generale (assieme alle opportune $n$ condizioni iniziali):

    $G\left(t,y(t),y'(t),...,y^{(n)}(t)\right)=0\qquad (1)$

la cui soluzione è rappresentata dalla funzione incognita $y(t):I\to J$, definita in un dominio $I$ contenuto o al più coincidente con l'insieme dei numeri reali, e a valori in un codominio $J$, anch'esso contenuto o al più coincidente con l'insieme dei numeri reali. Affinché $y(t)$ costituisca la soluzione della $(1)$, per ogni $t$ nell'insieme $I$, questa dovrà necessariamente appartenere allo spazio funzionale $C^{n}(I)$, definite in $I$, ed ivi continue e derivabili (fino all'ordine $n$).

Vogliamo ora interessarci a ricercare la soluzione $y(t)$ dell'equazione differenziale $(1)$ soltanto in un sottoinsieme del dominio $I$ di definizione della stessa, ovvero soltanto per determinati valori di $t$. Un caso tipico è quello per cui i punti in cui si desidera conoscere la $y(t)$ risultano equispaziati tra loro: ciò equivale ad effettuare una decomposizione $D$ periodica del dominio $I$ del tipo $D=(t_{1},t_{2},\dots ,t_{n},\dots )$, in cui si avrà $t_{i}-t_{i-1}=T$, per ogni $i=1,2,\dots ,n,\dots$ (dove $T$ è il periodo della decomposizione). Naturalmente, a seconda che $I$ sia limitato o meno, tale sarà anche $D$: nel primo caso, $D$ sarà un insieme (ordinato) finito, mentre nel secondo caso sarà un insieme (ordinato) numerabile. È chiaro allora che potremo omettere la variabile $t$ nell'indicare l'insieme $D$ e riferirci ai singoli elementi ordinati di quest'ultimo attraverso un indice intero $k$ che ne individui la posizione al suo interno, ovvero, supponendo di includere anche lo $0$ (a meno di semplici traslazioni) in $D$, si avrà $t_{k}=kT$. Abbiamo ottenuto dunque la sequenza $y[k]:D\to J$.

## Equazioni alle differenze finite del primo ordine

Il passo successivo consiste nel "tradurre" discretamente le operazioni di derivazione della funzione incognita $y(t)$ che compaiono nella $(1)$. A tal proposito, prendiamo come esempio la derivata prima e consideriamo il generico punto $t_{0}$ appartenente a $I$, che supponiamo essere anche un elemento di $D$. Essendo $y(t)$ appartenente a $C^{n}(I)$, essa sarà derivabile in $t_{0}$, in cui esisteranno, finite, derivata destra e sinistra, uguali naturalmente tra loro. Avremo pertanto:

    $y'(t_{0})=\lim _{t{\rightarrow }t_{0}^{+}}{\frac {y(t)-y(t_{0})}{t-t_{0}}}=\lim _{t{\rightarrow }t_{0}^{-}}{\frac {y(t_{0})-y(t)}{t_{0}-t}}\qquad (2)$

Approssimiamo la derivata prima con il suo rapporto incrementale destro, ovvero, se $t_{1}=t_{0}+T$, si potrà scrivere:

    $y'(t_{0})\approx {\frac {y(t_{1})-y(t_{0})}{T}}\qquad (3)$

Generalizzando la $(3)$ ad un punto qualsiasi $k$ appartenente a $D$ ($t=kT$), si potrà scrivere:

    $y'(t=kT)\approx {\frac {y[k+1]-y[k]}{T}}\qquad (4)$

La $(4)$ rappresenta la differenza finita prima in avanti di $y(t)$ in $t=kT$; l'errore che si commette nell'approssimare la derivata col rapporto incrementale sarà naturalmente tanto minore quanto più $T$ è piccolo. Per quanto detto a proposito della $(2)$, approssimiamo ora la derivata prima con il suo rapporto incrementale sinistro, ovvero, se $t_{1}=t_{0}-T$, si potrà scrivere:

    $y'(t_{0})\approx {\frac {y(t_{0})-y(t_{1})}{T}}\qquad (5)$

Generalizzando la $(5)$ ad un punto qualsiasi $k$ appartenente a $D$ ($t=kT$), si potrà scrivere:

    $y'(t=kT)\approx {\frac {y[k]-y[k-1]}{T}}\qquad (6)$

La $(6)$ rappresenta la differenza finita prima all'indietro di $y(t)$ in $t=kT$. Le formule $(4)$ e $(6)$ sono due modi possibili per discretizzare l'operazione di derivazione e che ci hanno condotto, in modo naturale, ad esprimere la $y'(t)$ in funzione del periodo di decomposizione $T$ dell'intervallo $I$, e dei valori di $y(t)$ nel punto di derivazione considerato $kT$ e nel punto immediatamente dopo (o immediatamente prima) che compare nell'ordinamento di $D$. Possiamo dunque, intanto, interpretare la $(1)$, nel caso in cui $n=1$, nel caso discreto come:

    $F(k,y[k],y[k+1])=0\qquad (7)$

se si è utilizzata la $(4)$ per approssimare la derivata prima, o come:

    $F(k,y[k],y[k-1])=0\qquad (8)$

se si è utilizzata la $(6)$ per approssimare la derivata prima. Le formule $(7)$ e $(8)$ rappresentano due equazioni alle differenze finite del primo ordine, in forma implicita. Si noti che, così come nel caso continuo abbiamo bisogno di una condizione iniziale $y(t_{0})$ per risolvere una generica ODE del primo ordine, allo stesso modo necessitiamo di una condizione iniziale $y[m]$ per risolvere una generica equazione alle differenze finite del primo ordine. Equazioni simili sono risolvibili in modo molto semplice; vediamone un esempio:

    $3y[n]+y[n-1]=0\ {\mbox{per}}\ n\geq 0\qquad (9)$

Si tratta di un'equazione alle differenze finite, del primo ordine, omogenea, lineare, a coefficienti costanti (questi infatti non dipendono dalla variabile indipendente $n$). Si noti che abbiamo esplicitamente limitato lo studio della soluzione della $(9)$ a quei valori di $n$ appartenenti a $\mathbb {N} _{0}^{+}$. Prendiamo come condizione iniziale $y[-1]=1$ e procediamo ricorsivamente:

    $3y[0]+y[-1]=0\Rightarrow y[0]=-{\frac {1}{3}}$
    $3y[1]+y[0]=0\Rightarrow y[1]={\frac {1}{9}}$
    $3y[2]+y[1]=0\Rightarrow y[2]=-{\frac {1}{27}}$
    $\dots$

Possiamo già ricavare una forma esplicita per la soluzione $y[n]$, ovvero $y[n]={\tfrac {1}{(-3)^{n+1}}}$, per $n\geq 0$, e che rappresenta di fatto una successione a termini di segno alterno e convergente a $0$ per $n$ tendente a $\infty$ . La radice dell'equazione algebrica omogenea associata alla $(9)$, ovvero $3h+1=0\Rightarrow h=-1/3$, viene tipicamente chiamata polo dell'equazione. Tutte le equazioni alle differenze, lineari e a coefficienti costanti, di ordine qualsiasi, una volte note tutte le condizioni iniziali, possono essere risolte ricorsivamente. Per completare la "discretizzazione" della $(1)$, dobbiamo discretizzare le restanti derivate: questo viene fatto in maniera analoga a quanto visto per l'approssimazione della derivata prima, utilizzando i rapporti incrementali $k$-esimi della derivata $(k-1)$-esima in modo progressivo. Si arriva dunque a poter scrivere la forma implicita di una generica equazione alle differenze finite, di ordine $n$:

    $F(k,y[k],y[k-1],...,y[k-n])=0\qquad (10)$

avendo supposto di utilizzare sempre le differenze finite prime all'indietro per l'approssimazione delle derivate.

Le equazioni alle differenze vengono introdotte in modo naturale nello studio di fenomeni che avvengono con cadenza periodica (o comunque approssimabili o descrivibili con modelli matematici di tipo discreto); in particolare, si prestano molto bene a descrivere il comportamento dinamico di sistemi tempo-discreti o spazio-discreti nell'elaborazione dei segnali e nei controlli automatici. Per questi casi, particolarmente importanti sono le equazioni alle differenze finite lineari e a coefficienti costanti: tali proprietà, infatti, riflettono in modo naturale le caratteristiche del sistema di cui l'equazione (o il sistema di equazioni) rappresenta il modello matematico. La linearità discende dall'omonima proprietà del sistema che si sta descrivendo, mentre la costanza dei coefficienti rispetto al tempo (discreto) discende dall'invariabilità dei parametri del sistema al variare del tempo. Un classico esempio è un simulatore discreto di un semplice ramo RC: scegliendo come grandezza d'ingresso la tensione applicata ai capi del ramo e come grandezza d'uscita la tensione misurata ai capi della capacità, è possibile discretizzare l'equazione differenziale che modellizza matematicamente questo semplice sistema, al fine di ridurla ad un'equazione alle differenze finite, del primo ordine, facilmente gestibile da un calcolatore, tipicamente un software CAD di progettazione di circuiti elettronici. Supposti i componenti R e C lineari e stazionari, tale sarà anche l'equazione differenziale descrittiva e quindi la diretta implementazione discreta di questa.

## Equazioni alle differenze finite, lineari ed omogenee

Una generica equazione alle differenze finite, di ordine $m$, lineare, a coefficienti costanti ed omogenea, si può scrivere in forma implicita come:

    $a_{m}y[n-m]+a_{m-1}y[n-(m-1)]+...+a_{1}y[n-1]+a_{0}y[n]=0\ {\mbox{per}}\ n\geq 0\qquad (11)$

e facilmente riducibile in forma normale (in funzione, cioè, solo di $y[n-m]$). Date le $m$ condizioni iniziali $\{y[-1],...,y[-m]\}$, è possibile risolvere la $(11)$ in maniera ricorsiva. Ad ogni modo, se l'ordine $m$ dell'equazione è maggiore o uguale a 4, è possibile ricavare la soluzione come somma di modi naturali, come somma cioè di opportune funzioni del tempo discreto $n$, attraverso la ricerca delle radici dell'equazione algebrica associata alla $(11)$; tale equazione si ricava sostituendo, alla generica differenza finita $y[n-p]$, un'incognita elevata alla potenza $m-p$, diciamo $h^{m-p}$. Facciamo un esempio:

    $2y[n-2]+3y[n-1]+y[n]=0\ {\mbox{con}}\ y[-1]=1,y[-2]=0\ {\mbox{per}}\ n\geq 0\qquad (12)$

L'equazione algebrica associata alla $(12)$ è:

    $h^{2}+3h+2=0$

le cui radici sono $-1$ e $-2$. La soluzione generica della $(12)$ è dunque:

    $y[n]=A(-1)^{n}+B(-2)^{n}$

dipendente da due costanti $A$ e $B$, che si ricavano a partire dalle condizioni iniziali:

    $y[-1]=-A-{\frac {1}{2}}B=1$
    $y[-2]=A+{\frac {1}{4}}B=0$

Risolvendo il sistema lineare in due equazioni e due incognite, si ricava $A=1$ e $B=-4$, per cui la soluzione è:

    $y[n]=(-1)^{n}-4(-2)^{n}\ {\mbox{per}}\ n\geq 0$

Osserviamo che i poli dell'equazione sono $-1$ e $-2$, e che i modi naturali aperiodici della soluzione sono rappresentati da $(-1)^{n}$ e da $(-2)^{n}$.

## Equazioni alle differenze finite, lineari e complete

Più complesso è il caso in cui la $(11)$ non sia omogenea, ma presenti un termine noto (in generale, una funzione $x[n]$, da sola o assieme alle sue differenze finite, nota su tutto l'intervallo in cui varia $n$):

    ${\begin{aligned}a_{m}y[n-m]+a_{m-1}&y[n-(m-1)]+\dots +a_{1}y[n-1]+a_{0}y[n]=\\&b_{r}x[n-r]+b_{r-1}x[n-(r-1)]+\dots +b_{1}x[n-1]+b_{0}x[n]=u[n]\ {\mbox{per}}\ n\geq 0\qquad (13)\end{aligned}}$

Il termine noto può dunque esprimersi, in generale, attraverso una sequenza $u[n]$, nota a priori. A seconda della forma del termine noto, si utilizzano tecniche analoghe a quelle utilizzate per la risoluzione delle ODE non omogenee. Si decompone, innanzitutto, la $(13)$ nella sua omogenea associata e si risolve quella, ottenendo una soluzione generale; a questa si somma una soluzione particolare dipendente dalla forma di $u[n]$. Nella teoria dei sistemi, tale decomposizione ha un ben preciso significato fisico: il termine noto rappresenta la sollecitazione (o ingresso) del sistema, e l'uscita è la sequenza incognita che occorre ricavare (risposta del sistema). L'equazione omogenea associata alla $(13)$ avrà, come soluzione, la cosiddetta evoluzione libera del sistema, la quale dipende solo dalle condizioni iniziali e dai poli (parametri intrinseci) di quest'ultimo. La soluzione particolare costituisce, invece, la risposta forzata del sistema, direttamente dipendente dall'ingresso $u[n]$. Con un esempio:

    $2y[n]+y[n-1]=2^{n}\ {\mbox{con}}\ y[-1]=1\ {\mbox{e per}}\ n\geq 0\qquad (14)$

Osserviamo che la restrizione dell'intervallo su cui varia $n$ per la ricerca della soluzione coincide con la proprietà del sistema (descritto dal modello matematico $(14)$) di causalità. La soluzione generale dell'equazione omogenea associata alla $(14)$ è:

    $y_{h}[n]=A(-{\frac {1}{2}})^{n}$

con $A$ costante da determinare in base alle condizioni iniziali. Per cercare una soluzione particolare della $(14)$, proviamo a cercarne una della forma:

    $y_{p}[n]=B(2)^{n}$

e sostituiamo nella $(14)$:

    $2B(2)^{n}+B(2)^{n-1}=2^{n}$

da cui $B={\tfrac {2}{5}}$ e poiché:

    $y[n]=y_{h}[n]+y_{p}[n]=A(-{\frac {1}{2}})^{n}+{\frac {2}{5}}(2)^{n}$

si avrà, utilizzando l'informazione sulla condizione iniziale $y[-1]=1$, che: $A=-{\tfrac {2}{5}}$ e quindi:

    $y[n]=-{\frac {2}{5}}(-{\frac {1}{2}})^{n}+{\frac {2}{5}}(2)^{n}\ {\mbox{per}}\ n\geq 0.$

Si nota l'evoluzione libera, dipendente dal polo $h=-{\tfrac {1}{2}}$ e dalla condizione iniziale, e l'evoluzione forzata, direttamente dipendente dall'ingresso $u[n]=2^{n}$.

Aggiungiamo che l'evoluzione libera di un sistema tempo-discreto LTI (lineare e tempo-invariante), il cui modello matematico consiste dunque in una (o più) equazioni alle differenze finite, lineari, a coefficienti costanti, tende asintoticamente a $0$ per $n$ che tende a infinito se e solo se tutti i suoi poli sono, in modulo, minori di $1$ (sistema asintoticamente stabile).

## Risoluzione di equazioni differenziali per mezzo della trasformata Z

Un altro modo per risolvere equazioni differenziali lineari, a coefficienti costanti, di ordine $m$ qualsiasi, è passare per la trasformata $Z[\cdot ]$. Ricordiamo che, in via del tutto generale, l'operazione di $Z$-trasformazione sulla sequenza $y[n]$, definita su un dominio $A\subseteq Z$, sia esso limitato (tra $N_{0}$ e $N_{1}$), sia esso illimitato, superiormente e/o inferiormente viene definita come:

    $Z[y[n]]=\sum _{k=N_{0}}^{N_{1}}y[k]z^{-k}=Y(z)\qquad (15)$

dove $Y(z)$ è una funzione che assume valori complessi, e definita su un sottoinsieme del campo complesso ($z=\rho e^{j\Omega }$; si può dimostrare che $\Omega$ è una quantità adimensionale ed è data dal prodotto $\Omega =\omega T$, con $\Omega$ pulsazione angolare e $T$ periodo di campionamento della funzione $y(t)$ da cui si ricava la sequenza $y[n]$). Per la proprietà derivativa della trasformata $Z$, si ha che, restringendo lo studio a sequenze causali:

    $Z[y[n-1]]=z^{-1}Y(z)+y[-1]\qquad (16)$

e, più in generale:

    $Z[y[n-m]]=z^{-m}Y(z)+z^{-(m-1)}y[-1]+\dots +y[-m]\qquad (17)$

Utilizzando la $(17)$ ed effettuando la $Z$-trasformazione di ambo i membri della $(13)$, si ottiene:

    $a_{m}z^{-m}Y(z)+a_{m-1}z^{-(m-1)}Y(z)+\dots +a_{1}z^{-1}Y(z)+a_{0}Y(z)-H(z)=U(z)\Rightarrow Y(z)={\frac {U(z)+H(z)}{D(z)}}\qquad (18)$

dove

    $D(z)=a_{m}z^{-m}+a_{m-1}z^{-(m-1)}+\dots +a_{1}z^{-1}+a_{0}$

ed $H(z)$ è una funzione polinomiale in $z$, con coefficienti dati da combinazioni lineari delle condizioni iniziali e $U(z)$ è la $Z$-trasformata della sequenza nota $u[n]$. In questa forma, la $Y(z)$ è rappresentata da un rapporto di polinomi in $z$ ed è facilmente riducibile (in modo analogo a ciò che si fa per ridurre in fratti semplici una trasformata di Laplace composta da un rapporto di polinomi in $s$) ad una somma di monomi semplici, le cui antitrasformate sono del tipo $Aa^{n}u_{-1}[n]$, con $u_{-1}[n]$ sequenza gradino unitario. Ridotta in fratti semplici la $(18)$, si antitrasforma $Y(z)$ per ottenere la sequenza soluzione $y[n]$.

Osserviamo come sia possibile scrivere la $(18)$ anche come:

    $Y(z)={\frac {H(z)}{D(z)}}+{\frac {U(z)}{D(z)}}={\frac {H(z)}{D(z)}}+{\frac {N_{u}(z)}{D(z)D_{u}(z)}}$

dove si è posto $U(z)={\tfrac {N_{u}(z)}{D_{u}(z)}}$, cioè come la somma di due contributi: il primo, dipendente dalle condizioni iniziali, che si evolve secondo i poli (modi naturali) del sistema (di cui la $(18)$ rappresenta il modello matematico; i poli sono proprio le radici dell'equazione caratteristica associata alla $(13)$, ovvero gli zeri del polinomio $D(z)$); il secondo, composto a sua volta da una parte che evolve secondo i modi naturali del sistema ("innescata" dall'applicazione dell'ingresso, fattore $D(z)$) e da un'altra componente che evolve esattamente come l'ingresso (cioè secondo i poli di $U(z)$, gli zeri cioè del polinomio a denominatore $D_{u}(z)$), che costituisce la risposta forzata del sistema. Se quest'ultimo è asintoticamente stabile, la parte della risposta complessiva che tende a $0$ al tendere di $n$ all'infinito è chiamata anche risposta transitoria, mentre quella parte di $y[n]$ che rimane comunque limitata al tendere di $n$ all'infinito (a seconda della limitatezza di $u[n]$) viene anche definita risposta forzata del sistema in esame.

