Namespaces
Variants
Actions

Difference between revisions of "Travelling-wave method"

From Encyclopedia of Mathematics
Jump to: navigation, search
(Importing text file)
 
m (tex encoded by computer)
Line 1: Line 1:
 +
<!--
 +
t0940501.png
 +
$#A+1 = 30 n = 0
 +
$#C+1 = 30 : ~/encyclopedia/old_files/data/T094/T.0904050 Travelling\AAhwave method
 +
Automatically converted into TeX, above some diagnostics.
 +
Please remove this comment and the {{TEX|auto}} line below,
 +
if TeX found to be correct.
 +
-->
 +
 +
{{TEX|auto}}
 +
{{TEX|done}}
 +
 
A direct method for the numerical solution of problems in the calculus of variations (cf. [[Variational calculus|Variational calculus]]). It is used for solving problems of [[Optimal control|optimal control]] of small dimensions, but with constraints on the phase coordinates and control functions. After discretization of the functional and the system of differential equations, the original problem is reduced to the minimization of a functional:
 
A direct method for the numerical solution of problems in the calculus of variations (cf. [[Variational calculus|Variational calculus]]). It is used for solving problems of [[Optimal control|optimal control]] of small dimensions, but with constraints on the phase coordinates and control functions. After discretization of the functional and the system of differential equations, the original problem is reduced to the minimization of a functional:
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940501.png" /></td> <td valign="top" style="width:5%;text-align:right;">(1)</td></tr></table>
+
$$ \tag{1 }
 +
= \sum _ { k= } 0 ^ { N- }  1
 +
F ^ { 0 } ( x _ {k} , x _ {k+} 1 , u _ {k} , t _ {k} , t _ {k+} 1 ),
 +
$$
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940502.png" /></td> <td valign="top" style="width:5%;text-align:right;">(2)</td></tr></table>
+
$$ \tag{2 }
 +
x _ {k+} 1 - x _ {k}  = \tau F( x _ {k} , x _ {k+} 1
 +
, u _ {k} , t _ {k} , t _ {k+} 1 ),
 +
$$
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940503.png" /></td> </tr></table>
+
$$
 +
= 0 \dots N - 1,
 +
$$
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940504.png" /></td> <td valign="top" style="width:5%;text-align:right;">(3)</td></tr></table>
+
$$ \tag{3 }
 +
( x _ {k} , u _ {k} )  \in  G _ {k} ,\  k = 0 \dots N .
 +
$$
  
Here <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940505.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940506.png" /> are vectors of the phase coordinates and the controls, respectively, at the node <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940507.png" /> (with respective dimensions <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940508.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t0940509.png" />), <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405010.png" /> being considered as constant on each interval <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405011.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405012.png" /> are given domains of an <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405013.png" />-dimensional space (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405014.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405015.png" /> describe the boundary conditions), and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405016.png" /> is the subdivision step of the initial interval <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405017.png" />.
+
Here $  x _ {k} $,  
 +
$  u _ {k} $
 +
are vectors of the phase coordinates and the controls, respectively, at the node t _ {k} $(
 +
with respective dimensions $  n $
 +
and $  m $),  
 +
$  u _ {k} $
 +
being considered as constant on each interval $  ( t _ {k} , t _ {k+ 1 }  ) $,
 +
$  G _ {k} $
 +
are given domains of an $  ( n + m) $-
 +
dimensional space ( $  G _ {0} $
 +
and $  G _ {N} $
 +
describe the boundary conditions), and $  \tau = ( T - t _ {0} )/N $
 +
is the subdivision step of the initial interval $  T - t _ {0} $.
  
The travelling-wave method is used in the case <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405018.png" />, which is typical of practical problems, and for which the use of other methods based on the variation of spatial states (cf. [[Travelling-tube method|Travelling-tube method]]; [[Local variations, method of|Local variations, method of]]) is difficult because of the labour involved in determining the control function.
+
The travelling-wave method is used in the case $  n \geq  m $,  
 +
which is typical of practical problems, and for which the use of other methods based on the variation of spatial states (cf. [[Travelling-tube method|Travelling-tube method]]; [[Local variations, method of|Local variations, method of]]) is difficult because of the labour involved in determining the control function.
  
The given initial approximation <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405019.png" />, which satisfies (2) and (3), is improved in the sense of criterion (1) on each segment between <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405020.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405021.png" /> (<img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405022.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405023.png" /> are given), and this segment is gradually shifted by one node each time from the beginning of the trajectory to its end and back again (hence the name  "travelling wave" ).
+
The given initial approximation $  ( x _ {0}  ^ {o} \dots x _ {N}  ^ {o} , u _ {0}  ^ {o} \dots u _ {N- 1 }  ^ {o} ) $,  
 +
which satisfies (2) and (3), is improved in the sense of criterion (1) on each segment between t _ {k} $
 +
and t _ {k+ p }  $(
 +
$  x _ {k} $,  
 +
$  x _ {k+ p }  $
 +
are given), and this segment is gradually shifted by one node each time from the beginning of the trajectory to its end and back again (hence the name  "travelling wave" ).
  
 
One obtains, for each wave, a non-linear programming problem: The minimization of
 
One obtains, for each wave, a non-linear programming problem: The minimization of
  
<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405024.png" /></td> <td valign="top" style="width:5%;text-align:right;">(4)</td></tr></table>
+
$$ \tag{4 }
 +
\Delta I _ {k}  = \
 +
\sum _ { l= } k ^ { k+ }  p- 1
 +
F ^ { 0 } ( x _ {l} , x _ {l+} 1 , u _ {l} , t _ {l} , t _ {l+} 1 )
 +
$$
  
with <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405025.png" /> constraints of the type (2) and under the conditions (3). In the practical realization of the travelling-wave method no attempt is made to solve problem (4); instead, the increments <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405026.png" /> are given to each one of the <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405027.png" /> free parameters and, if <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405028.png" /> (in equation (4)) decreases and the conditions (3) are fulfilled, a new trajectory is obtained. If the trajectory remains unchanged during a complete period of the wave, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405029.png" /> is replaced by a fraction of it.
+
with $  p $
 +
constraints of the type (2) and under the conditions (3). In the practical realization of the travelling-wave method no attempt is made to solve problem (4); instead, the increments $  \pm  h _ {i} $
 +
are given to each one of the $  r $
 +
free parameters and, if $  \Delta I _ {k} $(
 +
in equation (4)) decreases and the conditions (3) are fulfilled, a new trajectory is obtained. If the trajectory remains unchanged during a complete period of the wave, $  h _ {i} $
 +
is replaced by a fraction of it.
  
If <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/t/t094/t094050/t09405030.png" />, the travelling-wave method becomes identical with the method of local variations.
+
If $  m = n $,  
 +
the travelling-wave method becomes identical with the method of local variations.
  
 
====References====
 
====References====

Revision as of 08:26, 6 June 2020


A direct method for the numerical solution of problems in the calculus of variations (cf. Variational calculus). It is used for solving problems of optimal control of small dimensions, but with constraints on the phase coordinates and control functions. After discretization of the functional and the system of differential equations, the original problem is reduced to the minimization of a functional:

$$ \tag{1 } I = \sum _ { k= } 0 ^ { N- } 1 F ^ { 0 } ( x _ {k} , x _ {k+} 1 , u _ {k} , t _ {k} , t _ {k+} 1 ), $$

$$ \tag{2 } x _ {k+} 1 - x _ {k} = \tau F( x _ {k} , x _ {k+} 1 , u _ {k} , t _ {k} , t _ {k+} 1 ), $$

$$ k = 0 \dots N - 1, $$

$$ \tag{3 } ( x _ {k} , u _ {k} ) \in G _ {k} ,\ k = 0 \dots N . $$

Here $ x _ {k} $, $ u _ {k} $ are vectors of the phase coordinates and the controls, respectively, at the node $ t _ {k} $( with respective dimensions $ n $ and $ m $), $ u _ {k} $ being considered as constant on each interval $ ( t _ {k} , t _ {k+ 1 } ) $, $ G _ {k} $ are given domains of an $ ( n + m) $- dimensional space ( $ G _ {0} $ and $ G _ {N} $ describe the boundary conditions), and $ \tau = ( T - t _ {0} )/N $ is the subdivision step of the initial interval $ T - t _ {0} $.

The travelling-wave method is used in the case $ n \geq m $, which is typical of practical problems, and for which the use of other methods based on the variation of spatial states (cf. Travelling-tube method; Local variations, method of) is difficult because of the labour involved in determining the control function.

The given initial approximation $ ( x _ {0} ^ {o} \dots x _ {N} ^ {o} , u _ {0} ^ {o} \dots u _ {N- 1 } ^ {o} ) $, which satisfies (2) and (3), is improved in the sense of criterion (1) on each segment between $ t _ {k} $ and $ t _ {k+ p } $( $ x _ {k} $, $ x _ {k+ p } $ are given), and this segment is gradually shifted by one node each time from the beginning of the trajectory to its end and back again (hence the name "travelling wave" ).

One obtains, for each wave, a non-linear programming problem: The minimization of

$$ \tag{4 } \Delta I _ {k} = \ \sum _ { l= } k ^ { k+ } p- 1 F ^ { 0 } ( x _ {l} , x _ {l+} 1 , u _ {l} , t _ {l} , t _ {l+} 1 ) $$

with $ p $ constraints of the type (2) and under the conditions (3). In the practical realization of the travelling-wave method no attempt is made to solve problem (4); instead, the increments $ \pm h _ {i} $ are given to each one of the $ r $ free parameters and, if $ \Delta I _ {k} $( in equation (4)) decreases and the conditions (3) are fulfilled, a new trajectory is obtained. If the trajectory remains unchanged during a complete period of the wave, $ h _ {i} $ is replaced by a fraction of it.

If $ m = n $, the travelling-wave method becomes identical with the method of local variations.

References

[1] N.N. Moiseev, "Elements of the theory of optimal systems" , Moscow (1975) (In Russian)
[2] I.A. Vatel, A.F. Kononenko, "A numerical scheme for solving optimal control problems" USSR Comput. Math. Math. Phys. , 10 : 1 (1970) pp. 85–94 Zh. Vychisl. Mat. i Mat. Fiz. , 10 : 1 (1970) pp. 67–73
[3] I.A. Vatel', A.F. Kononenko, Algoritm. i Program. (Informatsion. Byull.) : 2 (1972) pp. 7

'

How to Cite This Entry:
Travelling-wave method. Encyclopedia of Mathematics. URL: http://encyclopediaofmath.org/index.php?title=Travelling-wave_method&oldid=13820
This article was adapted from an original article by I.B. VapnyarskiiI.A. Vatel (originator), which appeared in Encyclopedia of Mathematics - ISBN 1402006098. See original article