Controle e Servomecanismos em Engenharia
Controle e Servomecanismos em Engenharia
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
CONTROLE E SERVOMECANISMOS
Resp. Transitória Engenharia da Computação
Sist.
Realimentados
2024-2
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Sistema de controle
Introdução
Aula 01 Um conjunto de dispositivos de forma a
Aula 02
manter uma variável fı́sica dentro de
Modelagem
especificações de interesse.
Resp. Transitória
Sist. Aplicações
Realimentados
Praticamente toda a área tecnológica envolve
Lugar das Raı́zes
algum sistema de controle: fontes de
Resp. em
Frequência alimentação, equipamentos eletrônicos,
Controle PID projetos aeroespaciais, sistemas industriais para
Controle Moderno controle de pressão, temperatura, vazão, nı́vel,
etc.
Objetivos
Em geral, garantir o valor de alguma grandeza
fı́sica, na presença de alguma restrição ou
perturbação, com custo mı́nimo.
Uma Breve História da Teoria de Controle
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
1763-1775 James Watt desenvolve a máquina a vapor.
Introdução
Aula 01
c.1890 Aleksandr Lyapunov desenvolve a teoria da
Aula 02
estabilidade que leva seu nome.
Modelagem
Resp. Transitória
1932 Harry Nyquist desenvolve um procedimento para
Sist.
determinação de estabilidade.
Realimentados
1934 Harold L. Hazen discute o projeto de
Lugar das Raı́zes
servomecanismos a relé capazes de seguir uma
Resp. em
Frequência referência.
Controle PID
1938 Hendrik W. Bode desenvolve a análise pelo
Controle Moderno diagrama que leva seu nome.
1948 Walter Evans desenvolve o método do lugar das
raı́zes.
1960- Uso de métodos no domı́nio do tempo e variáveis
de estado.
1980- Desenvolvimento do controle robusto e H∞ .
A Máquina a Vapor de Watt
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Modelagem
Resp. Transitória
Sist.
Realimentados
Resp. em
Frequência
Controle PID
Controle Moderno
Terminologia
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Aula 01 d(t)
Aula 02
Controle Moderno
Sinais envolvidos no sistema: ▶ Erro: ε(t)
▶ Referência: r(t) ▶ Erro compensado: ũ(t)
▶ Saı́da: y(t) ▶ Perturbação: d(t)
▶ Medido (da saı́da): z(t) ▶ Atuação : u(t)
Um Sistema de Controle: Forno Elétrico
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Modelagem
Resp. Transitória
Sist.
Realimentados
Resp. em
Frequência
Controle PID
Controle Moderno
Exercı́cio
No esquema acima, quem é a planta? E o sensor? E o
compensador? E a entrada e a saı́da?
Desafios em um Projeto de Sistema de Controle
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Controle Moderno
algum outro parâmetro. Ex.: Overshoot e
tempo de resposta.
Digitalização do controle Ao envolver uma CPU para o
projeto do compensador, diversas outras
considerações para o projeto devem ser feitas.
Ex.: Erro de quantização.
Malha Aberta versus Malha Fechada
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02 Sistema em malha fechada Ou realimentado. O controle é
Modelagem
determinado de acordo com o erro presente na
Resp. Transitória
saı́da.
Sist.
Realimentados Sistema em malha aberta Não há leitura da saı́da para
Lugar das Raı́zes determinação da ação de controle.
Resp. em
Frequência
Importante!
Controle PID
Controle Moderno
Um sistema realimentado tende a ser menos sensı́vel a
perturbações externas. Por outro lado, um sistema em
malha aberta tende a ser mais simples e barato em sua
implementação.
Funções Complexas
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
▶ Uma variável complexa, s, é aquela que se escreve como
Aula 02
Modelagem s = σ + jω (1.1)
Resp. Transitória
Resp. em
Frequência F (s) = Fx (s) + jFy (s) (1.2)
Controle PID
Controle Moderno
Observação
A unidade imaginária, j, tem a seguinte propriedade:
j 2 = −1 (1.3)
Magnitude, Ângulo e Conjugado
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
▶ A magnitude do número complexo s é
Introdução
Aula 01 p
Aula 02
|s| = σ2 + ω2 (1.4)
Modelagem
Controle PID
▶ O conjugado do número complexo s é
Controle Moderno
s∗ = σ − jω (1.6)
Observação
As definições de magnitude, ângulo e conjugado se aplicam
de forma análoga a uma função complexa F (s).
Relação de Euler
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
▶ A partir das expansões de Taylor das funções seno e
Modelagem cosseno, pode-se verificar a relação de Euler:
Resp. Transitória
Lugar das Raı́zes ▶ Duas relações importantes podem, também, ser obtidas:
Resp. em
Frequência
Controle PID
ejθ + e−jθ
Controle Moderno cos θ = (1.8a)
2
ejθ − e−jθ
sen θ = (1.8b)
j2
Formas Cartesiana e Polar
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Im(s)
Se |s| = 1, então
Aula 01 j
Aula 02 s
s = cos φ +j sen φ = ejφ
ω
Modelagem | {z } | {z }
Resp. Transitória σ ω
Sist.
Realimentados φ Caso |s| =
̸ 1, basta fazer
−1 σ 1 Re(s)
Lugar das Raı́zes
s
Resp. em s = |s|
Frequência |s|
Controle PID
|{z}
Unitário!
Controle Moderno −j
Observação
Todo número complexo, s ̸= 0 pode ser escrito como
s = σ + jω = |s|ejφ (1.9)
Funções Analı́ticas
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução Definição
Uma função complexa, F , é dita analı́tica em Ω ⊂ C se, ∀s ∈ Ω,
Aula 01
Aula 02
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Definição
Modelagem
Seja F uma função complexa, p, z ∈ C.
Resp. Transitória i Os pontos de C, onde F não é analı́tica são ditos singulares.
Sist.
Realimentados ii Os pontos singulares onde F ou alguma de suas derivadas
Lugar das Raı́zes tende ao infinito são ditos polos.
Resp. em
Frequência iii Se p é um polo de F , tal que |(s − p)n F (s)| < ∞ para
Controle PID
algum n ∈ N. Então, p é dito um polo de ordem n.
Controle Moderno iv Se F (z) = 0, então z é dito um zero de F .
F (s)
v Se z é zero de F , tal que lim = 0 para
so →z(s − so )n−1 s=z
algum n ∈ N. Então, z é dito um zero de ordem n.
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Controle Moderno
▶ Multiplicando e dividindo pelo conjugado, chega-se a:
σ+1 −ω
F (σ + jω) = +j
(σ + 1)2 + ω 2 (σ + 1)2 + ω 2
| {z } | {z }
Fx Fy
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
▶ Assim, as derivadas parciais são
Introdução
Aula 01
Aula 02 ∂Fx ω 2 − (σ + 1)2 ∂Fy
= =
Modelagem ∂σ ((σ + 1)2 + ω 2 )2 ∂ω
Resp. Transitória
Controle Moderno
▶ Assim, pode-se calcular F ′ (s):
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Modelagem
Definição (Transformada de Laplace)
Resp. Transitória Seja f : R → R uma função. A Transformada de
Sist. Laplace (TL) de f , L[f (t)] é dada por
Realimentados
Observação
Sob quais condições uma função de variável real possui TL?
Existência da Transformada de Laplace
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01 Teorema (Existência da TL)
Aula 02
Modelagem
Seja f : R+ → R uma função. A integral (1.11) converge
Resp. Transitória
se:
Sist. i) f for seccionalmente contı́nua em seu domı́nio, e;
Realimentados
Controle Moderno
Observação
Caso as condições anteriores se verifiquem, então:
▶ f é dita de ordem exponencial;
▶ σ é dita abscissa de convergência.
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Sist. Note que esta expressão tenderá a zero se, e somente se:
Realimentados
Controle PID
Assim, a TL dessa função pode ser calculada e vale:
Controle Moderno Z ∞
L[f (t)] = e−αt e−st d t
0
Z ∞
1
= e−(s+α)t d t =
0 s+α
Atraso no Tempo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Sejam L[f (t)] = fˆ(s) e α ≥ 0. A TL de f (t − α) é dada por:
Aula 01
Aula 02
Z ∞
Modelagem L[f (t − α)] = f (t − α)e−st d t
0
Resp. Transitória
Portanto:
L[f (t − α)] = e−αs fˆ(s) (1.12)
Funções Degrau, Pulso e Impulso I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ A função degrau unitário e sua TL são, respectivamente
Aula 01
(
1 se t ≥ 0
Aula 02
Sist. Z ∞ Z ∞
Realimentados 1
L[1̃(t)] = 1̃(t)e−st d t = e−st d t = (1.14)
Lugar das Raı́zes 0 0 s
Resp. em
Frequência ▶ Dado, to > 0, define-se a função pulso como
Controle PID
1 se 0 ≤ t ≤ t
Controle Moderno o
pto (t) = to (1.15)
0 caso contrário
1
pto (t) = [1̃(t) − 1̃(t − to )] (1.16)
to
Funções Degrau, Pulso e Impulso II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
1 1 − e−to s
Modelagem
p̂to (s) = L[1̃(t)] − e−to s L[1̃(t)] = (1.17)
Resp. Transitória to to s
Sist.
Realimentados ▶ A função impulso ou delta de Dirac é definida com to → 0
Lugar das Raı́zes em (1.15): (
Resp. em ∞ se t = 0
Frequência δ(t) = (1.18)
Controle PID 0 caso contrário
Controle Moderno
▶ Fazendo o mesmo limite em (1.17):
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Teorema
Modelagem
Seja L[f (t)] = fˆ(s). Então:
Resp. Transitória
d f (t)
Sist.
Realimentados
L = sfˆ(s) − f (0) (1.20)
dt
Lugar das Raı́zes
Resp. em
Frequência
Observação
Controle PID
Controle Moderno
Para uma diferenciação de n-ésima ordem:
n n−1 <n−1−i>
d f (t)
= sn fˆ(s) − (1.21)
X
L si f (0)
d tn i=0
Teoremas do Valor Inicial e do Valor Final
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Teorema (Valor Final)
Modelagem Seja L[f (t)] = fˆ(s), com todos os seus polos no semi-plano
Resp. Transitória esquerdo e, no máximo, um polo na origem, então:
Sist.
Realimentados
lim f (t) = lim sfˆ(s) (1.22)
Lugar das Raı́zes t→∞ s→0
Resp. em
Frequência
Prof. Victor L.
Yoshimura
Introdução
Definição (Convolução)
Aula 01
Aula 02
Sejam f e g funções reais. Sua convolução é definida por
Modelagem Z t
Resp. Transitória f (t) ∗ g(t) = f (t − τ )g(τ ) d τ (1.24)
Sist. 0
Realimentados
Resp. em Observação
Frequência
Controle PID
A convolução é comutativa. De fato, basta fazer ξ = t − τ
Controle Moderno
em (1.24) para verificar tal proposição.
Teorema (Convolução)
Se as TLs de f e g existem, então:
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Definição (Transformada Inversa de Laplace)
Modelagem Sejam fˆ : C → C uma função complexa e c ∈ R, maior que
Resp. Transitória todas as partes reais das singularidades de fˆ. A
Sist. Transformada Inversa de Laplace (TIL) de fˆ(s), L−1 [fˆ(s)], é
Realimentados
Prof. Victor L.
Yoshimura
Controle Moderno
▶ Se n ≥ m, a função racional é dita própria;
▶ Se n > m, a função racional é dita estritamente própria;
▶ Para funções racionais estritamente próprias, a
Expansão em Frações Parciais (EFP) segue casos
distintos.
1º caso: Polos reais e distintos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
ˆ como
Neste caso, reescreve-se d(s)
Introdução ˆ = (s − p1 )(s − p2 ) . . . (s − pn )
d(s)
Aula 01
Aula 02
e (1.27) reescreve-se na forma
Modelagem
n
Resp. Transitória
X ri
fˆ(s) =
Sist. s − pi
Realimentados i=1
Prof. Victor L.
Yoshimura s+3
Calculemos a TIL de fˆ(s) = . Para tanto,
Introdução
s2 + 3s + 2
Aula 01 observe que:
Aula 02
Modelagem s+3 r1 r2
Resp. Transitória
fˆ(s) = = +
(s + 1)(s + 2) s+1 s+2
Sist.
Realimentados
Os resı́duos são calculados com (1.28):
Lugar das Raı́zes
Resp. em s+3
Frequência r1 = (s + 1) =2
Controle PID (s + 1)(s + 2) s=−1
Controle Moderno
s+3
r2 = (s + 2) = −1
(s + 1)(s + 2) s=−2
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Considere um par conjugado de polos da forma
Aula 02 p1,2 = α ± jω. Assim:
Modelagem
Resp. em ds + f ω s−α
Frequência fˆ(s) = = k1 + k2
Controle PID s2 + bs + c (s − α)2 + ω 2 (s − α)2 + ω 2
Controle Moderno
com k1 e k2 a determinar. Aplicando a TIL, obtém-se o
resultado desejado.
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02 2s + 12
Calculemos a TIL de fˆ(s) = . Seus polos são
Modelagem s2
+ 2s + 5
Resp. Transitória −1 ± j2. Reescrevendo a função, tem-se:
Sist.
Realimentados 2 s+1
fˆ(s) = k1 2
+ k2
Lugar das Raı́zes (s + 1) + 4 (s + 1)2 + 4
Resp. em
Frequência
Donde, comparando as formas, conclui-se que k1 = 5 e
Controle PID
k2 = 2. Usando (1.30), tem-se:
Controle Moderno
Prof. Victor L.
Yoshimura
Neste caso, assume-se o denominador de fˆ(s) sob a forma
Introdução
Aula 01
ˆ = (s − p)n
d(s)
Aula 02
Modelagem
Então, pode-se reescrever a função sob a forma
Resp. Transitória
n n
Sist.
X ri X
Realimentados fˆ(s) = ⇒ (s − p)n fˆ(s) = ri (s − p)n−i
(s − p)i
Lugar das Raı́zes i=1 i=1
Resp. em Pode-se verificar que, com o resultado acima, os resı́duos podem ser
Frequência
calculados com:
Controle PID
1 dn−i
Controle Moderno ri = [(s − p)n fˆ(s)]s=p (1.31)
(n − i)! d sn−i
Usando a tabela de TLs, tem-se que:
−1 ri ri
L = ti−1 ept (1.32)
(s − p)i (i − 1)!
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução s+1
Aula 01 Determine-se a TIL de fˆ(s) = . Calculando os
Aula 02 (s + 3)2
Modelagem resı́duos:
Resp. Transitória
Sist.
1 d2−2
Realimentados
r2 = [s + 1]s=−3 = −2
(2 − 2)! d s2−2
Lugar das Raı́zes
Resp. em
Frequência
1 d2−1
r1 = [s + 1]s=−3 = 1
Controle PID (2 − 1)! d s2−1
Controle Moderno
Assim, usando (1.32), tem-se
s+1
L−1 = (1 − 2t)e−3t
(s + 3)2
Solução de Equações Diferenciais Ordinárias
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Um Sistema Linear, Contı́nuo e Invariante no
Tempo (SLIT-C) é descrito pela EDO:
Introdução
Aula 01 n m
⟨i⟩ ⟨j⟩
(1.33)
Aula 02
X X
Modelagem
αi y = βj u
i=0 j=0
Resp. Transitória
Sist.
Realimentados
onde u é a entrada e y, a saı́da (e n, a ordem do sistema).
Lugar das Raı́zes
A Função Transferência (FT) desse SLIT-C é definida por:
Resp. em
ŷ(s) L[y(t)]
Frequência
ĝ(s) = = (1.34)
Controle PID û(s) L[u(t)] y(0)=ẏ(0)=...= y
⟨n−1⟩
(0)=0
Controle Moderno
Combinando (1.33) e (1.34):
m
X
βj sj
j=0
ĝ(s) = n (1.35)
X
i
αi s
i=0
Funções Transferência: Observações
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Observação
Modelagem
Resp. Transitória
A Função Transferência:
Sist. i É uma caracterı́stica do sistema, independendo da
Realimentados
entrada aplicada;
Lugar das Raı́zes
Prof. Victor L.
Yoshimura
Introdução
Aula 01
F
2
A EDO para a posição do satélite é:
Aula 02
Modelagem
Centro J θ̈ = T ⇒ J ÿ = u
Resp. Transitória + de
Massa
Sist. Aplicando a TL sob Condições
Realimentados
J → momento de inércia
FTs: Identificação de Sistemas via Impulso
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aula 01
Aula 02
Sist.
Realimentados
y(t) = L−1 [ĝ(s)] (1.36)
u=δ
Lugar das Raı́zes
Resp. em
Assim, para obter a FT de forma experimental, faz-se:
Frequência
1º Aplique um impulso (delta de Dirac) à entrada;
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Um modelo é um conjunto de equações que descrevem
Aula 03
Aula 04
com um certo grau de precisão um sistema.
Aula 05
Aula 06
Aula 07
Importante!
Resp. Transitória Um compromisso existe entre:
Sist. Simplicidade e Precisão
Realimentados
Controle PID
▶ Sistemas invariantes no tempo: Aqueles cuja estrutura
Controle Moderno não se altera ao longo do tempo;
▶ Os SLIT-Cs são descritos por EDOs lineares com
coeficientes constantes.
Diagramas de Blocos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Diagramas de Blocos (DBs): representações das
Aula 03
Aula 04
interações entre os elementos de um sistema.
Aula 05
Aula 06
▶ Um bloco é a representação para a operação
Aula 07
matemática sobre o sinal de entrada que produz a saı́da.
Resp. Transitória
Resp. em
Frequência Observação
Controle PID
Como as FTs, os DBs fornecem o comportamento dinâmico
Controle Moderno
do sistema, mas nada falam sobre sua estrutura interna.
Diferentes sistemas podem ter o mesmo diagrama de blocos.
Diagramas de Blocos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Também compõem um diagrama de blocos:
Modelagem ▶ Ponto de soma:
Aula 03
Aula 04
Aula 05
û
Aula 06
Aula 07
Resp. Transitória v̂ û + v̂ − ŵ
Σ
Sist.
Realimentados −
Lugar das Raı́zes
ŵ
Resp. em
Frequência
▶ Ponto de ramificação:
Controle PID
Controle Moderno
A(s) A(s)
A(s)
Exemplo: Sistema em Realimentação Unitária
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
r̂(s) ε̂(s) ŷ(s)
Modelagem Σ ĝ(s)
Aula 03
Aula 04 −
Aula 05
Aula 06
Aula 07
Resp. em
Frequência Assim:
Controle PID ŷ(s) = ĝ(s)r̂(s) − ĝ(s)ŷ(s)
Controle Moderno
ŷ(s) ĝ(s)
= (2.1)
r̂(s) 1 + ĝ(s)
Sistema em Malha Fechada
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
r̂(s) ε̂(s) ŷ(s)
Aula 03 Σ ĝ(s)
Aula 04
−
Aula 05
Aula 06
ẑ(s)
Aula 07 ĥ(s)
Resp. Transitória
Sist.
Realimentados ▶ Função Transferência de Malha Aberta (FTMA): É a
Lugar das Raı́zes razão entre o sinal do sensor e o de erro.
Resp. em
Frequência ẑ(s) = ĝ(s)ĥ(s)ε̂(s) ⇒ F T M A(s) = ĝ(s)ĥ(s) (2.2)
Controle PID
Controle Moderno
▶ Função Transferência de Feed Forward (FTFF): É a
razão entre o sinal de saı́da e o de erro.
ŷ(s) = ĝ(s)ε̂(s) ⇒ F T F F (s) = ĝ(s) (2.3)
Função Transferência de Malha Fechada
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
r̂(s) ε̂(s) ŷ(s)
Modelagem Σ ĝ(s)
Aula 03
Aula 04
−
Aula 05
ẑ(s)
Aula 06 ĥ(s)
Aula 07
Resp. Transitória
Resp. em
Frequência
ε̂(s) = r̂(s) − ĥ(s)ŷ(s) e que ŷ(s) = ĝ(s)ε̂(s)
Controle PID
Controle Moderno
Assim:
ŷ(s) ĝ(s)
F T M F (s) = = (2.4)
r̂(s) 1 + ĝ(s)ĥ(s)
Sistemas sob Perturbação
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem ˆ
d(s)
Aula 03
Aula 04
Aula 05 r̂(s) ε̂(s) û(s) ûd (s) ŷ(s)
Aula 06 Σ ĉ(s) Σ ĝ(s)
Aula 07
−
Resp. Transitória
Sist. ẑ(s)
Realimentados
ĥ(s)
Lugar das Raı́zes
Resp. em
Frequência ˆ ≡ 0 e com (2.4):
Como se tem um SLIT-C, fazemos d(s)
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Modelagem ˆ
d(s)
Aula 03
Aula 04
Aula 05 r̂(s) ε̂(s) û(s) ûd (s) ŷ(s)
Aula 06 Σ ĉ(s) Σ ĝ(s)
Aula 07
−
Resp. Transitória
Sist. ẑ(s)
Realimentados
ĥ(s)
Lugar das Raı́zes
Resp. em
Frequência
Fazendo r̂(s) ≡ 0 e, usando (2.4):
Controle PID
Prof. Victor L.
Yoshimura
Sist.
Realimentados ▶ Se |ĉ(s)ĝ(s)ĥ(s)|, |ĉ(s)ĥ(s)| ≫ 1, então a perturbação é
Lugar das Raı́zes suprimida;
Resp. em
Frequência r̂(s)
▶ Neste caso, tem-se ŷ(s) = .
Controle PID ĥ(s)
Controle Moderno Um sistema em malha fechada pode ser projetado para:
▶ rejeitar perturbações;
▶ igualar saı́da e entrada (seguir a referência).
Simplificação de Diagramas de Blocos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03 Observação
Aula 04
Aula 05 Blocos podem ser postos em série somente se a saı́da do
Aula 06
Aula 07 bloco antecedente não for afetada pelo bloco subsequente.
Resp. Transitória
Sist. Observação
Realimentados
Simplificar o diagrama de blocos implica em tornar mais
Lugar das Raı́zes
Resp. em
complexa a FT de cada bloco.
Frequência
Ao simplificar blocos deve-se:
Controle PID
Controle Moderno
▶ Manter o produto das FTs no caminho direto;
▶ Manter o produto das FTs em torno de cada laço.
Exemplo: Carga Afetando Bloco Anterior
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03 u R1 x R2 y
Aula 04
Aula 05
Aula 06
Aula 07
C1 C2
Resp. Transitória
Sist.
Realimentados
Resp. em
Frequência ▶ Note que o sinal x (tensão em C1 ) é afetado pela carga
Controle PID
(R2 e C2 );
Controle Moderno
▶ Assim, não é possı́vel fazer uma associação série de
blocos com FTs para cada par RC.
Operações com Blocos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
I Associação série
Aula 03
Aula 04
Aula 05
û(s) ŷ(s)
Aula 06 ĝ1 (s) ĝ2 (s)
Aula 07
Resp. Transitória
Lugar das Raı́zes ŷ(s) = ĝ1 (s)(ĝ2 (s)û(s)) = (ĝ1 (s)ĝ2 (s))û(s) (2.8)
Resp. em
Frequência
Controle PID
Então este diagrama é equivalente a:
Controle Moderno
û(s) ŷ(s)
ĝ1 (s)ĝ2 (s)
Operações com Blocos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Controle PID
r̂(s) ĝ(s) ŷ(s)
Controle Moderno
1 + ĝ(s)ĥ(s)
Operações com Blocos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
III Deslocamento de bloco a jusante de ponto de soma:
Modelagem
Aula 03
v̂(s)
Aula 04
Aula 05
Aula 06
û(s) − ŷ(s)
Aula 07
ĝ(s) Σ
Resp. Transitória
Sist.
Realimentados
Notando a relação
Lugar das Raı́zes
Resp. em
Frequência ŷ(s) = ĝ(s)û(s) − v̂(s) = ĝ(s)(û(s) − [ĝ(s)]−1 v̂(s)) (2.9)
Controle PID
Prof. Victor L.
Yoshimura
IV Deslocamento de bloco a jusante de ponto de
Introdução
Modelagem
ramificação:
Aula 03
Aula 04
Aula 05
û(s) ŷ(s)
Aula 06
ĝ(s)
Aula 07
Resp. Transitória
ẑ(s)
Sist.
Realimentados
Notando a relação
Lugar das Raı́zes
Resp. em
Frequência ŷ(s) = ẑ(s) = ĝ(s)û(s) (2.10)
Controle PID
Prof. Victor L.
Yoshimura
Sist.
Realimentados
ĥ(s)
Lugar das Raı́zes
Resp. em
Frequência Após usar (2.9):
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Resp. Transitória
−
Sist. r̂ ŷ
Realimentados Σ Σ ĝ1 Σ ĝ2 ĝ3
Lugar das Raı́zes
−
Resp. em
Frequência
Prof. Victor L.
Yoshimura
Movendo ĝ1 a jusante do 3º ponto de soma:
Introdução
Modelagem
ĥ2
Aula 03
ĝ1
Aula 04
Aula 05
Aula 06
Aula 07 −
r̂ ŷ
Resp. Transitória Σ Σ Σ ĝ1 ĝ2 ĝ3
Sist.
Realimentados
−
Lugar das Raı́zes
Resp. em ĥ1
Frequência
Controle PID
Prof. Victor L.
Yoshimura
Introdução
ĥ2
Modelagem
Aula 03
−
Aula 04
r̂ ŷ
Aula 05
Aula 06
Σ Σ ĝ1 Σ ĝ2 ĝ3
Aula 07
−
Resp. Transitória
Sist.
Realimentados ĥ1
Lugar das Raı́zes
Resp. em
Frequência
Controle PID
Controle Moderno
Observação
Prof. Victor L.
Yoshimura
Resp. Transitória
t ≥ to , determina completamente o comportamento do
Sist.
sistema para t ≥ to ;
Realimentados
▶ Variáveis de estado: Cada uma das variáveis que
Lugar das Raı́zes
compõem o estado;
Resp. em
Frequência ▶ Espaço de estado: O espaço Rn , cujos eixos
Controle PID
representam cada uma das variáveis de estado;
Controle Moderno
▶ Equações do espaço de estado: São equações que
relacionam as entradas, u, as saı́das, y, e as variáveis de
estado, x, com as derivadas destas, ẋ.
Espaço de Estado: Introdução
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Em forma matricial, representa-se um SLIT-C no espaço de
Modelagem
Aula 03
estado como:
Aula 04
Aula 05
Aula 06 ẋ = Ax + Bu (2.12a)
Aula 07
Prof. Victor L.
Yoshimura
Resp. Transitória u ẋ R x y
B Σ C Σ
Sist.
Realimentados
Resp. em
Frequência A
Controle PID
Controle Moderno
Observação
As matrizes A, B, C, D serão funções do tempo, caso o
sistema seja variante no tempo.
Exemplo: Sistema Massa-Mola-Amortecedor
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução y
Modelagem u
Aula 03
Aula 04 k b
Aula 05 m
Aula 06
Aula 07
Resp. Transitória
Sist.
Realimentados
Famort = −bẏ
Exemplo: Sistema Massa-Mola-Amortecedor
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
y
Introdução
u
Modelagem
Aula 03 k b
Aula 04 m
Aula 05
Aula 06
Aula 07
Resp. Transitória
Sist.
Realimentados Pela lei de Newton:
Lugar das Raı́zes
b k 1
Resp. em
Frequência
u + Fmola + Famort = mÿ ⇒ ÿ + ẏ + y = u
m m m
Controle PID
Prof. Victor L.
Yoshimura
Introdução y
Modelagem
Aula 03 u
Aula 04
Aula 05 k b
Aula 06
m
Aula 07
Resp. Transitória
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Controle PID
Admitindo que sI − A seja invertı́vel:
Controle Moderno
ŷ(s)
ĝ(s) = = C(sI − A)−1 B + D (2.13)
û(s)
Funções Transferência e Espaço de Estado
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Lembrando que:
Introdução
Modelagem 1
Aula 03 (sI − A)−1 = adj(sI − A)
Aula 04
Aula 05
det(sI − A)
Aula 06
Aula 07
Tem-se:
Resp. Transitória
Sist.
1
Realimentados ĝ(s) = (C adj(sI − A)B + det(sI − A)D)
det(sI − A)
Lugar das Raı́zes
(2.14)
Resp. em
Frequência
Observação
Controle PID
Prof. Victor L.
Yoshimura
Obtenhamos uma representação no espaço de estado para o
Introdução
sistema:
Modelagem
Aula 03
Aula 04 r̂(s) 1
ŷ(s)
1
Aula 05 Σ s s+a
Aula 06
Aula 07 −
compensador planta
Resp. Transitória
Sist. 1
Realimentados s+b
Lugar das Raı́zes
sensor
Resp. em
Frequência
Controle PID
Observação
Controle Moderno
A solução a ser apresentada não é a única! Fica como
exercı́cio encontrar uma forma de mostrar todas as
representações possı́veis no espaço de estado para o
diagrama de blocos apresentado.
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Escolhendo as saı́das dos blocos como variáveis de estado:
Introdução
r̂(s) x̂1 (s) 1
ŷ(s) = x̂2 (s)
Modelagem 1
Aula 03
Σ s s+a
Aula 04 −
Aula 05 compensador planta
Aula 06
Aula 07 x̂3 (s) 1
Resp. Transitória s+b
Sist. sensor
Realimentados
Resp. em
Donde:
Frequência
1
Controle PID x̂1 (s) = [r̂(s) − x̂3 (s)] sx̂1 (s) = −x̂3 (s) + r̂(s)
s
Controle Moderno
1
x̂2 (s) = x̂1 (s) ⇒ sx̂2 (s) = x̂1 (s) − ax̂2 (s)
s+a
1 sx̂3 (s) = x̂2 (s) − bx̂3 (s)
x̂3 (s) = x̂2 (s)
s+b
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Modelagem
Aula 03
r̂(s) 1
x̂1 (s) 1
ŷ(s) = x̂2 (s)
Aula 04
Σ s s+a
Aula 05
−
Aula 06 compensador planta
Aula 07
Resp. Transitória
x̂3 (s) 1
s+b
Sist.
Realimentados
sensor
Lugar das Raı́zes
Resp. em
Frequência
Aplicando a TIL:
Controle PID
0 0 −1 1
Controle Moderno
ẋ = 1 −a 0 x + 0 r
0 1 −b 0
y= 0 1 0 x
Representação de SLIT-Cs no Espaço de Estado
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03
Considere um SLIT-C-SISO, como em (1.33):
Aula 04
Aula 05
n m
Aula 06
X ⟨i⟩ X ⟨j⟩
Aula 07 αi y = βj u
Resp. Transitória i=0 j=0
Sist.
Realimentados
Controle Moderno
▶ O sistema não envolve derivadas da entrada: m = 0;
▶ O sistema envolve derivadas da entrada: 0 < m ≤ n.
Representação de Sistemas: caso m = 0
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. Transitória
Uma escolha conveniente de variáveis de estado é:
Sist.
Realimentados
⟨i−1⟩ d ⟨i−2⟩
Lugar das Raı́zes xi = y = y = ẋi−1 , i = 1, 2, . . . , n (2.15)
Resp. em
dt
Frequência
Controle PID
Aplicando (2.15) a (1.33), com m = 0:
Controle Moderno
n−1
αi βo
(2.16)
X
ẋn = − xi+1 + u
i=0
αn αn
Representação de Sistemas: caso m = 0
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Com (2.15)-(2.16), escreve-se:
Introdução
0 1 0 ... 0 0 0
Modelagem
Aula 03
0 0 1 ... 0 0 0
.. .. .. .. .. .. ..
Aula 04
. . . . . . .
Aula 05
ẋ = x + u
Aula 06
Aula 07 0 0 0 ... 1 0 0
Resp. Transitória 0 0 0 ... 0 1 0
Sist. − ααno − ααn1 − ααn2 ... − ααn−2
n
− ααn−1
n
βo
αn
Realimentados
(2.17a)
Lugar das Raı́zes
(2.17b)
Resp. em y= 1 0 0 ... 0 0 x
Frequência
Controle PID
É possı́vel aplicar a TL em (1.33), originando
Controle Moderno
βo
ĝ(s) = n (2.18)
X
i
αi s
i=0
Representação de Sistemas: caso m = 0
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
u βo
R R R x1 = y
Modelagem
Aula 03 αn Σ
Aula 04
Aula 05
−
Aula 06
αn−1 αn−2 αo
Aula 07
αn αn αn
Resp. Transitória
Sist.
Realimentados
Σ Σ
Lugar das Raı́zes
Resp. em
Frequência ⟨i−1⟩ d ⟨i−2⟩
xi = y = y = ẋi−1 , i = 1, 2, . . . , n
Controle PID dt
Controle Moderno n−1
X αi βo
ẋn = − xi+1 + u
i=0
αn αn
Representação de Sistemas: caso 0 < m ≤ n
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
É possı́vel que (1.33) esteja completa (m = n).
Introdução
n n
X ⟨i⟩ X ⟨j⟩
Modelagem
αi y = βj u
Aula 03 i=0 j=0
Aula 04
Aula 05
Aula 06
Aula 07 ▶ Se não estiver completa, faça alguns (mas, não todos!)
Resp. Transitória βj = 0;
Sist.
Realimentados
▶ Não é possı́vel proceder com (2.15) (solução
Lugar das Raı́zes inconsistente, não-única).
Resp. em
Frequência
Uma abordagem é definir
Controle PID
ẋ1 = −αo y + βo u (2.19)
Controle Moderno
Prof. Victor L.
Yoshimura
Introdução
A partir de x2 , faz-se (analogamente)
Modelagem
Aula 03 ẋ2 = x1 − α1 y + β1 u (2.20)
Aula 04
Aula 05
Aula 06
Aula 07
Novamente, substituindo no resultado anterior
Resp. Transitória
n−2 n−2
X ⟨i⟩ X ⟨j⟩
Sist.
Realimentados
αi+2 y − βj+2 u = x2
Lugar das Raı́zes
i=0 j=0
Resp. em
Frequência Continuando e obtendo ẋn :
Controle PID
Tem-se
αn y − βn u = xn (2.22)
Representação de Sistemas: caso 0 < m ≤ n
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Com (2.19)-(2.22), escreve-se:
Introdução
βo − αo αβnn
0 0 0 ... 0 − ααno
Modelagem
Aula 03
1 0 0 . . . 0 − ααn1 β − α βn
Aula 04 1 1 αn
Aula 05
α2 βn
Aula 06 0 1 0 . . . 0 − αn β2 − α2 αn
Aula 07
ẋ =
.. .. .. . . .. ..
x +
. u
Resp. Transitória ..
. . . . . .
Sist.
Realimentados
αn−2 βn
0 0 0 . . . 0 − αn β − αn−2 αn
n−2
Lugar das Raı́zes
Resp. em 0 0 0 . . . 1 − ααn−1
n βn−1 − αn−1 αβnn
Frequência
Controle PID
(2.23a)
h i h i
Controle Moderno
y = 0 0 0 ... 0 1
αn x+ βn
αn u (2.23b)
Observação
Em SLIT-C-SISOs, haverá transmissão direta se, e somente
se, n = m!
Representação de Sistemas: caso 0 < m ≤ n
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03 É possı́vel aplicar a TL em (1.33), originando
Aula 04
Aula 05
Aula 06 n
X
Aula 07
βj sj
Resp. Transitória
j=0
Sist. ĝ(s) = n (2.24)
Realimentados X
i
αi s
Lugar das Raı́zes
i=0
Resp. em
Frequência
Prof. Victor L.
Yoshimura
Resp. em
Frequência Devemos modelar os três elementos principais dos sistemas
Controle PID mecânicos:
Controle Moderno ▶ Massa;
▶ Mola linear;
▶ Amortecedor viscoso.
Sistemas Mecânicos: Massa
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
x
Modelagem F
Aula 03
Aula 04
Aula 05
m
Aula 06
Aula 07
Resp. Transitória
Sist.
Realimentados
A relação força/posição é dada pela 2ª lei de Newton:
Lugar das Raı́zes
Resp. em
d(mv)
Frequência
F = = mv̇ = mẍ (2.25)
Controle PID dt
Controle Moderno
Ainda, a massa retém energia cinética, segundo a relação:
1
Ec = m(ẋ)2 (2.26)
2
Sistemas Mecânicos: Mola Linear
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03
x
Aula 04
k
Aula 05
Aula 06
F
Aula 07
Mola
linear!
Resp. Transitória
Sist.
Realimentados A relação força/posição é dada pela lei de Hooke:
Lugar das Raı́zes
Resp. em F = kx (2.27)
Frequência
Controle PID
Ainda, a mola retém energia potencial elástica, segundo a
Controle Moderno
relação:
1
Ep = kx2 (2.28)
2
Sistemas Mecânicos: Amortecedor
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
b x
Modelagem
Aula 03
Aula 04
Aula 05
F
Aula 06
Aula 07
Resp. Transitória
Sist.
Realimentados ▶ O amortecedor linear é constituı́do por um êmbolo
Lugar das Raı́zes envolto por óleo: atrito viscoso!
Resp. em
Frequência
▶ A relação força/posição é dada por
Controle PID F = bẋ (2.29)
Controle Moderno
Observação
O amortecedor, ao contrário dos elementos anteriores,
dissipa energia!
Exemplo: Suspensão Automotiva
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
y
Introdução
Forças sobre o bloco
Modelagem m
Aula 03
Aula 04 Fmola = k(u − y)
Aula 05
Aula 06
k b
Famort = b(u̇ − ẏ)
Aula 07 u
Resp. Transitória
Sist.
Realimentados
Resp. em
Frequência Com a 2ª lei de Newton:
Controle PID
Controle Moderno
k(u − y) + b(u̇ − ẏ) = mÿ
Logo:
b k b k
ÿ + ẏ + y = u̇ + u
m m m m
Exemplo: Suspensão Automotiva
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
y
Introdução Ao aplicar a TL:
m
Modelagem
Aula 03
b k
Aula 04
s+
Aula 05
ĝ(s) = m m
Aula 06
Aula 07
k b u b k
s2 + s +
Resp. Transitória m m
Sist.
Realimentados
Controle PID
k k
Controle Moderno
0 −
m x + m u
ẋ =
b b
1 −
h im m
y= 0 1 x
Molas e Amortecedores Torcionais
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03 ▶ Na rotação, tem-se molas com relação entre torque e
Aula 04
Aula 05 ângulo de deformação:
Aula 06
Aula 07
Prof. Victor L.
Yoshimura
θ2
Introdução
k
Modelagem
k1 , b1 k2 , b2
Aula 03
Aula 04
θ1
Aula 05
Aula 06
Aula 07
Resp. Transitória
Sist.
Para cada cilindro, tem-se:
Realimentados
Introdução
y y′ xg = x + l sen θ
m Centro
Modelagem yg = l cos θ
2l
de
Aula 03 θ ×
Aula 04
massa
Aula 05 Rotação de (xg , yg ) em relação
à polia:
Aula 06
u x
Aula 07
M
Resp. Transitória
Sist.
J θ̈ = Fy l sen θ − Fx l cos θ
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
y y′ Para θ ≈ 0:
Aula 03
m Centro
2l
Aula 04 de
Aula 05
θ ×
massa
J θ̈ = Fy lθ − Fx l
Aula 06
Aula 07 Fx = mẍ + mlθ̈
u x
Resp. Transitória Fy = mg
M
Sist.
Realimentados M ẍ = u − Fx
Lugar das Raı́zes
Resp. em
Frequência
Donde se obtém:
Controle PID
Controle Moderno
(M + m)ẍ + mlθ̈ = u
(J + ml2 )θ̈ − mglθ + mlẍ = 0
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03
Em Sistemas Elétricos, em geral, relaciona-se correntes e
Aula 04
Aula 05
tensões;
Aula 06
Aula 07 Observação
Resp. Transitória
Para evitar erro de modelagem em Sistemas Elétricos, use as
Sist.
Realimentados leis de Kirchhoff!
Lugar das Raı́zes
São três os elementos principais dos Sistemas Elétricos:
Resp. em
Frequência ▶ Resistor;
Controle PID ▶ Indutor;
Controle Moderno
▶ Capacitor.
Sistemas Elétricos: Resistor
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
vr
Introdução
Modelagem
Aula 03
Aula 04 R ir
Aula 05
Aula 06
Resp. Transitória
Sist.
Realimentados vr (t) = Rir (t) (2.32)
Lugar das Raı́zes
Controle Moderno
Er = vr (t)ir (t) d t (2.33)
to
Observação
Teve alguma ideia para calcular a energia no amortecedor?
Sistemas Elétricos: Indutor
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
vl
Aula 03
Aula 04
Aula 05
Aula 06
L il
Aula 07
Resp. Transitória
A relação tensão/corrente é dada por
Sist.
Realimentados
d il (t)
Lugar das Raı́zes
vl (t) = L (2.34)
Resp. em dt
Frequência
Prof. Victor L.
Yoshimura
Introdução
vc
Modelagem
Aula 03
Aula 04
Aula 05 ic
Aula 06 C
Aula 07
Resp. Transitória
Sist.
A relação tensão/corrente é dada por
Realimentados
Prof. Victor L.
Yoshimura
y − x2
Controle PID
Controle Moderno
y = R x1 −
rc
Neste caso, obter a equação de saı́da primeiro facilita a
modelagem.
Exemplo: Um Circuito Elétrico
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
rl L x1
Introdução
Modelagem
Aula 03
rc
Aula 04
Aula 05
Aula 06
u R y
Aula 07
x2 C
Resp. Transitória
Sist.
Realimentados
Controle Moderno
− L(R + rc ) L(R + r )
c x + L u
ẋ =
R 1
− 0
(R + rc )C (R + rc )C
Rrc R
y= x
R + rc R + rc
Exemplo: Um Circuito Elétrico
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura rl L x1
Introdução
Modelagem rc
Aula 03
Aula 04
Aula 05
u R y
Aula 06
Aula 07 x2 C
Resp. Transitória
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Introdução
▶ Procedimento prático para circuitos elétricos:
Modelagem
Aula 03 i) Escolha como variáveis de estado: tensões de
Aula 04
Aula 05
capacitores e correntes de indutores;
Aula 06 ii) Escolha as correntes de malha e expresse-as como
funções das variáveis de estado e suas derivadas;
Aula 07
Resp. Transitória
iii) Escreva as equações de tensão de malha e elimine as
Sist.
Realimentados variáveis que não sejam de estado.
Lugar das Raı́zes ▶ Outro procedimento prático para circuitos elétricos:
Resp. em i) Escolha como variáveis de estado: tensões de
Frequência
capacitores e correntes de indutores;
Controle PID
ii) Substitua cada capacitor (indutor) por uma fonte de
Controle Moderno
tensão (corrente);
iii) Obtenha a corrente (tensão) em cada capacitor
(indutor) e iguale a C v̇c (Li̇l ).
O Amplificador Operacional
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03 vp 0A +
Aula 04
0V − vo
Aula 05 vn
Aula 06
Aula 07
0A
Resp. Transitória
Introdução
Modelagem C
Aula 03
R1 i R2 i
Aula 04 u y
Aula 05
Aula 06
Aula 07 −
+
Resp. Transitória
Sist.
Realimentados
Prof. Victor L.
Yoshimura
R L i(t)
Introdução
Modelagem
J, b
Aula 03
vin (t) + va (t)
Aula 04 − ω
Aula 05
Aula 06
Aula 07
Prof. Victor L.
Yoshimura
Introdução R L i(t)
Modelagem
Aula 03
Aula 04 J, b
vin (t) + va (t)
Aula 05
− ω
Aula 06
Aula 07
Resp. Transitória
Para o circuito elétrico:
Sist.
Realimentados
d i(t)
Lugar das Raı́zes vin (t) − Ri(t) − L − Kω(t) = 0
dt | {z }
Resp. em
va (t)
Frequência
Controle PID
Para a parte mecânica:
Controle Moderno
d ω(t)
J + bω(t) = Ki(t)
dt | {z }
T (t)
O Motor de Corrente Contı́nua
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução R L i(t)
Modelagem
Aula 03
J, b
Aula 04
vin (t) + va (t)
Aula 05 − ω
Aula 06
Aula 07
Resp. Transitória
Fazendo x1 = i, x2 = ω e u = vin :
Sist.
Realimentados R
−K
1
−
Lugar das Raı́zes ẋ = KL L x+ L u
−b 0
Resp. em J J
Frequência y= 0 1 x
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 03
▶ O escoamento de fluidos depende do número de
Aula 04
Aula 05
Reynolds
Aula 06 ρvD
Aula 07 R= (2.39)
Resp. Transitória
µ
Sist.
Realimentados ▶ ρ é a massa especı́fica do fluido (em kg m−3 );
Lugar das Raı́zes ▶ v é a velocidade de escoamento;
Resp. em ▶ D é o diâmetro da tubulação, e;
Frequência ▶ µ é a viscosidade dinâmica (em Pa · s).
Controle PID
▶ R ≤ 2000: escoamento laminar, e;
Controle Moderno
▶ R > 2000: escoamento turbulento (não linear).
Sistemas de Controle de Nı́vel: Resistências
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
R = K −1 (2.42)
h
Resp. Transitória
Q
Sist.
Realimentados Escoamento Turbulento:
Lugar das Raı́zes √
Resp. em A resistência da válvula é dada Q=K h (2.43)
Frequência
por:
Controle PID
dh 2√
Controle Moderno R= (2.40) R= h (2.44)
dQ K
Sistemas de Controle de Nı́vel: Capacitâncias
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
A capacitância do tanque é
Introdução dada por:
Modelagem
Aula 03
h dV
Aula 04
Q C= (2.45)
Aula 05
Aula 06
dh
Aula 07
Resp. Transitória
Para secção transversal constante, V = Ah, logo
Sist.
Realimentados
Controle PID
Controle Moderno
Observação
Para modelar sistemas de nı́vel, use a equação da vazão:
dV d V d h (2.45)
Q= = = C ḣ (2.47)
dt dh dt
Exemplo: Suprimento de Fluido a um Tanque
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Com (2.41)
Introdução
Modelagem
1
qo = h
Aula 03
Qin + qin R
Aula 04
Aula 05
Aula 06 Usando (2.47)
Aula 07
Resp. Transitória C
C ḣ = V̇ = qi − qo
Sist.
H +h
Realimentados
Controle Moderno
Logo:
ĥ(s) R
=
q̂i (s) RCs + 1
Exemplo: Tanques Comunicantes
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
q
Modelagem
Aula 03
Aula 04
Aula 05 C1 C2
Aula 06
Aula 07
h1
Resp. Transitória
R1 R2
h2
Sist.
Realimentados
Controle PID
Para o registro 1, de (2.41):
Controle Moderno
ĥ1 (s) 1
q̂1 (s)
Σ R1
h1 − h2 = R1 q1 −
ĥ2 (s)
Exemplo: Tanques Comunicantes
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
q
Modelagem
Aula 03
Aula 04
Aula 05
Aula 06
Aula 07
C1 C2
Resp. Transitória
h1
Sist. R1 R2
Realimentados
h2
Lugar das Raı́zes
Resp. em
q1 q2
Frequência
Controle PID
Para o registro 2, de (2.41):
Controle Moderno
h2 = R2 q2 ĥ2 (s) 1
q̂2 (s)
R2
Exemplo: Tanques Comunicantes
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
q
Modelagem
Aula 03
Aula 04
Aula 05 C1 C2
Aula 06
Aula 07
h1
Resp. Transitória
R1 R2
h2
Sist.
Realimentados
Controle PID
Para o tanque 1, de (2.47):
Controle Moderno
q̂(s) 1
ĥ1 (s)
Σ sC1
C1 ḣ1 = q − q1 −
q̂1 (s)
Exemplo: Tanques Comunicantes
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
q
Modelagem
Aula 03
Aula 04
Aula 05 C1 C2
Aula 06
Aula 07
h1
Resp. Transitória
R1 R2
h2
Sist.
Realimentados
Controle PID
Para o tanque 2, de (2.47):
Controle Moderno
q̂1 (s) 1
ĥ2 (s)
Σ sC2
C2 ḣ2 = q1 − q2 −
q̂2 (s)
Exemplo: Tanques Comunicantes
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução q
Modelagem
Aula 03
Aula 04
Aula 05
Aula 06
C1 C2
Aula 07
Resp. Transitória
h1
R1 R2
h2
Sist.
Realimentados
Prof. Victor L.
Yoshimura
▶ Sistemas térmicos: transferência de calor;
Introdução
▶ A rigor são sistemas distribuı́dos (n = ∞);
Modelagem
Aula 03 ▶ A transferência de calor se dá por:
Aula 04
Aula 05 ▶ Condução;
Aula 06
Aula 07
▶ Convecção (envolve transferência de massa);
Resp. Transitória
▶ Radiação (altas temperaturas).
Sist. ▶ O modelo similar ao de Sistemas de Nı́vel.
Realimentados
▶ Nı́vel (h) ↔ Temperatura (θ);
Lugar das Raı́zes
▶ Vazão (q) ↔ Fluxo de Calor (h);
Resp. em
Frequência
Controle PID
Observação
Controle Moderno Para modelar Sistemas Térmicos:
θ = Rh (2.48)
h = C θ̇ (2.49)
Exemplo: Aquecedor de Fluido
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Relação entre o tanque e a saı́da:
Introdução
Modelagem (2.48) TL
Aula 03
Fluido quente θ = Rho ⇒ θ̂(s) = Rĥo (s)
Aula 04
Resp. Transitória
(2.49) hi − ho
Aquecedor θ̇ =
Sist. C
Realimentados
Prof. Victor L.
Yoshimura
Fluido quente
Introdução Usando a linearidade
Modelagem (hi = 0):
Aula 03
Aula 04
Aula 05 Aquecedor 1
Aula 06 C θ̇ = htotal = (θi − θ)
Aula 07 R
Resp. Transitória
Fluido frio
Sist.
Aplicando TL:
Realimentados
Prof. Victor L.
Yoshimura
Elementos podem apresentar comportamento não-linear, da forma:
Introdução
Prof. Victor L.
Yoshimura
Introdução
u
Modelagem
Aula 03
Note que:
Aula 04
y
Aula 05
1 √ dy 1
Aula 06
y=√ F ⇒ = √
Aula 07
m k dF 2 kF
Resp. Transitória
Sist. Aplicando
r (2.52) em
Realimentados
mg
Lugar das Raı́zes
k mg mg, :
Resp. em
k
Frequência r
Controle PID mg 1
y− ≈ √ (u − mg)
Controle Moderno
Mola não linear: k 2 mgk
Fmola = ky 2
Sistemas MIMO Não Lineares
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Considere a i-ésima saı́da de um elemento não linear, da
Modelagem
Aula 03
forma
Aula 04
Aula 05
yi = fi (u1 , u2 , . . . , um ) (2.53)
Aula 06
Aula 07 Tome um ponto de operação (ũ1 , ũ2 , . . . , ũm , ỹi ) e,
Resp. Transitória similarmente a (2.52):
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Aplicando a cada saı́da, tem-se:
Introdução
∂f1 ∂f1 ∂f1
Modelagem
Aula 03
···
∂u1 ∂u2 ∂um
Aula 04 y1 − ỹ1 u1 − ũ1
Aula 05 ∂f2 ∂f2 ∂f2
y2 − ỹ2 ··· u2 − ũ2
Aula 06
∂u ∂u2 ∂um
Aula 07 .. ≈
1 ..
. .. .
Resp. Transitória
.
··· ··· ···
Sist. yn − ỹn
∂fn ∂fn
u − ũ
∂fn m m
Realimentados
···
Lugar das Raı́zes ∂u1 ∂u2 ∂um
| {z }
Resp. em J
Frequência
(2.55)
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Dinâmica não linear:
Modelagem
Aula 03
Aula 04 T = −mgl sen θ
Aula 05
Aula 06
Aula 07
No ponto (θ̃, T̃ ) = (0, 0):
Resp. Transitória
l
Sist. ∂T
Realimentados T − T̃ = (θ − θ̃)
Lugar das Raı́zes ∂θ θ=θ̃
Resp. em
Frequência m Ou seja:
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
▶ A resposta no tempo de SLIT-Cs é composta de:
Introdução
Modelagem
▶ Resposta transitória: para t < ∞;
▶ Resposta em regime permanente: para t → ∞.
Resp. Transitória
Aula 08
Aula 09
▶ Entre outros critérios, deve-se verificar:
Sist.
▶ Estabilidade: Na perturbação, o sistema
Realimentados ▶ Retorna ao ponto de operação?
Lugar das Raı́zes ▶ Apresenta oscilações sustentadas?
Resp. em ▶ Diverge do ponto de operação?
Frequência
▶ Erro em regime permanente
Controle PID
▶ O sistema segue a referência dada?
Controle Moderno ▶ Se sim, quanto tempo leva para o fazer (estabilidade
relativa)?
Observação
Iremos estudar a resposta transitória de sistemas de 1ª e de
2ª ordens. E isto é suficiente, como será visto mais adiante!
Sistemas de 1ª Ordem
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Descrição geral saı́da/referência de um sistema de controle
Aula 08 de 1ª ordem:
Aula 09
ŷ(s) 1
Sist. = (3.1)
Realimentados r̂(s) τs + 1
Lugar das Raı́zes
Escrevendo sob a forma de realimentação unitária (2.1)
Resp. em
Frequência
1
Controle PID
ŷ(s)
Controle Moderno = τs (3.2)
r̂(s) 1
1+
τs
Sistemas de 1ª Ordem: Resposta ao Degrau
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Sist.
Realimentados
Observação
Verifique que A = −B = 1 em (3.3).
Sistemas de 1ª Ordem: Resposta ao Degrau
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Introdução y, ε
Modelagem
1,0
Resp. Transitória
0,9
Aula 08
Aula 09
0,8
Sist.
Realimentados 0,7
Resp. em 0,5
Frequência
0,4
Controle PID
0,2
0,1
0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 t
τ
Sistemas de 1ª Ordem: Resposta à Rampa
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Retornando a (3.1) e aplicando r(t) = t, t ≥ 0:
Resp. Transitória
1 1 A B C
Aula 08
Aula 09
ŷ(s) = = 2+ + (3.6)
τ s + 1 s2 s s τs + 1
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Introdução y ε
rampa, ,
τ τ
Modelagem
5,0
Resp. Transitória
Aula 08 4,5
Aula 09
4,0
Sist.
Realimentados 3,5
1,0
0,5
0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 t
τ
Sistemas de 1ª Ordem: Resposta ao Impulso
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Retornando a (3.1) e aplicando r(t) = δ(t):
Resp. Transitória
Aula 08
1
Aula 09
ŷ(s) = (3.9)
Sist. τs + 1
Realimentados
Observação
Note que ε = −y, ∀t > 0.
Sistemas de 1ª Ordem: Resposta ao Impulso
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução τy
Modelagem
1,0
Resp. Transitória
0,9
Aula 08
Aula 09
0,8
Sist.
Realimentados 0,7
Resp. em 0,5
Frequência
0,4
Controle PID
0,2
0,1
0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 t
τ
Um Lema Importante!
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Sist.
Realimentados Lema
Lugar das Raı́zes Seja u a entrada de um SLIT-C, que produz a resposta
Resp. em forçada y. Então, ao submeter o sistema a u̇ em sua
Frequência
Controle PID
entrada, a resposta forçada será ẏ.
Controle Moderno
Demonstração.
Exercı́cio!
Sistemas de 2ª Ordem
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
A forma geral de um sistema de 2ª ordem é:
Resp. Transitória
Aula 08
Aula 09 ŷ(s) ωn2
= 2 (3.11)
Sist.
Realimentados
r̂(s) s + 2ζωn s + ωn2
Lugar das Raı́zes
onde:
Resp. em
Frequência ωn → frequência natural não-amortecida;
Controle PID ζ → amortecimento.
Controle Moderno Um parâmetro importante para estes sistemas é a
atenuação, dada por:
σ = ζωn (3.12)
DB em Realimentação Unitária
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Aula 08 r̂(s) ε̂(s) 2
ωn
ŷ(s)
Aula 09 Σ s(s+2ζωn )
Sist. −
Realimentados
Resp. em
Frequência Aplicando a regra da realimentação unitária (2.1) a (3.11):
Controle PID
Controle Moderno
ŷ(s) ωn2 ĝ(s)
= 2 2
=
r̂(s) s + 2ζωn s + ωn 1 + ĝ(s)
Sistemas de 2ª Ordem: Resposta ao Degrau
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ É a principal resposta dos SLIT-Cs;
Resp. Transitória
Aula 08 ▶ Para projetos, normalmente as plantas são reduzidas a
Aula 09
Sist.
um par de polos;
Realimentados
▶ Fazendo r = 1̃:
Lugar das Raı́zes " #
ω 2
Resp. em
Frequência
y(t) = L−1 n
(3.13)
s(s2 + 2ζωn s + ωn2 )
Controle PID
Controle Moderno
Observação
O comportamento de sistemas de 2ª ordem é altamente
dependente de ζ!
Caso Subamortecido (0 < ζ < 1)
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Neste caso, escreve-se:
Resp. Transitória
Aula 08
A Bs + C
Aula 09
ŷ(s) = + 2 (3.14)
Sist. s s + 2ζωn s + ωn2
Realimentados
Observação
Mostre que A = −B = 1 e que C = −2ζωn .
Caso Subamortecido (0 < ζ < 1)
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aplicando a técnica de completar quadrados:
Resp. Transitória
Aula 08 s2 + 2ζωn s + ζ 2 ωn2 + ωd2 = (s + ζωn )2 + ωd2
Aula 09 | {z }
2
ωn
Sist.
Realimentados
(3.16)
Caso Subamortecido (0 < ζ < 1)
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Aplicando a TIL a (3.16):
Introdução
ζ
Modelagem y(t) = 1̃(t) − e−ζωn t cos ωd t − p e−ζωn t sen ωd t
Resp. Transitória
1 − ζ2
Aula 08
Aula 09
e−ζωn t
p
Sist.
1 − ζ2
= 1̃(t) − p sen ωd t + arctg (3.17)
Realimentados
1 − ζ2 ζ
Lugar das Raı́zes | {z }
Resp. em arccos ζ
Frequência
Observação
(3.17)-(3.18) são as respostas mais importantes dos
SLIT-C-SISOs!
Caso Criticamente Amortecido (ζ = 1)
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Neste caso, (3.13) toma a forma mais simples:
Resp. Transitória
Aula 08 " #
Aula 09
ωn2
Sist. y(t) = L−1 (3.19)
Realimentados s(s + ωn )2
Lugar das Raı́zes A B1 B2
Resp. em = L−1 + + (3.20)
Frequência s s + ωn (s + ωn )2
Controle PID
Controle Moderno
Aplicando a TIL a (3.20):
Prof. Victor L.
Yoshimura As raı́zes de s2 + 2ζωn s + ωn2 são reais e distintas:
Introdução q
Modelagem s = −ζωn ± ωn ζ 2 − 1 (3.22)
Resp. Transitória
Aula 08 Aplicando a EFP a (3.13):
Aula 09
Sist. "
Realimentados
−1 A B
Lugar das Raı́zes y(t) = L + p +
s s + ζωn − ωn ζ 2 − 1
Resp. em
Frequência
#
C
Controle PID + p (3.23)
Controle Moderno
s + ζωn + ωn ζ 2 − 1
Observação
Mostre que A = 1,
−1 1
B= p p eC= p p .
2 ζ 2 − 1(ζ − ζ 2 − 1) 2 ζ 2 − 1(ζ + ζ 2 − 1)
Caso Sobreamortecido (ζ > 1)
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aplicando a TIL:
Modelagem
Resp. Transitória
√ √
Aula 08 1 e−(ζ+ ζ 2 −1)ωn t e−(ζ− ζ 2 −1)ωn t
Aula 09 y(t) = 1 + p 2 p − p
Sist. 2 ζ −1 ζ + ζ2 − 1 ζ − ζ2 − 1
Realimentados
(3.24)
Lugar das Raı́zes
Resp. em Observação
Frequência
Se ζ ≫
p 1, o polo em (ζ − ζ 2 − 1)ωn domina o polo em
p
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Aula 08
y(ζ=0,2) , y(ζ=0,5) , y(ζ=0,707) , y(ζ=1)
Aula 09
1,6
Sist. 1,5
Realimentados 1,4
1,3
1,2
Lugar das Raı́zes 1,1
1,0
Resp. em 0,9
0,8
Frequência 0,7
0,6
Controle PID 0,5
0,4
0,3
Controle Moderno 0,2
0,1
0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 t
ω
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura r = 1̃, y
1 + Mp
Introdução
Modelagem
Resp. Transitória
Aula 08
Aula 09
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Mp , ωn tp /10, β/100[◦ ]
Introdução 1,0
Modelagem 0,9
Controle Moderno ζ
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0
Observação
Muito do projeto depende do amortecimento, ζ. Se for baixo, o
sistema responde rapidamente (os tempos envolvidos serão
pequenos), mas apresentará overshoot excessivo e vice-versa.
Muitos projetos tentam fazer 0,4 < ζ < 0,8.
Resposta de Segunda Ordem com Zero
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Considere uma FT de 2ª ordem subamortecida, com zero,
Resp. Transitória
Aula 08
normalizada:
Aula 09
s
Sist. γζωn +1
Realimentados ĝ(s) = 2 , γ>0 (3.30)
s
Lugar das Raı́zes
ωn + 2ζ ωsn + 1
Resp. em
Frequência
Considerando (1.20) e (3.17) para a resposta ao degrau:
Controle PID
Controle Moderno
e−ζωn t 1
y(t) = 1̃(t) − p sen (ωd t + arccos ζ) − sen ωd t
1−ζ 2 γζ
(3.31)
Resposta de Segunda Ordem com Zero
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Para ζ = 0,5:
Modelagem y(γ=1) , y(γ=2) , y(γ=4) , y(γ=∞)
1,8
Resp. Transitória 1,6
Aula 08 1,4
Aula 09 1,2
1,0
Sist. 0,8
Realimentados 0,6
0,4
Lugar das Raı́zes 0,2
1,0 2,0 3,0 4,0 5,0 6,0 7,0 8,0
Resp. em ωn t
Frequência
Controle PID
Prof. Victor L.
Yoshimura
Para recalcular o overshoot, derive-se (3.31):
Introdução
e−ζωn t 1
Modelagem
Resp. Transitória
ẏ = p ζωn sen (ωd t + arccos ζ) − sen ωd t −
Aula 08
1−ζ 2 γζ
1
Aula 09
Sist. −ωd cos (ωd t + arccos ζ) − cos ωd t
Realimentados γζ
Lugar das Raı́zes
Resp. em
Os pontos crı́ticos, (tcr , ycr ), são obtidos com ẏ = 0, o que
Frequência leva a
Controle PID p
Controle Moderno 1 − ζ2
(γ − 1) sen ωd tcr + cos ωd tcr = 0. (3.32)
ζ
Do qual se conclui que
π
ωn tcr = p ⇐⇒ γ = 1. (3.33)
2 1 − ζ2
Resposta de Segunda Ordem com Zero
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Se γ ̸= 1 e cos ωd tcr ̸= 0, as soluções de (3.32) são os
valores positivos de
Introdução
p
Modelagem 1 1 − ζ2
Resp. Transitória
ωn tcr =p arctg . (3.34)
Aula 08
1 − ζ2 ζ(1 − γ)
Aula 09
7
Sist.
Realimentados 6.5 γ=1
γ=2
Lugar das Raı́zes 6 γ=3
Resp. em 5.5 γ=4
Frequência γ→∞
5
Controle PID
4.5
ωn tcr
Controle Moderno 4
3.5
3
2.5
2
1.5
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
ζ
Resposta de Segunda Ordem com Zero
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Substituindo (3.33) e (3.34) em (3.31), obtem-se
Introdução
!
Modelagem
−πζ
Resp. Transitória − exp p
2 1 − ζ2
π 1
Aula 08
Aula 09 Mp = p sen + arccos ζ − (3.35)
1 − ζ2 2 ζ
Sist.
Realimentados
e
Lugar das Raı́zes
Resp. em
Frequência p !
ζ 1 − ζ2
Controle PID − exp arctg "
− γ)
p
ζ(1
p
Controle Moderno
1−ζ 2
1 − ζ2
Mp = sen arctg +
ζ(1 − γ)
p
1− ζ2
! p #
1 1 − ζ2
+ arccos ζ − sen arctg , (3.36)
γζ ζ(1 − γ)
respectivamente.
Resposta de Segunda Ordem com Zero
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
1.2
Introdução
1.1 γ=1
Modelagem γ=2
1 γ=3
Resp. Transitória
Aula 08
0.9 γ=4
Aula 09
γ→∞
Sist. 0.8
Realimentados
0.7
Lugar das Raı́zes
Mp
0.6
Resp. em
Frequência 0.5
Controle PID
0.4
Controle Moderno
0.3
0.2
0.1
0
0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8
ζ
Redução de Modelos: Considerações Iniciais
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ Muitos SLIT-Cs são de ordem elevada;
Modelagem ▶ Sua representação por FTs é simples, porém o projeto
Resp. Transitória
Aula 08
pode ser bastante complexo;
Aula 09 ▶ Como os polos dominantes possuem resposta mais
Sist.
Realimentados
persistente:
Lugar das Raı́zes
▶ É possı́vel aproximar ĝ(s) por outra FT, ĝr (s), de
Resp. em
ordem menor?
Frequência ▶ Como fazer isto?
Controle PID ▶ Qual o erro introduzido?
Controle Moderno ▶ Assunções:
▶ ĝr tem ordem ñ < n (não faz sentido o caso ñ = n,
evidentemente);
▶ ĝr pode apresentar m̃ zeros, mas m̃ ≤ ñ;
▶ Pares de polos complexos conjugados devem ser retidos
e com constantes multiplicativas conjugadas.
Procedimento via Polos Dominantes I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. em ñ
bi
(3.38)
X
Frequência
ĝr (s) =
Controle PID
i=1
s − pi
Controle Moderno
Onde os ñ polos dominantes foram retidos no modelo.
Observação
Em geral, não é interessante fazer bi = ai no modelo
reduzido, pois o erro introduzido torna-se maior.
Procedimento via Polos Dominantes II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem Observação
Resp. Transitória
Aula 08
Proposta: A escolha dos coeficientes bi tratará do
Aula 09 comportamento em regime permanente de ĝr para t → ∞.
Sist.
Realimentados
Justificativa: Os polos retidos já tratam do comportamento
Lugar das Raı́zes
transitório.
Resp. em Considere as entradas u = ti 1̃, i = 0, 1, . . . , ñ. Para cada
Frequência
Controle PID
entrada, pode-se determinar o regime permanente com:
Controle Moderno i
ĝ(s) X 1 di−j ĝ(0)
= + r̂i (s) (3.39)
si+1 j=1 sj d si−j
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Modelagem 150(s + 1)
ĝ(s) =
Resp. Transitória (s + 2)(s2 + 2s + 2)(s2 + 6s + 18)
Aula 08
Aula 09
Sist. ▶ Polos: −1 ± j, −2 e −3 ± j3
Realimentados
▶ Derivadas de ĝ: ĝ(0) = 2,0833 e ĝ ′ (0) = −1,7361
Lugar das Raı́zes
▶ Usando (3.40):
Resp. em
Frequência
−1 −1
h i h i
Controle PID −1+j −1−j b1 2,0833
−1 −1 = ⇒ b1,2 = 0,3472 ∓ j1,7361
(−1+j)2 (−1−j)2
b2 −1,7361
Controle Moderno
Prof. Victor L.
Yoshimura
Introdução
Sist. −
Realimentados
Aula 10
z(t) Controlador Automático
Aula 11
Aula 12
Aula 13
Sensor
Lugar das Raı́zes
Resp. em
Frequência
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Introdução
se ε(t) > 0
(
Aula 12
Uon
Aula 13
u(t) = (4.1)
Lugar das Raı́zes Uoff se ε(t) < 0
Resp. em
Frequência
Implementação eletrônica:
Controle PID
Controle Moderno +
referência
− erro
sensor
Compensador On/Off II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Pode haver a inserção, acidental ou não, de um laço de
Introdução histerese em (4.1):
Modelagem
Resp. Transitória
Uon
se ε(t) > εl
Sist.
Realimentados
u(t) = Uoff se ε(t) < −εl (4.2)
Aula 10
caso contrário
u(t )
Aula 11 −
Aula 12
Aula 13
Controle Moderno
ε u
−εl εl ε
Compensadores Proporcional e Integral
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura ▶ O compensador On/Off é não-linear!
Introdução ▶ Os demais compensadores aqui estudados serão lineares.
Modelagem ▶ O compensador P é o mais simples e, ainda, veremos
Resp. Transitória
que é de estudo importante.
Sist.
Realimentados ▶ Sua compensação é proporcional ao erro lido:
Aula 10
Aula 11
Aula 12 u(t) = kp ε(t) (4.3)
Aula 13
▶ Ou, no domı́nio s:
û(s) ki
= (4.5)
ε̂(s) s
Compensador Proporcional mais Integral
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura É a combinação dos dois anteriores:
Introdução Z t
kp
Modelagem u(t) = kp ε(t) + ε(t) d t (4.6)
Resp. Transitória
Ti 0
Sist.
Realimentados kp
Aula 10
Observe que ki = . Aplicando a TL:
Aula 11 Ti
Aula 12
Aula 13
û(s) kp Ti s + kp
Lugar das Raı́zes = (4.7)
Resp. em
ε̂(s) Ti s
Frequência
Ti t
Compensador Proporcional mais Derivativo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
É uma tentativa de inserir um efeito antecipativo no
Modelagem
controle:
Resp. Transitória u(t) = kp ε(t) + kp Td ε̇(t) (4.8)
Sist.
Realimentados Aplicando a TL:
Aula 10
Aula 11
û(s)
Aula 12
Aula 13
= kp + kp Td s (4.9)
ε̂(s)
Lugar das Raı́zes
Resp. em
Frequência u(t)
Controle PID
PD
Controle Moderno
kp Td P
Td t
Controle e
Servomecanismos Compensador Proporcional mais Integral mais Derivativo
Prof. Victor L.
Yoshimura
Introdução
Modelagem
É a tentativa de combinar as propriedades dos
Resp. Transitória
compensadores anteriores.
Sist. Z t
Realimentados kp
Aula 10 u(t) = kp ε(t) + ε(t) d t + kp Td ε̇(t) (4.10)
Aula 11 Ti 0
Aula 12
Aula 13
Aplicando a TL:
Lugar das Raı́zes
Resp. em
û(s) kp Td Ti s2 + kp Ti s + kp
Frequência
= (4.11)
Controle PID ε̂(s) Ti s
Controle Moderno
Observação
Os compensadores PD e PID não são causais!
Efeito do Sensor na Dinâmica do Sistema
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Controle Moderno
Observação
Analise o comportamento dinâmico do LM35, conforme seu
datasheet.
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Considere o sistema em realimentação unitária:
Introdução r̂(s) ε̂(s) ŷ(s)
Modelagem Σ ĉ(s) ĝ(s)
Resp. Transitória −
compensador planta
Sist.
Realimentados
Aula 10
Aula 11 Defina as FTs do compensador e da planta:
Aula 12
Aula 13
n̂g (s)
Lugar das Raı́zes ĝ(s) = (4.12)
Resp. em dˆg (s)
Frequência
Prof. Victor L.
Yoshimura
Observação
Se a entrada for um degrau, a FTMA deverá ter um polo na
origem para erro nulo!
Compensadores Integral e Derivativo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ A principal função do compensador I é eliminar o erro
Modelagem
Resp. Transitória
em regime em sistemas sob entrada degrau.
Sist.
ki
Realimentados
ĉint (s) = (4.16)
Aula 10
Aula 11
s
Aula 12
Aula 13 ▶ O compensador D “aumenta a estabilidade” do sistema;
Lugar das Raı́zes
Controle PID
▶ O compensador D apresenta problemas com
Controle Moderno
amplificação de ruı́dos;
▶ Pior ainda, é anti-causal! (Como implementar seu
efeito?)
Exemplo: Controle de Nı́vel
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Modelo da planta (vide
Introdução aula 07):
Modelagem
Aula 13
Controle Moderno
Calculemos o erro em regime com (4.15)
RCs + 1 1 1
ε(∞) = lim s = >0
s→0 RCs + 1 + kp R s 1 + kp R
Exemplo: Servomotor DC
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução ˆ
d(s)
Modelagem
r̂(s) ε̂(s) 1
ŷ(s)
Resp. Transitória ĉ(s)
Σ Σ s(Js+b)
Sist. −
Realimentados
Aula 10
Aula 11
Aula 12
Aula 13
▶ Planta com polo na origem: erro nulo para a referência;
Lugar das Raı́zes
Resp. em
▶ Porém, não para a perturbação (não está em
Frequência
realimentação unitária);
Controle PID
▶ Considerando compensador P:
Controle Moderno
1
ε̂(s) s(Js + b) 1
=− =− 2
ˆ
d(s) kp Js + bs + kp
1+
s(Js + b)
Exemplo: Servomotor DC
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
ˆ
d(s)
Introdução
Modelagem
r̂(s) ε̂(s) 1
ŷ(s)
Σ ĉ(s) Σ s(Js+b)
Resp. Transitória
−
Sist.
Realimentados
Aula 10
Prof. Victor L.
Yoshimura
Introdução
Modelagem ˆ
d(s)
Resp. Transitória
Resp. em
▶ Porém... Polinômios incompletos são sempre instáveis!
Frequência
▶ Os problemas anteriores podem ser resolvidos com um
Controle PID
PI bem projetado.
Controle Moderno
Observação
Refaça este exemplo, com o modelo obtido na Aula 06.
Exemplo: Posicionamento de Satélite
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
r̂(s) ε̂(s) θ̂(s)
1
Σ ĉ(s) Js2
Resp. Transitória
−
Sist.
Realimentados
Aula 10
Controle Moderno
Para compensador P, existirão oscilações sustentadas, pois
os polos serão: s
kp
p1,2 = ±j
J
Exemplo: Posicionamento de Satélite
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Determinemos a relação entre o erro e a FTMA na
realimentação unitária:
Introdução
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Este é um caso mais simples de (4.14), originando:
Sist.
s
Realimentados
ε(∞) = lim ε(t) = lim sε̂(s) = lim r̂(s) (4.20)
Aula 10
t→∞ s→0 s→0 1 + ĝ(s)
Aula 11
Aula 12
Aula 13
Para três sinais de entrada, determinaremos o erro estático:
Lugar das Raı́zes
▶ Entrada degrau: origina o erro estático de posição;
Resp. em
Frequência
▶ Entrada rampa: origina o erro estático de velocidade;
Controle PID
Controle Moderno
▶ Entrada quadrática: origina o erro estático de
aceleração.
Erro Estático de Posição
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Importante!
Para anular o erro estático de posição, é necessário que o
sistema seja, ao menos, de tipo 1.
Erro Estático de Velocidade
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Sendo a referência uma rampa, tem-se:
Introdução
s 1 1
Modelagem
ε(∞) = lim 2
= lim (4.23)
Resp. Transitória s→0 1 + ĝ(s) s s→0 sĝ(s)
Sist.
Realimentados Define-se a constante de erro estático de velocidade e,
Aula 10
Aula 11 referindo-se a (4.19):
Aula 12
Aula 13
Lugar das Raı́zes
se N = 0
0
Resp. em
Frequência
Kv = lim sĝ(s) = K se N = 1 (4.24)
s→0
∞ se N > 1
Controle PID
Controle Moderno
Importante!
Se o sistema for de tipo 0, o erro estático de velocidade
torna-se infinito. Para anulá-lo, é necessário que o sistema
seja, ao menos, de tipo 2.
Erro Estático de Aceleração
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Sendo a referência r(t) = t2 1̃(t), tem-se:
Introdução
s 1 1
Modelagem
ε(∞) = lim 3
= lim 2 (4.25)
Resp. Transitória s→0 1 + ĝ(s) s s→0 s ĝ(s)
Sist.
Realimentados Define-se a constante de erro estático de aceleração e,
Aula 10
Aula 11 referindo-se a (4.19):
Aula 12
Aula 13
Lugar das Raı́zes se N ≤ 1
0
Resp. em
Frequência
2
Ka = lim s ĝ(s) = K se N = 2 (4.26)
s→0
Controle PID
∞ se N > 2
Controle Moderno
Importante!
Se o sistema for de tipo 0 ou 1, o erro estático de aceleração
é infinito. Para anulá-lo, é necessário que o sistema seja, ao
menos, de tipo 3.
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Considerando
Introdução r̂(s) ŷ(s) n̂g (s)
Σ ĝ(s) ĝ(s) = (4.27a)
Modelagem
− dˆg (s)
Resp. Transitória
n̂h (s)
Sist. ĥ(s) ĥ(s) = (4.27b)
Realimentados dˆh (s)
Aula 10
Aula 11 Assim, observe que são possı́veis as formas:
Aula 12
Aula 13
ŷ(s) n̂g (s)dˆh (s)
Lugar das Raı́zes
= (4.28a)
Resp. em r̂(s) dˆg (s)dˆh (s) + n̂g (s)n̂h (s)
Frequência
m
, n
Controle PID
(4.28b)
X X
j
= βj s αi si , αi , βj ∈ R
Controle Moderno
j=0 i=0
m
, n
(s − pi ), K ∈ R e pi , zj ∈ C
Y Y
=K (s − zj )
j=1 i=1
(4.28c)
1º Caso: Apenas Polos Reais
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. em n
(4.30)
X
Frequência y(t) = Ao + Ai epi t
Controle PID i=1
Controle Moderno
Observação
Quando um sistema de ordem qualquer com polos todos
reais será estável?
1º Caso: Apenas Polos Reais
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Observação
Modelagem
Controle Moderno
relação aos demais), sua exponencial reduz-se muito
mais rapidamente;
v O efeito deste polo pode, então, ser desprezado;
vi Efeitos análogos ocorrem para polos múltiplos.
2º Caso: Presença de Polos Não Reais
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Aplicando a TIL a (4.31b)
Introdução
n n2
Modelagem
e−ζk ωk t (Bk cos ωdk t+Ck sen ωdk t)
X X
y(t) = Ao + Ai epi t +
Resp. Transitória
i=1 k=1
Sist.
Realimentados
(4.32)
Aula 10
Aula 11 Observação
Aula 12
Aula 13 Quando um sistema de ordem qualquer será estável?
Lugar das Raı́zes
Resp. em Observação
Frequência
Prof. Victor L.
Yoshimura
Introdução
Modelagem Importante!
Resp. Transitória
Polo dominante: Se p1 e p2 são polos com parte real
Sist.
Realimentados negativa e Re(p1 ) ≫ Re(p2 ), então os efeitos de p2 são
Aula 10
Aula 11
suprimidos pelos de p1 . Diz-se que p1 domina p2 .
Aula 12
Aula 13
Importante!
Lugar das Raı́zes
Resp. em
Estabilidade: Para garantir que y(t) não cresça
Frequência indefinidamente (para entrada limitada), devemos assegurar
Controle PID que todos os polos da FTMF estejam no semi-plano
Controle Moderno
complexo esquerdo. Neste caso, diz-se que o sistema é
estável. A estabilidade de um sistema é uma caracterı́stica
inerente, independendo da entrada aplicada.
Polinômios e Estabilidade
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
▶ Os polos da FTMF devem estar no semi-plano esquerdo;
Sist.
Realimentados
Aula 10
▶ Fatorar polinômios é uma tarefa complexa;
Aula 11
Aula 12
▶ Acima do 5º grau, é impossı́vel através de operações
Aula 13
elementares em seus radicais (teorema de Abel-Ruffini);
Lugar das Raı́zes
Resp. em
Frequência
Observação
Controle PID Objetivo do Critério de Routh-Hurwitz (CRH): determinar se
Controle Moderno um sistema é estável sem fatorar seus polinômios.
O Critério de Routh-Hurwitz
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Considere a FTMF do sistema em malha fechada como em
Resp. Transitória
(2.4):
Sist.
Realimentados ŷ(s) ĝ(s) n̂(s)
Aula 10
= = (4.33)
Aula 11
r̂(s) 1 + ĝ(s)ĥ(s) ˆ
d(s)
Aula 12
Aula 13
Referindo ao denominador de (4.33), considere o polinômio:
Lugar das Raı́zes
Controle PID
Prof. Victor L.
Yoshimura
Com coeficientes todos positivos, escreva o arranjo de Routh:
Introdução
Modelagem sn ao a2 a4 a6 ...
Resp. Transitória sn−1 a1 a3 a5 a7 ...
Sist.
Realimentados
sn−2 b1 b2 b3 b4 ...
Aula 10
Aula 11
Aula 12
onde ai = αn−i e :
Aula 13
Prof. Victor L.
Yoshimura
Introdução
Modelagem
ˆ = 5s4 + 4s3 + 3s2 + 2s + 1 é estável?
O polinômio d(s)
Resp. Transitória
Sist.
Realimentados s4 5 3 1
Aula 10
Aula 11
s3 4 2 Dividindo esta linha por 2...
Aula 12
Aula 13
s3 2 1
Lugar das Raı́zes
s2 0,5 1
Resp. em s1 −3
Frequência
s0 1
Controle PID
Prof. Victor L.
Yoshimura
Resp. Transitória
ˆ = s3 + s2 + s + k
d(s)
Sist.
Realimentados
Aula 10
é estável?
Aula 11
Aula 12 s3 1 1
Aula 13
s2 1 k
Lugar das Raı́zes
Resp. em
s1 1−k
Frequência s0 k
Controle PID
0<k<1
para estabilidade.
Elemento Nulo na 1ª Coluna
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
▶ De acordo com (4.32) (e (4.30)) a parte real da raiz
Sist.
Realimentados determina a rapidez do transitório;
Aula 10
Aula 11 ▶ Assim, é interessante não só garantir a estabilidade de
Aula 12
Aula 13 um sistema, mas esta queda;
Lugar das Raı́zes
▶ Ou seja, Re(pi ) < σmin , i = 1, 2, . . . , n;
Resp. em
Frequência ▶ Isto equivale a fazer, no arranjo de Routh:
Controle PID
s = ŝ + σ (4.38)
Controle Moderno
Exemplo 3
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
ˆ = s3 + s2 + s + k é estável para
Sabemos que d(s)
Resp. Transitória
Sist.
0 < k < 1. Determinemos a faixa de valores para a qual os
Realimentados
Aula 10
polos terão parte real menor que −0,5.
Aula 11 Para tanto, com (4.38):
Aula 12
Aula 13
Controle Moderno Donde concluı́mos que não existe k para este decaimento
mı́nimo. (Por quê?)
Exemplo 4
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Considere o sistema de controle
Introdução
Modelagem r̂(s) 1
ŷ(s)
Resp. Transitória
Σ ĉ(s) (s−1)(s+2)
Sist.
−
Realimentados
Aula 10
Aula 11
Aula 12
Aula 13
Determinemos os compensadores PI capazes de estabilizar
Lugar das Raı́zes
este sistema. Usando (4.7):
Resp. em
Frequência kp s + ki
ĉ(s) =
Controle PID s
Controle Moderno
Assim, a FTMF torna-se:
ŷ(s) kp s + ki
= 3 2
r̂(s) s + s + (kp − 2)s + ki
Exemplo 4
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
ki
Resp. Transitória
3
Sist. s 1 kp − 2
Realimentados
s2 1 ki
Aula 10
Região factı́vel
Aula 11 s1 kp − 2 − ki
Aula 12
Aula 13
s0 ki
Lugar das Raı́zes
ki = 0
Ou seja:
Resp. em kp
Frequência
Controle PID
ki > 0
kp > ki + 2
Controle Moderno
kp = ki + 2
Sistemas Incertos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ Sistema Incerto: É aquele que tem ao menos um
Modelagem parâmetro cujo valor não é precisamente conhecido;
Resp. Transitória ▶ Fatores que conduzem a incertezas: limitações de
Sist.
Realimentados
fabricação, envelhecimento, variações na operação,
Aula 10
Aula 11
entre outros;
Aula 12
Aula 13
▶ Considere a FT de um sistema incerto:
Lugar das Raı́zes m
X
Resp. em bj sj
Frequência
(
j=0 0 < ai ≤ ai ≤ āi
Controle PID ĝ(s) = n , com ¯ (4.39)
Controle Moderno
X
i bj ≤ bj ≤ b̄j
ai s ¯
i=0
Prof. Victor L.
Yoshimura
Controle Moderno
KD3 (s) = a0 + ā1 s + ā2 s2 + a3 s3 + a4 s4 + ā5 s5 + . . .
¯ ¯ ¯
(4.40c)
KD4 (s) = ā0 + a1 s + a2 s2 + ā3 s3 + ā4 s4 + a5 s5 + . . .
¯ ¯ ¯
(4.40d)
Exemplo 5 (Prova P2 2020-1)
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. em
estabilidade robusta.
Frequência
Controle PID
ŷ(s) ĉ(s)ĝ(s) ki
= = 4
Controle Moderno r̂(s) 1 + ĉ(s)ĝ(s) s + cs + bs2 + as + ki
3
Exemplo 5 (Prova P2 2020-1)
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
s4 1 b ki
Modelagem
Resp. Transitória
s3 c a
Sist.
s2 bc − a cki
Realimentados s1 a(bc − a) − c2 ki
Aula 10
Aula 11 s0 ki
Aula 12
Aula 13
Resp. em
▶ Para a linha s2 , a condição mais desfavorável ocorre em bc − ā = 1;
Frequência
¯¯
a(bc − a)
▶ Para a linha s1 , deve-se ter ki < , o que leva a
Controle PID c2
Controle Moderno
a(b̄c̄ − a) ā(bc − ā) ā(b̄c − ā) a(bc̄ − a)
ki < min ¯ ¯ , ¯¯ 2 , ¯2 ,¯ ¯ 2 ¯ = 0,75
c̄2 c c c̄
¯ ¯
▶ A linha s0 apenas exige ki > 0.
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
▶ A partir das condições (3.26)-(3.29), é possı́vel escolher
Sist.
Realimentados posições de polos adequadas para a FTMF;
Aula 10
Aula 11
Aula 12
▶ O processo descrito é chamado alocação de polos;
Aula 13
Prof. Victor L.
Yoshimura
r̂(s) ε̂(s) ŷ(s)
Introdução Σ ĝ(s)
Modelagem −
Resp. Transitória
Sist.
Realimentados
Aula 10
Escrevamos as FTs da planta e do compensador como:
Aula 11
Aula 12
Pm
g j
n̂g (s) j=0 βj s
(4.41)
Aula 13
ĝ(s) = = Png
Lugar das Raı́zes dˆg (s) i=0 αi s
i
Resp. em
Frequência Pm
c j
Controle PID n̂c (s) j=0 bj s
ĉ(s) = = Pnc (4.42)
Controle Moderno dˆc (s) i=0 ai s
i
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Resp. Transitória n̂c (s)n̂g (s) + dˆc (s)dˆg (s) = q̂(s) (4.45)
Sist.
Realimentados
Aula 10
que é chamada equação diofantina.
Aula 11
Aula 12
Aula 13
Observação
Lugar das Raı́zes ▶ Note que n̂c (s) e dˆc (s) são as incógnitas de (4.45);
Resp. em
Frequência ▶ Sob quais circunstâncias existe solução para (4.45)?
Controle PID ▶ A solução, se existe, é única?
Controle Moderno
▶ Como encontrar a solução?
▶ Como não há produto de incógnitas, a dependência
entre estas e os coeficientes de dˆg e n̂g é linear.
Solução da Equação Diofantina
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Observando (4.41) e (4.42), note que:
Introdução
Observação
Como αi , βj são dados e ai , bj são as incógnitas, a Equação
Diofantina é um problema linear!
Solução da Equação Diofantina
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ Dois polinômios são iguais se seus respectivos
Modelagem coeficientes o são.
Resp. Transitória ▶ Para os termos de grau zero de (4.45):
Sist.
Realimentados
Aula 10 α0 a0 + β0 b0 = q0
Aula 11
Aula 12
Aula 13 ▶ Para o termo de primeiro grau:
Lugar das Raı́zes
Resp. em α1 a0 + α0 a1 + β1 b0 + β0 b1 = q1
Frequência
Controle PID
▶ Para o termo de i-ésimo grau:
Controle Moderno
i
X
(αi−k ak + βi−k bk ) = qi
k=0
A Matriz de Sylvester
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
▶ De (4.45), surge um sistema de equações lineares;
Introdução
Modelagem
▶ Convém utilizar a matriz de Sylvester.
Resp. Transitória ▶ Defina o bloco de Sylvester de dˆg :
Sist.
Realimentados
Aula 10
Aula 11
α0 0 0 ... 0 0 ... 0
Aula 12
α1 α0 0 ... 0 0 ... 0
Aula 13
α2 α1 α0 ... 0 0 ... 0
Lugar das Raı́zes
.. .. .. .. .. .. .. ..
Resp. em . . . . . . . .
Frequência SD =
α αng −1 αng −2 . . . α0 0 . . . 0
Controle PID ng
Controle Moderno 0
αng αng −1 . . . α1 α0 . . . 0
. .. .. . .. . . ..
. ..
. ..
. . . . . .
0 0 0 . . . 0 0 . . . αng
∈ R(nc +ng +1)×(nc +1) (4.46)
Solução da Equação Diofantina
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura ▶ O bloco do numerador é feito similarmente e com
Introdução dimensão SN ∈ R(nc +ng +1)×(mc +1) ;
Modelagem ▶ Escrevamos:
Resp. Transitória
h i′
Sist. a = a0 a1 . . . anc
Realimentados
Aula 10
Aula 11 h i′
Aula 12
Aula 13
b = b0 b1 . . . bmc
Lugar das Raı́zes h i′
Resp. em q = q0 q1 . . . qnc +ng
Frequência
h i
Controle PID
S = SD SN
Controle Moderno
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
▶ Não se garante que (4.50) atenderá às especificações de
Introdução
projeto;
Modelagem
Resp. Transitória
▶ A solução encontrada deve ser testada;
Sist. ▶ É comum ter que se projetar um controlador que
Realimentados
Aula 10 imponha um tipo mı́nimo;
Aula 11
Aula 12 ▶ Procedimento simples:
Aula 13
ĝ
Lugar das Raı́zes i Troque ĝ por g̃ˆ = k , onde k leva ao tipo desejado;
s
Resp. em
Frequência
ii Solucione a Equação Diofantina, obtendo c̃ˆ;
c̃ˆ
Controle PID iii A solução verdadeira será k .
Controle Moderno s
Observação
Principal desvantagem da Equação Diofantina: os polos
desejados devem ser exatamente especificados, o que leva a
pouca flexibilidade do método.
Exemplo: Tanques Comunicantes
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura Devemos inserir um polo na origem. Logo, usaremos:
Introdução
ĝ(s) 0,4
Modelagem = 3
s 4s + 4,8s2 + s
Resp. Transitória
Sist.
Realimentados
Para garantir solução da Equação Diofantina, façamos
Aula 10 mc = 2. Assim:
Aula 11
Aula 12
Aula 13
b2 s2 + b1 s + b0
Lugar das Raı́zes sĉ(s) = (o polo extra está em ĝ/s)
a1 s + a0
Resp. em
Frequência
Os blocos de Sylvester de ĝ/s são:
Controle PID
Controle Moderno
0 0 0,4 0 0
1 0 0 0,4 0
SD = 4,8 1 SN = 0 0 0,4
4 4,8 0 0 0
0 4 0 0 0
Exemplo: Tanques Comunicantes
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
▶ Para Mp ≤ 0,1, deve-se ter ζ ≥ 0,6;
Introdução
▶ Para ts < 10s, deve-se ter ζωn > 0,4 (critério 2%);
Modelagem
Resp. Transitória
▶ Façamos ζ = 0,65 e ωn = 0,7.
Sist. ▶ Polos dominantes: −0,455 ± j0,532
Realimentados
Aula 10
Aula 11
▶ Os demais polos serão reais, com σ dez vezes maior.
Aula 12
Aula 13
Assim:
Lugar das Raı́zes
Controle PID
= s4 + 10,9s3 + 34,6s2 + 27,7s + 12,3
Controle Moderno
▶ Aplicando (4.47):
" #
a h i′
= 2,425 0,25 30,75 63,19 56,78
b
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Introdução
Modelagem Observação
Resp. Transitória Os polos de malha fechada são chamados raı́zes do sistema.
Sist.
Realimentados ▶ A localização das raı́zes determina:
Lugar das Raı́zes ▶ a estabilidade (aula 12), e;
Aula 14
Aula 15
▶ o desempenho do sistema (aulas 08 e 09).
Resp. em ▶ O cálculo das raı́zes pode ser feito numericamente;
Frequência
Prof. Victor L.
Yoshimura
Introdução
Considere o sistema de controle e (2.4):
Modelagem
Resp. Transitória
r̂(s) ŷ(s)
Σ ĝ(s) ŷ(s) ĝ(s)
Sist.
Realimentados − =
r̂(s) 1 + ĝ(s)ĥ(s)
Lugar das Raı́zes ĥ(s)
Aula 14
Aula 15
▶ Considere, em ĝ(s), a inclusão do compensador.
Resp. em
Frequência
▶ Para encontrar o LR, deve-se resolver:
Controle PID
Controle Moderno
ĝ(s)ĥ(s) = −1 (5.1)
Importante!
O método LR, determina os polos da FTMF, a partir de
informações da FTMA. Logo, não é necessário determinar a
FTMF!
Conceitos Preliminares II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Como as raı́zes são complexas, pode-se reescrever (5.1):
Introdução
Modelagem
|ĝ(s)ĥ(s)| = 1 (5.2a)
Resp. Transitória
Sist.
arg(ĝ(s)ĥ(s)) = 180◦ (2k + 1), k ∈ N (5.2b)
Realimentados
Resp. em m m
Frequência
Y Y
(s − zj ) |s − zj |ejψj
Controle PID
j=1 j=1
Controle Moderno ĝ(s)ĥ(s) = K n =K n (5.3)
Y Y
jφi
(s − pi ) |s − pi |e
i=1 i=1
Importante!
Como determinar as raı́zes à medida em que K ∈ (0, ∞) é
alterado?
Conceitos Preliminares III
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura As condições (5.2) são reescritas com (5.3):
Introdução
n
, m
Modelagem
(5.4a)
Y Y
K= |s − pi | |s − zj |
Resp. Transitória
i=1 j=1
Sist.
Realimentados
m n
Lugar das Raı́zes
φi = 180◦ (2k + 1), k ∈ N (5.4b)
X X
Aula 14
ψj −
Aula 15 j=1 i=1
Resp. em
Frequência
Controle PID
Observação
Controle Moderno ▶ Com a forma racional para a FTMA, se reescreve (5.1):
ˆ + K n̂(s) = 0
d(s) (5.5)
ˆ = n ≥ m = gr(n̂);
▶ É evidente que gr(d)
▶ Cada polo gera uma raiz e seu traçado é dito um ramo
do LR.
1ª Regra: Simetria, Inı́cio e Fim do Traçado
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Lema (Simetria)
Modelagem
Resp. Transitória
O eixo real é um eixo de simetria do LR
Sist.
Realimentados Demonstração.
Lugar das Raı́zes Basta notar que a equação caracterı́stica tem coeficientes
Aula 14
Aula 15 todos reais, logo, as raı́zes complexas aparecem em pares
Resp. em conjugados.
Frequência
Controle PID
Lema (Inı́cio e Fim do Traçado)
Controle Moderno
Os ramos do LR iniciam-se em seu polo e terminam em um
zero da FTMA.
Demonstração.
Basta fazer k → 0 e k → ∞ em (5.5).
2ª Regra: LR sobre o Eixo Real
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Introdução Observação
Modelagem
Cada ramo LR inicia em um polo e termina em um zero. E
Resp. Transitória
se n > m?
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Introdução Observação
Modelagem Se existe LR no eixo real e fora dele, de onde partem tais
Resp. Transitória
ramos?
Sist.
Realimentados
Controle Moderno
Demonstração.
Note que sb é raiz múltipla de (5.5), para um K > 0 bem
escolhido. Logo, é raiz da derivada em relação a s de (5.5),
para o mesmo K. Combinando estes resultados, chega-se a
(5.7).
Sumário das Regras para o Traçado do LR
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura Assı́ntotas com (5.6):
K
Introdução ĝ(s)ĥ(s) = 0−1−2
s(s + 1)(s + 2) σa = = −1
Modelagem 3−0
Im(s)
Resp. Transitória K →∞ θak = (2k − 1)60◦
Sist.
Realimentados K=6 Quebras com (5.7) e (5.4a):
s = j1,41
Lugar das Raı́zes
sb1 = −0,42
Aula 14
Aula 15
3s2 +6s+2 = 0
sb2 = −1,58
Resp. em
Frequência
K →∞
σa
Controle PID K = |sb1 (sb1 + 1)(sb1 + 2)|
Controle Moderno Re(s) = 0,38
Cruzamento de Im(s):
K = 0,38
s = −0,42 s3 1 2
√
s2 3 K (s = ±j 2)
s1 6−K (K = 6)
s0 K
K →∞
A Região Ω de Desempenho Garantido
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura ▶ Overshoot, de (3.28):
Introdução
!
ζ
Modelagem ζ = 0,5 exp −π p < Mmax
Resp. Transitória ζ = 0,71 Im(s)
1 − ζ2
(5.8)
Sist.
Realimentados
Observação
ζ = 0,71
Desempenho garantido: o
ζ = 0,5
desempenho real é, pelo menos,
o das especificações!
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Introdução
▶ O projeto de sistemas de controle deve atender a certas
Modelagem
especificações:
Resp. Transitória
▶ Precisão (erro);
Sist.
Realimentados
▶ Estabilidade relativa;
Lugar das Raı́zes
▶ Velocidade de resposta.
Aula 14
Aula 15
▶ Em muitos casos, um compensador P é suficiente;
Resp. em
Frequência
▶ Em outros, o traçado do LR deve ser alterado;
Controle PID ▶ A planta é considerada dada. Esta alteração é
Controle Moderno responsabilidade do compensador!
Observação
O conhecimento do traçado do LR facilita o projeto do
compensador!
Avanço e Atraso de Fase
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
▶ O compensador pode assumir qualquer estrutura: P, PI,
Introdução
PD, etc.
Modelagem
▶ Outra estrutura compensadora: avanço (atraso) de fase.
Resp. Transitória
Sist.
Realimentados C1 C2
Lugar das Raı́zes
Aula 14
Aula 15
Resp. em ε R1 R2 R3 R4 u
Frequência
Controle PID
− −
Controle Moderno
+ +
s−z
ĉ(s) = kc (5.10)
s−p
Prof. Victor L.
Yoshimura
Sist.
“mais estável”);
Realimentados
▶ ĉ(0) < kc , os erros estáticos tendem a aumentar (vide
Lugar das Raı́zes
Aula 14
aula 11);
Aula 15
▶ O compensador atraso aproxima-se do PI;
Resp. em
Frequência
▶ Faz o contrário das caracterı́sticas do PD;
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Reescreva (5.10) na forma avanço (0 < α < 1):
Modelagem
Resp. Transitória 1
Sist.
s+
Realimentados ĉ(s) = kc T (5.11)
1
Lugar das Raı́zes s+
Aula 14 αT
Aula 15
Prof. Victor L.
Yoshimura
r̂(s) 4
ŷ(s)
Introdução Σ ĉ(s) s(s+2)
Modelagem −
Resp. Transitória
Sist.
Realimentados
▶ Especificações: εp (∞) = 0, ts ≤ 2s (ε ≤ 2%) e
Lugar das Raı́zes
Aula 14 Mp ≤ 16%;
Aula 15
Resp. em
▶ A especificação de erro já foi atendida pela planta;
Frequência
▶ Não é possı́vel obter as raı́zes para as especificações
Controle PID
Controle Moderno
dadas com compensador P (Por quê?);
▶ Escolha das raı́zes desejadas:
(3.29)
ts ≤ 2s (ε ≤ 2%) =⇒ σ ≥ 2
(3.28)
Mp = 16% =⇒ ζ = 0,5 =⇒ β = 60◦
√
sd = −2 ± j2 tg 60◦ = −2 ± j2 3
Exemplo de Projeto
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Critério angular (15.4b):
Introdução
4
Modelagem
Resp. Transitória
arg = − arg(sd ) − arg(sd + 2) = −210◦
s(s + 2) s=sd
Sist.
Realimentados
Resp. em s + 2,5
Frequência arg kc = 30◦ ⇒ p = −4,7
Controle PID s−p s=sd
Controle Moderno
Para que a raiz seja escolhida, deve-se ajustar o ganho com
(5.4a):
|sd + 4,7|
K = 4kc = |sd ||sd + 2| = 17,4
|sd + 2,5|
Compensação em Atraso de Fase
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Reescreva (5.10) na forma avanço (α > 1, constantes de
Resp. Transitória tempo):
Ts + 1
Sist.
Realimentados
ĉ(s) = αkc (5.12)
αT s + 1
Lugar das Raı́zes
Aula 14 Procedimento para projeto:
Aula 15
Resp. em
1. Estabeleça a localização das raı́zes dominantes;
Frequência
2. Determine o valor a ser adicionado à constante de erro;
Controle PID
Controle Moderno 3. Escolha o polo e o zero para esta adição, mas próximos
entre si e da origem (assim, não alterarão o LR);
4. Desenhe o novo LR;
5. Determine kc , com (5.4a), para as raı́zes desejadas.
Exemplo de Projeto
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
r̂(s) 1,06
ŷ(s)
Modelagem Σ ĉ(s) s(s+1)(s+2)
Resp. Transitória −
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Introdução
▶ As raı́zes dominantes atuais, sd = −0,3307 ± j0,5864,
Modelagem
já atendem as exigências de acomodação e overshoot.
Resp. Transitória
Sist.
▶ A forma (5.12) é útil neste caso, pois já se tem
Realimentados
αkc = 10;
Lugar das Raı́zes
1 1
Aula 14
▶ Escolhe-se = 0,05 e = 0,005
Aula 15
T αT
Resp. em
Frequência
▶ Para o ajuste de ganho, usa-se (5.4a):
Controle PID 1 s(s + 0,005)(s + 1)(s + 2)
K= = 0,9656
Controle Moderno 1,06 s + 0,05 s=sd
Observação
Qual o preço pago pela melhoria do erro estático de
velocidade?
Compensação em Avanço e Atraso de Fase
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Introdução
Modelagem r̂(s) 4
ŷ(s)
Resp. Transitória Σ ĉ(s) s(s+0,5)
Sist.
−
Realimentados
Controle PID
▶ A constante de erro de velocidade, tampouco.
Controle Moderno
Importante!
Para atender às especificações, pode ser interessante
posicionar as raı́zes desejadas “mais para dentro” da região
Ω. Aqui, faremos sd = −3 ± j4.
Exemplo de Projeto
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ Etapa avanço:
Modelagem Posicionando o zero do compensador em −1, note que a
Resp. Transitória deficiência angular é de 48◦ . Logo:
Sist.
Realimentados
arg(sd − p) = 48◦ ⇒ p = −6,6
Lugar das Raı́zes
Aula 14
Aula 15
Do critério de magnitude (5.4a):
Resp. em
Frequência s+1 4
kc = 1 ⇒ kc = 7,1
Controle PID s + 6,6 s(s + 0,5)
Controle Moderno
▶ Etapa atraso:
Note que Kv = 8,6 com a inserção do avanço. Usando (5.12)
(sem kc ), devemos fazer α ≈ 9,3. Usando a recomendação de que
o polo e o zero devem ser próximos entre si e da origem,
escolhe-se T = 10.
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Aula 19
Tópico 6 - Resposta em Frequência
Aula 20 Aula 16: Diagramas de Bode
Aula 17: O Critério de Nyquist
Aula 18: Parâmetros de Desempenho
Controle PID Aula 19: Resposta em Malha Fechada
Aula 20: Projeto de Compensadores
Controle Moderno
Tópico 7 - Controle PID
Aula 21: Controle PID
Aula 19
Tópico 6 - Resposta em Frequência
Aula 20 Aula 16: Diagramas de Bode
Aula 17: O Critério de Nyquist
Aula 18: Parâmetros de Desempenho
Controle PID Aula 19: Resposta em Malha Fechada
Aula 20: Projeto de Compensadores
Controle Moderno
Tópico 7 - Controle PID
Aula 21: Controle PID
Prof. Victor L.
Yoshimura
Considere um SLIT-C-SISO estável, submetido a uma
Introdução excitação da forma:
Modelagem
Resp. em n̂(s)
Frequência
Aula 16
ĝ(s) = n (6.2)
Y
Aula 17
Aula 18
s − si
Aula 19 i=1
Aula 20
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Resp. Transitória
ω Umax
a = Umax ĝ(s) (s + jω) =− ĝ(−jω)
Sist. s2 + ω2 s=−jω j2
Realimentados
(6.5)
Lugar das Raı́zes
Como a é o conjugado de a:
Resp. em
Frequência
Umax
(6.6)
Aula 16
Aula 17 a= ĝ(jω)
Aula 18 j2
Aula 19
Defina-se
Aula 20
Controle PID
Controle Moderno
ĝ(jω) = |ĝ(jω)|ejφ (6.7a)
Im(ĝ(jω))
φ = arctg (6.7b)
Re(ĝ(jω))
Avanço e Atraso na Resposta Senoidal IV
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Aplicando (6.5), (6.6) e (6.7) a (6.4), com
Modelagem
t ≫ max {− Re(si )−1 }:
Resp. Transitória i=1:n
Sist.
Realimentados Umax Umax
Lugar das Raı́zes
y(t) = |ĝ(jω)|ejφ ejωt − |ĝ(jω)|e−jφ e−jωt
j2 j2
Resp. em
Frequência ej(ωt+φ) − e−j(ωt+φ)
Aula 16 = Umax |ĝ(jω)|
Aula 17 j2
Aula 18
Aula 19 = Umax |ĝ(jω)| sen(ωt + φ) (6.8)
Aula 20
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Sist.
▶ Resposta em Frequência: Resposta em regime a uma
Realimentados
entrada senoidal.
Lugar das Raı́zes
▶ A representação gráfica da resposta em frequência,
Resp. em
Frequência ĝ(jω), pode ser feita por:
Aula 16
Aula 17
▶ Diagramas de Bode (logarı́tmicos);
Aula 18 ▶ Diagramas de Nyquist (polares);
Aula 19
Aula 20
▶ Diagramas de Black-Nichols (log magnitude versus
Controle PID fase).
Controle Moderno
Diagramas de Bode
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
A resposta em frequência pode ser reescrita como:
Introdução
Prof. Victor L.
Yoshimura
Introdução Observação
Modelagem
Diagramas de Bode: Dois gráficos, ||ĝ|| × log ω e
Resp. Transitória
arg(ĝ) × log ω.
Sist.
Realimentados O traçado dos diagramas de Bode necessita de papel
Lugar das Raı́zes mono-log!
Resp. em
Frequência
▶ Como a FT é uma função racional, determinamos
Aula 16 apenas 4 termos:
Aula 17
Aula 18 ▶ Termo constante;
Fator integrador ou derivador;
Aula 19
Aula 20
▶
Controle PID
▶ Fator de 1ª ordem;
Controle Moderno
▶ Fator de 2ª ordem.
▶ Diagramas exatos: Difı́cil construção, não apresenta
erro;
▶ Diagramas assintóticos: Simples, porém com erros.
Termos Constante e Integrador/Derivador
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Ganho constante (K > 0):
Introdução
Controle PID
Controle Moderno
Observação
Termo constante: ganho constante e fase nula para
ω ∈ [0, ∞).
Termo integrador/derivador: variação de ±20dB por década,
desde a origem e fase constante em ±90◦ .
!±1
Controle e ω
Servomecanismos
Termo de 1ª Ordem 1 + j
Prof. Victor L.
Yoshimura
ωn
Introdução Determinação de ganho e de fase:
Modelagem
±1
ω ω
Resp. Transitória
Sist.
|| · || = 20 log 1+j = ±20 log 1 + j
Realimentados
ωn ωn
Lugar das Raı́zes
2 !
ω
Resp. em = ±10 log 1 + (6.15a)
Frequência ωn
Aula 16
Aula 17
ω
Aula 18
Aula 19 arg(·) = ± arctg (6.15b)
Aula 20 ωn
Controle PID
Regras assintóticas:
Controle Moderno
▶ Ganho: Nulo até ωn , a partir de então ±20dB por
década.
▶ Fase: Nula até uma década antes de ωn . ±90◦ uma
década após ωn . Variação linear entre estes valores.
!±1
Controle e ω
Servomecanismos
Termo de 1ª Ordem 1 + j
Prof. Victor L.
Yoshimura
ωn
Introdução
Modelagem
Resp. Transitória
ω
Sist. log ||ĝ||[dB] arg(ĝ)[◦ ]
Realimentados ωn
Lugar das Raı́zes
−3 ±4,34 · 10−6 ±0,057
Resp. em
Frequência −2 ±4,34 · 10−4 ±0,57
Aula 16
Aula 17
−1 ±4,32 · 10−2 ±5,7
Aula 18
Aula 19
0 ±3,01 ±45
Aula 20
1 ±20,04 ±84,29
Controle PID
2 ±40 ±89,43
Controle Moderno
3 ±60 ±89,94
!±1
Controle e ω
Servomecanismos
Termo de 1ª Ordem 1 + j
Prof. Victor L.
Yoshimura
ωn
Introdução 0
Modelagem
−10
Magnitude (dB)
Resp. Transitória
−20
Sist.
Realimentados
−30
Lugar das Raı́zes
−40
Resp. em −3 −2 −1 0 1
10 10 10 10 10
Frequência
Frequency (Hz)
Aula 16
Aula 17
Aula 18
0
Aula 19
Aula 20 −20
Phase (degree)
Controle PID
−40
Controle Moderno
−60
−80
−3 −2 −1 0 1
10 10 10 10 10
Frequency (Hz)
Controle e
2 !±1
ω ω
Servomecanismos
Termo de 2ª Ordem 1 + 2ζ j + j
Prof. Victor L. ωn ωn
Yoshimura
Introdução
Modelagem
Resp. Transitória
Observação
Sist. Interessa apenas o caso subamortecido. Por quê?
Realimentados
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Considere a FT senoidal:
Introdução 1 1
ĝ(jω) = = 2 (6.19)
Modelagem ˆ
d(jω)
ω
ω
Resp. Transitória
1 + 2ζ j + j
ωn ωn
Sist.
Realimentados
Controle Moderno
q −1
Mr = ||ĝ(jωr )|| = 2ζ 1 − ζ 2 (6.20b)
Observação
1
Não há ressonância para ζ > √ ≈ 0,707.
2
ζ = 0,1, ζ = 0,3, ζ = 0,5 e ζ = 0,707
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução 10
Modelagem 0
Magnitude (dB)
Resp. Transitória −10
Sist. −20
Realimentados
−30
Lugar das Raı́zes
−2 −2 −1 −1 0
10 5.0x10 10 5.0x10 10
Resp. em
Frequência Frequency (Hz)
Aula 16
Aula 17
0
Aula 18
Aula 19
Aula 20 −50
Phase (degree)
Controle PID
−100
Controle Moderno
−150
−2 −2 −1 −1 0
10 5.0x10 10 5.0x10 10
Frequency (Hz)
Sistemas de Fase Mı́nima
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Importante!
Resp. Transitória Considere um SLIT-C estável. Diz-se que este é de fase
Sist. mı́nima se seus zeros estiverem no semi-plano complexo
Realimentados
esquerdo.
Lugar das Raı́zes
Resp. em
Frequência Importante!
Aula 16
Aula 17
Sejam ĝ(s) e ĥ(s) SLIT-Cs estáveis, ĝ(s) de fase mı́nima, e
Aula 18
Aula 19
||ĝ(jω)|| = ||ĥ(jω)||, ∀ω > 0. Então a excursão de fase de ĝ
Aula 20
é menor ou igual à de ĥ.
Controle PID
Prof. Victor L.
Yoshimura
Considere o sistema:
Introdução
Observação
Note que, com o TVF, s → 0 ⇒ ω → 0.
Curva de Ganho e Erros Estáticos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. em Kv
Frequência ||ĝ(jω)|| = 20 log , N =1 (6.23)
Aula 16 jω
Aula 17
Aula 18
Aula 19 Em (4.23), ||ĝ(jKv )|| = 0dB. Partindo de (4.26):
Aula 20
Controle PID Ka
Controle Moderno ||ĝ(jω)|| = 20 log , N =2 (6.24)
(jω)2
√
Similarmente, para ||ĝ(j Ka )|| = 0, obtém-se a constante
de erro.
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Aula 19
Tópico 6 - Resposta em Frequência
Aula 20 Aula 16: Diagramas de Bode
Aula 17: O Critério de Nyquist
Aula 18: Parâmetros de Desempenho
Controle PID Aula 19: Resposta em Malha Fechada
Aula 20: Projeto de Compensadores
Controle Moderno
Tópico 7 - Controle PID
Aula 21: Controle PID
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Controle PID
▶ A tı́tulo ilustrativo, faremos o estudo dos termos
Controle Moderno possı́veis em uma FT.
Termo Constante e Termo Integrador/Derivador
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução Im(ĝ)
Modelagem ω→∞ Termo constante:
Resp. Transitória
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
O atraso de transporte, no domı́nio da frequência é escrito
Modelagem
como:
Resp. Transitória
Sist.
Realimentados ĝ(jω) = e−jωT = cos ωT − j sen ωT (6.27)
Lugar das Raı́zes
Resp. em
Frequência
Observação
Aula 16
Aula 17
O DN do atraso de transporte é a circunferência unitária
Aula 18
Aula 19
centrada na origem, girando no sentido horário.
Aula 20
Resp. Transitória 1
Sist. Im(ĝ) ĝ(jω) = −j e−jωT
Realimentados
ω
π
Lugar das Raı́zes
1 −j ωT + 2
Resp. em = e
Frequência ω
Aula 16
Aula 17
Aula 18
Assim:
Aula 19
Re(ĝ)
Aula 20 ▶ |ĝ(j0)| → ∞;
Controle PID
▶ arg(ĝ(j0)) → −90◦ ;
Controle Moderno ωT = 2π
▶ |ĝ(j∞)| → 0
ωT ≈ 0
▶ O diagrama “gira” no
sentido horário.
Regras Gerais do Diagrama de Nyquist
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
dos eixos;
Aula 19
Aula 20
Prof. Victor L.
Yoshimura
Seja F : C → C uma função complexa. Então:
Introdução
1º Se l ⊂ D(F ) é uma curva fechada que não contém
Modelagem
singularidades de F , então F (l) ⊂ CD(F ) é fechada;
Resp. Transitória
N =Z −P (6.28)
Funções Complexas e Seus Mapeamentos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Im(s) Im(F (s))
Resp. Transitória
Sist.
Realimentados
Resp. em
Frequência l F (l)
F
Aula 16
Aula 17 Re(s) =⇒ Re(F (s))
Aula 18
Aula 19
Aula 20
Controle PID
Controle Moderno
Mapeamento e Estabilidade de SLIT-Cs
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Considere o sistema realimentado, (2.4), então:
Introdução F (s) = 1 + ĝ(s)ĥ(s) (6.29)
Modelagem
Prof. Victor L.
Yoshimura
Teorema (Critério de Estabilidade de Nyquist)
Introdução
Se a FTMA tem k polos no semi-plano direito (no domı́nio) e
Modelagem
(6.30) vale, então o DN dessa FTMA, com ω ∈ (−∞, ∞), deverá
Resp. Transitória
circundar k vezes o ponto −1 + j0 no sentido anti-horário para
Sist.
Realimentados
que o sistema seja estável.
Lugar das Raı́zes
Importante!
Resp. em
Frequência Em outras palavras, adaptando (6.28):
Aula 16
(6.31)
Aula 17
Aula 18 Pf = N + Pa
Aula 19
Aula 20
▶ N é contado positivo no sentido horário e em torno de
Controle PID
−1 + j0;
Controle Moderno
▶ O critério vale mesmo com polos ou zeros da FTMA sobre a
origem.
▶ Se o DN passar por −1 + j0, há raiz sobre o eixo imaginário
(E daı́?);
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Aula 19
Tópico 6 - Resposta em Frequência
Aula 20 Aula 16: Diagramas de Bode
Aula 17: O Critério de Nyquist
Aula 18: Parâmetros de Desempenho
Controle PID Aula 19: Resposta em Malha Fechada
Aula 20: Projeto de Compensadores
Controle Moderno
Tópico 7 - Controle PID
Aula 21: Controle PID
Prof. Victor L.
Yoshimura ▶ Pode ser realizada com a mudança de variável proposta
Introdução em (4.38);
Modelagem ▶ Para realimentação unitária tome a FTMA estável ĝ(s);
Resp. Transitória
▶ Considere o mapeamento:
Sist.
Realimentados
Im(s) Im(ĝ)
Lugar das Raı́zes
Resp. em
Frequência
Aula 16
Aula 17
Aula 18
G
Aula 19
Re(s) =⇒ Re(ĝ)
Aula 20 −1 + j0
Controle PID
Controle Moderno
Importante!
O diagrama de Nyquist não pode abraçar o ponto −1 + j0!
Margens de Fase e de Ganho
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Observação
Para estabilidade absoluta, deve-se ter M F, M G > 0.
Margens de Fase e de Ganho
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. Transitória j
Sist. ωf
Realimentados 0
ω
Lugar das Raı́zes −M G
MG
Resp. em arg(ĝ) −1 10 20 1
Frequência
Aula 16 Re(ĝ)
Aula 17
Aula 18 −180 + M F MF
Aula 19
Aula 20 −180
ωg ω
Controle PID −j
Controle Moderno
Observação
O critério de estabilidade com M F e M G supõe sistema de
fase mı́nima!
Correlação entre M F , Mr e ζ
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura M F/100, Mr /5
1,0
Introdução
Modelagem 0,9
Sist. 0,7
Realimentados
0,6
Lugar das Raı́zes
0,5
Resp. em
0,4
Frequência
Aula 16 0,3
Aula 17
Aula 18 0,2
Aula 19
Aula 20
0,1
Controle PID ζ
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0
Controle Moderno
Prof. Victor L.
Yoshimura
M F/100, Mr /5
Introdução
1,0
Modelagem
0,9
Resp. Transitória
0,8
Sist.
0,7
Realimentados
Resp. em 0,5
Frequência 0,4
Aula 16
Aula 17 0,3
Aula 18
Aula 19 0,2
Aula 20
0,1
Controle PID
ζ
Controle Moderno 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0
M F ≈ 100ζ (6.34)
Frequência de Corte, Banda Passante e Dinâmica
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
▶ A frequência de corte é a frequência ωb , tal que se
Introdução
ω > ωb , então:
Modelagem
Resp. Transitória
||ĝ(jω)|| − ||ĝ(j0)|| ≤ −3dB (6.35)
Sist.
Realimentados
Exemplo
Aula 19
Aula 20
Controle PID
Faça a análise da planta em realimentação unitária
Controle Moderno
k
ĝ(s) =
s(s + 1)(s + 5)
para k = 10 e k = 100.
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Aula 19
Tópico 6 - Resposta em Frequência
Aula 20 Aula 16: Diagramas de Bode
Aula 17: O Critério de Nyquist
Aula 18: Parâmetros de Desempenho
Controle PID Aula 19: Resposta em Malha Fechada
Aula 20: Projeto de Compensadores
Controle Moderno
Tópico 7 - Controle PID
Aula 21: Controle PID
Prof. Victor L.
Yoshimura Considere sistema em malha fechada e realimentação unitária:
→
Introdução Observe que OA = ĝ(jωa ) e que
→
Modelagem
P A = 1 + ĝ(jωa ). Logo:
Resp. Transitória Im(ĝ)
→
Sist.
Realimentados ||OA|| ĝ(jωa )
→ = (6.36a)
Lugar das Raı́zes
||P A|| 1 + ĝ(jωa )
Resp. em P
Frequência
Re(ĝ)
φ ĝ(jωa )
Aula 16
−1 + j0 θ arg =φ−θ
Aula 17
Aula 18
1 + ĝ(jωa )
Aula 19 (6.36b)
Aula 20
Defina:
Controle PID A ω
a
Controle Moderno ŷ(jω)
= M (jω)ejα(jω) (6.37)
r̂(jω)
Observação
Deseja-se determinar os lugares geométricos de ganho e de fase
constantes, para a FTMF.
Lugares de Magnitude Constante I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Escreva ĝ(jω) = X(jω) + jY (jω) = X + jY . Tendo em
Resp. Transitória
vista (6.36a), a magnitude, M , é dada por:
Sist.
Realimentados |X + jY |2 X2 + Y 2
M2 = =
Lugar das Raı́zes |1 + X + jY |2 (1 + X)2 + Y 2
Resp. em
Frequência
Aula 16
ou, ainda:
Aula 17
Aula 18
Aula 19 (1 − M 2 )X 2 − 2M 2 X − M 2 + (1 − M 2 )Y 2 = 0 (6.38)
Aula 20
Controle PID
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Resp. em tg A − tg B
Frequência Definindo N = tg α e usando tg(A − B) = :
Aula 16 1 + tg A tg B
Aula 17
Aula 18
Aula 19 Y Y
Aula 20 −
Controle PID N= X 1 + X ⇒ X2 + X + Y 2 − 1 Y = 0
Y Y N
Controle Moderno 1+
X 1+X
Lugares de Fase Constante II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem 1 1
Somando + a ambos os membros:
Resp. Transitória 4 (2N )2
Sist.
Realimentados 2 2
1 1 1 1
Lugar das Raı́zes X+ + Y − = + (6.40)
Resp. em
2 2N 4 (2N )2
Frequência
Aula 16
Aula 17
Aula 18
Observação
Aula 19
Aula 20 i O N -cı́rculo de α e o de α ± 180◦ são o mesmo.
Controle PID
ii Os M -cı́rculos e N -cı́rculos, no diagrama de Nyquist,
Controle Moderno
dão origem à carta de Hall.
Diagrama de Black e a Carta de Nichols
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura ▶ Diagrama de Black: É a representação gráfica de
Introdução ||ĝ|| × arg(ĝ);
Modelagem ▶ Frequência implı́cita, mas MF e MG são facilmente
Resp. Transitória
reconhecı́veis;
Sist.
Realimentados ▶ Carta de Nichols: Plotagem dos M -cı́rculos e
Lugar das Raı́zes N -cı́rculos no diagrama de Black;
Resp. em
Frequência ▶ Esses cı́rculos não aparecerão como tais no diagrama de
Aula 16
Aula 17 Black. (Por quê?)
Aula 18
Exemplo
Aula 19
Aula 20
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Em muitos casos, é difı́cil estabelecer um modelo
Resp. Transitória
analı́tico para sistemas;
Sist.
Realimentados ▶ Um modelo aproximado pode ser obtido por testes na
Lugar das Raı́zes
resposta em frequência;
Resp. em
Frequência ▶ Procedimento:
Aula 16
Aula 17 1º Excitar a entrada do sistema com um sinal senoidal:
Aula 18
Aula 19
amplitude, frequência e fase conhecidas;
Aula 20 2º Coletar o sinal (senoidal) de saı́da: amplitude e fase;
Controle PID 3º Repetir para diversas frequências;
Controle Moderno 4º Traçar os diagramas de Bode e obter uma FT
aproximada;
Sistemas de Fase Mı́nima
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Importante!
Introdução
Ao observar o diagrama de fase em ω → ∞:
Modelagem
Resp. Transitória
arg(ĝ(j∞)) = −90◦ (n − m) ⇔ fase mı́nima (Por quê?)
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
▶ O caso mais comum ocorre quando há uma variação
Introdução
constante da fase com a frequência;
Modelagem
Resp. Transitória
▶ Assim, pode-se assumir que há um atraso de transporte
Sist. no sistema, da forma:
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aula 19
Tópico 6 - Resposta em Frequência
Aula 20 Aula 16: Diagramas de Bode
Aula 17: O Critério de Nyquist
Aula 18: Parâmetros de Desempenho
Controle PID Aula 19: Resposta em Malha Fechada
Aula 20: Projeto de Compensadores
Controle Moderno
Tópico 7 - Controle PID
Aula 21: Controle PID
Prof. Victor L.
Yoshimura
Introdução
▶ Se a compensação introduzida não for tipo P, o DN
Modelagem
Resp. Transitória
altera seu formato;
Sist. ▶ Assim, é mais conveniente o uso de diagramas de Bode
Realimentados
para projetar compensadores;
Lugar das Raı́zes
Resp. em
▶ Caracterı́sticas desejadas na resposta em malha aberta:
Frequência ▶ alto ganho na região de baixa frequência;
Aula 16
Aula 17 ▶ cruzamento de ganho com inclinação de −20dB por
Aula 18
Aula 19
década;
Aula 20 ▶ atenuação significativa na região de alta frequência.
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Este compensador foi apresentado em (5.11) e será reescrito
Modelagem
como:
Resp. Transitória Ts + 1
Sist.
ĉ(s) = αkc , 0<α<1 (6.45)
Realimentados
αT s + 1
Lugar das Raı́zes Observe que ĉ(j0) = αkc e que ĉ(j∞) = kc . Para o restante
Resp. em de seu DN:
Frequência
Aula 16
Aula 17 kc (1 + α) jωT + 1 1+α
Aula 18 ĉ(jω) − = kc α −
Aula 19 2 jαωT + 1 2
Aula 20
kc (1 − α) jαωT − 1
Controle PID =
Controle Moderno
2 jαωT + 1
1−α
= kc (6.46)
2
Compensação em Avanço de Fase II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
φm [◦ ]
Introdução 90
Modelagem 80
Resp. Transitória 70
Im(ĉ) 60
Sist.
Realimentados 50
Lugar das Raı́zes 40
Resp. em 30
Frequência φm αkc kc (1+α) kc
Aula 16
2 Re(ĉ) 20
Aula 17
10
Aula 18
Aula 19
Aula 20 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 α
Controle PID Avanço máximo de fase:
Controle Moderno
1−α
φm = arcsen (6.47)
1+α
Frequência onde φm ocorre:
1
ωm = √ (6.48)
T α
Compensação em Avanço de Fase - Projeto
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
O procedimento para o projeto segue os passos:
i Determine o ganho (αkc ) necessário para a exigência de
Introdução
Modelagem
erro estático;
Resp. Transitória
ii Plote os diagramas de Bode e avalie a margem de fase;
Sist. iii Determine α a partir de (6.47);
Realimentados
iv Escolha a nova frequência de cruzamento de ganho
Lugar das Raı́zes
com:
Resp. em
Frequência ||ĝ(j ω̃g )|| = 10 log α (6.49)
Aula 16
Aula 17
Aula 18
v Use (6.48) para determinar T ;
Aula 19
Aula 20
vi Determine kc com os resultados anteriores.
Controle PID
Exemplo
Controle Moderno
Projete um compensador avanço para Kv ≥ 20, M F ≥ 50◦
e M G ≥ 10dB, para a planta:
4
ĝ(s) =
s(s + 2)
Compensação em Atraso de Fase I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Este compensador foi apresentado em (5.12) e será reescrito
Modelagem como:
Ts + 1
Resp. Transitória ĉ(s) = βkc , β>1 (6.50)
Sist.
βT s + 1
Realimentados
Controle PID
Prof. Victor L.
Yoshimura
▶ O compensador atraso afeta a MG. Considere
▶ ωf ≫ T −1 > (βT )−1 , e
Introdução
▶ Ke e K̃e as constantes de erro atual e desejada,
Modelagem
respectivamente.
Resp. Transitória
▶ Aplicando TVF a ĉĝ, tem-se
Sist.
Realimentados
Lugar das Raı́zes lim sĉĝ = lim ĉ lim sĝ = βkc Ke = K̃e
s→0+ s→0+ s→0+
Resp. em
Frequência
Aula 16
▶ Por outro lado, na frequência ωf tem-se
Aula 17
Aula 18
Aula 19 jωf T + 1 K̃e
Aula 20 ĉ(jωf ) = βkc ≈
Controle PID
jβωf T + 1 βKe
Controle Moderno ▶ Assim, o fator a ser adicionado à MG atual é o ganho
nessa frequência, ou seja
!
K̃e
∥ĉ(jωf )∥ ≈ 20 log (6.51)
βKe
Compensação em Atraso de Fase - Projeto
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
O procedimento para o projeto segue os passos:
i Determine o ganho (βkc ) necessário para a exigência de erro
Introdução
estático;
Modelagem
Resp. Transitória
ii Plote os diagramas de Bode e determine a frequência ωgdes
Sist.
que teria a M F desejada (com uma folga);
Realimentados
iii Escolha o zero do compensador em valores bem abaixo de
Lugar das Raı́zes
ωgdes (pelo menos, uma década);
Resp. em
Frequência iv Determine β para impor ganho de 0dB em ωgdes com:
Aula 16
Aula 17
Aula 18
||ĝ(jωgdes )|| = 20 log β (6.52)
Aula 19
Aula 20
v Determine o polo e kc .
Controle PID
Prof. Victor L.
Yoshimura
Sua formulação pode ser escrita como a combinação de
Introdução
(6.45) e (6.50):
Modelagem
Resp. Transitória T1 s + 1 T2 s + 1
Sist.
ĉ(s) = αβkc , 0<α<1eβ>1
Realimentados
αT1 s + 1 βT2 s + 1
(6.53)
Lugar das Raı́zes
Resp. em
Para αβ = 1:
Frequência
Aula 16
Im(ĉ)
Aula 17
Aula 18 com:
Aula 19
Aula 20
ωl 1
Controle PID
ω→∞
ωl = √ (6.54)
ω=0
Re(ĉ) T1 T2
Controle Moderno
Observação
Devido às exigências de projeto, deve-se fazer T2 ≫ T1 .
Comentários sobre os Compensadores
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura i O compensador avanço realiza a compensação pelo seu
Introdução avanço de fase, aumentando a M F .
Modelagem ii O compensador atraso o faz pela atenuação em altas
Resp. Transitória
frequências, aumentando a M G.
Sist.
Realimentados iii O compensador avanço tende a aumentar a largura de
Lugar das Raı́zes banda, reduzindo o tempo de acomodação, porém,
Resp. em
Frequência
pode introduzir ruı́do no canal de controle;
Aula 16
Aula 17
iv O compensador avanço requer maior ganho para corrigir
Aula 18
Aula 19
o erro em regime.
Aula 20
Prof. Victor L.
Yoshimura
Introdução
▶ A estrutura do PID foi apresentada em (4.10)-(4.11):
Modelagem
Resp. Transitória kp
Sist.
ĉpid (s) = kp + + kp Td s (7.1)
Realimentados Ti s
Lugar das Raı́zes
▶ Amplamente utilizada na indústria, devido à natureza
Resp. em
Frequência incerta das plantas;
Controle PID ▶ Além disto, é comum haver apenas dados experimentais
Aula 21
Prof. Victor L.
Yoshimura
▶ Supõe planta sem integradores, nem polos conjugados
Introdução dominantes;
Modelagem
▶ Forma assumida para a planta:
Resp. Transitória
Sist. K
Realimentados ĝ(s) = e−Ls (7.2)
Ts + 1
Lugar das Raı́zes
Resp. em
Frequência
▶ Aplique um degrau à entrada da planta e verifique sua saı́da.
Controle PID Sintonia proposta:
Aula 21
Controle Moderno y
T
K
se P
L
T T
ĉ(s) = 0,9 + 0,27 2 se PI
L L s
−1 2
0,6T (s + L )
se PID
s
L T +L t (7.3)
2º Método de Ziegler-Nichols
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ Monte um sistema em realimentação unitária com
Modelagem
compensador P;
Resp. Transitória
Resp. em
▶ De posse de kcr e Tcr (perı́odo da oscilação
Frequência sustentada), a sintonia proposta é
Controle PID
se P
Aula 21
0,5kcr
Controle Moderno
kcr
ĉ(s) = 0,45kcr + 0,54 se PI (7.4)
Tcr s
−1 2
0,075k T (s + 4Tcr )
se PID
cr cr
s
Comentários aos Métodos de ZN
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
1º As regras de ZNs são muito usadas em plantas não
Introdução
precisamente conhecidas (controle robusto de sistemas
Modelagem
incertos);
Resp. Transitória
Sist.
2º essas regras podem ser utilizadas em plantas
Realimentados precisamente conhecidas;
Lugar das Raı́zes
3º Se a planta apresentar integrador, as regras podem não
Resp. em
Frequência ser aplicáveis;
Controle PID
Aula 21
4º Ambos os métodos almejam overshoot inferior a 25%
Controle Moderno para entrada degrau.
Exemplo
Projete compensadores PID para as plantas:
(s + 2)(s + 3) 1
ĝ(s) = e ĝ(s) =
s(s + 1)(s + 5) s(s + 1)(s + 5)
Controle PID Convencional
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução ˆ
d(s)
Modelagem 1
Resp. Transitória
r̂(s) û(s) ŷ(s)
Sist. Σ 1/Ti s Σ kp Σ ĝ(s)
Realimentados
−
Lugar das Raı́zes
Td s
Resp. em n̂(s)
Frequência
Σ
Controle PID
Aula 21
Td s
ĉd (s) = , γ ≈ 0,1 (7.5)
1 + γTd s
Controle PI-D
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Troca-se o impulso por um pulso abrupto (“set point
Resp. Transitória kick”);
Sist. ▶ A fim de evitá-lo, usa-se o PI-D:
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Existem aplicações onde não se deseja um sinal degrau
Resp. Transitória na atuação;
Sist. ▶ Para evitá-lo, usa-se o I-PD:
Realimentados
Prof. Victor L.
Yoshimura
As FTMFs para os diferentes PIDs são:
Introdução
Modelagem ŷ(s) ˜
ĉ(s)ĝ(s)
Resp. Transitória = (7.6)
r̂(s) 1 + ĉpid (s)ĝ(s)
Sist.
Realimentados
Observação
ˆ não se altera nas estruturas apresentadas. (Sim,
ŷ(s)/d(s)
provar este fato e (7.6)-(7.7) é um exercı́cio!)
Controle com Dois Graus de Liberdade I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Considere o sistema de controle
Introdução
Modelagem ˆ
d(s)
Resp. Transitória
Sist.
Realimentados
r̂(s) ŷ(s)
Σ ĉ(s) Σ ĝ(s)
Lugar das Raı́zes
−
Resp. em n̂(s)
Frequência Σ
Controle PID
Aula 21
Controle Moderno
ĝ(s) h i r̂(s)
(7.8)
ˆ
ŷ(s) = ĉ(s) 1 −ĉ(s) d(s)
1 + ĉ(s)ĝ(s)
n̂(s)
Observação
Como a planta é dada, ao projetar ĉ(s) para o canal
referência-saı́da, os demais ficam fixos.
Controle com Dois Graus de Liberdade II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. Transitória ˆ
d(s)
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Observação
Estes parâmetros dependem de uma função a ser
estabelecida: a sensibilidade.
Sensibilidade
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. em Importante!
Frequência
Controle PID
A sensibilidade é o quociente entre a variação relativa da
Aula 21
saı́da para uma variação relativa da FTMF. Ou seja:
Controle Moderno
∆ŷ(s)
ŷ(s) d ŷ(s) ĝ(s) d ln ŷ(s)
S(s) = lim = = (7.10)
∆ĝ(s)→0 ∆ĝ(s) d ĝ(s) ŷ(s) d ln ŷ(s)
ĝ(s)
Sensibilidade: Sistemas em Malha Fechada
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Considere o sistema realimentado com um grau de
Resp. Transitória
Sist.
liberdade;
Realimentados
▶ A planta, ĝ(s), é incerta, mas não o compensador ĉ(s).
Lugar das Raı́zes
Resp. em
▶ Assim, tem-se com a definição (7.10):
Frequência
d ŷ ĝ (3.4) 1 + ĉĝ d ĉĝ 1
Controle PID S(s) = = = (7.11)
Aula 21 d ĝ ŷ ĝ7→ĉĝ ĉ d ĝ 1 + ĉĝ 1 + ĉĝ
Controle Moderno
Observação
Note que a sensibilidade coincide com a FT da referência
para o erro para sistemas sem perturbação e ruı́do.
Sensibilidade: Sistemas sob Perturbação e Ruı́do
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Neste caso, lembrando que ε = r − y e, aplicando (7.11) a (7.8),
Introdução tem-se:
Modelagem 1 ĉ(s)ĝ(s)
ε̂(s) = ˆ
(r̂(s) − ĝ(s)d(s)) + n̂(s) (7.12)
Resp. Transitória
1 + ĉ(s)ĝ(s) 1 + ĉ(s)ĝ(s)
Sist. | {z } | {z }
Realimentados S(s) T (s)
Lugar das Raı́zes
Resp. em
Frequência
Importante!
Controle PID ▶ T (s) é chamada sensibilidade complementar e corresponde à
Aula 21
Controle Moderno
FT do ruı́do para o erro.
▶ Note que S + T = 1. Logo, são parâmetros conflitantes!
▶ Obedecer a todos os critérios estabelecidos para o controle
robusto é impossı́vel! Há “trade off”.
Exercı́cio
Estude para o caso com dois graus de liberdade.
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
▶ Foi mostrado em (2.14) que os autovalores de A são os
Introdução polos da FT;
Modelagem
▶ A representação no espaço de estado não é única;
Resp. Transitória
Sist.
▶ Tome, então uma transformação linear bijetora
Realimentados
(isomorfismo), T ∈ Rn×n , e um vetor x = Tz.
Lugar das Raı́zes
Aplicando a (2.12):
Resp. em
Frequência
Observação
Aula 25
Aula 26
Aula 27
Aula 28 Ou seja, com um isomorfismo adequadamente escolhido,
pode-se trocar as matrizes de estado A 7→ T−1 AT,
B 7→ T−1 B, C 7→ CT e D 7→ D e o vetor de estado x 7→ z.
O que acontece com os autovalores com a aplicação de T?
Autovalores, Polos e Representação II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Considere a EDO (2.12a), com u = 0:
Resp. Transitória
ẋ = Ax, x(0) = xo ∈ Rn
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Substituindo (8.2) na EDO:
Introdução
Modelagem
∞
X ∞
X
Resp. Transitória
ifi ti−1 = A fi ti
Sist.
i=1 i=0
Realimentados
Prof. Victor L.
Yoshimura Substituindo em (8.2):
Introdução ∞
Ai ti
(8.3)
X
Modelagem x(t) = xo
Resp. Transitória i=0
i!
| {z }
Sist.
Realimentados eAt
Lugar das Raı́zes
Resp. em Importante!
Frequência
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Lema (Propriedades da Matriz Exponencial)
Resp. Transitória Sejam A, B ∈ Rn×n , e a, b ∈ R. Valem, então, as
Sist. propriedades:
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aplicando a TL a (2.12a), com u = 0:
Resp. Transitória
L{ẋ} = L{Ax} ⇒ sx̂ − xo = Ax̂ ⇒ x̂ = (sI − A)−1 xo
Sist.
Realimentados
Controle Moderno
Aula 22 Observação
Aula 23
Aula 24
Aula 25
Comparando (8.3) e (8.6):
Aula 26
Aula 27
Aula 28
L−1 {(sI − A)−1 } = eAt (8.7)
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
" #
0 1
Introdução
Solucionemos o sistema ẋ = x, x(0) = xo .
Modelagem −2 −3
Resp. Transitória Observe que: " #
Sist. s −1
Realimentados sI − A =
2 s+3
Lugar das Raı́zes
Resp. em
Frequência
Para inverter esta matriz, usamos:
Controle PID " #
Controle Moderno −1 adj(sI − A) 1 s+3 1
(sI − A) = = 2
Aula 22
Aula 23
det(sI − A) s + 3s + 2 −2 s
Aula 24
s+3 1
Aula 25
Aula 26 (s + 1)(s + 2) (s + 1)(s + 2)
Aula 27
Aula 28
= −2 s
(s + 1)(s + 2) (s + 1)(s + 2)
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Considere a EDO (2.12a):
Introdução
Modelagem
ẋ = Ax + Bu, x(0) = xo ∈ Rn
Resp. Transitória
Sist.
Realimentados Rearranjando a EDO e pré-multiplicando por e−At :
Lugar das Raı́zes
Resp. em
e−At [ẋ − Ax] = e−At Bu
Frequência
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Aplicando a TL a (2.12a):
Resp. Transitória
Sist.
L{ẋ} = L{Ax + Bu} ⇒ sx̂ − xo = Ax̂ + Bû
Realimentados
⇒ x̂ = (sI − A)−1 (xo + Bû)
Lugar das Raı́zes
Resp. em
Frequência Aplicando a TIL ao resultado anterior, escreve-se:
Controle PID
Controle Moderno x = L−1 {(sI − A)−1 }xo + L−1 {(sI − A)−1 Bû} (8.9)
Aula 22
Aula 23
Observação
Aula 24
Aula 25
Aula 26
Aula 27 (8.8) é obtida de (8.9) com a aplicação do Teorema da
Convolução.
Aula 28
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Solucionemos o sistema
Introdução 0 1 0
Modelagem
ẋ = x+ u, x(0) = xo , u = 1̃
−2 −3 1
Resp. Transitória
Resp. em
x = eAt xo + eA(t−τ ) Bu d τ
Frequência 0
2e−t − e−2t e−t − e−2t
Controle PID
= x +
Controle Moderno
−2e−t + 2e−2t −e−t + 2e−2t o
Aula 22
Z t −(t−τ )
− e−2(t−τ )
Aula 23 e
+ dτ
Aula 24
Aula 25 0 −e−(t−τ ) + 2e−2(t−τ )
Aula 26
Aula 27
Aula 28
Portanto:
"1 1 #
2e−t − e−2t e−t − e−2t
− e−t + e−2t
x= x + 2 2
−2e−t + 2e−2t −e−t + 2e−2t o e−t − e−2t
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
▶ Até agora, vimos esquemas de realimentação pela
Sist. leitura da saı́da;
Realimentados
Controle Moderno
▶ Objetivo da realimentação de estado: transferir o
Aula 22
Aula 23
estado de um sistema xo 7→ xf no tempo to 7→ tf ;
Aula 24
Aula 25
▶ Pergunta imediata: “Sob quais condições isto é
Aula 26
Aula 27
possı́vel?”
Aula 28
Definição de Controlabilidade
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
À solução de (2.12a), (8.8), com a definição acima, to = 0 e
xf = 0:
Introdução
Z tf
Modelagem
x(tf ) = xf = 0 = e Atf
xo + eA(tf −τ ) Bu d τ
Resp. Transitória
0
Sist.
Realimentados Assim: Z tf
Lugar das Raı́zes
xo = − e−Aτ Bu d τ
Resp. em 0
Frequência
Controle PID
Usando a fórmula de interpolação de Sylvester:
Controle Moderno n−1 Z tf
Aula 22
X
Aula 23 xo = Ak B −αk (τ )u d τ
Aula 24
Aula 25
k=0 |0 {z }
Aula 26 βk
Aula 27
Aula 28 β0
β1
An−1 B (8.10)
= B AB ... ..
.
βn−1
Matriz de Controlabilidade II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Considere o sistema massa-mola-amortecedor (aula 04):
Introdução
Modelagem
0 1 0
Resp. Transitória
ẋ = b k x + 1 u
Sist.
Realimentados − −
h mi m m
Lugar das Raı́zes
Resp. em
y= 1 0 x
Frequência
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Introdução
Calculemos a matriz de controlabilidade de (8.13):
Modelagem
Resp. Transitória
0 0
Sist. .. .. Im
. .
Realimentados −αn−1 Im
B = 0 AB =
0
. . . An−1 B =
(α2n−1 − αn−2 )Im
Lugar das Raı́zes
0
I0 ..
Resp. em
.
0
Frequência m
Im −αn−1 Im
Controle PID
Controle Moderno
Aula 22
que possui sempre n colunas LI!
Aula 23
Aula 24 Lema
Aula 25
Aula 26
Aula 27
i A FCC é controlável.
Aula 28
ii Todo SLIT-C similar à FCC é controlável.
Alguns Comentários
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ No item ii do lema anterior, vale a recı́proca!
Resp. Transitória
Sist.
▶ A controlabilidade no plano s ocorre se não ocorrer
Realimentados
cancelamentos na FT;
Lugar das Raı́zes
Resp. em
▶ Pode ser necessário determinar a controlabilidade de
Frequência saı́da;
Controle PID
▶ Neste caso, é possı́vel mostrar que a matriz:
Controle Moderno
Aula 22 h i
Aula 23
Aula 24
Cy = CB CAB . . . CAn−1 B D (8.14)
Aula 25
Aula 26
Aula 27 deve ser tal que rank(Cy ) = p, para tal controlabilidade.
Aula 28
Diagrama de Blocos em Malha Fechada
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
D
Resp. Transitória
Sist.
Realimentados r u ẋ R x y
Lugar das Raı́zes
M Σ K B Σ C Σ
Resp. em
Frequência
A
Controle PID
Controle Moderno
Aula 22
Aula 23 ▶ O sinal de controle é:
Aula 24
Aula 25
Aula 26
Aula 27
u = K(Mr + x) (8.15)
Aula 28
Prof. Victor L.
Yoshimura
Introdução
Modelagem
D
Resp. Transitória
Sist.
Realimentados r u ẋ R x y
Lugar das Raı́zes
M Σ K B Σ C Σ
Resp. em
Frequência
A
Controle PID
Controle Moderno
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Tome um SLIT-C-SISO;
Resp. Transitória
Prof. Victor L.
Yoshimura Dada a EDO (2.12), se pode escrever sua FCC com:
Introdução
0 1 0 ... 0 0 0
Modelagem
Resp. Transitória
0 0 1 ... 0 0 0
.. .. .. .. .. .. .
. .
Sist.
Realimentados
ẋ = . . . . . x + . u
0 0 0 ... 1 0 0
Lugar das Raı́zes
Resp. em
0 0 0 ... 0 1 0
αn−2
Frequência
− ααno − ααn1 − ααn2 . . . − αn − ααn−1
n
1
Controle PID
(8.18a)
Controle Moderno h i h i
Aula 22
Aula 23
y = βo β1 β2 . . . βn−2 βn−1 x + βn u (8.18b)
Aula 24
Aula 25
Aula 26
Aula 27 Observação
Aula 28
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Sist.
▶ Seja um SLIT-C-SISO com realimentação de estado,
Realimentados dado por (2.12), (8.15) e r = 0.
Lugar das Raı́zes
▶ Considere, ainda, uma condição inicial, x(0) = xo ,
Resp. em
Frequência devida a uma perturbação.
Controle PID
▶ Sua solução, a partir de (8.3) torna-se:
Controle Moderno
Aula 22
Aula 23 x = e(A+BK)t xo (8.19)
Aula 24
Aula 25
Aula 26
Aula 27
Aula 28
Alocação de Polos: Condição Suficiente I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Lema
Introdução Sejam (A, B) ∈ Rn×n × Rn×1 um par controlável,
Modelagem (Ã, B̃) ∈ Rn×n × Rn×1 o par canônico controlável e C e C˜ suas
Resp. Transitória respectivas matrizes de controlabilidade. Então
Sist. −1
Realimentados T = C C˜ (8.20)
Lugar das Raı́zes
Prof. Victor L.
Yoshimura
Sejam λ1 , λ2 , . . . , λn ∈ C os autovalores escolhidos (possivelmente
em pares conjugados). Então, pode-se escrever o polinômio
Introdução caracterı́stico desejado:
Modelagem n n−1
Y X
Resp. Transitória P (s) = (s − λi ) = sn + αi+1 si , α1 , α2 , . . . , αn ∈ R
Sist. i=1 i=0
Realimentados
(8.21)
Lugar das Raı́zes
Seja T a similaridade que leva (Ã, B̃) em (A, B). Logo, tem-se:
Resp. em
(8.22)
Frequência u = Kx = KTx̃ = k̃1 k̃2 . . . k̃n x̃
Controle PID
Controle Moderno
Aula 22 Assim, referindo ao sistema em FCC, tem-se, em (2.12):
Aula 23
Aula 24 x̃˙ = Ãx̃ + B̃u = (T−1 AT + T−1 BKT)x̃
Aula 25
Aula 26
Aula 27
Aula 28 Deseja-se fazer (lembre da invariância de autovalores)
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
▶ As condições (8.23) sempre podem ser atendidas, se o
Sist. sistema for controlável;
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
1º Verifique se o sistema dado é controlável;
Resp. Transitória 2º Escreva o polinômio caracterı́stico de A:
Sist.
Realimentados n−1
(8.24)
X
Lugar das Raı́zes det(sI − A) = sn + ai+1 si
Resp. em i=0
Frequência
Prof. Victor L.
Yoshimura
Posicione os autovalores em −2 ± j4 e −10 para o SLIT-C
com matrizes:
Introdução
Modelagem 0 1 −1 0
Resp. Transitória A = −1 −5 0 B = 1
Sist.
Realimentados
−1 −5 −1 1
Lugar das Raı́zes
1º passo: Matriz de controlabilidade
Resp. em
Frequência
Controle PID 0 0 1
det C = det 1 −5 25 = −1 ⇒ Sistema controlável
Controle Moderno
Aula 22
Aula 23
1 −6 31
Aula 24
Aula 25
Aula 26
Aula 27
2º passo: Polinômio caracterı́stico
Aula 28
s −1 1
det(sI − A) = det 1 s + 5 0 = s3 + 6s2 + 5s + 1
1 5 s+1
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura 3º passo: Transformação de similaridade
Introdução
Com o polinômio caracterı́stico do passo anterior e
Modelagem
observando (8.18a):
Resp. Transitória
0 0 1 5 6 1 1 0 0
Sist. −1 (8.20)
Realimentados
C˜ = 0 1 −6 ⇒ C˜ = 6 1 0 ⇒ T = 0 1 1
Lugar das Raı́zes
1 −6 31 1 0 0 0 0 1
Resp. em
Frequência 4º passo: Polinômio desejado
Controle PID
h i h i
KT = −199 −55 −8 ⇒ K = −199 −55 47
Fórmula de Ackermann I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Resp. em
Frequência
Observe as identidades:
Controle PID
Controle Moderno
A 0 = I A 1 = A + BK A 2 = A2 + ABK + BKA ...
Aula 22
Aula 23
Aula 24
Ou seja:
Aula 25 k−1
A k = Ak + (8.26)
X
Aula 26
Aula 27
Ai−1 BKA k−i
Aula 28
i=1
Fórmula de Ackermann II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Tome os autovalores desejados, λ1 , λ2 , . . . , λn ∈ C, gerando
Introdução (8.21). Então:
Modelagem
Resp. Transitória
n n−1
αj+1 sj = det(sI − A )
Y X
Sist.
P (s) = (s − λj ) = sn +
Realimentados j=1 j=0
Lugar das Raı́zes
Prof. Victor L.
Yoshimura
Introdução
Em notação matricial:
Modelagem
Resp. Transitória
KA n−1 + αn KA n−2 + . . . + α2 K
KA n−2 + αn KA n−3 + . . . + α3 K
Sist.
Realimentados P (A) = −C
..
.
Lugar das Raı́zes
Resp. em
K
Frequência
Controle PID
Controle Moderno
Lema (Fórmula de Ackermann)
Aula 22
Aula 23
Nas condições do Teorema da Alocação de Polos, a
Aula 24
Aula 25
expressão
Aula 26
h i
K = 0 0 . . . 0 −1 C −1 P (A) (8.27)
Aula 27
Aula 28
Prof. Victor L.
Yoshimura
Para que a saı́da do sistema siga uma referência degrau, a
Introdução
FT (8.17) deve apresentar valor unitário para s = 0. Assim:
Modelagem
Resp. Transitória
[D − (C + DK)(A + BK)−1 B]KM = 1 (8.28)
Sist.
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
▶ Estratégia da realimentação de estado: medir todas as
Sist. variáveis de estado e usar (8.15);
Realimentados
Controle Moderno
▶ Realimentação dinâmica da saı́da: Usar o estado
Aula 22
Aula 23
estimado no lugar do real em (8.15);
Aula 24
Aula 25
▶ Pergunta imediata: “Sob quais condições isto é
Aula 26
Aula 27
possı́vel?”
Aula 28
Definição de Observabilidade
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Definição (Observabilidade)
Resp. Transitória
Sist.
O SLIT-C homogêneo (2.12) é dito observável se, dada a
Realimentados
saı́da do sistema em um intervalo [to , tf ], então é possı́vel
Lugar das Raı́zes
determinar x(to ).
Resp. em
Frequência
▶ A observabilidade responde à questão: “É possı́vel
Controle PID
determinar o estado a partir da leitura da saı́da?”
Controle Moderno
Aula 22 ▶ Se a observabilidade for verificada para um dado
Aula 23
Aula 24 sistema, diz-se que o par (A, C) é observável;
Aula 25
Aula 26 ▶ Como verificar a observabilidade de um dado par
Aula 27
Aula 28 (A, C)?
Matriz de Observabilidade I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Resp. Transitória
y = Cx = CeAt xo
Sist.
Realimentados Aplicando novamente a fórmula de interpolação de Sylvester:
Lugar das Raı́zes
Resp. em
n−1
X n−1
X
k
Frequência y=C αk (t)A xo = αk (t)CAk xo
Controle PID k=0 k=0
Controle Moderno n−1
X
Aula 22
Aula 23
= CAk αk (t)xo
Aula 24
k=0
Aula 25
Aula 26
Aula 27 Sem perda de generalidade, com xo = 0, deve-se resolver o
Aula 28
sistema:
CAk xo = 0, k = 0, 1, . . . , n − 1
Matriz de Observabilidade II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Assim, a transformação linear
Introdução
Modelagem
C
Resp. Transitória
CA
(8.30)
Sist.
Realimentados
O= ..
.
Lugar das Raı́zes
CAn−1
Resp. em
Frequência
(8.31)
Aula 26
Aula 27 rank(O) = n
Aula 28
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Considere o sistema massa-mola-amortecedor (aula 04):
Resp. Transitória
0 1 0
Sist.
Realimentados ẋ = b x + 1 u
k
Lugar das Raı́zes − −
h mi m m
Resp. em
Frequência y= 1 0 x
Controle PID
Controle Moderno
Aula 22
Determinemos a matriz de observabilidade:
Aula 23 " #
Aula 24
C
Aula 25
O= = I2
Aula 26
Aula 27
CA
Aula 28
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Controle Moderno
Aula 22
▶ Desta forma, obtém-se a Forma Canônica
Aula 23
Aula 24
Observável (FCO);
Aula 25
Aula 26
▶ O vetor de estado não é o mesmo!
Aula 27
Aula 28
Diagrama de Blocos
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
D
Modelagem
Resp. Transitória r u ẋ R x y
M Σ K B Σ C Σ
Sist.
Realimentados
Controle PID D
Controle Moderno
Aula 22
Aula 23 A
Aula 24
Aula 25
Aula 26 x̂˙ R
Aula 27 B Σ C Σ
Aula 28 x̂ −
ŷ − y
L
Observador de Estado
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura ▶ Nem todas as variáveis de estado possuem sensor;
Introdução
▶ Assim, x não está completo e deve-se estimá-lo a partir
Modelagem de y;
Resp. Transitória ▶ Estratégia: construir um “clone” do sistema original, da
Sist.
Realimentados
forma
Lugar das Raı́zes
Resp. em
x̂˙ = Ax̂ + Bu + L(ŷ − y) (8.34a)
Frequência
ŷ = Cx̂ + Du (8.34b)
Controle PID
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ Sob quais condições tem-se x̂ = x?
Resp. Transitória
Sist.
▶ Sob quais condições o sistema mostrado é estável?
Realimentados
▶ A inserção do observador altera a dinâmica do sistema?
Lugar das Raı́zes
De que forma?
Resp. em
Frequência ▶ O projeto do observador é afetado pela escolha de K e
Controle PID
de M?
Controle Moderno
Aula 22 ▶ Existe alguma relação entre a dinâmica do observador e
Aula 23
Aula 24 da planta?
Aula 25
Aula 26 ▶ O problema a ser resolvido é chamado estimação de
Aula 27
Aula 28 estado.
Erro de Observação
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ Devemos, primeiro determinar a “eficiência” do
Modelagem
Resp. Transitória
observador (2.12);
Sist. ▶ Será que ele consegue fazer x̂ = x?
Realimentados
Controle PID
Controle Moderno
ε̇ = ẋ − x̂˙ = Ax + Bu − (Ax̂ + Bu + L(ŷ − y))
Aula 22
Aula 23
= Aε + L(y − ŷ)
Aula 24
(2.12b)
Aula 25
(8.34b)
Aula 26
Aula 27
= Aε + L(Cx + Du − (Cx̂ + Du))
Aula 28
= (A + LC)ε (8.36)
Dualidade e Projeto de Observadores
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Lema
Introdução Seja ε o erro de estimação do SLIT-C-SISO (2.12)-(8.34a).
Modelagem Então ε → 0, para t → ∞ se, e somente se, ∃L ∈ Rn×1 , tal
Resp. Transitória que A + LC tenha seus autovalores no semi-plano complexo
Sist. esquerdo.
Realimentados
ν = L′ χ (8.38)
Observadores de Estado: Projeto
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura 1º O projeto é feito a partir do sistema dual. Verifique se este é
controlável, ou seja, se o primal é observável;
Introdução
Observação
Aula 25
Aula 26
Aula 27
Aula 28
▶ Adapte os passos 3 a 5: use a matriz de controlabilidade dual
e o polinômio alocador para os autovalores do observador.
▶ É possı́vel aplicar a fórmula de Ackermann de forma
adequada.
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Posicione os autovalores do observador em −2 ± j0,2 para o
Introdução SLIT-C com matrizes:
Modelagem " #
Resp. Transitória −1 1 h i
A= C= 1 0
Sist. 0 5
Realimentados
Prof. Victor L.
Yoshimura
3º passo: Transformação de similaridade
Introdução Com o polinômio caracterı́stico do passo anterior e
Modelagem observando (8.18a):
Resp. Transitória
0 1 −1 −4 1 (8.20) −5 1
Sist. C˜ = ⇒ C˜ = ⇒ T=
Realimentados 5 4 1 0 1 0
Lugar das Raı́zes
Resp. em
4º passo: Polinômio desejado
Frequência
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Introdução Lema
Modelagem Sejam A, B, C ∈ (Rn×n , Rn×m , Rm×m ). Então:
" #!
Resp. Transitória A B
Sist.
i λ∈Λ ⇔ λ ∈ Λ(A) ou λ ∈ Λ(C), e;
Realimentados
0 C
" #
Lugar das Raı́zes
A B
Resp. em ii A multiplicidade de λ em é a soma de suas
Frequência 0 C
Controle PID multiplicidades em A e em C.
Controle Moderno
Teorema (Separação)
Aula 22
Aula 23
Aula 24
Aula 25 O projeto das matrizes de realimentação de estado e de
Aula 26
Aula 27 ganho da saı́da podem ser realizados independentemente e
alocando os polos de A + BK e de A + LC em posições
Aula 28
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. em
▶ Logo, o projeto da matriz de ganho da referência, M,
Frequência
pode ser feito com (8.28);
Controle PID
▶ Pode-se usar (8.29) para a norma euclidiana possı́vel;
Controle Moderno
Aula 22
Aula 23
▶ É interessante posicionar os autovalores do observador
Aula 24
Aula 25
de 2 a 5 vezes mais à esquerda do último autovalor da
Aula 26
Aula 27
planta (realimentada).
Aula 28
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
▶ Para SLIT-Cs, muitos critérios de estabilidade estão
Introdução
disponı́veis:
Modelagem
▶ Critério de Routh-Hurwitz;
Resp. Transitória
▶ Critério de Nyquist.
Sist.
Realimentados ▶ Se o sistema não for linear ou não for invariante, esses
Lugar das Raı́zes
critérios são inválidos.
Resp. em
Frequência ▶ A Teoria de Lyapunov (formulada em 1892!) é aplicável
Controle PID a tais situações.
Controle Moderno
Aula 22
▶ Para SLIT-Cs, a Teoria de Lyapunov proporciona uma
Aula 23
Aula 24
análise e sı́ntese de controladores mais versátil.
Aula 25
Aula 26
Aula 27
Importante!
A Teoria de Lyapunov não invalida o que estudamos até
Aula 28
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Considere o sistema:
" # " #
Introdução −1 2 −x1 + 2x2
Modelagem
ẋ = x=
0 −5 −5x2
Resp. Transitória
Sist.
Realimentados 10
Resp. em 6
Frequência
4
Controle PID
Controle Moderno 2
Aula 22
x2
0
Aula 23
Aula 24
−2
Aula 25
Aula 26 −4
Aula 27
Aula 28 −6
−8
−10
−12 −10 −8 −6 −4 −2 0 2 4 6 8 10 12
x1
Definições de Estabilidade I
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Observação
Introdução
Modelagem
Se xe ̸= 0, pode-se fazer ze = 0, onde z é obtida por uma
Resp. Transitória mudança adequada de variáveis de estado (translação) em x.
Sist.
Realimentados Definição (Estabilidade no Sentido de Lyapunov)
Lugar das Raı́zes
O ponto de equilı́brio xe é dito (localmente) estável no
Resp. em
Frequência sentido de Lyapunov (ou Lyapunov-estável) se, dado ε > 0,
Controle PID então ∃δ > 0, tal que:
Controle Moderno
Aula 22
Aula 23
||xo − xe || < δ ⇒ ||ϕ(t; xo , to ) − xe || < ε
Aula 24
Aula 25
Aula 26
Aula 27
Aula 28
Definição (Estabilidade Assintótica)
O ponto de equilı́brio xe é dito (localmente)
assintoticamente estável se, além de Lyapunov-estável,
verifica-se, para t → ∞, ϕ(t; xo , to ) → xe .
Definições de Estabilidade II
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Controle PID
ε xe δ → ∞, então a
Controle Moderno
estabilidade é dita global.
Aula 22
Aula 23
▶ Se dado ε > 0 não existir
Aula 24
Aula 25
δ (mesmo muito pequeno)
Aula 26 para as definições dadas, o
Aula 27 Lyapunov-estável
Aula 28 ponto de equilı́brio é dito
Assintoticamente estável
instável.
Instável
Positividade
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Considere P ∈ Rn×n , P = P′ . Uma importante classe de
funções Rn → R é dada pela forma quadrática:
Introdução
Sist.
Da Álgebra Linear, tem-se:
Realimentados ▶ V é positiva definida se, e somente se, todos os
Lugar das Raı́zes
autovalores de P forem positivos (P é uma matriz
Resp. em
Frequência positiva-definida);
Controle PID ▶ V é negativa definida se, e somente se, todos os
Controle Moderno autovalores de P forem negativos (P é uma matriz
Aula 22
Aula 23 negativa-definida);
Aula 24
Aula 25
Aula 26
Lema (Critério de Sylvester)
Aula 27
Aula 28 P ∈ Rn×n , P = P′ , é positiva-definida se, e somente se:
p11 ⋆
p11 > 0, > 0, . . . , |P| > 0
p21 p22
O Teorema da Estabilidade de Lyapunov
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Teorema (Lyapunov)
Introdução
Seja (8.41) um SC tendo a origem como ponto de equilı́brio.
Modelagem
Se houver uma função V : Rn → R com derivadas parciais
Resp. Transitória
de primeira ordem contı́nuas, tal que:
Sist.
Realimentados
i V é positiva-definida, e;
Lugar das Raı́zes
Resp. em
ii V̇ é negativa-definida.
Frequência
então o ponto de equilı́brio é assintoticamente estável.
Controle PID
Adicionalmente, se V → ∞ para ||x|| → ∞, a estabilidade é
Controle Moderno
Aula 22 global.
Aula 23
Aula 24
Aula 25
▶ Sistemas mecânicos (e elétricos) são estáveis se a
Aula 26
Aula 27
energia (que é positiva-definida) decai (derivada
Aula 28 negativa-definida);
▶ Em outros sistemas, pode ser difı́cil definir “energia”;
▶ Ideia de Lyapunov: definir uma função energia fictı́cia
no sistema.
Exemplo: Função “Energia”
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Para o exemplo anterior, seja V (x) = x′ Ix = ||x||2 = x21 + x22 .
Modelagem
Calculando sua derivada:
Resp. Transitória ′
−x1 + 2x2 2 2
V̇ = ∇V · ẋ = 2x1 2x2 = −2x1 + 4x1 x2 − 10x2
Sist. −5x2
Realimentados
= −2[(x1 − x2 )2 + 4x22 ]
Lugar das Raı́zes
Resp. em
Frequência
que é negativa-definida. Por outro lado, também tem-se:
Controle PID
′ 2 2 ′ 2 −2
Controle Moderno V̇ = ∇V · ẋ = −2x1 + 4x1 x2 − 10x2 = −x x
Aula 22
−2 10
Aula 23
Aula 24
Aula 25
cujos autovalores são −1,5 e −10,5. Ainda, pelo critério de
Aula 26 Sylvester:
Aula 27
2 −2
Aula 28
det[2] = 2 det = 16
−2 10
Logo, a origem deste sistema é global e assintoticamente estável.
Exemplo: Função “Energia”
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
10
Introdução
8
Modelagem
6
Resp. Transitória
Sist. 4
Realimentados
2
Lugar das Raı́zes
x2
0
Resp. em
Frequência −2
Controle PID −4
Controle Moderno −6
Aula 22
Aula 23 −8
Aula 24
Aula 25 −10
Aula 26 −12 −10 −8 −6 −4 −2 0 2 4 6 8 10 12
Aula 27 x1
Aula 28
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Teorema (Lyapunov Recı́proco)
Sist.
Realimentados Seja (8.41) um SC tendo a origem como ponto de equilı́brio.
Lugar das Raı́zes Se houver uma função W : Rn → R com derivadas parciais
Resp. em
Frequência
de primeira ordem contı́nuas, tal que:
Controle PID i W é positiva-definida em algum Ω ⊆ Rn ;
Controle Moderno
Aula 22
ii Ẇ é positiva-definida em Ω, e;
Aula 23
Aula 24
iii 0 ∈ Ω.
Aula 25
Aula 26 então o ponto de equilı́brio é instável.
Aula 27
Aula 28
Alguns Fatos sobre Estabilidade de SLIT-Cs
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução Lema
Modelagem Considere o SLIT-C homogêneo (2.12). Se A for invertı́vel,
Resp. Transitória então a origem é o único ponto de equilı́brio do sistema.
Sist.
Realimentados Lema
Lugar das Raı́zes Em um SLIT-C, as seguintes afirmativas são equivalentes:
Resp. em
Frequência i A origem é global e assintoticamente estável;
Controle PID
ii Os autovalores de A estão no semi-plano complexo
Controle Moderno
Aula 22
esquerdo;
iii Existem c > 0 e λ < 0, tais que:
Aula 23
Aula 24
Aula 25
Aula 26
Aula 27 ||ϕ(t; xo , 0)|| = ||eAt xo || ≤ c · eλt (8.43)
Aula 28
Prof. Victor L.
Yoshimura
Considere (8.42) como uma função candidata a Lyapunov.
Introdução Derivemo-la:
Modelagem
Resp. Transitória
V̇ (x) = ẋ′ Px + x′ Pẋ = (Ax)′ Px + x′ PAx
Sist.
Realimentados = x′ (A′ P + PA)x (8.44)
Lugar das Raı́zes
A′ P + PA = −Q (8.45)
Exemplo
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura Considere o SLIT-C homogêneo:
" #
Introdução
0 1
Modelagem ẋ = x
−1 −1
Resp. Transitória
Sist.
Realimentados Façamos Q = I e apliquemos (8.45):
Lugar das Raı́zes " #
Resp. em p ⋆
Frequência P= 1
p2 p3
Controle PID
O que origina:
Aula 26
Aula 27
Aula 28
−2p2 = −1
" #
1,5 0,5
P=
p1 − p2 − p3 = 0 0,5 1
2p2 − 2p3 = −1
Controle e
Servomecanismos Sumário
Prof. Victor L.
Tópico 1 - Introdução
Yoshimura Aula 01: Introdução aos Sistemas de Controle
Aula 02: Transformada de Laplace e Funções Transferência
Prof. Victor L.
Yoshimura
▶ Na Teoria de Lyapunov aplicada a para SLIT-Cs, vimos
Introdução
Modelagem
que todo sistema assintoticamente estável possui uma
Resp. Transitória
função de Lyapunov quadrática;
Sist. ▶ Como usar este fato para proporcionar estabilidade a
Realimentados
Controle Moderno
▶ Como, ainda, projetar este K para tornar algum ı́ndice
Aula 22
Aula 23
de desempenho o melhor possı́vel?
Aula 24
Aula 25
▶ Os resultados a serem apresentados aqui dependem de
Aula 26
Aula 27
estudos de um tópico da Matemática conhecido como
Aula 28
Cálculo de Variações;
▶ Assim sendo, as demonstrações serão omitidas e as
afirmações serão tomadas como verdadeiras.
Índices de Desempenho
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
▶ O ı́ndice de desempenho de um sistema é dado por uma
Resp. Transitória forma integral a ser minimizada:
Sist. Z ∞
Realimentados
Importante!
Aula 22
Aula 23
Aula 24
h i′
Aula 25 Fato: Se L for uma forma quadrática em x′ u′ , então é
Aula 26
Aula 27 possı́vel minimizar (8.46) com uma lei linear da forma (8.15)
Aula 28
(com r = 0).
O Regulador Linear Quadrático
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Importante!
Modelagem Considerando o ı́ndice geral (8.46), um problema de
Resp. Transitória otimização de grande importância é o do Regulador Linear
Sist. Quadrático (LQR), dado por:
Realimentados
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Controle Moderno
J= [(ε(t))2 + k(ε̇(t))2 ] d t, k>0
Aula 22
0
Aula 23
Aula 24 Observe que sua EDO é:
Aula 25
Aula 26
Aula 27 ε=r−y
Aula 28 ÿ + 2ζ ẏ + y = r =⇒ ε̈ + 2ζ ε̇ + ε = r̈ + 2ζ ṙ = 0, ∀t > 0
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Resp. em
▶ Note que, se ρ for grande, o critério irá reduzir o sinal
Frequência
de controle, u;
Controle PID
▶ Por outro lado, se Q for “muito positiva-definida”, o
Controle Moderno
Aula 22 critério irá fazer as trajetórias próximas à origem;
Aula 23
Aula 24 ▶ Estas caracterı́sticas são conflitantes!
Aula 25
Aula 26
Aula 27
Aula 28
Observação
O que é Q “muito positiva-definida”? Tome α ≫ 0 e faça
Q − αI ≻ 0.
A Equação Algébrica de Riccati
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Teorema
Modelagem
Seja (4.1) um SLIT-C-SISO, Q ∈ Rn×n , Q ⪰ 0 e ρ > 0. Se
Resp. Transitória
existir P ∈ Rn×n , P ⪰ 0, tal que (8.50) valha, então:
Sist.
Realimentados
i A realimentação de estado (8.51) minimiza (8.49) em
Lugar das Raı́zes
V (xo );
Resp. em
Frequência ii A origem é global e assintoticamente estável;
Controle PID
A′ P + PA − ρ−1 PBB′ P + Q = 0 (8.50)
Controle Moderno
−1 ′
Aula 22
Aula 23
u = −ρ B Px (8.51)
Aula 24
Aula 25
Observação
Aula 26
Aula 27
Aula 28
Prof. Victor L.
Yoshimura
Teorema
Introdução Sejam V ∈ R1×n , com Q = V′ V, para o critério quadrático
Modelagem (8.47) e a FT:
Resp. Transitória
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura
Considere uma planta sujeita a duas origens de ruı́do:
▶ w(t) ∈ Rn , um ruı́do aditivo presente nas variáveis de
Introdução
Modelagem
estado, e;
Resp. Transitória ▶ v(t) ∈ R, um ruı́do aditivo na medição da saı́da.
Sist.
Realimentados Assim, a representação da planta fica
Lugar das Raı́zes
Resp. em
ẋ = Ax + Bu + w (8.55a)
Frequência
y = Cx + Du + v. (8.55b)
Controle PID
Controle Moderno
Aula 22
Retornando à Aula 25, calculando o erro, obtém-se
Aula 23
Aula 24
Aula 25 ε̇ = (A + LC)ε + Lv + w (8.56)
Aula 26
Aula 27
Aula 28
Observação
Um observador com um ganho considerável minimiza a
dinâmica do erro, mas amplifica ruı́dos!
Ruı́do Branco Determinı́stico
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
▶ Um ruı́do é dito branco quando constante no eixo da
Modelagem frequência;
Resp. Transitória ▶ Escrevamos os ruı́dos presentes no controle como
Sist.
Realimentados
w = Wδ(t), W ∈ Rn (8.57a)
Lugar das Raı́zes
√
Resp. em v = µδ(t) (8.57b)
Frequência
{z W} + µe
| {z }
εw εv
Desempenho do Observador
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Teorema
Resp. Transitória
Seja (8.55) um SLIT-C-SISO, (8.56) a equação que descreve
Sist.
Realimentados o erro de observação, W ∈ Rn e µ > 0. Se existir
Lugar das Raı́zes P ∈ Rn×n , P ≻ 0, que satisfaça a
Resp. em
Frequência
AP + PA′ − µ−1 PCC′ P + WW′ = 0, (8.60)
Controle PID
Controle Moderno
Aula 22
então a matriz de ganho do observador
Aula 23
Aula 24
Aula 25 L = −µ−1 PC′ (8.61)
Aula 26
Aula 27
Aula 28
assegura que min J(L) = Tr(P).
Observador Ótimo via Lugar das Raı́zes Simétrico
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Resp. Transitória
Teorema
Sist.
Realimentados Considere a FT para o sistema (8.55)-(8.56)
Lugar das Raı́zes
Resp. em
φ(s) = C(sI − A)−1 W. (8.62)
Frequência
Prof. Victor L.
Yoshimura
Introdução
Prof. Victor L.
Yoshimura ▶ Em sistemas de 2ª ordem, pode-se plotar o gráfico
Introdução
(x1 , x2 ).
Modelagem
dω dω
Resp. Transitória ω = θ̇ =⇒ θ̈ = ω̇ = θ̇ = ω
dθ dθ
Sist.
Realimentados
▶ Substituindo nas EDs não-linear e linearizada, tem-se
Lugar das Raı́zes
Resp. em
R
g g
Frequência
ω d ω + sen(θ) d θ = 0 ⇒ ω 2 − 2 cos(θ) = c
Controle PID l l
Controle Moderno R
Aula 22 g g
Aula 23 ω d ω + θ d θ = 0 ⇒ ω2 + θ2 = d
Aula 24 l l
Aula 25
Aula 26 ▶ Ainda, com x1 = θ e x2 = ω, tem-se:
Aula 27
Aula 28
ẋ1 = x2 x2e = 0
g
ẋ2 = − sen(x1 ) x1e = kπ, k ∈ Z
l
Plano de Fase
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura 8
Introdução
6
Modelagem
Resp. Transitória
4
Sist.
Realimentados
Frequência
0
Controle PID
Controle Moderno −2
Aula 22
Aula 23
Aula 24
Aula 25
−4
Aula 26
Aula 27
Aula 28 −6
−8
−6 −4 −2 0 2 4 6
θ
Sistemas Tipo Lur’e
Controle e
Servomecanismos
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Um sistema tipo Lur’e (lê-se “Lur i”) tem a forma
Resp. Transitória
Sist.
ẋ = Ax + Bu (8.64a)
Realimentados
y = Cx + Du (8.64b)
Lugar das Raı́zes
Resp. em
u = −ϕ(y), (8.64c)
Frequência
Prof. Victor L.
Yoshimura
Prof. Victor L.
Yoshimura ED não linear:
Introdução
b g
Modelagem θ̈ + 2
θ̇ + sen(θ) = 0
l ml l
Resp. Transitória
Prof. Victor L.
Yoshimura
Resp. Transitória
b g
θ̈ + 2
θ̇ + sen(θ) = 0
Sist. l ml l
Realimentados
Prof. Victor L.
Yoshimura
Introdução
Modelagem
Teorema
Resp. Transitória
Considere o sistema tipo Persidiskii (8.66). Se existir
Sist.
P ∈ Rn×n , diagonal e positiva-definida, tal que A′ P + PA
Realimentados seja negativa-definida, então a origem do sistema é global e
Lugar das Raı́zes
assintoticamente estável.
Resp. em
Frequência A prova deste teorema se apoia na construção da função de
Controle PID Lyapunov
Controle Moderno
n Z xi
(8.67)
X
Aula 22 V (x) = 2 pi ψi (ξ) d ξ.
Aula 23
i=1 0
Aula 24
Aula 25
Exercı́cio
Aula 26
Aula 27