I am working on sub-Riemannian geometry and try to understand what are the tools to find the equations of a sub-Riemannian problem. Here is an example:


Let us consider the system defined by a lagrangian:


\begin{equation}
L=\frac{1}{2}m(\overset{\cdot}{x}^2+\overset{\cdot}{y}^2)+\frac{1}{2}I\overset{\cdot}{\theta}^2
\end{equation}


And take the following non-holonomic constraint


\begin{equation}
\mathscr{D}=span\left\{\cos\theta\frac{\partial}{\partial x}+\sin\theta\frac{\partial}{\partial y},\frac{\partial}{\partial\theta}\right\}
\end{equation}
So we have,
\begin{equation}
\mathscr{D}^\perp=span\{\omega\}
\end{equation}
with $\omega=-\sin\theta dx+\cos\theta dy$


If I understand clearly the sub-Riemannian problem is to minimize
\begin{equation}
\int\_0^TLdt
\end{equation}
under the constaint $\dot{z}(t)\in\mathscr{D}(z(t))$. Now, how do I find these equations? So far I am not asking how to solve these equations but just to understand how to find them in the first place.


Thanks in advance for your help.


Edit:


I quite understand your answer, thank you for the clear explanation. I am though not really familiar with canonical coordinates. I agree with the expression of the Hamiltonian.


Now, I am right if I say that the Hamilton's equations are given by:


\begin{equation}
 \begin{cases}
 \overset{\cdot}{x}(t)&=\frac{\partial H}{\partial p\_x}=p\_x\cos^2\theta +p\_y\cos\theta\cdot\sin\theta \\
 \overset{\cdot}{y}(t)&=\frac{\partial H}{\partial p\_y}=p\_y\sin^2\theta +p\_x\cos\theta\cdot\sin\theta \\
 \overset{\cdot}{\theta}(t)&=\frac{\partial H}{\partial p\_\theta}=p\_\theta\\
 \overset{\cdot}{p\_x}(t)&=-\frac{\partial H}{\partial x}= 0 \\
 \overset{\cdot}{p\_y}(t)&=-\frac{\partial H}{\partial y}= 0 \\
 \overset{\cdot}{p\_\theta}(t)&=-\frac{\partial H}{\partial \theta}= (p\_x^2-p\_y^2)\sin(\theta)cos(\theta)+p\_xp\_y(\sin^2(\theta)-\cos^2(\theta)) \\
 \end{cases}
\end{equation}


I'm not sure to understand how to interpret them.
Are these the equations of the sub-Riemannian problem?



How to find equations of a sub-Riemannian problem
You can reformulate your problem in the language of geometric control (where your dynamical system is sometimes called Dubin's car). Let
$$
X\_1=\frac{1}{\sqrt{m}}(\cos\theta\partial\_x + \sin\theta\partial\_y),\qquad X\_2=\frac{1}{\sqrt{I}}\partial\_\theta$$
be a basis of your distribution. Horizontal curves are trajectories $z :[0,T]\to M$ (here $M=\mathbb{R}^2\times S^1$), that are tangent to the distribution, that is there exist $u\_1,u\_2:[0,T]\to \mathbb{R}$ (let us forget about regularity issues), such that
$$
\dot z(t) = u\_1(t) X\_1|\_{\gamma(t)} + u\_2(t)X\_2|\_{\gamma(t)}, \qquad t\in [0,1].
$$
You look for horizontal trajectories between fixed endpoints, that minimize the cost functional:
$$
\int L(\gamma(t))dt = \frac{1}{2}\int\_0^T \left(u\_1(t)^2 + u\_2(t)^2\right)dt.
$$
(this should clarify the normalization of $X\_1,X\_2$).


This is a very simple optimal control problem. An application of the Pontryagin Maximum Principle yields first order necessary conditions for solutions of this problem. In this very simple all solutions correspond to projections on $M$ of integral trajectories of an Hamiltonian system on $T^\*M$, with Hamiltonian given by


$$
H(\lambda) = \frac{1}{2}\sum\_{i=1}^2 \lambda(X\_i)^2, \qquad \lambda \in T^\*M.
$$


In canonical coordinates $(x,y,\theta,p\_x,p\_y,p\_\theta)$ on $T^\*M$, your Hamiltonian reads explicitly


$$
H=\frac{1}{2}(p\_x \cos\theta + p\_y \sin\theta)^2 + \frac{1}{2}p\_\theta^2.
$$


Then your "geodesic equations" correspond to the Hamilton's equations for $H$.


A good reference to this general kind of optimal control problems can be found in the following book (chapter 12, Theorem 12.10, case (b) is vacuous in your case, so only case (a) remains):


*Agrachev, Andrei A.; Sachkov, Yuri L.*, Control theory from the geometric viewpoint., Encyclopaedia of Mathematical Sciences 87. Control Theory and Optimization II. Berlin: Springer (ISBN 3-540-21019-9/hbk). xiv, 412 p. (2004). [ZBL1062.93001](https://zbmath.org/?q=an:1062.93001).


For a specifically sub-Riemannian reference, see Chapter 4.3 in


*Agrachev, Andrei; Barilari, Davide; Boscain, Ugo*, [A comprehensive introduction to sub-Riemannian geometry. From the Hamiltonian viewpoint.](http://dx.doi.org/10.1017/9781108677325), [ZBL07073879](https://zbmath.org/?q=an:07073879).


The references proposed by Robert Bryant are also very good, but the two aforementioned one are closer (I think) to the spirit of your question.


