10-Realimentação de Estados
10-Realimentação de Estados
10-Realimentação de Estados
e observadores de estados
Sistemas Realimentados
Conteúdo
! Justificativa
! Modelos em variáveis de estados
! Decomposição de funções de transferência
! Relações entre modelos no espaço de
estados e funções de transferência
! Controlabilidade
! Controle via realimentação de estados
! Observadores de estados
! Realimentação integral de estados
Justificativa
! . Equação diferencial
Esse reator possui duas malhas de controle que interagem entre si.
Como controlar usando apenas controladores PID?
Obtendo G(s) do modelo
em variáveis de estados
! A FT é obtida aplicando Laplace
Simulando os dois modelos
no Matlab
Simulando os dois modelos
no Matlab
! Simulação de G(s):
! Para entradas degrau,
rampa, etc:
! Y=step(g);
! Y=impulse(g)
! Y=lsim(g,u,t);
Simulando os dois modelos
no Matlab
! Simulação do modelo em
variáveis de estado
! O mesmo de G(s), mas
também pode incluir
estados iniciais.
! y=lsim(s1,u,t,[1;2]); Resposta para
estados iniciais
! [y,t,x]=lsim(s1);
Obtendo o modelo em variáveis
de estado a partir de G(s)
! Rearranjando
Decomposição direta
Decomposição direta
! Escolhendo os estados
! Vem
Decomposição direta
Na forma matricial,
Forma
canônica
controlável
Decomposição direta
! Saída,
Diagrama completo
Observações
! .
Decomposição paralela
! .
Decomposição em cascata
! .
Decomposições
! .
Autovetores
O sistema
é controlável se
! Seja
Número de entradas
! .
Forma canônica controlável
• Estabilidade
• Erro em regime
• Sobreelevação
• Tempo de estabelecimento
• IAE
Y(s) ωn2
= M(s) = 2
R(s) s + (2ζωn + K2)s + K1
Aplicação ao protótipo de
segunda ordem
A escolha da sobreelevação UP e do tempo de estabelecimento
Ts permitem definir ζ̄, ω̄n e com isso K1 e K2.
a2 4
ζ̄ = a = log(UP/100) ω̄n =
π2 + a2 Tsζ̄
Y(s) ωn2
= M(s) = 2
R(s) s + (2ζωn + K2)s + K1)
Aplicação ao protótipo de
segunda ordem
Y(s) ωn2
= M(s) = 2
R(s) s + (2ζωn + K2)s + K1)
! Seja o sistema
x· = Ax + Bu
y = Cx + Du
! Equações,
x· = Ax + B(r − Kx)
y = Cx + D(r − Kx)
x· = (A − BK )x + Br
y = (C − DK )x + Dr
Sistema em malha fechada
! .
Exemplo
Ou
Alocação de polos no
Matlab
c=[1 0]
M=ss(a-b*k, b,c,0) Sistema em malha fechada
! Limitação da rotina place: os polos não podem ser
repetidos
FT de MA e MF
2 ! Apenas o
Gma = 2 denominador da FT
s +s−2 muda após a
realimentação de
estados, pois os polos
são alterados pela
realimentação.
2
Gmf = 2 ! O numerador somente
s + 3s + 2 muda se um polo for
alocado para a
mesma posição de um
zero.
Vantagens e desvantagens da
realimentação de estados
! A principal vantagem é a possibilidade de escolher
quaisquer polos de malha fechada desejados.
! No caso abaixo, um controlador PID não permite
escolher todos os coeficientes do denominador!
!
Vantagens e desvantagens da
realimentação de estados
! Uma desvantagem é que o erro em regime
não é considerado no projeto. Além disto, o ganho
da FT de MF fica menor quando se esocolhe polos
que dêem respostas mais rápidas.
Vantagens e desvantagens da
realimentação de estados
! smf=ss(a-b*K,b,c,0)
! gmf=tf(smf)
! step(g,gmf)
Vantagens e desvantagens da
realimentação de estados
0.8
0.4
Malha aberta: g
0.2
0
0 1 2 3 4 5 6 7 8 9 10
! smf=ss(a-b*K,b,c,0) Amplitude
0.04
0.03
! gmf=tf(smf) 0.02
Malha fechada: gmf
0.01
0
0 1 2 3 4 5 6 7 8 9 10
Time (seconds)
! step(g,gmf)
Vantagens e desvantagens da
realimentação de estados
M=ss(a-b*K,b,c,0)
[y,t,x]=step(M);
u=1-K*(x’);
plot(t,u);
Controle via realimentação integral
de estados
Tipo 1
polo 1/s
ki
FT de Malha aberta: G(s)
s
Os polos de G(s) foram alocados no SPE pela matriz K
Realimentação integral de estados
! Em MF,
Realimentação integral de estados
! Sendo
ζ = 0.707
4
ωn = ≥ 2.8
Ts * ζ
Sistema aumentado:
Exemplo:
[y,t,x]=step(M);
Polos de MF em {-2-2j, -2+2j, -10}
plot(t,x);
1.4
x1
1.2 x2
x3
1
0.8
0.6
0.4
0.2
-0.2
0 0.5 1 1.5 2 2.5 3 3.5
Tempo(s)
Exemplo
[y,t,x]=step(M); Polos de MF em {-2-2j, -2+2j, -1}
plot(t,x); Escolheu-se o terceiro polo em -1,
mais lento que os outros 2.
1.5
x1
x2
x3
0.5
0
0 1 2 3 4 5 6 7
Tempo(s)
Exemplo
Sinal de controle u(t)
0
s=-10
s=-1
-1
-3
-5
-6
0 0.5 1 1.5 2 2.5 3 3.5
Tempo(s)
Exemplo
k3
k1, k2
Realimentação com
observadores
! .
Planta real
Cópia da planta
(Modelo)
Equação do observador
Como calcular L?
Erro:
! .
L=place(A’,C',polos)
L=L’ (transposta)
x· = Ax + Bu
u(t) = − Kx(t)
∞
∫0
J= (x*Qx + u*RU )dt Controle ótimo!
K=lqr(A,B,Q,R)
m1=ss(a-b*k1,b,[1 0],0)
m2=ss(a-b*k2,b,[1 0],0)
Controle ótimo
0.5
m1
0.45
m2
0.4
0.35
Saídas para as duas escolhas
0.3
0.25
0.2
∫0
0.1
J= (x*Qx + u*RU )dt
0.05
0
0 1 2 3 4 5 6 7 8 9 10
Tempo(s)
Controle ótimo
Sinal de controle u(t)
0.5
m1
0.45
m2
0.4
0.35
0.3
Respostas mais rápidas requerem
0.25 sinal de controle de maior.
Um maior valor para R em relação a
0.2
Q reduz o sinal de controle, e deixa o
0.15 sistema mais lento.
∞
∫0
0.1
J= (x*Qx + u*RU )dt
0.05
0
0 1 2 3 4 5 6 7 8 9 10
Tempo(s)
Comentários finais sobre controle
moderno: