2014 dxdy logo

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

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




 
 Посчитать интеграл
Сообщение02.10.2015, 15:27 
Здравствуйте!

Есть интеграл:
$   \int_{0}^{x} z^b \cdot e^{a z} \cdot dz  $
Надо его найти. Параметры действительные $a>0, b>0, 0<x<1$, то есть могут быть и дробными.

Сделал чтобы он был похож на неполную гамма-функцию:

$  \frac{1}{(-a)^{1+b}}  \int_{0}^{x} (-az)^b \cdot e^{a z} \cdot d(-az)     $

После замены переменной $y=-az$ происходит следующее:

$  \frac{1}{(-a)^{1+b}}  \int_{0}^{-ax} y^b \cdot e^{-y} \cdot dy     $

То есть получается отрицательное число в верхнем пределе интегрирования. Мат. пакеты отказываются такое считать (при попытке посчитать неполную гамма-функцию с такими параметрами: $\Gamma(b+1,-ax)$), хотя метод трапеций для $z^b \cdot e^{a z}$ работает нормально. Что делать? Как точно вычислить интеграл?

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 16:11 
Аватара пользователя
Zazaqa в сообщении #1058439 писал(а):
Сделал чтобы он был похож на неполную гамма-функцию:

$  \frac{1}{(-a)^{1+b}}  \int_{0}^{x} (-az)^b \cdot e^{a z} \cdot d(-az)     $

Вас не смущает, что под "произвольную" степень попало отрицательное число:$(-az)^b$ ?

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 16:22 
Zazaqa в сообщении #1058439 писал(а):
Мат. пакеты отказываются такое считать

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

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 16:39 
Pphantom в сообщении #1058457 писал(а):
Не совсем профильное для раздела "решение"... но все же Maxima выполняет и приведение к неполным гамма-функциям, и вычисление значений, что называется, "из коробки". Поэтому, если надо просто посчитать, то см.выше, а если надо именно так, то распишите, как Вы "сделал чтобы", похоже, что там ошибка.

Нет ошибки (вроде), так как $ \frac{1}{(-a)^1} $ идет под дифференциал $d$, а $ \frac{1}{(-a)^b} $ домножается к $ (z)^b $.

Brukvalub в сообщении #1058452 писал(а):
Вас не смущает, что под "произвольную" степень попало отрицательное число:$(-az)^b$ ?

А что с этим можно сделать?

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 16:40 
Аватара пользователя
Zazaqa в сообщении #1058464 писал(а):
Brukvalub в сообщении #1058452

писал(а):
Вас не смущает, что под "произвольную" степень попало отрицательное число:$(-az)^b$ ?
А что с этим можно сделать?

Просто можно нужно так не делать!

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 17:07 
Brukvalub в сообщении #1058466 писал(а):
Просто можно нужно так не делать!

Вы намекаете на то, что я неправильно привожу? Или что приводить надо к чему-то другому?
Если к неполной гамме приводить, то другого варианта не вижу.

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 17:10 
Аватара пользователя
Zazaqa в сообщении #1058475 писал(а):
Если к неполной гамме приводить, то другого варианта не вижу.

Это беспощадный аргумент! Если нельзя, но очень хочется, то можно! :D

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 21:15 
Аватара пользователя
Раскладывайте с другого конца.

 
 
 
 Re: Посчитать интеграл
Сообщение02.10.2015, 22:15 
ИСН в сообщении #1058591 писал(а):
Раскладывайте с другого конца.

Что это означает?

 
 
 
 Re: Посчитать интеграл
Сообщение03.10.2015, 13:17 
Аватара пользователя
$\int\limits_0^xze^z\;dz=-\int\limits_x^0(x-t)e^{x-t}\;dt=e^x\int\limits_0^x(x-t)e^{-t}\;dt$, а этот Вы предположительно знаете, как выразить через неполную гамму.

 
 
 
 Re: Посчитать интеграл
Сообщение03.10.2015, 20:33 
ИСН в сообщении #1058780 писал(а):
$\int\limits_0^xze^z\;dz=-\int\limits_x^0(x-t)e^{x-t}\;dt=e^x\int\limits_0^x(x-t)e^{-t}\;dt$, а этот Вы предположительно знаете, как выразить через неполную гамму.

Я понял что Вы хотите сделать, что так можно менять знаки, оставляя пределы интегрирования теми же. Но степень $b$ мешает и от нее не избавиться (так как дробная).

 
 
 
 Re: Посчитать интеграл
Сообщение04.10.2015, 13:28 
Аватара пользователя
Ах, дробная. Хм. Тогда таки плохо.

 
 
 
 Re: Посчитать интеграл
Сообщение04.10.2015, 15:22 
ИСН в сообщении #1059046 писал(а):
Ах, дробная. Хм. Тогда таки плохо.

В первом посте у меня получалась вот такая штука:

$ \frac{1}{(-a)^{1+b}} \int_{0}^{-ax} y^b \cdot e^{-y} \cdot dy = \frac{1}{(-a)^{1+b}} \Gamma(1+b, -ax)$

Во-первых, $ \Gamma(1+b, -ax)$ в мат. пакетах возвращает комплексное число $\gamma =  \Gamma(1+b, -ax)$.
Во-вторых, интегрируя численно обнаружил, что $ \int_{0}^{x} z^b \cdot e^{a z} \cdot dz = \frac{1}{a^{1+b}} |\gamma|$, где $| \cdot |$ есть модуль числа (длина вектора на комплексной плоскости).

Осталось найти тот самый промежуточный шаг решения. Я пока не знаю, как это сделать.

 
 
 [ Сообщений: 13 ] 


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