Compartmental OED

Compartmental OED
State dimension:	1
Differential states:	12
Discrete control functions:	1

The Compartmental OED problem looks for an optimal measurement strategy to determine a single parameter in a one-dimensional ODE model, where we can directly measure the single state. The following description is taken from TODO.

Here we consider the open one-compartment model with first-order absorption input fitted by Button (unpublished Ph.D. thesis, Texas A&M University, 1979) to the results of an experiment in which six horses each received 15 mg/kg of theophylline as aminophylline by intragastric administration.

The optimal integer control functions shows bang bang behavior.

Mathematical formulation

For a single parameter $p$ the original initial value problem is given by $\dot{x}(t)=:f(t,p)=0.2+0.8\cdot t+0.3\cdot (\sin (p\cdot t)+\cos (p\cdot t)\cdot p\cdot t)-2.5\cdot \sin (50\cdot t),x(0)=x_0.$

We assume both $x_0$ and $t_f$ to be fixed and are only interested in when to measure, with an upper bound $M$ on the measuring time. We can measure the state directly, i.e. $h(x(t))=x(t)$.

Now we formulate the OED problem:

$\begin{array}{lll}\underset{y,G,F,z,w}{\min }& & \text{trace}(F^{-1}(t_f))\\ \text{subject to}\\ \dot{y}(t)& =& f(t,p)\\ \dot{G}(t)& =& f_p(y(t),p)\\ \dot{F}(t)& =& w(t)(h_y(y(t))G(t){)}^T(h_y(y(t))G(t))\\ \dot{z}(t)& =& w(t),\\ y(0)& =& y_0\\ G(0)& =& 0\\ F(0)& =& 0,\\ z(0)& =& 0\\ w(t)& \in & {𝒲}\\ z(t_f)& \le & M\end{array}$

Parameters

These fixed values are used within the model:

$x_0=0.1;t_f=2;{𝒲}=[0,1];M=0.2;p=15$

Reference Solutions

Here is one local solution to the above control problem.

• Reference solution plots
• 
  States and measurement control for $\theta =(0.05884,4.298,21.80)$.

Miscellaneous and Further Reading

This problem was introduced by Sebastian Sager.