Esercizio studio sistema tempo contintinuo

claudio_p88
Dato il sistema con $A=((0,-2,1),(8,0,4),(0,0,-1))$ $B=((0),(1),(1))$ $C=(0,0,1)$ e $D=1$
a)Studiare eccitailità e osservabilità;
b)Calcolare $phi(t), H(t),psi(t),omega(t)$;
c)Si calcoli l'insieme degli stati iniziali, per cui la risosta in evoluzione libera nello stato tende asintoticamente a 0;
d)si calcoli la risposta forzata e, se esiste, permanente all'ingresso $u(t)=t$;
Mi calcolo matrice di raggiungibilità $(B, AB, A^2B)=((0,-1,-9),(1,4,-12),(1,-1,1))$ che ha $det(B, AB, A^2B)=58ne0$, dunque il sistema dovrebbe essere completamente raggiungibile, e quindi eccitabile.
Ci calcoliamo ora la matrice di osservabilità $((C),(CA),(CA^2))=((0,0,1),(0,0,-1),(0,0,1))$, che ha rango = 1, quindi solo uno dei modi è osservabile, la matrice A ha autovalori $lambda_1=-1,lambda_{2,3}=4i$, inoltre la matrice A dovrebbe non essere diagonalizzabile, in quanto la molteplicità algebrica dell'autovalore 4i non corrisponde con quella geometrica. Gli autovettori risultanti dovrebbero essere per $lambda_1$ $u_1=(k frac{9}{17},k frac{-4}{17},k)^T=(9k,-4k,17k)^T$ per $lambda_{2,3}$ $u_{2,3}=(1,2i,0)^T$, ora mi chiedo se A è non diagonalizzabile allora la matrice tranzione nello stato $phi(t)$ non esiste?

Risposte
D4lF4zZI0
E' vero la matrice $A$ non è diagonalizzabile, ma puoi arrivare alla matrice di transizione per altra via; riflettici: hai tutti gli strumenti ;)

claudio_p88
dovrebbe essere $ (SI-A)^-1=frac{1}{(s+1)(s^2+16)}((s^2+s,-2-2s,-8+s),(8+8s,s+s^2,4s+8),(0,0,s^2+16)) $ allora avrò che
$ frac{s^2+s}{(s+1)(s^2+16)}=\frac{s}{s^2+16}$ antitrasformando \(\displaystyle cos(4t) \)
$ frac{-2-2s}{(s+1)(s^2+16)}=\frac{-2}{s^2+16}$ antitrasformando $-frac{1}{2}sin(4t)$
$ frac{-8+s}{(s+1)(s^2+16)}=\frac{-9}{17(s+1)}+frac{9s+8}{17(s^2+16)}$ antitrasformando =$frac{1}{17}(-9e^{-t}+9cos(4t)+2sin(4t))$;
$ frac{8+8s}{(s+1)(s^2+16)}=\frac{8}{s^2+16}$ antitrasformando $2sin(4t)$;
$ frac{4s+8}{(s+1)(s^2+16)}=\frac{-5}{17(s+1)}+frac{5(s-16)}{17(s^2+16)}$ antitrasformando $frac{4}{17}(e^{-t}-cos(4t)+frac{9}{2}sin(4t))$;
$ frac{s^2+16}{(s+1)(s^2+16)}=\frac{1}{s+1}$ antitrasformando \(\displaystyle e^{-t}\).
Quindi dovrebbe essere $phi(t)=((cos(4t),-frac{1}{2}sin(4t),frac{1}{17}(-9e^{-t}+9cos(4t)+2sin(4t))),(2sin(4t),cos(4(t)),frac{4}{17}(e^{-t}-cos(4t)+frac{9}{2}sin(4t))),(0,0,e^{-t}))$.
E $H(t)=((-frac{1}{2}sin(4t)+frac{1}{17}(-9e^{-t}+9cos(4t)+2sin(4t))),(cos(4t)+frac{4}{17}(e^{-t}-cos(4t)+frac{9}{2}sin(4t))),(e^{-t}))$ e
$psi(t)=((0,0,e^{-t}))$, dunque solo modo \(\displaystyle e^{-t} \) relativo all'autovalore -1 è osservabile.

D4lF4zZI0
Dovrebbe esserci un errore con l'inversa, quella che mi trovo io ( controllata con Matlab ) è la seguente:
$ (sI-A)^-1=( ( s/(s^2+16) , -2/(s^2+16) , (s-8)/((s^2+16)(s+1)) ),( 8/(s^2+16) , s/(s^2+16) , (4(s+2))/((s^2+16)(s+1)) ),( 0 , 0 , 1/(s+1) ) ) $
ricontrolla il termine di posto $2-3$

claudio_p88
Ho modificato e provato a scrivere la soluzione ;-)

claudio_p88
Inoltre $W(s)=frac{CB}{sI-A}+D= frac{1}{(s+1)(s^2+16)}+2$?

