Difference between revisions of "Double Tank"
FelixMueller (Talk | contribs) (→Parameters) |
FelixMueller (Talk | contribs) (→Reference Solution) |
||
| Line 56: | Line 56: | ||
== Reference Solution == | == Reference Solution == | ||
The problem was solved in Matlab (version 2014b) with [[:Category: switch | Switch]]. The switch codes rely on [http://homepages.laas.fr/henrion/software/gloptipoly3/ GloptiPoly 3] and the solver [http://sedumi.ie.lehigh.edu/ SeDuMi 1.3] (optimization over symmetric cones). | The problem was solved in Matlab (version 2014b) with [[:Category: switch | Switch]]. The switch codes rely on [http://homepages.laas.fr/henrion/software/gloptipoly3/ GloptiPoly 3] and the solver [http://sedumi.ie.lehigh.edu/ SeDuMi 1.3] (optimization over symmetric cones). | ||
| − | By introducing the lifts <math>l_i=\sqrt{x_i}</math>, algebraically constrained as <math>l_i^2=x_i, \; l_i\geq 0,</math> the problem is recast with polynomial data. In this way way switch in connection with GloptiPoly3 can be applied. GloptiPoly 3 is a Matlab package developed by Didier Henrion for generalized problems of moments and needs polynomial input data. The exact code used to solve the problem can be found under [[Double Tank | + | By introducing the lifts <math>l_i=\sqrt{x_i}</math>, algebraically constrained as <math>l_i^2=x_i, \; l_i\geq 0,</math> the problem is recast with polynomial data. In this way way switch in connection with GloptiPoly3 can be applied. GloptiPoly 3 is a Matlab package developed by Didier Henrion for generalized problems of moments and needs polynomial input data. The exact code used to solve the problem can be found under [[Double Tank (switch)]]. The calculated objective is 4.7296 with the following trajectories of states and controls: |
<gallery caption="Reference solution plots" widths="180px" heights="140px" perrow="2"> | <gallery caption="Reference solution plots" widths="180px" heights="140px" perrow="2"> | ||
Image:State_sequence_Tank.jpg| State trajectory over time. | Image:State_sequence_Tank.jpg| State trajectory over time. | ||
Revision as of 14:12, 17 February 2016
| Double Tank | |
|---|---|
| State dimension: | 1 |
| Differential states: | 2 |
| Discrete control functions: | 1 |
| Interior point equalities: | 2 |
The double tank problem is a basic example for a switching system. It contains the dynamics of an upper and a lower tank, connected to each other with a pipe. The goal is to minimize the deviation of a certain fluid level 
in the lower tank. The problem was introduced and discussed in a variety of publications for the optimal control of constrained switched systems, e.g. Henrion et al. and Vasudevan et al.
Mathematical formulation
![\begin{array}{llll}
\displaystyle \min_{\sigma} & \displaystyle \int_{0}^{T} & k_1(x_2-k_2)^2 \; \text{d}t\\[1.5ex]
\mbox{s.t.} & \dot{x}_1 & = c_{\sigma(t)}-\sqrt{x_1}, \\[1.5ex]
& \dot{x}_2 & = \sqrt{x_1}-\sqrt{x_2}, \\[1.5ex]
& x(0) & = (2,2)^T, \\[1.5ex]
& \sigma(t) & \in \{1,2\} \qquad &\text{for } t\in[0,T].\\[1.5ex]
\end{array}](https://mintoc.de/images/math/e/5/1/e510dbcc83fc45b2984d452d9298c738.png)
The two states of the system correspond to the fluid levels of an upper and a lower tank.
The output of the upper tank flows into the lower tank, the output of the lower
tank exits the system, and the flow into the upper tank is restricted to either 
[lt/s]
or 
[lt/s]. The dynamics in each mode are then derived using Torricelli’s law, as
shown in constraints 1 and 2. The objective of the optimization is to have the fluid level in the
lower tank equal to [m], as reflected in the cost function.
Parameters
In an exemplary test, the parameters were chosen to be:
| Symbol | Value |
|
|
|
|
|
|
|
|
|
|
|
|
|
Reference Solution
The problem was solved in Matlab (version 2014b) with Switch. The switch codes rely on GloptiPoly 3 and the solver SeDuMi 1.3 (optimization over symmetric cones).
By introducing the lifts 
, algebraically constrained as 
the problem is recast with polynomial data. In this way way switch in connection with GloptiPoly3 can be applied. GloptiPoly 3 is a Matlab package developed by Didier Henrion for generalized problems of moments and needs polynomial input data. The exact code used to solve the problem can be found under Double Tank (switch). The calculated objective is 4.7296 with the following trajectories of states and controls:
- Reference solution plots
State trajectory over time.
Control sequence over time.
Source Code
Model descriptions are available in
References
There were no citations found in the article.
CLAEYS, Mathieu; DAAFOUZ Jamal; HENRION Didier Modal occupation measures and LMI relaxations for nonlinear switched systems control. arXiv preprint arXiv:1404.4699 (2014)
VASUDEVAN, Ramanarayan, et al. Consistent Approximations for the Optimal Control of Constrained Switched Systems---Part 2: An Implementable Algorithm. SIAM Journal on Control and Optimization, 2013, 51. Jg., Nr. 6, S. 4484-4503.
