# Monique CHYBA, Gautier PICOT Graduate course on ... gautier/Picot June 18 2015... Monique CHYBA,...

date post

08-Jul-2020Category

## Documents

view

8download

0

Embed Size (px)

### Transcript of Monique CHYBA, Gautier PICOT Graduate course on ... gautier/Picot June 18 2015... Monique CHYBA,...

Introduction to control theory and applications

Monique CHYBA, Gautier PICOT

Department of Mathematics, University of Hawai'i at Manoa

Graduate course on Optimal control University of Fukuoka

18/06/2015

1/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

What is control theory ?

I Control theory is a branch of mathematics that studies the properties of control systems i.e dynamical systems whose behavior can be modi�ed by a command

I General Mathematical formalism of a control system :

ẋ(t) = f (t, x(t), u(t))

where I t ∈ [t0tf ] is the time variable, I x is the state variable de�ned on [t0tf ] and valued in a smooth variable M, I u is a measurable bounded function de�ned on [t0tf ], valued in a smooth

variable U, called the control variable, I f : R×M × U → TM is a smooth application.

I Goal : Bring the state variable from a given initial condition to a given �nal condition i.e solve a boundary value problem{

ẋ(t) = f (t, x(t), u(t)) x(t0) = x0 ∈ M, x(tf ) = xf ∈ M.

2/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Schematic diagram

Figure: A dynamical system controlled by a feedback loop. We call the error the di�erence between the reference (the desired output) and the the measured output. This error is used by the controller to design a control on the system so that the measured output gets closer to the reference.

Source : Wikipedia

3/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Questions that arise...

Controllability of a control system ? Is it possible to bring the state variable from any initial condition to any �nal condition in a �nite time ?

I linear systems { ẋ(t) = A(t)x(t) + B(t)u(t) + r(t)

x(t0) = x0

where x(t) Rn, u(t) Rn, A(t) ∈Mn(R), B(t) ∈Mn,m(R) and r(t) ∈Mn,1(R) for all t ∈ [t0, tf ].

→ Kalman condition constraints on the control ?

I nonlinear systems

→ way more di�cult → Poincaré's Reccurence theorem, Poisson-stability, linearization, local controllability

4/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Questions that arise...

Stabilization of a control system ? How can we make a control system insensitive to perturbations ?

Example : if (xe , ue) is an equilibrium point of the autonomous control system

ẋ(t) = f (x(t), u(t))

i.e f (xe , ue) = 0.

Does it exist a control u such that, for all � > 0, there exists η > 0 such that, for all x0 ∈ B(xe , η) and all t ≥ 0, the solution to the system{

ẋ(t) = f (t, x(t), u(t)) x(t0) = x0

satis�es ||x(t)− xe || ≤ � ? I linear systems → controlability I nonlinear systems → Lyapunov functions

5/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Questions that arise...

Optimal control ? Can we determine the optimal solutions to a control system for a given optimization criterion ?

→ Find the solution to the boundary value problem{ ẋ(t) = f (t, x(t), u(t))

x(t0) = x0 ∈ M0, x(tf ) = xf ∈ M1

which minimizes the cost

min u(.)∈U

∫ tf t0

f0(t, x(t), u(t))dt + g(tf , xf )

where f0 : R×M × U → R is smooth and g : R×M → R is continuous.∫ tf t0

f0(t, x(t), u(t))dt : Lagrange cost

g(tf , xf ) : Mayer cost

see

6/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Origins of optimal control theory : the Brachistochrone

In 1696, Johan Bernoulli challenged his contemporary with the following problem :

Consider two points A and B such that A is above B. Assume that a object is located a the point A with no initial velocity and is only subjected to the gravity. What is the curve between A and B so that the object travels from A to B in minimal time ?

Remark : We know that the straight line is the shortest way between two points. Is it the fastest way ?

NO !

The fastest way is a cycloid arc whose tangent line at the point A is vertical.

7/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Origins of optimal control theory : the Brachistochrone

A short movie which illustrates this result.

8/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Applications of optimal control theory : Treatment of HIV

Model of the interaction of HIV and T-cells in the immune system

→ dT (t) dt

= s1 − s2V (t)

B1 + V (t) − µT (t)− kV (t)T (t) + u1(t)T (t)

dV (t)

dt =

g(1− u2(t))V (t) B2 + V (t)

− cV (t)T (t).

where I T (t) : Concentration of una�ected T cells I V (t) : Concentration of HIV particles I (u1, u2) : Control terms, action of the treatment I s1 − s2V (t)B1+V (t) : proliferation of una�ected T cells I µT (t) : natural loss of una�ected T cells I kV (t)T (t) : loss by infection I g(1−u2(t))V (t)

B2+V (t) : proliferation of virus

I cV (t)T (t) : viral loss

9/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Applications of optimal control theory : Treatment of HIV

Objective : maximizing the e�ciency of the treatment

→ max (u1,u2)

∫ tf 0

T (t)− (A1u21(t) + A2u22(t))dt.

i.e

I maximizing the number of una�ected T cells during the treatment I and minimizing the systemic cost of the treatment

A1, A2 : 2 constants/weights A1u

2

1 (t) + A2u

2

2 (t) : severity of side e�ects of the treatment.

Results : The optimal control (u1, u2) can be written as a feeback control i.e function of T and V . The optimal synthesis can be simulated numerically.

10/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Applications of optimal control theory : a 2-sectors economy model

Context : An economy consisting of 2 sectors. The sector 1 produces �nancial goods and the sector 2 produces consumption goods. Denote x1(t) and x2(t) the productions in sectors 1 and 2 and u(t) the proportion of investment allocated to sector 1.

Assumption : Increase in production in each sector is proportional to the investment allocated to each sector.

Problem : Maximizing the total consumption over interval of time [0,T ].

→ Optimal control problem ẋ1(t) = αu(t)x1(t)

ẋ2(t) = α(1− u(t))x2(t) maxu(.)∈U

∫ T 0 x2(t)dt

x1(0) = a1, x2(0) = a2

where α is some constant of proportionality and a1 and a2 are the initial productions in sectors 1 and 2.

11/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Applications of optimal control theory : a 2-sectors economy model

This problem can be solved analytically. The optimal solution (u∗, x∗1 , u ∗ 2 )

is

u∗(t) =

{ 1 if 0 ≤ t ≤ T − α

2 0 if T − α

2 < t ≤ T

x∗1 (t) =

{ a1eαt if 0 ≤ t ≤ T − α2

a1eαT−2 if T − α2 < t ≤ T

x∗2 (t) =

{ a2 if 0 ≤ t ≤ T − α2

a2e(αt−αT+2)e αT−2

if T − α 2 < t ≤ T

12/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Identi�cation of the Fragmentation Role in the Amyloid Assembling Processes and Application to their Optimization

Context : Apply techniques from geometric control to a kinetic model of amyloid formation which will take into account the contribution of fragmentation to the de novo creation of templating interfaces to design optimal strategies for accelerating the current ampli�cation protocols, such as the Protein Misfolding Cyclic Ampli�cation (PMCA). The objective is to reduce the time needed to diagnose many neurodegenerative diseases.

Fibril fragmentation : Fibril fragmentation has been reported to enhance the polymerization process underlying the behavior of some speci�c prions. There is a signi�cant lack of knowledge concerning the fragmentation process and the de novo generation of templating interfaces, both in the mechanisms of its occurrence and its contribution to the acceleration of the pathology.

13/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Protein Misfolding

I Transmissible spongiform encephalopathies (TSEs)

I Aggregation-Fragmentation

14/32 Monique CHYBA, Gautier PICOT Optimal control and application to space transfers

Compartmental model of amyloid formation

We denote by x li (t), l = 1, · · · , ki , the density of polymers of size i in compartment l at a given time t. The corresponding rate of change due to elongation is then described as follows :

r(u(t)) [ ki−1∑

s=1

τ l,si−1x s i−1 −

ki+1∑ r=1

τ r,li x l i

] .

The parameter τ l,si−1 is the growth rate of polymers of size i − 1 in compartment s that grow in compartment l of polymers of size i , and τ r,li is the growth rate of polymers of size i in compartment l that grow into polymers of size i + 1 (in compartment r).

The c

Recommended

*View more*