D4lF4zZI0
Ricontrolla $W(s)=C(sI-A)^-1B+D=(s^3+2s^2+16s+32)/((s+1)(s^2+16))$

claudio_p88
sì, hai ragione scusa... Per il resto lo studio dell'osservabilità e dell'eccitabilità è svolto bene? Per il punto c, mi servirebbe un suggerimento che non riesco a capire da dove cominciare...

claudio_p88
per il punto d, invece si ha che $Y(s)=W(s)U(s)$ e abbiamo che $u(t)=t$ allora $U(s)=frac{1}{s^2}$,
$y(t)=frac{-1}{s}+frac{1}{s^2}+frac{1}{s+1}$, per quando riguarda la risposta a regime, non dovrebbe esistere in quanto c'è un polo di molteplicità algebrica due nell'origine.

D4lF4zZI0
Allora la traccia di chiede di trovare l'insieme di stati iniziali per il quale l'uscita in evoluzione libera tende a $0$.
La risposta in evoluzione libera vale:
$ y_l(t)=Ce^(At)x(0)== ( 0 \ \ 0 \ \ 1 )( ( cos(4t) , -1/2sen(4t) , 1/17(-9e^(-t)+9cos(4t)+2sen(4t)) ),(2sin(4t) ,cos(4t) ,4/17(e^(-t)-cos(4t)+9/2sin(4t)) ),( 0 , 0 , e^(-t) ) )( ( x_1(0) ),( x_2(0) ),( x_3(0) ) ) $
Quindi:
$ y_l(t)=( ( 0 ),( 0 ),( e^(-t) ) ) ( ( x_1(0) ),( x_2(0) ),( x_3(0) ) )=x_3(0)e^(-t) $
Come si vede, l'evoluzione libera del sistema tende asintoticamente a zero per qualunque terna di valori degli stati iniziali ( ed è proprio nulla quando $x_3(0)=0$ )

D4lF4zZI0
Per il punto d c'è qualcosa che non va:
$ Y(s)=U(s)W(s)=1/s^2(s^3+2s^2+16s+32)/((s+1)(s^2+16))=1/s^2(1+16/((s+1)(s^2+16)))=1/s^2+16/(s^2(s+1)(s^2+16)) $
A te l'antitrasformata

claudio_p88
$ y (t)=2t+e ^{-t}-1$ la risposta a regime permanente non esiste in quanto il limite va a infinito. Una domanda l'evoluzione forzata vale $ e^{-At} B $?

D4lF4zZI0
La risposta forzata che hai calcolato è sbagliata; ricontrollala
Per quanto riguarda l'altra domanda, l'uscita di un sistema vale:
$ y(t)=Ce^(At)x(0)+Cint_(0)^(t) e^(A(t-tau))Bu(tau) d tau $
quindi l'evoluzione forzata è data dall'integrale di convoluzione:
$ y_f(t)=Cint_(0)^(t) e^(A(t-tau))Bu(tau) d tau $

claudio_p88
Ho modificato

D4lF4zZI0
Allora
$ Y(s)=1/s^2 (s^3+2s^2+16s+32)/((s+1)(s^2+16))=1/s^2 ((s+1)(s^2+16)+s^2+16)/((s+1)(s^2+16))=
1/s^2 (1+1/(s+1))=1/s^2-1/s+1/s^2+1/(s+1)=2/s^2-1/s+1/(s+1) $
Da cui:
$ y(t)=2t-1+e^(-t) $

claudio_p88
$ Y(s)=U(s)W(s)=1/s^2(s^3+2s^2+16s+32)/((s+1)(s^2+16))=frac { s+2}{s^2(s+1)} $ come fá a venirti$ Y(s)=U(s)W(s)=1/s^2(1+16/((s+1)(s^2+16)))=1/s^2+16/(s^2(s+1)(s^2+16)) $?

D4lF4zZI0
Hai ragione, avevo semplificato male, ora ho modificato ;)

claudio_p88
oddio, stavo impazzendo,hahaha.

claudio_p88
Ho un ulteriore punto che mi chiede: è possibile stabilizzare asintoticamente il sistema con reazione dello stato? è possibile stabilizzare asintoticamente il sistema con reazione dell'uscita? Per quanto riguarda la reazione dell'uscita, dovrebbe essere impossibile, in quanto il sistema non è del tutto osservabile, per la reazione dello stato, dovrebbe essere possibile, anche se non ho ben presente come fare, tento di abbozzare una soluzione:
calcolo $A+BK$ e impongo k tale che ... E da qui non saprei come andare avanti.

D4lF4zZI0
Con la retroazione dell'uscita riesci a stabilizzare il sistema; per quanto riguarda invece la retroazione dello stato, avendo un sistema non osservabile, come potresti retroazionare lo stato?

claudio_p88
Sí, hai ragione

Rispondi
Per rispondere a questa discussione devi prima effettuare il login.