2014 dxdy logo

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

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




На страницу 1, 2  След.
 
 Определённый интеграл с синками
Сообщение20.12.2015, 02:57 
Аватара пользователя
Подскажите, пожалуйста, каким путём можно посчитать такой интеграл:

$$\[\int\limits_{-\infty }^{+\infty }{\exp \left( -\frac{{{\left( x+{{y}_{3}} \right)}^{2}}}{{{a}^{2}}} \right)\operatorname{sinc}\left( \frac{{{\left( x-{{y}_{0}} \right)}^{2}}}{{{b}^{2}}} \right)\operatorname{sinc}\left( \frac{{{\left( x+{{y}_{0}} \right)}^{2}}}{{{b}^{2}}} \right)dx}\]$$

У самого соображений нет, пытался скормить онлайн интегратору, но он выдал первообразную, которая обращается в нули на бесконечности. Вызывает сомнение вообще полученный результат.

Буду очень благодарен за подсказки.

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 03:10 
B@R5uk в сообщении #1083722 писал(а):
он выдал первообразную, которая обращается в нули на бесконечности.

Как-то навскидку это не очевидно. Вы проверяли?

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 03:15 
Аватара пользователя
Там на вкладке Traditional Form не правильная картинка. Если расшифровать код с вкладки Input Form, то получится такая штука:

$$\[\begin{matrix}
  \int{\exp \left( -\frac{{{\left( x+{{y}_{3}} \right)}^{2}}}{{{a}^{2}}} \right)\operatorname{sinc}\left( \frac{{{\left( x-{{y}_{0}} \right)}^{2}}}{{{b}^{2}}} \right)\operatorname{sinc}\left( \frac{{{\left( x+{{y}_{0}} \right)}^{2}}}{{{b}^{2}}} \right)dx}=\frac{a}{8{{b}^{4}}}\exp \left( -\frac{{{\left( x+{{y}_{3}} \right)}^{2}}}{{{a}^{2}}} \right)\times  \\ 
  \times {{\operatorname{sinc}}^{2}}\left( 2a\left( {{a}^{2}}\left( -3x+5{{y}_{3}} \right)-2\left( x-{{y}_{3}} \right)\left( {{x}^{2}}-2y_{0}^{2}+3y_{3}^{2} \right) \right)+\sqrt{\pi }\operatorname{erf}\left( \frac{x+{{y}_{3}}}{a} \right)\times  \right. \\ 
  \times \left. \left( 3{{a}^{4}}-4{{a}^{2}}\left( y_{0}^{2}-3y_{3}^{2} \right)+4{{\left( y_{0}^{2}-y_{3}^{2} \right)}^{2}} \right)\exp \left( \frac{{{\left( x+{{y}_{3}} \right)}^{2}}}{{{a}^{2}}} \right) \right) \\ 
\end{matrix}\]$$

Синк ограничен, а экспонента перед ним стремится к нулю, при стремлении икса к любой из бесконечности.

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 03:23 
А зачем Вы экспоненту из аргумента вытащили? Нормальная там картинка, можете с Output Form свериться.

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 03:26 
Аватара пользователя
Otta в сообщении #1083731 писал(а):
можете с Output Form свериться
Там тоже в знаменателе
$$\[8{{b}^{4}}\exp \left( \frac{{{\left( x+{{y}_{3}} \right)}^{2}}}{{{a}^{2}}} \right)\]$$

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 03:35 
Слушайте, не, он вообще какую-то дурь выдает.
Он не знает этой функции. И воспринимает запись sinc как постоянную величину $\sin c$.

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 03:39 
Аватара пользователя
Otta в сообщении #1083734 писал(а):
воспринимает запись sinc как постоянную
Так вот в чём дело было! Жаль, что не знает такой функции. А если расписать через отношение синуса к аргументу, то он отказывается интегрировать.

Может как-нибудь можно этот интеграл руками взять?

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 04:31 
B@R5uk в сообщении #1083736 писал(а):
А если расписать через отношение синуса к аргументу, то он отказывается интегрировать.

Ничего удивительного, когда каждый множитель в отдельности не имеет первообразной в элементарных функциях. Можно попробовать считать
$$\int\limits_{-\infty }^{+\infty }\exp\left(- (x+y_3)^2 \right)\frac{\sin b (x-y_0)^2}{(x-y_0)^2}\frac{\sin c(x+y_0)^2}{(x+y_0)^2} \, dx$$
дифференцируя по параметрам. По крайней мере с одним синусом он должен считаться, ну а с двумя может выйдет, если тригонометрию припахать. Не первообразные считать, конечно, а сам несобственный интеграл.

