Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y...
Transcript of Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y...
![Page 1: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/1.jpg)
Aula 3
13 Ago 2019
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 1 / 39
![Page 2: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/2.jpg)
Resumo da aula passada - DMC
Informacao requerida sobre a planta:
Resposta a degrau g(n), n = 1, 2, . . . ,Ns
(assume-se g(0) = 0 e g(n) = g(Ns),∀n ≥ Ns).
Parametros de projeto:
Peso do controle ρ
Horizonte de predicao N
Horizonte de controle M
Inicializacao:
Fazer G =
g(1) 0 · · · 0g(2) g(1) · · · 0
......
. . ....
g(N) g(N − 1) · · · g(N −M + 1)
Calcular KMPC = [ 1 0 · · · 0 ](GTG + ρIM)−1GT
Fazer k = 0, u(−1)=0 e ∆u(−1)=∆u(−2)= · · ·=∆u(−Ns + 1)=0
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 2 / 39
![Page 3: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/3.jpg)
Rotina principal:
1 Ler y(k) (saıda da planta) e yref (valor de referencia)
2 Fazer r = [yref ]N3 Calcular
f (k + i |k) = y(k) +Ns−1∑n=1
[g(n + i)− g(n)]∆u(k − n), i = 1, 2, . . . ,N
4 Fazer f = [f (k + 1|k) f (k + 2|k) · · · f (k + N|k)]T
5 Calcular o incremento no controle: ∆u(k) = KMPC (r − f)
6 Atualizar o controle aplicado a planta: u(k) = u(k − 1) + ∆u(k)
7 Fazer k = k + 1
8 Aguardar o proximo instante de amostragem e retornar ao passo 1
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 3 / 39
![Page 4: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/4.jpg)
Topicos da aula de hoje
DMC: Sintonia de parametros (Perıodo de amostragem, M, N, ρ)
Implementacao em Matlab/Simulink
Exemplos
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 4 / 39
![Page 5: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/5.jpg)
Orientacoes para sintonia dos parametros
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 5 / 39
![Page 6: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/6.jpg)
DMC: Orientacoes para sintonia dos parametros
Referencias:
1 SHRIDHAR, R.; COOPER, D.J. A tuning strategy for unconstrainedSISO model predictive control. Industrial and Engineering ChemistryResearch, v. 36, n. 3, p. 729-746, 1997.
2 SHRIDHAR, R.; COOPER, D.J. A tuning strategy for unconstrainedmultivariable model predictive control. Industrial and EngineeringChemistry Research, v. 37, n. 10, p. 4003-4016, 1998.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 6 / 39
![Page 7: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/7.jpg)
Consideracoes gerais
O perıodo de amostragem T deve ser escolhido de modo a se ter umarepresentacao parcimoniosa da resposta a degrau da planta.
O horizonte de predicao N deve ser aproximadamente igual ao tempode acomodacao da resposta a degrau.
O horizonte de controle M usualmente e escolhido entre 1 e 6.
Valores menores de M tendem a deixar o controlador menos “agressivo”(resposta mais lenta, menor sensibilidade a ruıdo).
A diferenca N −M deve ser suficientemente grande (da ordem dotempo de acomodacao) para que o efeito das ultimas variacoes nocontrole se manifeste dentro do horizonte de predicao.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 7 / 39
![Page 8: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/8.jpg)
Orientacoes especıficas
Suponha que a dinamica da planta possa ser aproximada por um modelo deprimeira ordem com atraso descrito pela seguinte funcao de transferencia:
Y (s)
U(s)=
Ke−θs
τs + 1
Os parametros do DMC podem ser escolhidos da seguinte forma:
Perıodo de amostragem: T = max(0,1τ ; 0,5θ)
Horizontes de predicao e controle:
N = round
[5τ + θ
T+ 1
], M = round
[τ + θ
T+ 1
]
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 8 / 39
![Page 9: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/9.jpg)
“Regra de Bryson” para escolha do peso ρ
J(y,∆u) =N∑i=1
[y(k + i |k)− yref ]2︸ ︷︷ ︸(1)
+ ρ
M∑i=1
[∆u(k + i − 1|k)]2︸ ︷︷ ︸(2)
Ideia: Escolher ρ de modo que as parcelas (1) e (2) tenhamaproximadamente a mesma magnitude.
Dificuldade: Como saber, de antemao, os valores de y(k + i |k) e∆u(k + i − 1|k) a serem considerados para uso dessa regra ?
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 9 / 39
![Page 10: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/10.jpg)
J(y,∆u) =N∑i=1
[y(k + i |k)− yref ]2︸ ︷︷ ︸(1)
+ ρ
M∑i=1
[∆u(k + i − 1|k)]2︸ ︷︷ ︸(2)
Sugestao: Considerar a aplicacao de um degrau na entrada, comamplitude adequada para conduzir a saıda ate yref = 1.
Seja gss o valor de regime da resposta a degrau unitario. Para que a saıdaseja conduzida ate yref = 1, deve-se aplicar uma entrada degrau comamplitude 1/gss . Nesse caso, as parcelas (1) e (2) serao dadas por
(1) =N∑i=1
[g(i)
gss− 1
]2
, (2) = ρ
(1
gss
)2
Logo, sugere-se tomar
ρ = (gss)2N∑i=1
[g(i)
gss− 1
]2
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 10 / 39
![Page 11: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/11.jpg)
Observacao sobre a escolha de ρ
De maneira geral, valores maiores para o peso de controle ρ tendem adeixar o controlador menos “agressivo”.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 11 / 39
![Page 12: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/12.jpg)
Exemplo: Aquecedor de agua
ÁguaFria
ÁguaQuente
y
Aumento da temperaturaem relação ao ambienteu
Exemplo adaptado de CAMACHO, E. F.; BORDONS, C. Model PredictiveControl. London: Springer-Verlag, 1999.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 12 / 39
![Page 13: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/13.jpg)
Funcao de transferencia adotada na simulacao
H(s) =0,3
s + 0,2e−2s =
1,5
5s + 1e−2s =
K
τs + 1e−θs
K = 1,5 τ = 5 θ = 2
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 13 / 39
![Page 14: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/14.jpg)
Escolha dos parametros do controlador DMC
K = 1,5 τ = 5 θ = 2
Perıodo de amostragem:
T = max(0,1τ ; 0,5θ) = max(0,5; 1) = 1
Horizonte de predicao:
N = round
[5τ + θ
T+ 1
]= round
[27
1+ 1
]= 28
Horizonte de controle:
M = round
[τ + θ
T+ 1
]= round
[7
1+ 1
]= 8
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 14 / 39
![Page 15: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/15.jpg)
Uso da “Regra de Bryson”
ρ = (gss)2N∑i=1
[g(i)
gss− 1
]2
Diagrama de simulacao:
resposta degrau aquecedor.mdl
gss = 1.5
N = 28
rho = (gss^2)*sum((g(1:N)/gss - 1).^2)
Resultado: rho = 9.1
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 15 / 39
![Page 16: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/16.jpg)
Implementacao em Matlab/Simulink
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 16 / 39
![Page 17: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/17.jpg)
Arquivos Matlab a serem empregados
dmc.m
aquecedor dmc.mdl
resposta degrau aquecedor.mdl
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 17 / 39
![Page 18: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/18.jpg)
DMC: Implementacao em Matlab/Simulink (S-function)
Controlador PreditivoDMC
yref
y(k)u(k)
Interior damáscara
u(k)1
Zero-OrderHold
S-function
Parâmetros:KMPC, N, g, Ns, T
dmc
Armazena u(k-1)z1
y(k)2yref
1u(k)
u(k-1)
du(k)
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 18 / 39
![Page 19: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/19.jpg)
DMC: Detalhes da mascara
Parametros:
Peso do controle (escalar ρ)
Horizonte de predicao (escalar N)
Horizonte de controle (escalar M)
Resposta a degrau da planta (array g)
Perıodo de amostragem (escalar T )
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 19 / 39
![Page 20: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/20.jpg)
Comandos executados na inicializacao:
Ns = length(g);
gaux = g;
if N > Ns
gaux(Ns+1:N) = g(Ns);
end
col = gaux(1:N);
row = [gaux(1) zeros(1,N-1)];
G = toeplitz(col,row);
G = G(:,1:M);
GAIN = inv(G’*G + rho*eye(M))*G’;
KMPC = GAIN(1,:);
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 20 / 39
![Page 21: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/21.jpg)
Alguns comentarios sobre S-functions
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 21 / 39
![Page 22: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/22.jpg)
DMC: Codigo da S-function
Arquivo dmc.m
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 22 / 39
![Page 23: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/23.jpg)
Exemplo 1: Aquecedor de agua
ÁguaFria
ÁguaQuente
y
Aumento da temperaturaem relação ao ambienteu
Adaptado de CAMACHO, E.F.; BORDONS, C. Model Predictive Control.London: Springer-Verlag, 1999.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 23 / 39
![Page 24: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/24.jpg)
Roteiro para apresentacao do exemplo
1 Definir perıodo de amostragem no workspace.
2 Obter a resposta a degrau unitario em malha aberta:resposta degrau aquecedor.mdl
3 Abrir diagrama aquecedor dmc.mdl, inserir parametrosρ = 9.1,N = 28,M = 8 na mascara e executar a simulacao.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 24 / 39
![Page 25: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/25.jpg)
Simulacoes a serem realizadas:
Inserir perturbacao de saıda tipo degrau.
Inserir perturbacao de entrada tipo degrau.
Alterar ganho da planta.
Inserir perturbacao de saıda tipo rampa.
Variar ρ.
Inserir ruıdo de medida e variar ρ. Monitorar y na saıda da planta,antes do ponto de insercao do ruıdo.
Remover o ruıdo de medida. Fixar ρ = 0.91. Fixar N = 28 e variar Mde 8 ate 1.(Interpretacao: Calcular y = G∆u + f para k = 0 com M = 1).
Fixar M = 2 e variar N de 28 ate 2.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 25 / 39
![Page 26: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/26.jpg)
Exemplo 2: Pendulo de haste rıgida
m
l
up
m senyp
m
yp
yp
Entrada: Torque up
Saıda: Angulo yp
Parametros:
Massa m
Aceleracao da gravidade γ
Comprimento l
Coeficiente de atrito viscoso b
ml2yp = up −mlγsenyp − bl2yp ⇒ yp =upml2− γ
lsenyp −
b
myp
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 26 / 39
![Page 27: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/27.jpg)
yp =upml2− γ
lsenyp −
b
myp
No equilıbrio:
upml2− γ
lsenyp = 0
up = mγlsenyp
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 27 / 39
![Page 28: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/28.jpg)
Obtencao da resposta a degrau g
Plantau up yp
+
+
+
- y
up yp
Sistema inicialmente em repouso em k = 0 com yp = yp e up = up.
u(k) =
{Au, k ≥ 00, k < 0
,
Resposta a degrau para uso no DMC:
g(k) =y(k)
Au=
yp(k)− ypAu
, k > 0
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 28 / 39
![Page 29: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/29.jpg)
Malha de controle
ControladorDMC
Plantau
yref
y
up
up yp
yp
+
+
+
-yp,ref
+
-
yp
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 29 / 39
![Page 30: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/30.jpg)
Arquivos a serem usados
1 dmc.m
2 parametros pendulo.m
3 resposta degrau pendulo.mdl
4 pendulo dmc.mdl
Parametros do pendulo (unidades SI): m = 1, l = 1, γ = 10, b = 2
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 30 / 39
![Page 31: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/31.jpg)
Resposta a degrau
yp = 5π/180 rad
0 5 10 15 20 25 30 35 40 45 500
0.02
0.04
0.06
0.08
0.1
0.12
0.14g(n)
n
yp = 60π/180 rad
0 5 10 15 20 25 30 35 40 45 500
0.05
0.1
0.15
0.2
0.25
0.3
0.35g(n)
n
Perıodo de amostragem: T = 0,1s
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 31 / 39
![Page 32: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/32.jpg)
Horizontes adotados: N = 50, M = 10.
Peso do controle (com base na resposta a degrau partindo deyp = 60π/180 rad):
ρ = (gss)2N∑i=1
[g(i)
gss− 1
]2
>> gss = g(end)
>> N = 50
>> rho = (gss^2)*sum((g(1:N)/gss - 1).^2)
rho =
0.3190
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 32 / 39
![Page 33: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/33.jpg)
Remocao da constante de polarizacao na saıda da planta eno valor da referencia
ControladorDMC
Plantau
yref
y
up
up yp
yp
+
+
+
-yp,ref
+
-
yp
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 33 / 39
![Page 34: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/34.jpg)
Lei de controle: ∆u(k) = KMPC (r − f)
r =
yrefyref
...yref
=
yp,ref − ypyp,ref − yp
...yp,ref − yp
f =
f (k + 1|k)f (k + 2|k)
...f (k + N|k)
f (k + i |k) =
yp(k)−yp︷︸︸︷y(k) +
Ns∑n=1
[g(n + i)− g(n)]∆u(k − n)
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 34 / 39
![Page 35: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/35.jpg)
Lei de controle: ∆u(k) = KMPC (r − f)
r − f =
yp,ref − ypyp,ref − yp
...yp,ref − yp
−
yp(k)− yp +∑Ns
n=1 · · ·yp(k)− yp +
∑Nsn=1 · · ·
...
yp(k)− yp +∑Ns
n=1 · · ·
=
yp,refyp,ref
...yp,ref
−
yp(k) +∑Ns
n=1 · · ·yp(k) +
∑Nsn=1 · · ·
...
yp(k) +∑Ns
n=1 · · ·
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 35 / 39
![Page 36: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/36.jpg)
ControladorDMC
Plantau
yref
yup yp
+
+
+
-yp,ref
+
-
00up
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 36 / 39
![Page 37: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/37.jpg)
Remocao da constante de polarizacao na entrada da planta
ControladorDMC
Plantau
yup yp
+
+
+
-+
--
+
ControladorDMC
Plantau
yref
yup yp
+
+
+
-yp,ref
+
-
00
00
yp,ref yref
up
upup
A ausencia da constante de polarizacao na entrada da planta ecompensada pela acao integral do controlador DMC.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 37 / 39
![Page 38: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/38.jpg)
Resumo da aula de hoje
DMC: Sintonia de parametros
Implementacao em Matlab/Simulink
Exemplo 1: Aquecedor de agua
Exemplo 2: Pendulo de haste rıgida (modelo de simulacao nao linear)
Consideracoes sobre as constantes de polarizacao na malha decontrole DMC.
Vale salientar que a abordagem DMC requer convergencia da resposta adegrau em malha aberta.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 38 / 39
![Page 39: Aula 3 - Instituto Tecnológico de Aeronáuticakawakami/ee254/EE254_2019_2oSem_Aula3...Quente y Aumento da temperatura u em relação ao ambiente Adaptado de CAMACHO, E.F.; BORDONS,](https://reader036.fdocument.pub/reader036/viewer/2022090611/607bcd8c99e6e02565262174/html5/thumbnails/39.jpg)
Topico da proxima aula
Uso de funcoes de transferencia.
EE-254 (Controle Preditivo) Aula 3 13 Ago 2019 39 / 39