2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 Восстановить Гамильтониан по траекториям.
Сообщение08.11.2010, 09:33 
Заслуженный участник
Аватара пользователя


30/10/10
1481
Ереван(3-й участок)
Возник такой вопрос,

на симплекическом многообразии $M$, параметризованном координатами $y$ задано семейство гладких кривых $y^i=y^i(t,C_1,\ldots,C_{2n})$, таких, что при разных параметрах $C_1,\ldots,C_{2n}$ эти кривые не пересекаются. Можно ли нарисовать функцию $H(y)$(Гамильтониан), такую, что эти кривые были решением системы уравнений

$\dot{y}^i=\left\{H,y^i\right\}$,

и главное, как это сделать? Можно ли это сделать неявно?

Насколько я понимаю, таких Гамильтонианов много, но может возможно представить какой-то класс?

(Оффтоп)

На всякий случай: Симплектическая 2-форма локально записывается как $\omega^{(2)}=\omega_{ij}dy^i\wedge dy^j$. Если определить $\omega^{ij}\omega_{jk}=\delta^i_k$, то $\left\{H,y^i\right\}\equiv \omega^{ij}\frac{\partial H}{\partial y^j}$
$\delta^i_k$- симвло Кронекера. По повторяющимся индексам подразумевается суммирование.


Заранее спасибо.

-- Пн ноя 08, 2010 10:52:13 --

Да кстати, по моему эти кривые дожны еще либо не пересекаться сами с собой, либо если пересеклись один раз, то... вобщем не знаю как это словами сформулировать. На языке формул
Если для какого-нибудь $t$ имеет место $y(t,C_1,\ldots,C_{2n})=y(t+T,C_1,\ldots,C_{2n})$, то это так и для любого $t^\prime$.
Т.е. частица всегда должна знать куда ей идти.

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение08.11.2010, 10:03 


02/10/10
376
Если заданы не просто кривые а поток, то вопрос о (локальной) гамильтоновости простой: проверьте сохраняет ли этот поток симплектическую структуру, если речь идет только о кривых, то являются ли эти кривые интегральными кривыми некоторого гамильтонового потока, это вопрос очень сложный, вообще говоря

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение08.11.2010, 15:47 
Заслуженный участник
Аватара пользователя


30/10/10
1481
Ереван(3-й участок)
moscwicz
Спасибо за ответ. А до меня этим вопросом никто не занимался? Может книжка есть какая или сатья? Вы не знаете?

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение12.11.2010, 13:54 


02/10/10
376
Всетаки вопрос вполне обозримый. Во всяком случае локально. Раз семейство кривыхзадано, имеем векторное поле $v(x)$, которое определено с точностью до умножение каждого вектора на число: $\lambda(x)v(x)$.
Если кривые не имеют особенностей, то $v\ne 0$. Пусть $\omega$ -- симплектическая форма, тогда надо написать уравнение $L_{\lambda v}\omega=0$. ($L$ -- производная Ли.) Это система урчп на функцию $\lambda$ Дальше все стандартно: условия разрешимости этой системы писать надо.
В случае двумерного многообразия эта система превращается в единственное уравнение, которое разрешимо в окрестности любой точки где $v\ne 0$

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение12.11.2010, 17:58 
Заслуженный участник
Аватара пользователя


30/10/10
1481
Ереван(3-й участок)
moscwicz в сообщении #373956 писал(а):
имеем векторное поле $v(x)$, которое определено с точностью до умножения каждого вектора на число: $\lambda(x)v(x)$.

А что такое $v(x)$в терминах первого поста?

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение12.11.2010, 18:10 


02/10/10
376
v(x) это поле касатеольных векторов к кривым

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение12.11.2010, 18:18 
Заслуженный участник
Аватара пользователя


30/10/10
1481
Ереван(3-й участок)
moscwicz в сообщении #374103 писал(а):
v(x) это поле касатеольных векторов к кривым

А понятно. т.е. $v(y)$. Но как это поле восстановить? Можно на примере? Имеем

$y_1=C_2\sin{(\omega t+C_1)},\quad y_2=C_2\omega\cos{(\omega t+C_1)}$
Как из них построить это векторное поле?

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение12.11.2010, 22:07 
Заслуженный участник
Аватара пользователя


30/10/10
1481
Ереван(3-й участок)
Наверно так:
Выражаем $C_1,\ldots,C_{2n}$, через $y$ при $t=0$.
Таким образом, кривая, проходящя через точку $y_0$(заметьте индекс снизу) есть
$y=y(C_i(y_0),t)$.
Далее, векторное поле в точке $y_0$ есть $v(y_0)=\left.\frac{\partial y^i(C_k(y_0))}{\partial t}\right|_{t=0}\frac{\partial}{\partial y^i}$.

Теперь, можем написать
$\omega^{(2)}(v(y))=dH,$
где $\omega^{(2)}$-симплектическая 2-форма.

Вроде так.

 Профиль  
                  
 
 Re: Восстановить Гамильтониан по траекториям.
Сообщение12.11.2010, 23:40 
Заслуженный участник
Аватара пользователя


30/10/10
1481
Ереван(3-й участок)
Пример:
$y_1=C_2\sin{(\omega t+C_1)},\quad y_2=C_2\omega\cos{(\omega t+C_1)}$

Для любой точки $y_0$:
$C_1(y_0)=\arctan{\left(\omega \frac{y^1_0}{y^2_0}}\right),\quad C_2=\sqrt{(y^1_0)^2+(y^2_0/\omega)^2}$
Соответственно, имеем

$y^1=\sqrt{(y^1_0)^2+(y^2_0/\omega)^2}\sin{\left(\omega t+\arctan{\left(\omega \frac{y^1_0}{y^2_0}}\right)\right)},\quad y^2=\omega\sqrt{(y^1_0)^2+(y^2_0/\omega)^2}\cos{\left(\omega t+\arctan{\left( \omega\frac{y^1_0}{y^2_0}\right)}\right)}$
$

И для векторного поля имеем

v(y_0)=\omega\sqrt{(y^1_0)^2+(y^2_0/\omega)^2}\cos{\left(\arctan{\left(\omega\frac{y^1_0}{y^2_0}}\right)\right)}\frac{\partial}{\partial y^1}-\omega^2\sqrt{(y^1_0)^2+(y^2_0/\omega)^2}\sin{\left(\arctan{\left(\omega\frac{y^1_0}{y^2_0}}\right)\right)}\frac{\partial}{\partial y^2}$

Теперь, для симплектичексой формы берем обычную антисимметричную форму
$dy_1 dy_2-dy_2 dy_1$
и из уравнения$ \omega^{(2)}(v(y))=dH$
получаем

$\omega\sqrt{(y^1)^2+(y^2/\omega)^2}\cos{\left(\arctan{\left(\omega\frac{y^1}{y^2}}\right)\right)}=y^2=\frac{\partial H}{\partial y^2}$

$\omega^2\sqrt{(y^1)^2+(y^2/\omega)^2}\sin{\left(\arctan{\left(\omega\frac{y^1}{y^2}}\right)\right)}=\omega y^1=\frac{\partial H}{\partial y^1}$

Вот, видимо, и все. Решение этой системы уравнений есть

$H=\frac{(y^2)^2+\omega^2 (y^1)^2}{2}$

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 9 ] 

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group