Я не пробовала, не факт, что выйдет. Шансов было бы больше если бы $y_i$ были нулевыми.

(Оффтоп)

И зачем Вы ставите столько фигурных скобок в формулах? В итоге проще набрать заново, чем чуть-чуть поправить.

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 10:11 
Аватара пользователя
Otta в сообщении #1083745 писал(а):
когда каждый множитель в отдельности не имеет первообразной в элементарных функциях
Ну, это ещё ничего не значит. У меня есть несколько вдохновляющих результатов:$$\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\sin \left( {{b}^{2}}{{x}^{2}} \right)dx}=\frac{\sqrt{\pi }}{\sqrt[4]{{{a}^{4}}+{{b}^{4}}}}\sin \left( \frac{1}{2} \arctg\frac{{{b}^{2}}}{{{a}^{2}}} \right)=$$$$=\sqrt{\frac{\pi }{2}}\sqrt{\frac{-{{a}^{2}}+\sqrt{{{a}^{4}}+{{b}^{4}}}}{{{a}^{4}}+{{b}^{4}}}}=\frac{\sqrt{\pi }}{2}\left( \frac{1}{\sqrt{-{{a}^{2}}+i{{b}^{2}}}}+\frac{1}{\sqrt{-{{a}^{2}}-i{{b}^{2}}}} \right)$$
$$\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\cos \left( {{b}^{2}}{{x}^{2}} \right)dx}=\frac{\sqrt{\pi }}{\sqrt[4]{{{a}^{4}}+{{b}^{4}}}}\cos \left( \frac{1}{2} \arctg\frac{{{b}^{2}}}{{{a}^{2}}} \right)=$$$$=\sqrt{\frac{\pi }{2}}\sqrt{\frac{{{a}^{2}}+\sqrt{{{a}^{4}}+{{b}^{4}}}}{{{a}^{4}}+{{b}^{4}}}}=\frac{\sqrt{\pi }}{2}\left( \frac{1}{\sqrt{{{a}^{2}}-i{{b}^{2}}}}+\frac{1}{\sqrt{{{a}^{2}}+i{{b}^{2}}}} \right)$$
$$\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\operatorname{sinc}\left( {{b}^{2}}{{x}^{2}} \right)dx}=\sqrt{\frac{2\pi }{{{a}^{2}}+\sqrt{{{a}^{4}}+{{b}^{4}}}}}=\frac{2\sqrt{\pi }}{\sqrt{{{a}^{2}}-i{{b}^{2}}}+\sqrt{{{a}^{2}}+i{{b}^{2}}}}$$Впрочем, элементарные функции не нужны, главное, чтобы полученные функции считались в матпакетах.

Otta в сообщении #1083745 писал(а):
зачем Вы ставите столько фигурных скобок в формулах?
Дело в том, что я вывожу формулы в Ворде, а не на бумаге. Это быстрее и удобнее, меньше шанс ошибиться при переписывании (копи-паст не ошибается, только копирует уже сделанные ошибки). Команда автоматической конвертации формул в $\TeX$ ставит фигурные скобки везде, где только можно.

Otta в сообщении #1083745 писал(а):
Можно попробовать считать дифференцируя по параметрам.
Я тоже про это думал. Сейчас пытаюсь сообразить случай, когда под интегралом только один синк, смещённый относительно гауссианы. Там есть проблема связанная как раз с этими ненулевыми игреками. Тот же онлайн-интегратор выдаёт для такого случая такую первообразную:$$\int{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\sin \left( {{b}^{2}}{{\left( x-y \right)}^{2}} \right)dx}=\frac{i\sqrt{\pi }\exp \left( -{{a}^{2}}{{y}^{2}} \right)}{4\sqrt{{{a}^{4}}+{{b}^{4}}}}\times $$$$\times \left( \sqrt{{{a}^{2}}-i{{b}^{2}}}\exp \left( \frac{{{a}^{4}}{{y}^{2}}}{{{a}^{2}}+i{{b}^{2}}} \right)\operatorname{erf}\left( \frac{\left( {{a}^{2}}+i{{b}^{2}} \right)x-i{{b}^{2}}y}{\sqrt{{{a}^{2}}+i{{b}^{2}}}} \right)- \right.$$$$\left. -\sqrt{{{a}^{2}}+i{{b}^{2}}}\exp \left( \frac{{{a}^{4}}{{y}^{2}}}{{{a}^{2}}-i{{b}^{2}}} \right)\operatorname{erf}\left( \frac{\left( {{a}^{2}}-i{{b}^{2}} \right)x+i{{b}^{2}}y}{\sqrt{{{a}^{2}}-i{{b}^{2}}}} \right) \right)$$Но даже в случае, когда $y=0$, и выражение для первообразной значительно упрощается, я не представляю, как подставить бесконечные пределы интегрирования, ведь под функцией ошибок аргумент имеет комплексный множитель. Другой интегратор (из которого я путём эксперимента выудил формулы выше) каким-то образом эту проблему с комплексным коэффициентом решает. Но это скрыто в недрах системы. К сожалению я не знаком с функцией ошибок на столько, чтобы решить эту задачу самостоятельно. Подскажите, пожалуйста, что тут можно сделать?

 
 
 
 Re: Определённый интеграл с синками
