2014 dxdy logo

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

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




Начать новую тему Ответить на тему
 
 Малые колебания систем (ТеорМех)
Сообщение10.12.2008, 16:43 
Аватара пользователя


13/11/07
56
Привет! Не понимаю как перейти к нормальным координатам в след задаче:
Для заданного потенциала нужно найти частоты малых колебаний, перейдя к нормальным координатам.
\[
U = \frac{k}{2} \cdot (x^2  + y^2 ) + \frac{{U_0 a^2 }}{{a^2  + xy}}
\]

Взял первые производные:
\[
\begin{array}{l}
 \frac{{\partial U}}{{\partial x}} = kx - \frac{{U_0 a^2 y}}{{(a^2  + xy)^2 }} \\
 \frac{{\partial U}}{{\partial y}} = ky - \frac{{U_0 a^2 x}}{{(a^2  + xy)^2 }} \\ 
 \end{array}
\]
Получил (?): \[
x = 0,y = 0
\]
Теперь рассматривая потенциал вблизи (0;0) нужно как-то перейти к нормальным координатам пользуясь условием малости "х" и "y". А как пользоваться условием не пойму.
Пытался разложить \[
(1 + \frac{{xy}}{{a^2 }})^{ - 1} 
\] но пришел к тому откуда вышел. Все частные производные до 2 в (0;0) равны 0 кроме смешанной она в ней \[
{a^{ - 2} }
\] и умножается на "ху"

 Профиль  
                  
 
 
Сообщение10.12.2008, 18:38 


10/03/07

473
Москва
Shpilev в сообщении #166469 писал(а):
Получил (?): \[ x = 0,y = 0 \]
Правильно. Первым делом найти положение равновесия.


Shpilev в сообщении #166469 писал(а):
Пытался разложить \[ (1 + \frac{{xy}}{{a^2 }})^{ - 1} \]
Правильно. Вторым делом разложить потенциальную энергию в окрестности равновесия до квадратичных членов, то есть

$$
\frac{{U_0 a^2 }}{{a^2 + xy}}=U_0-U_0xy/a^2+\dots
$$

(первое слагаемое у Вас и так квадратичное). А теперь приводим квадратичную форму

$$
\frac{k}{2} \cdot (x^2 + y^2 ) -U_0xy/a^2
$$

к диагональному виду. Только тут одна тонкость есть: нужно еще кинетическую энергию знать. Если она $\sim\dot x^2+\dot y^2$, то просто диагонализуете потенциальную, а если нет... тогда еще интереснее.

 Профиль  
                  
 
 Re: Малые колебания систем (ТеорМех)
Сообщение10.12.2008, 22:13 
Заслуженный участник
Аватара пользователя


30/01/06
72407
Shpilev писал(а):
Все частные производные до 2 в (0;0) равны 0 кроме смешанной она в ней \[{a^{ - 2} }\] и умножается на "ху"

Это то что надо. Теперь матрицу этих коэффициентов
$\left(\begin{array}{cc}
0 & a^{-2}/2\\
a^{-2}/2 & 0\end{array}\right)$
надо повернуть так, чтобы она приняла диагональный вид. Это и будет приведение к нормальным координатам. О кинетической энергии можете не заморачиваться, если ничего сложного не оговорено, она такая, какая нужно.

 Профиль  
                  
 
 
Сообщение14.12.2008, 16:10 
Аватара пользователя


13/11/07
56
Я решил, спасибо!
Однако преобразовывал методом Лагранжа.
А где можно прочитать про преобразование квадратичной формы к нормальному виду при помощи матрицы (так что бы с практическим примером)?

 Профиль  
                  
 
 
Сообщение14.12.2008, 17:12 


10/03/07

473
Москва
Shpilev в сообщении #167529 писал(а):
Однако преобразовывал методом Лагранжа.
А кинетическая энергия при этом не запортится? Вообще-то предполагалось, что приводить будете ортогональным преобразованием.


Shpilev в сообщении #167529 писал(а):
А где можно прочитать про преобразование квадратичной формы к нормальному виду при помощи матрицы
А чего там читать? Решаете задачу на собственные значения/векторы, из векторов составляете матрицу преобразования, из значений --- преобразованную матрицу квадратичной формы.

 Профиль  
                  
 
 
Сообщение15.12.2008, 00:08 
Аватара пользователя


13/11/07
56
Запортилась ))) Я потом поправлял :)

 Профиль  
                  
 
 
Сообщение17.12.2008, 09:47 
Заслуженный участник
Аватара пользователя


30/01/06
72407
Shpilev в сообщении #167529 писал(а):
А где можно прочитать про преобразование квадратичной формы к нормальному виду при помощи матрицы (так что бы с практическим примером)?

Это в курсе линейной алгебры и аналитической геометрии должно быть.

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

Модераторы: photon, whiterussian, profrotter, Jnrty, Aer, Парджеттер, Eule_A, Супермодераторы



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

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


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

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