Dielectrophoretic Particle OED

Dielectrophoretic Particle OED
State dimension:	1
Differential states:	13
Discrete control functions:	3

The Dielectrophoretic Particle OED problem is a variation of the Dielectrophoretic Particle problem. It looks for optimal time intervals to measure the two states in order to minimize the uncertainty of a follow-up parameter estimation problem for the two unknown parameters.

The mathematical equations form a small-scale ODE model. It also includes state sensitivities, the Fisher information matrix entries and integrated sampling states.

Mathematical formulation

We are interested in estimating the parameters $\alpha$ and $c$ of the initial value problem

$\begin{array}{rclll}\dot{x_1}(t)& =& x_2(t)\cdot u(t)+\alpha \cdot u(t{)}^2,& & t\in [0,t_f],x_1(0)=1,\\ \dot{x_2}(t)& =& -c\cdot x_2(t)+u(t),& & t\in [0,t_f],x_2(0)=0.\end{array}$

The initial values and $t_f=10$ are fixed. We are interested in how to choose the control $u$ and when to measure, with an upper bound $M$ on the measuring time. We can measure the states directly, $h^1(x(t))=x_1(t)$ and $h^2(x(t))=x_2(t)$. We use two different sampling functions, $w^1(\cdot )$ and $w^2(\cdot )$ in the same experimental setting. This can be seen either as a two-dimensional measurement function $h(x(t))$, or as a special case of a multiple experiment, in which $u(\cdot ),x(\cdot )$, and $G(\cdot )$ are identical.

Now we formulate the OED problem:

$\begin{array}{lll}\underset{x,G,F,z,w,u}{\min }& & \text{trace}(F^{-1}(t_f))\\ \text{subject to}\\ \dot{x}(t)& =& f(x(t),\theta )\\ \dot{G}(t)& =& f_x(x(t),\theta )G(t)+f_{\theta }(x(t),\theta )\\ \dot{F}(t)& =& \sum _{i=1}^{n_o}{w}_i(t)(h_x^i(x(t))G(t){)}^T(h_x^i(x(t))G(t))\\ \dot{z}(t)& =& w(t),\\ x(0)& =& x_0\\ G(0)& =& \frac{\partial x(0)}{\partial \theta }\\ F(0)& =& 0,\\ z(0)& =& 0\\ u(t)& \in & {𝒰}\\ w(t)& \in & {𝒲}\\ z_i(t_f)& \le & M_i\end{array}$

The evolution of the symmetric matrix $F:[0,t_f]\to {\mathbb{ℝ}}^{2\times 2}$ is given by the weighted sum of observability Gramians $h_y^i(y(t))G(t),i=1,2$ for each observed function of states.

Parameters

We use $t_f=8$, $\alpha =-0.75$, and $c=1$. The upper bound on the measurement time intervals is chosen as $M_1=M_2=2$.

Reference Solutions

Here is one local solution to the above control problem.

• Reference solution plots
• 
  States, control, and sampling functions for a local optimum.

References

There were no citations found in the article.