Сообщение20.12.2015, 13:50 
Аватара пользователя
Покопавшись в справочнике Рыжика и Градштейна обнаружил вот такой интеграл с экспонентами и синусами от квадратичных полиномов:$$\begin{matrix}
  \int\limits_{-\infty }^{+\infty }{\exp \left( -\left( a{{x}^{2}}+2bx+c \right) \right)\cos \left( p{{x}^{2}}+2qx+r \right)dx}= \\ 
  =\frac{\sqrt{\pi }}{\sqrt[4]{{{a}^{2}}+{{p}^{2}}}}\exp \left( \frac{a\left( {{b}^{2}}-ac \right)-\left( a{{q}^{2}}-2bpq+c{{p}^{2}} \right)}{{{a}^{2}}+{{p}^{2}}} \right)\times  \\ 
  \times \cos \left( \frac{1}{2}\arctg\frac{p}{a}-\frac{p\left( {{q}^{2}}-p r \right)-\left( {{b}^{2}}p-2abq+{{a}^{2}}r \right)}{{{a}^{2}}+{{p}^{2}}} \right) \\ 
\end{matrix}$$После небольших махинаций с параметрами вот этот интеграл$$\[I\left( y,a,b \right)=\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\operatorname{sinc}\left( {{b}^{2}}{{\left( x-y \right)}^{2}} \right)dx}\]$$свёлся к такому:
$$\[I\left( y,a,b \right)=\frac{2\sqrt{\pi }}{{{b}^{2}}}\int\limits_{0}^{b}{\frac{\beta }{\sqrt[4]{{{a}^{4}}+{{\beta }^{4}}}}\exp \left( -\frac{{{a}^{2}}{{\beta }^{4}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)\cos \left( \frac{1}{2}\arctg\frac{{{\beta }^{2}}}{{{a}^{2}}}+\frac{{{a}^{4}}{{\beta }^{2}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)d\beta }\]$$Даже и не знаю, как к нему подступиться. Косинус можно раскрыть, синус и косинус половинного арктангенса превратятся в радикалы, при этом корень четвёртой степени в знаменателе станет корнем второй степени. Типа вот так:$$\[\begin{matrix}
  \frac{2}{\sqrt[4]{{{a}^{4}}+{{b}^{4}}}}\cos \left( \frac{1}{2}\arctg\frac{{{b}^{2}}}{{{a}^{2}}}+\frac{{{a}^{4}}{{b}^{2}}{{y}^{2}}}{{{a}^{4}}+{{b}^{4}}} \right)= \\ 
  =\left( \frac{1}{\sqrt{{{a}^{2}}+i{{b}^{2}}}}+\frac{1}{\sqrt{{{a}^{2}}-i{{b}^{2}}}} \right)\cos \left( \frac{{{a}^{4}}{{b}^{2}}{{y}^{2}}}{{{a}^{4}}+{{b}^{4}}} \right)-\left( \frac{1}{\sqrt{-{{a}^{2}}+i{{b}^{2}}}}+\frac{1}{\sqrt{-{{a}^{2}}-i{{b}^{2}}}} \right)\sin \left( \frac{{{a}^{4}}{{b}^{2}}{{y}^{2}}}{{{a}^{4}}+{{b}^{4}}} \right) \\ 
\end{matrix}\]$$
Но вот что делать потом? Что можно сделать, например, с таким интегралом:$$\[\int\limits_{0}^{b}{\frac{\beta }{\sqrt{{{a}^{2}}+i{{\beta }^{2}}}}\exp \left( -\frac{{{a}^{2}}{{\beta }^{4}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)\cos \left( \frac{{{a}^{4}}{{\beta }^{2}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)d\beta }\]$$

 
 
 
 Re: Определённый интеграл с синками
Сообщение26.12.2015, 18:49 
Аватара пользователя
А если я путём численного эксперимента угадаю аналитический вид функции, задаваемой интегралом с параметрами, например, такой:$$\[F\left( {{x}_{1}},{{x}_{2}},b \right)=\int\limits_{-\infty }^{+\infty }{\exp \left( -{{x}^{2}} \right)\operatorname{sinc}\left( {{b}^{2}}{{\left( x-{{x}_{1}} \right)}^{2}} \right)\operatorname{sinc}\left( {{b}^{2}}{{\left( x-{{x}_{2}} \right)}^{2}} \right)dx}\]$$То как мне аналитически проверить, что угаданная функция действительно является интегралом?

Вот, например, если бы я интеграл с переменным верхним пределом, то мне было бы достаточно дифференцировать и вычесть. Здесь же такое не прокатывает. Существует ли метод сделать похожую проверку в моём случае?

 
 
 
 Re: Определённый интеграл с синками
Сообщение26.12.2015, 19:21 
У вас целых три свободные переменные. :-) Всяческие подстановки, дифференцирования по ним и тому подобное могут дать контрпример (а могут и не дать, конечно).

 
 
 
 Re: Определённый интеграл с синками
Сообщение26.12.2015, 20:20 
Аватара пользователя
Контрпример меня не очень интересует. На последнем этапе сравнения хотелось бы получить что-то вроде метода неопределённых коэффициентов.

 
 
 
 Re: Определённый интеграл с синками
Сообщение27.12.2015, 13:56 
Аватара пользователя
Проверьте, пожалуйста, мои выкладки. Не могу понять где ошибка.

У меня есть целевой интеграл, который я хочу посчитать:$$\[I\left( y,a,b \right)=\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\operatorname{sinc}\left( {{b}^{2}}{{\left( x-y \right)}^{2}} \right)dx}\]$$Я преобразую его следующим образом:$$\[{{b}^{2}}I\left( y,a,b \right)=\int\limits_{-\infty }^{+\infty }{\frac{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\sin \left( {{b}^{2}}{{\left( x-y \right)}^{2}} \right)}{{{\left( x-y \right)}^{2}}}dx}\]$$Затем дифференцирую по параметру $b$:$$\[\frac{\partial }{\partial b}\left( {{b}^{2}}I\left( y,a,b \right) \right)=2b\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\cos \left( {{b}^{2}}{{\left( x-y \right)}^{2}} \right)dx}\]$$Интеграл справа у меня посчитан и проверен (это слегка модернизированная формула 3.923 из Рыжика и Градштейна, страница 499), его значение:$$\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\cos \left( {{b}^{2}}{{\left( x-y \right)}^{2}} \right)dx}=\frac{\sqrt{\pi }}{\sqrt[4]{{{a}^{4}}+{{b}^{4}}}}\exp \left( -\frac{{{a}^{2}}{{b}^{4}}{{y}^{2}}}{{{a}^{4}}+{{b}^{4}}} \right)\cos \left( \frac{1}{2} \arctg\frac{{{b}^{2}}}{{{a}^{2}}}+\frac{{{a}^{4}}{{b}^{2}}{{y}^{2}}}{{{a}^{4}}+{{b}^{4}}} \right)$$С учётом этого у меня получается следующее:$$\[\frac{\partial }{\partial b}\left( {{b}^{2}}I\left( y,a,b \right) \right)=\frac{2b\sqrt{\pi }}{\sqrt[4]{{{a}^{4}}+{{b}^{4}}}}\exp \left( -\frac{{{a}^{2}}{{b}^{4}}{{y}^{2}}}{{{a}^{4}}+{{b}^{4}}} \right)\cos \left( \frac{1}{2}\arctg\frac{{{b}^{2}}}{{{a}^{2}}}+\frac{{{a}^{4}}{{b}^{2}}{{y}^{2}}}{{{a}^{4}}+{{b}^{4}}} \right)\]$$Теперь я интегрирую правую и левую часть по переменной $b$ от $0$ до $b$. Слева интеграл берётся сразу, нулевая подстановка даёт нуль, поэтому я получаю:$$\[I\left( y,a,b \right)=\frac{2\sqrt{\pi }}{{{b}^{2}}}\int\limits_{0}^{b}{\frac{\beta }{\sqrt[4]{{{a}^{4}}+{{\beta }^{4}}}}\exp \left( -\frac{{{a}^{2}}{{\beta }^{4}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)\cos \left( \frac{1}{2}\arctg\frac{{{\beta }^{2}}}{{{a}^{2}}}+\frac{{{a}^{4}}{{\beta }^{2}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)d\beta }\]$$Последний определённый интеграл с конечными пределами я тоже посчитал (хотя выражение получилось с комплексными числами, даёт оно чисто действительный результат):$$\[\begin{matrix}
  I\left( y,a,b \right)=\frac{i\sqrt{\pi }}{{{b}^{2}}}\left( \sqrt{{{a}^{2}}-i{{b}^{2}}}\exp \left( \frac{i{{a}^{2}}{{b}^{2}}{{y}^{2}}}{{{a}^{2}}-i{{b}^{2}}} \right)-\sqrt{{{a}^{2}}+i{{b}^{2}}}\exp \left( -\frac{i{{a}^{2}}{{b}^{2}}{{y}^{2}}}{{{a}^{2}}+i{{b}^{2}}} \right) \right)+ \\ 
  +\frac{\pi {{a}^{2}}y}{{{b}^{2}}}\exp \left( -{{a}^{2}}{{y}^{2}} \right)\left( \operatorname{erf}\left( \frac{i{{a}^{2}}y}{\sqrt{{{a}^{2}}-i{{b}^{2}}}} \right)-\operatorname{erf}\left( \frac{i{{a}^{2}}y}{\sqrt{{{a}^{2}}+i{{b}^{2}}}} \right) \right) \\ 
\end{matrix}\]$$

Теперь я хочу эти вычисления проверить. Для этого я задаюсь какими-нибудь случайными значениями параметров $a$, $b$ и $y$ и рассчитываю следующие три величины:$$\[{{V}_{1}}=\int\limits_{-\infty }^{+\infty }{\exp \left( -{{a}^{2}}{{x}^{2}} \right)\operatorname{sinc}\left( {{b}^{2}}{{\left( x-y \right)}^{2}} \right)dx}\]$$$$\[{{V}_{2}}=\frac{2\sqrt{\pi }}{{{b}^{2}}}\int\limits_{0}^{b}{\frac{\beta }{\sqrt[4]{{{a}^{4}}+{{\beta }^{4}}}}\exp \left( -\frac{{{a}^{2}}{{\beta }^{4}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)\cos \left( \frac{1}{2}\arctg\frac{{{\beta }^{2}}}{{{a}^{2}}}+\frac{{{a}^{4}}{{\beta }^{2}}{{y}^{2}}}{{{a}^{4}}+{{\beta }^{4}}} \right)d\beta }\]$$$$\[\begin{matrix}
  {{V}_{3}}=\frac{i\sqrt{\pi }}{{{b}^{2}}}\left( \sqrt{{{a}^{2}}-i{{b}^{2}}}\exp \left( \frac{i{{a}^{2}}{{b}^{2}}{{y}^{2}}}{{{a}^{2}}-i{{b}^{2}}} \right)-\sqrt{{{a}^{2}}+i{{b}^{2}}}\exp \left( -\frac{i{{a}^{2}}{{b}^{2}}{{y}^{2}}}{{{a}^{2}}+i{{b}^{2}}} \right) \right)+ \\ 
  +\frac{\pi {{a}^{2}}y}{{{b}^{2}}}\exp \left( -{{a}^{2}}{{y}^{2}} \right)\left( \operatorname{erf}\left( \frac{i{{a}^{2}}y}{\sqrt{{{a}^{2}}-i{{b}^{2}}}} \right)-\operatorname{erf}\left( \frac{i{{a}^{2}}y}{\sqrt{{{a}^{2}}+i{{b}^{2}}}} \right) \right) \\ 
\end{matrix}\]$$Первые две рассчитываю путём численного интегрирования, а третью — просто по формуле (правда пришлось найти реализацию для MATLAB функции ошибок комплексного аргумента). При этом я обнаруживаю, что вторая и третья величины совпадают вплоть до самых последних значащих цифр в формате числа (точность интегрирования поставил $10^{-12}$), но совершенно не совпадают с первой величиной. Отсюда я делаю вывод, что где-то допустил ошибку. Однако, как бы я не старался её найти ничего подозрительного не вижу. Нарекания мог бы вызывать промежуточный интеграл от экспоненты и косинуса с бесконечными пределами, однако он удачно прошёл точно такую же численную валидацию, какую я применяю к первому и последнему выражениям.

Очень надеюсь на вашу помощь.

 
 
 
 Re: Определённый интеграл с синками
Сообщение27.12.2015, 14:14 
Аватара пользователя
А не кошерно ли будет попросить Коши? В смысле припасть к теореме о вычетах?

 
 
 [ Сообщений: 18 ]  На страницу 1, 2  След.


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group