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.
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
Sistemi a Tempo Discreto (Dominio Z)
Differenze con Tempo Continuo
- Abbiamo Equazioni alle Differenze (no differenziali)
- In pratica se una funzione è differenziabile posso studiare f(x) in un piccolissimo intervallo (dx)
- Alle differenze è noto il valore di f(x) solo in punti presi K
- K ∈ Z, 1, 2, 3...
Sistema
- Nel Tempo Continuo e/o Stazionario - Proprio
- Nel Tempo Discreto - Stazionario - Proprio
Risolviamo Le Espressioni
- Iu: am y(k-m) + am-1 y(k-m-1) + ... + ao y(k) = bo(k)+ ... + bu(k)
- x(k+1) = A x(k) + B u(k)
- y(k) ...
Ricordo
.se Statico:
y(k) dipende solo da u(k)
Ia
- x(k+1) = β u(k) + A x(k)
- y(k) = D u(k)
NON TROVO x(k) NELLA RISPOSTA DEL SISTEMA
nella forma compatta
lineare & stazionario (A(x), B(x), C(x), D(x))
Dimostrazione Lagrange
- Prendi l'espressione (siso) del Iuv in Ki:
ay(k) + a y(k-1) + ... + am y(k-m) = bm u(k-m) + ... + b0 u(k)
- Isolo yk e applico la trasformata z
z[y(k)] = z[bmu(k-m) + ... + u(k)] - a1z[y(k)] - ... - amz[y(k-m)]
Y(z) =
- Applico la trasformata elementare f(k-m) = -k/m z-m F(z) = k/m2 -1 z-1
bm/a0 z{-m}[u(z) + ... + u(m-1)z-1 + ...]
b0/a0 z{-1}[z-1Y(z)] + ...
- Isolo e raggruppo Y(z) e trovo
yc(z) = Yiet(z) + Yforza(z) = ...
Hanno lo stesso denominatore
3) Trovo sostituendo k = k - k0
xf(k+1) = Ak-k0Bu(k0) + Ak-k0Bu(k0 + 1) + ... + ABu(k0 + (k-k0 - 1)) + Bu(k)
4) Ho trovato
xf(k+1) = Ak-k0Bu(k0) + Ak-k0Bu(k0 + 1) + ... + Bu(k1)
Ho trovato la tesi => Dimostrato
Stessa dimostrazione dello xe(k)
Poichè se sostituendo k = k + 2
xf(k) = A(k-k0)-2Bu(k0) + A(k-k0)-2Bu(k0 + 2) + ...
Ovvero
xf(k) = Σμ=k0k-2 A(k-μ)-1Bu(μ)
Stabilità Interna
Studio le variabili interne del sistema Σ (X(t))
Definizioni
Stato d'equilibrio: uno stato può essere definito "stato d'equilibrio" Xeq se vale la relazione:
X(t) = Xeq ∀ t ≥ t0 ⇒ X(t) ≡ Xeq
Vale solo per la risposta libera Xe (t) quindi
X• (t) = A Xe (t) + B u(t)
In pratica dopo t0 X sarà sempre fermo sul suo punto di equilibrio
Formule
- Dato che studiamo solo risposta libera u(t) = 0
- Esprimiamo nell condizioni X(t) ≡ Xeq
- A dimostra di Xeq può dipendere da t
- Perché non dipende da t
- Osserva X(t) - Xeq [0] = -1, se t ≥ 0
- Alcunché sostituendo X• (t) = X• eq = 0 y X(t) = Xeq
A Xeq = 0
GAS (Globlalmente Asintoticamente Stabile)
- Stabile → ∀ ε>0, ∃ δ(t0, ε) : ||X(t0) - Xe|| < δ(t0, ε) → ||X(t) - Xe|| < ε
- Attrattiva → ∀ X(t0) ∃ε ∀xe ∈ X → lim (t → ∞) ||X(t) - xe|| = 0
Criterio di stabilità
Ipotizzando di conoscere un certo polinomio caratteristico p(λ) da:
det (λI - A) = p(λ) → ∀posto = 0 fornisce 1 real 2 λ complessi e coniugati
Tempo Continuo
- Asintoticamente Stabile: Reλk < 0 ∀λ
- Semplicemente stabile: Reλe = 0
- Instabile: Reλ1 > 0 V = 1
Tempo Discreto
- Asintoticamente Stabile: ||xe|| < 1
- Semplicemente stabile: |λe| = 1 V = 1
- Instabile: ||λ1|| > 1 V > 1
CRITERIO di Riuth
dote il generico polinomio p(l)=amlm+am-1lm-1+...+a1l+a2.
è possibile costruire la tabella di Routh e si costruisca la tabella (T.d.Routh)
- m 3 3 am am-2
- m 2 2 am-1 am-3
- m 1 1 d1
- m 0 d2
- 2 1 d3
- 1 0 d4
Le prime due le prendi da p(l) ha m+1 righe (pardi vado fino a 0)
Si trova come legge
Con l che va da 2, 3, 4, ...
Dove:
b, c, d, e, ... = -1 det primo termine colonna prec.
es:
bm-4 = -1/am-1
det [ am am-2 ]
le due righe dello colonne successce qualla deolo det quella studiata
det am-6 am-7
det am am-2
det am-2 bm-2
Criterio di Jury
b ha lo stesso segno di quello di Ruth in K.
- Condizione necessaria e sufficiente affinché |λi| < 1
Ipotesi:
- dato il polinomio p(λ) = amλm + am-1λm-1 + ... + a2λ2 + a1λ + a0
Se m = 2
k = a2λ2 + a1 λ + a0
Ipotesi
- devono valere contemporaneamente:
- am = 1
- p(1) > 0
Es. 3λ2 - 2λ + 1 con λ = 1
3 - 2 + 1 = 2 > 0
- am = 1
sostituisco λ = -1
(-1)m p(-1) > 0
|am| > |a0|
Se m > 2
Ipotesi
- devono valere le 3 condizioni precedenti
- deve essere soddisfatte m-2 disuguaglianze tra il primo elemento della riga 1 e l'ultimo elemento
Tabella di Jury
b0
am
a0
- le m-1 coppie di righe (r e b c ... come nel criterio di Routh)
Sistema (lineale e stazionario) Controllabile
se:
- X0 (f*) fissato
- ∃ X* ∈ ℰt0, +∞ (sistema di termini finito t*
- ∃ U*(t) definita fra t0, t*
e possibile o mantenere del stato iniziale X(t0)
trasporre lo stato del sistema a X* attraverso l’ingresso U*(t)
Sotto-spazio di Non Raggiungibilita XNR
- e il sotto-spazio dello spazio di stato X
- Tutto cio che non e raggiungibile
e il complementare dell’insieme X∞
max X∞
tstart toc XNR
Sotto-spazio di Non Controllibilità XNC
- e il sotto-spazio dello spazio di stato X
- Tutto cio che non e controllabile
e il complementare di XC
max XC
(sotto-spazio di controllabilità) XNC
Disegno
(per facilità assuniamo X(t0) di dimensione 2x1 X(t) = [ ])
- In pratica se t = t0
- x1 = y ⇒ X(t) ∈ &overline;X(t*)
- X*(t0) sono i punti o chi combacia X(y,z)
fissato
tutti possibili u