[ Teoria dei sistemi e Controlli automatici ]esercizio su calcolo matrice di transizione

Messaggioda polid » 10/12/2018, 11:49

Buongiorno . Ho eseguito questo tipo di esercizio per la prima volta .
Sono riuscito ad arrivare alla matrice modale pero' non riesco a calcolare gli autovettori .
Mi servirebbe qualcuno che gli dia un'occhiata e verifica se il ragionamento eseguito e' giusto .

Per favore , datemi un'aiuto ..



Immagine


Svolgimento esercizio:

$ dot(x) =( ( 1 , 1 , 1 ),( 0 , -1 , 1 ),( 0 , 0 ,2 ) ) x+ ( ( 0 ),( 0 ),( 1 ) ) *u $

$ y= ( 1 \ \ 0 \ \ 1 ) x$

(questo è un sistema strettamente proprio perché’ non compare u(t) , non c’è dipendenza esplicita dell’uscita y(t) dell’ingresso u(t)

La matrice di stato è la seguente :

$ F = ( ( 1 , 1 , 2 ),(0 , -1 , 1 ),( 0 , 0 , 2 ) ) $

L’obiettivo è quello di calcolare la matrice esponenziale di F , e^Ft è definita nel seguente modo :

$ e^Ft=I+Ft+1/(2!)(Ft)^2+1/(3!)(Ft)^3...$

Mi calcolo gli auto valori della matrice .

Il determinate $ det(F-lambda I) = det(F)= -2$

Il polinomio caratteristico e’ : $ -lambda ^3+2lambda ^2+lambda -2=0 $

Le radici sono : $ lambda _1=2 , lambda _2=-1,lamda_3=1$

Mi calcolo gli auto valori per ognuna delle radici :
( Non scrivo tutti passaggi per questioni di spazio )

$ nu _1=(7,1,3) , nu _2=(-1,2,0),nu _3=(1,0,0) $

Finalmente mi posso calcolare la matrice esponenziale e^FT usando la forma canonica di Jordan e la matrice modale .

$ F = ( ( 1 , 1 , 2 ),(0 , -1 , 1 ),( 0 , 0 , 2 ) ) $

Conosciamo gli auto valori di questa matrice e sono reali e distinti . Noti gli auto valori , possiamo andare a determinare la forma canonica di Jordan equivalente ad F : essendo F una matrice di ordine 3 e avendo trovato 3 auto valori reali e distinti , la matrice di Jordan sarà una matrice di ordine 3 avente mini blocchi coincidenti con i tre auto valori .

$ J=( ( 2 , 0 , 0 ),( 0 , -1 , 0 ),( 0 , 0, 1 ) ) $

A questo punto mi calcolo l’esponenziale di questa matrice :

$e^(jt)=( ( e^(2t) , 0 , 0 ),(0 , e^(-t) , 0 ),( 0 ,0 , e^(t) ) )$

A questo punto l’obiettivo e’ l’applicazione della relazione $ e^(Ft) = Me^(jt)*M^(-1)$ , per cui devo individuare la matrice Modale e successivamente la sua inversa .

A questo scopo , dobbiamo trovare gli auto vettori associati ai due auto vettori della matrice F.

L’auto valore associato all’auto valore :

$ (lambda I- F)x=0$

Applicando questa formula , la forma estesa è la seguente :

$2*( ( 1 , 0 , 0 ),( 0 , 1 , 0 ),( 0 , 0 , 1 ) ) -( ( 1 , 1 , 2),( 0 , -1 , 1 ),( 0 , 0 , 2 ) ) *( ( x_1 ),( x_2 ),( x_3 ) ) =( ( 0 ),( 0 ),( 0 ) ) =
| ( x_1 , -x_2-2 , x_3 ),( , 3x_2 ,-x_3 ),( , 0 , ) |$

Lo stesso discorso per $lamda_2=-1$ e risulta :

$| ( -2x_1 , -x_2 , -2x_3 ),( , -x_3 , ),( , -3x_3 , ) |$

Per $ lambda_3 = 1 $ , vale : $| ( -x_2 , , -2x_3 ),( 2x_2 , ,-x_3 ),( , -x_3 , ) | $

A questo punto NON riesco ad andare avanti perché non riesco a calcolarmi gli

Auto vettori associati della matrice modale :


MI aiutate per favore ?
polid
Junior Member
Junior Member
 
Messaggio: 98 di 344
Iscritto il: 04/12/2015, 19:34

Re: [ Teoria dei sistemi e Controlli automatici ]esercizio su calcolo matrice di transizione

Messaggioda Magma » 10/12/2018, 11:59

Scusami, è un po' lungo a faccio fatica a leggere tutto quanto; qual è la matrice di cui di serve cercare autovalori e autovettori?

Se hai una matrice triangolare:


$ F = ( ( 1 , 1 , 2 ),(0 , -1 , 1 ),( 0 , 0 , 2 ) ) $

puoi notare che gli autovalori sono gli elementi della diagonale principale.

Testo nascosto, fai click qui per vederlo
Infatti presa una generica matrice triangolare superiore

$M:= ( ( mu_1 , a , b ),(0 ,mu_2 ,c ),( 0 , 0 , mu_3 ) ) $

ne cerchiamo il polinomio caratteristico

$det(M-lambdaI)=det ( ( mu_1-lambda , a , b ),(0 , mu_2-lambda ,c ),( 0 , 0 , mu_3-lambda ) ) $

applicando l'algoritmo di Sarrus otteniamo:

$(mu_1-lambda)*(mu_2-lambda)*(mu_3-lambda)=0$

da cui si ottiene che gli autovalori sono proprio

$lambda_1=mu_1, qquad lambda_2=mu_2, qquad lambda_3=mu_3$


Per cercare, ad esempio, l'autovettore relativo all'autovalore $lambda=2$, si risolve il seguente sistema lineare omogeneo:

$ (F-2I) barv=bar0$

$hArr = ( ( -1 , 1 , 2 ),(0 , -3 , 1 ),( 0 , 0 , 0 ) )((x),(y),(z))=((0),(0),(0))$


scegliendo come incognita libera $y$, si ottiene

${ ( y in RR ),( z=3y ),(x=y+2z=y+6y=7y):} hArr { ( x=7y ),( y in RR ),(z=3y):}$

Quindi un vettore soluzione del sistema lineare omogeneo è

$ ((x),(y),(z))= ((7y),(y),(3y))=y((7),(1),(3))$

quindi un autovettore per l'autovalore $lambda=2$ è (ad esempio per $y=1$)

$((7),(1),(3))$
Ultima modifica di Magma il 11/12/2018, 10:42, modificato 3 volte in totale.
Magma
Cannot live without
Cannot live without
 
Messaggio: 1355 di 3000
Iscritto il: 03/09/2015, 13:15

Re: [ Teoria dei sistemi e Controlli automatici ]esercizio su calcolo matrice di transizione

Messaggioda polid » 10/12/2018, 14:15

Ok. grazie per aver risposto.

Ci sono quasi , mi serve ancora un piccolo aiuto.

faccio un esempio :

$ ( ( 1 , 0 , 0 ),( 0 , 1 , 0 ),( 0 , 0 , 1 ) ) -( ( 1 , 1 , 2),( 0 , -1 , 1 ),( 0 , 0 , 2 ) ) *( ( x ),( y ),( z ) ) =( ( 0 ),( 0 ),( 0 ) ) $

corrisponde a :

$ ( ( 0 , -1 , -2),( 0 , 2 , -1 ),(0,0,-1)) *( ( x ),( y ),( z ) ) = $

Il risultato sono tre sistemi in tre incognite :


$ { ( -y-2z=0 ),( 2y-z=0 ),( -z=0 ):} $

A questo punto mi perdo . Da questo sistema devo trovare l'autovettore .

Devo utilizzare il metodo della sostituzione ?

Magma , ti sembra corretto un autovettore del tipo .. $ nu _1=(0 ,0 ,0 ) .... $ :?:

(E' un problema di coefficienti ....non capisco quali sono i coefficienti giusti per determinare il vettore .)
polid
Junior Member
Junior Member
 
Messaggio: 99 di 344
Iscritto il: 04/12/2015, 19:34

Re: [ Teoria dei sistemi e Controlli automatici ]esercizio su calcolo matrice di transizione

Messaggioda Magma » 10/12/2018, 14:46

polid ha scritto:
[autovettori per $lambda=1$]
$(1-F)= [( ( 1 , 0 , 0 ),( 0 , 1 , 0 ),( 0 , 0 , 1 ) ) -( ( 1 , 1 , 2),( 0 , -1 , 1 ),( 0 , 0 , 2 ) ) ]( ( x_1 ),( x_2 ),( x_3 ) ) =( ( 0 ),( 0 ),( 0 ) )$

corrisponde a :
$ ( ( 0 , -1 , -2),( 0 , 2 , -1 ),(0,0,-1)) ( ( x),( y ),( z ) ) = ((0),(0),(0))$

A questo punto mi perdo

I metodi per risolvere un sistema lineare sono vari; io, in genere, uso l'algoritmo di Gauss:

Testo nascosto, fai click qui per vederlo
$ ( ( 0 , -1 , -2),( 0 , 2 , -1 ),(0,0,-1))$

$R_1 -> R_1+2R_3$
$R_2 -> R_2+R_3$
$( ( 0 , -1 , 0),( 0 , 2 , 0 ),(0,0,-1)) $

$R_2 -> R_2+2R_1$
$( ( 0 , -1 , 0),( 0 , 0 , 0 ),(0,0,-1)) $


e il sistema di partenza è equivalente a quello avente la matrice ridotta per righe

$ ( ( 0 , -1 ,0),( 0 ,0 , 0 ),(0,0,-1)) ( ( x ),( y ),( z ) ) = ((0),(0),(0))$

le cui soluzioni sono
${ ( x in RR ),( y=0 ),( z=0 ):} hArr ((x),(0),(0))=x((1),(0),(0))$


infatti dal teorema di Kronecker-Rouché-Capelli sappiamo che un sistema in $2$ equazioni (indipendenti) e $3$ incognite, ha $oo^(3-2)=oo^1$ soluzioni.

Infine, un1 candidato autovettore è
$((1),(0),(0))$


polid ha scritto:
Magma , ti sembra corretto un autovettore del tipo .. $ nu _1=(0 ,0 ,0 ) .... $ :?:

Mai! Un autovettore è non nullo per definizione.

Note

  1. Tutti i multipli di $((1),(0),(0))$ sono autovettori relativi al medesimo autovalore!
Ultima modifica di Magma il 10/12/2018, 21:54, modificato 3 volte in totale.
Magma
Cannot live without
Cannot live without
 
Messaggio: 1356 di 3000
Iscritto il: 03/09/2015, 13:15

Re: [ Teoria dei sistemi e Controlli automatici ]esercizio su calcolo matrice di transizione

Messaggioda polid » 10/12/2018, 15:26

Grazie Magma

Devo imparare anche io l'algoritmo di Gauss .

Non deve essere difficile .
polid
Junior Member
Junior Member
 
Messaggio: 100 di 344
Iscritto il: 04/12/2015, 19:34

Re: [ Teoria dei sistemi e Controlli automatici ]esercizio su calcolo matrice di transizione

Messaggioda Magma » 10/12/2018, 18:12

Prego, comunque ho fatto delle correzioni al post precedente. :wink:
Magma
Cannot live without
Cannot live without
 
Messaggio: 1360 di 3000
Iscritto il: 03/09/2015, 13:15


Torna a Ingegneria

Chi c’è in linea

Visitano il forum: Nessuno e 1 ospite