Вычисление тройного интеграла с использованием замены





Здравствуйте. Никак не могу понять как правильно в данном примере сделать замену. Дан интеграл, где область G ограничена 5-ю поверхностями: x^2+y^2<=a^2 y^2+z^2<=a^2 x>=0 z>=0 y>=0 Пробовал переходить к сферическим и цилиндрическим системам координат, все собственно показано на фотографии. Но интегралы там не очень радужные выходят. Что WolframAlpha не справляется с вычислением, что wxMaxima. Что собственно не так? UPD. Попытка с цилиндрической системой координат. http://savepic.ru/12506402.jpg

  • Автор сообщения: Петр

Всего: 9 комментариев.
Добавить комментарий

Александр Емелин    03.12.2016 в 07:17
Здравствуйте, Пётр. Почему интегралы не радужные? Они таки довольно простые - решаем через цилиндрические координаты, и во внутреннем интеграле проводим замену:
t = a^2+r^2+z^2
Тогда dt=2zdz (поскольку а - константа, и r тоже считается константой), откуда: zdz=dt/2
Находим новые пределы интегрирования и берём внутренний интеграл. Более подробно метод замены в определённом интеграле можно посмотреть здесь:
http://mathprofi.ru/opredelennye_integraly_primery_reshenij.html

Петр    03.12.2016 в 12:20
Я, собственно, так и пытался. Но именно после вычисления самого внутреннего интеграла там и выходит не очень радужный)
Добавил ссылку на новую фотографию с тем, что выходит.

Александр Емелин    03.12.2016 в 12:41
Придётся ломать интегралы по частям,
обозначая u = r^2
Решение громоздкое, но вполне реальное.

Но вообще такое чувство, что где-то в условии опечатка, возможно, здесь будет целесообразно проконсультироваться с преподавателем.

Однако, интересно, что этот интеграл всё же берётся. Существует ещё одна гипотеза решения, а именно переход к нестандартной (криволинейной) системе координат с вычислением Якобиана перехода. Но по первой оглядке что-то ничего в голову не приходит(

Петр    03.12.2016 в 13:18
Условие все же выходит верное, скачал Wolfram Mathematica, выводит как раз то, что указано в ответах. (ln(3)-1)/16. Но вычисляет он это где-то секунд 10.

Александр Емелин    03.12.2016 в 13:31
Придётся возиться с интегралами, и решение тут удобно разделить на две части - сначала один интеграл, затем другой, и, наконец, разность.

Петр    03.12.2016 в 16:32
Вот что в итоге выходит. Второй решается легко, а вот первый...) Поседеешь, пока решать будешь.
http://savepic.ru/12503379.jpg

Александр Емелин    04.12.2016 в 07:39
откуда у вас синус в 4-й степени в знаменателе появился?
И куда делся параметр а? Интегрируется же от 0 до а (я так понял, вы его единицей заменили)

Петр    04.12.2016 в 20:55
http://savepic.ru/12474418.jpg 1 часть интеграла

http://savepic.ru/12461106.jpg
http://savepic.ru/12474421.jpg 2 часть
а сокращается, если б он на данном этапе не сократился, он бы был в ответе

Александр Емелин    05.12.2016 в 07:02
Как-то странно вы интегрируете по частям, тут лучше использовать стандартный алгоритм:
обозначаем dv = r / (...)^2, находим v = - 1/(2*(...)) и используем обычную формулу, только с пределами интегрирования:
http://mathprofi.ru/opredelennye_integraly_primery_reshenij.html
(последний параграф)

После чего нужно взять два интеграла по "фи" - один с тригонометрическими функциями, другой с ними же и логарифмом.
Последний интеграл не подарочный (замена + снова интегрирование частям), но вполне решабельный.