Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
vuoi
o PayPal
tutte le volte che vuoi
Notazioni
Sia T l'insieme dei tempi:
- T = ℝ tempo continuo
- T = ℤ tempo discreto
Notazione omificata:
Δt(t) x(t) t ∈ ℝ
xk t ∈ ℤ
x ∈ ℝn, u ∈ ℝp
SISTEMI:
Δx(t) = Ax(t) + Bu(t)
x (t+1) = Ax(t) + Bu(t)
Differenze tra tempo continuo e discreto
- Un sistema a tempo discreto, al contrario di quello continuo, a risposta libera partendo da un valore diverso da zero, potrebbe andare a zero in tempo finito
- A tempo discreto posso avere oscillazioni senza autovalori immaginari, cosa che non è possibile a tempo continuo
Raggiungibilità
Sia x(t) = φ(t, x0, u(·)) e x(0) = 0 2i può sempre scomporre come
φ(t, 0, u(·)) + φ(t, x0, 0)
RISPOSTA LIBERA
RISPOSTA FORZATA
RISP. FORZATA
Xe(t)
RISPOSTA FORZATA
Xp(t)
abbandonato con adatto comportamento
Def
Uno stato x ∈ ℝn è RAGGIUNGIBILE se ∃ T t ∈ T, x(0) = x
se partendo da condizioni iniziali nulle riesco a raggiungere x0
Se sistema è RAGGIUNGIBILE ∀ x ∈ ℝn è RAGGIUNGIBILE
Def
Uno stato x ∈ ℝn è CONTROLLABILE a(0) se ∃ T t ∈ T, x(t; u) = 0
5.1.2
è controllabile e applicando un controllo, riesco a farlo tornare a 0
Il sistema è CONTROLABILE se ∀x ∈ ℝn, x è CONTROLLABILE
Proprietà
a partire da x
- x ∈ ℝm raggiungibile ⇔ ∃φ | φ(t1,x0) = x
- x ∈ ℝm controllabile ⇔ ∀x̄ : x̄ = φ(t1,x0) per t > 0 ⇒ anche x̄ è controllabile
Sottoinsieme del ℝm
Xr: insieme degli stati raggiungibili del sistema S
Xc: controllabili
- Xc = X̄sc
Proprietà
Xr è sottoinsieme di ℝm
- ∀xa, xb ∈ Xr ∃ λ ∈ [0,1] tale che λxa+βxb ∈ Xr
- ∀t ∃ φ(tb,0,ub(t)): esendo xa raggiungibile
- XS = φ(tb,0,ub(t))
2) ker (P!') è l'insieme dei vettori x : P!'x = 0 , cioè xP0 , x[j⋯i] Am-1B = 0
ovvero xTAKB = 0 ∀ i = 0, ⋯, m .
Il lemma sarà vero se la (16) è vera ⇒ vale la (17)
Se x ∈ Rm soddisfa la (16), derivandola un numero arbitrario di volte rispetto al tempo
si ricava xTAi è AiB = 0 ∀t ∈ [0, ti] i = 0, ___, ________ che per i = 0 implica la (14)
Se x ∈ Rm soddisfa la (17) per C.H.: Δm' = ∑ m-1 A'' i > 0, perciò se per un qualunque j : j⊥:
vettore soddisfa la xTAit'B = 0 , per moltiplicando la (19) per xi e postmoltiplicando B si
ricava xTAmt'B = 0
(17) x… 2 A… Bi = 0
(17) + C.H sette studik A ∀x > m ⇒ xm B = 0 ∀ k > 0
16) ε ∞ Ai = ∑ AskB
(17) xT AkB = x∑ σ… B = ∑ xA
(17) + C.H ⇒ (16) ⇒ x ∈ ker(P) ⇒ xi ε ker(C(t))
P̂=[ Br Brr ] =
Arr Br [.] 0 0 0 ] 0Quando cambiamo base
La forma di Kalman non è unica per la scelta di T-1
Ma alcune cose non cambiano cioè:
det( λΙ- ) = det( λΙ-rr ) det( λΙ-uu ) = Paut(λ) non dip. dalla base scelta
6 (Αrr) autovalor sotto sys.[facilmente] scindibile
6 (Πilr) non facile
Sottosistema
La parte non cagg. influisce nel sist. cogg. ingibile
xr xr ∑uMATRICI DI TRASFER. APPELLO SCOLTO A QUELLO
Matrice di trasferimento (Su non contribuisce)
W(λ)= B·( λΙ--1 )B + D = ( λΙ-rr)-1 x Br 0Algoritmo che c'è solo sull'iscritto
Assegniamo autovalori
Δx(t) = Ax(t) + βu(t)
u(t) = Fx(t) + γ(t)
Lavoriamo su σ(A+BF)
A serve per cambiare gli autovalori (i contorni)
G aspettiamo di non cambiare gli autovalori di Aut
θ(Au) ⊂ θ(A+BF) ∀F ∈ Rm×n
λ = TλT
μ = Tμ
λ = Tx
x = Tx
F = Fr - Fr1
σ(A+BF) = θ(A+BF)
Con thm C
θ(TΔT-1 + TBFT-1) = θ(T(A+BF)T-1)
Costruiamo F = [Fr | Fu]
A + BF = [Arr Aru] + Bf [Frr Fru] = [Arr+BFrr Aru+BrFu]Aru 0 0 | 0 Aru Bru 0 0 0 0
Disegno omino θ(Rru) ⊂ θ(A+BF)
θ(A+BF) = θ(Rru) ⋄ θ- (Arr+BrFru)
Regolarità
Se s ∈ regg ⇒ Sc ∉ regg cioè se (A,B) regg ⇒ ∃F ∈ Rm×n (A+BF,B) ∈ regg
ƴ(t) = -Fx(t) + ω(t)
- w
- N
- u
È vero che rank[A+BF - λI ; B] = m ∀λ ∈ σ Sì
[A+BF - λI ; B] = [A+BF - λI - BF] [B∞] = |λ I - λ I| = λ I - λ I
-F Ip
Aem = a0ea
Aen = en - a0e1
Def.A1 = Amn-2B3 onAmn-2B2 ... am,m-1 An,m-1B
As = esAm-n+1 es
AeB3
Eqs.:
EB = θ(ABtqnmB)
1 omn : ok on-1
0
a1 o3
0 1
∅ on on-m-2
1 am
0 1
Proiezione F = FT-1 − F.Ft
F = fm n-2
A+BF = an − fm ... a-p0
1 0
0 0
-oo − fo = o
fa − ao
Det(An-A) = on − am
1 am-2 on-3
a-1 o0 θ
Ai
π(x) =(x-λ1)(x-λ2)...(x-λm)
λn − a1 λn−1 ... + a2 λo ∅ Δoo
Osservabilità (verificare il fatto usato)
S.1 Δx(t) = A x(t) + B u(t) S.2 y(t) = C x(t) + Du(t)
x ∈ ℝm, u ∈ ℝp, y ∈ ℝp
Problema
Sono noti A, B, C, D, t ∈ [0, tf] t.c. pos. misurando l’uscita etcc...
È possibile ricavare univocamente X(0) o x(tf)?
Def.
S osservabile se è possibile ricavare x(0)
Proprietà
Osservabile e determinabile
È possibile calcolare esattamente ye(τ 0 S non è osservabile
Proprietà
x in osservabile → Ψ(t)x = 0, Φ(t)x è inosservabile
Φ(t)x eA
Dim. V2 = 0
Φ(τ)x = 0 → P(t)Φ(t)x = 0
C Aτfx
UNO STATO INOSSERVABILE IN SOLUZIONE LIBERA SARÀ ANCORA INOSSERVABILE
∀ τ ≥ 0 Y(t): Y (t) = ( eAT*) - CAxTr x =
Y(t) = CeA(t-τ)x = C eA(t-tf)
t = 0 (A(x(0) + Bu(t)))
CA*Tr(t)Bu(t) → Bu(t)
∀t ≥ 0 ∫ Bu(τ) e u(τ) dτ ∀ t ≥ 0
x(t) = eA(t-γη)(B u(γ) Bu(γ)) → Y(τ)