2014 dxdy logo

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

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




Начать новую тему Ответить на тему
 
 вычисление континуального интеграла
Сообщение24.10.2010, 20:49 


17/09/09
226
Помогите разобраться, кто знает. Нужно вычислить вот такое среднее, где $x_n(k)$ - фурье образ $n-$ой компоненты векторного (Бозе) поля $\Bold{X}(\Bold{r})$.
$<x_n(p)x_m(q)>=\frac{\int\prod_{k,i}dx_i(k)x_n(p)x_m(q)e^{-\sum x_i(-k)A_{ij}x_j(k)}}{\int\prod_{k,i}dx_i(k)e^{-\sum x_i(-k)A_{ij}x_j(k)}}$. Интуитивно понимаю, что должно получится что-то типа $\delta(p+q)(A_{ij})^{-1}_{nm}$. Кто-нибудь может рассказать детально как провести вычисления? Желательно, без диагонализации матрицы $A$, а просто в явном виде?

 Профиль  
                  
 
 Re: вычисление континуального интеграла
Сообщение25.10.2010, 00:14 
Заслуженный участник


21/08/10
2462
Kamaz в сообщении #365835 писал(а):
Помогите разобраться, кто знает. Нужно вычислить вот такое среднее, где $x_n(k)$ - фурье образ $n-$ой компоненты векторного (Бозе) поля $\Bold{X}(\Bold{r})$.
$<x_n(p)x_m(q)>=\frac{\int\prod_{k,i}dx_i(k)x_n(p)x_m(q)e^{-\sum x_i(-k)A_{ij}x_j(k)}}{\int\prod_{k,i}dx_i(k)e^{-\sum x_i(-k)A_{ij}x_j(k)}}$. Интуитивно понимаю, что должно получится что-то типа $\delta(p+q)(A_{ij})^{-1}_{nm}$. Кто-нибудь может рассказать детально как провести вычисления? Желательно, без диагонализации матрицы $A$, а просто в явном виде?


Это не сложно. Во-первых, поскольку фурье-образы (а значит комплексные величины), надо расписать через действительную и мнимую часть (так что интеграл "удвоится"). Далее надо выкинуть
$x_n(p)x_m(q) $ из под верхнего интеграла а вместо этого вставить $e^{\sum x_i J_i}$ Т.е. прейти к производящему функционалу средних заданного типа. Нужные вам средние получаются из этого функционала дифференцированием по J и последующей (!) подстановкой J=0. Производящий функционал легко и просто вычисляется сдвигом переменной интегрирования (сдвиг надо сделать так, чтобы исчез линейный по переменным интегрирования член в экспоненте). Вот, собственно, и все. Но пожалуй было бы проще сделать то же самое в координатном представление а к фурье-образам перейти уже в ответе.

 Профиль  
                  
 
 Re: вычисление континуального интеграла
Сообщение25.10.2010, 06:32 


17/09/09
226
Спасибо за совет. Итак, следуя вашему совету у меня получилось:
$<x_n(r')x_m(r'')>=\frac{\delta ^2}{\delta J_n(r') \delta J_m(r'')}e^{\int dr J_i(r)(A^{-1})_{ij}J_j(r)}=
(A^{-1})_{nm}(r')\delta (r'-r'')+(A^{-1})_{mn}(r')\delta (r'-r'')$

В последнем равенстве положено уже $J=0$. Переходя в фурье представление получаем:
$ (A^{-1})_{nm}(p+p')+(A^{-1})_{mn}(p+p')$

Странный ответ...что-то тут не то...

-- Пн окт 25, 2010 10:48:02 --

А вот если не переходить в координатное пространство, то ответ получается такой $[(A^{-1})_{nm}(p)+(A^{-1})_{mn}(p)]\delta (p+p')$

 Профиль  
                  
 
 Re: вычисление континуального интеграла
Сообщение25.10.2010, 12:38 
Заслуженный участник


21/08/10
2462
Kamaz в сообщении #365926 писал(а):
Странный ответ...что-то тут не то...


Где-то вы наврали (разбирайтесь сами). Куда дельта-функция делась? Подсказка: она возникает из второго преобразования Фурье (два преобразования по двум координатам, одно снимается координатной дельта-функцией а второе дает дельта-функцию сохранения импульса).

Если бы был нелокальный случай, то вместо координатной дельта-функции была бы просто некая функция от разности. Делаем замену переменных: переходим к разности и сумме координат. Фурье-преобразование по сумме все равно дает сохранение импульса.

 Профиль  
                  
 
 Re: вычисление континуального интеграла
Сообщение25.10.2010, 14:51 


17/09/09
226
Вы все правильно говорите, и так действительно было бы, но только если бы $A^{-1}$ не зависила от координат. первое интегрирование снимается дейсвительно координатной дельта функцией, и получаем $\int dr' A^{-1}(r')e^{i(p+p')r}=A^{-1}(p+p')$

 Профиль  
                  
 
 Re: вычисление континуального интеграла
Сообщение25.10.2010, 20:50 
Заслуженный участник


21/08/10
2462
Kamaz в сообщении #366034 писал(а):
Вы все правильно говорите, и так действительно было бы, но только если бы $A^{-1}$ не зависила от координат. первое интегрирование снимается дейсвительно координатной дельта функцией, и получаем $\int dr' A^{-1}(r')e^{i(p+p')r}=A^{-1}(p+p')$


Ах у вас нелокальный случай... Тогда наврано в производящем функционале. Должно быть $e^{\int\int J_i(r')K_{ij}(r'-r)J_j(r)dr'dr}$. Никакой координатной дельта-функции в этом случае не будет. И просто к обратной матрице ядро не сводится.

Если же K зависит не от разности координат а более общая зависимость (например вида $F(r) \delta (r-r')$), то закона сохранения импульса и быть не может -- нет трансляционной инвариантности.

В общем проверьте все внимательно с самого начала. Простая небрежность в вычислениях. Тщательнее нужно:-)

 Профиль  
                  
 
 Re: вычисление континуального интеграла
Сообщение26.10.2010, 05:21 


17/09/09
226
Действительно, небрежность! Я уже разобрался, спасибо за помощь!!

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

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



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

Сейчас этот форум просматривают: Geen


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

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