2014 dxdy logo

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

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




 
 Преобразовать ряд к конечной сумме
Сообщение12.12.2017, 20:03 
При вычислении интеграла $I=\int\limits_0^1\frac{a\cdot t^n}{b+a\cdot t}dt$ для $\left|\frac{a}{b}\right|\ge1$ получил, что $I=\left(\frac{-b}{a}\right)^n\left[\ln\left(\frac{a+b}{b}\right)+\sum\limits_{k=1}^n\left(\frac{1}{k}\left(\frac{-a}{b}\right)^k\right)\right]$. Здесь $n=0,1,2,3...$. В случае, когда $\left|\frac{a}{b}\right|<1$ (особенно для больших $n$) $I=-\sum\limits_{k=0}^{\infty}\left[\frac{1}{k+n+1}\left(\frac{-a}{b}\right)^{k+1}\right]$ (это равенство можно задать и через гипергеометрическую функцию). В последнем выражении не устраивает бесконечная сумма. Подскажите, пожалуйста, как вычислить $I$ для $\left|\frac{a}{b}\right|<1$ не прибегая к бесконечным суммам или гипергеометрической функции.

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение12.12.2017, 20:32 
van341 в сообщении #1274429 писал(а):
$I=\left(\frac{-b}{a}\right)^n\left[\ln\left(\frac{a+b}{b}\right)+\sum\limits_{k=1}^n\left(\frac{1}{k}\left(\frac{-a}{b}\right)^k\right)\right]$...
$I=-\sum\limits_{k=0}^{\infty}\left[\frac{1}{k+n+1}\left(\frac{-a}{b}\right)^{k+1}\right]$
Подскажите, пожалуйста, как вычислить $I$ для $\left|\frac{a}{b}\right|<1$ не прибегая к бесконечным суммам
Слишком сильное требование.
Гармоническое число - оно и есть гармоническое число, асимптотика его известна. Если гармонический ряд начать с другого места - будет тоже асимптотика, но другая. Но даже асимптотический ряд - это же ряд, его нельзя свернуть в что-то более простое.
А вот как Вы вычисляете логарифм без бесконечных рядов - я не понимаю.

(Оффтоп)

$I(x)=-\sum\limits_{k=1}^{\infty}\frac{x^k}{n+k}$

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение12.12.2017, 20:48 
van341
А причём тут величина $\[\frac{a}{b}\]$? (за исключением случаев, когда интеграл сходится)?
Разложение на простые дроби всегда одно
$$\[\frac{{a{t^n}}}{{b + at}} = {( - 1)^n}\frac{{{b^n}}}{{{a^{n - 1}}}}\frac{1}{{b + at}} + \sum\limits_{k = 0}^{n - 1} {{{( - 1)}^{n - k - 1}}{{(\frac{b}{a})}^{n - k - 1}}{t^k}} \]$$
Откуда
$$\[\int\limits_0^1 {\frac{{a{t^n}}}{{b + at}}dt}  = {( - 1)^n}\frac{{{b^n}}}{{{a^n}}}\ln (1 + \frac{a}{b}) + \sum\limits_{k = 0}^{n - 1} {\frac{{{{( - 1)}^{n - k - 1}}}}{{k + 1}}{{(\frac{b}{a})}^{n - k - 1}}} \]$$
P.S.Присоединяюсь к предыдущему ответу. Если хотите что-то считать, проще это делать численно.

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение12.12.2017, 20:55 
Ms-dos4
К сожалению для $\left|\frac{a}{b}\right|<1$ и больших $n$ это равенство не выполняется.

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение12.12.2017, 21:19 
van341
Попробую угадать. Вы считаете сумму в какой нибудь системе компьютерной алгебры? Тогда просто смотрите на точность. Равенство верное, просто величина суммы и первого слагаемого очень велики и когда вы считаете сумму (по факту разность), то реальный результат задавлен ошибкой округления. Вот например, в Mathematicа, вычисляем при $\[a = 1\]$, $\[b = 7\]$ и $\[n = 20\]$. Интеграл равен
Код:
NIntegrate[(a*t^n)/(b + a*t), {t, 0, 1}]

$\[\int\limits_0^1 {\frac{{a{t^n}}}{{b + at}}dt}  \approx 0.00598657477282...\]$
Считаем сумму
Код:
N[(-1)^n*(b/a)^n*Log[1 + a/b] +
  Sum[(-1)^(n - k - 1)/(k + 1)*(b/a)^(n - k - 1), {k, 0, n - 1}]]

$$\[{( - 1)^n}\frac{{{b^n}}}{{{a^n}}}\ln (1 + \frac{a}{b}) + \sum\limits_{k = 0}^{n - 1} {\frac{{{{( - 1)}^{n - k - 1}}}}{{k + 1}}{{(\frac{b}{a})}^{n - k - 1}}}  =  - 4\]$$
Оой, как так то. Ан-нет, всё нормально, если указать системе считать с точностью, например 1000 знаков, то
Код:
N[(-1)^n*(b/a)^n*Log[1 + a/b] +
  Sum[(-1)^(n - k - 1)/(k + 1)*(b/a)^(n - k - 1), {k, 0,
    n - 1}], 1000]

$$\[{( - 1)^n}\frac{{{b^n}}}{{{a^n}}}\ln (1 + \frac{a}{b}) + \sum\limits_{k = 0}^{n - 1} {\frac{{{{( - 1)}^{n - k - 1}}}}{{k + 1}}{{(\frac{b}{a})}^{n - k - 1}}}  = 0.00598657477282...\]$$
И чем больше зададите, тем лучше точность (обратите внимание, она НЕ равна заданным 1000 знакам).

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение12.12.2017, 21:25 
Ms-dos4
Ой ой)) Все верно. Начал мудрить там, где не надо. Спасибо))

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение12.12.2017, 21:26 
van341
Только вот хотел бы повторить, это как раз указывает на то, что НЕ НУЖНО вычислять рядами. Это по меньшей мере глупо. Намного эффективнее просто считать интеграл численно.

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение13.12.2017, 21:58 
Аватара пользователя
Мучает меня смутное сомнение, что если взять сумму геометрической прогрессии и проинтегрировать...

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение14.12.2017, 12:57 
Для больших $n$ может оказаться выгоднее вычислять с помощью ряда. $I=\frac a{(n+1)(a+b)}+o(\frac 1{n})$ (явно выписан первый член ряда).

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение14.12.2017, 15:53 
Аватара пользователя

(Оффтоп)

Или я чего-то не понимаю, или одно из двух...

$I=-\sum\limits_{k=0}^{\infty}\left[\frac{1}{k+n+1}\left(\frac{-a}{b}\right)^{k+1}\right]=-\sum\limits_{k=0}^{\infty}\left[\frac{1}{k+n+1}r^{k+1}\right]$
где
$r=\frac{-a}{b}$
I=-r^{-n}\sum\limits_{k=0}^{\infty}\left[\frac{1}{k+n+1}r^{k+n+1}\right]=-r^{-n}J
где
$J=\sum\limits_{k=0}^{\infty}\left[\frac{1}{k+n+1}r^{k+n+1}\right]$
$\frac {dJ}{dr}=\sum\limits_{k=0}^{\infty}\left[r^{k+n+1}\right]=\frac{r^{n+1}}{1-r}$

(Оффтоп)

intégre qui pourra

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение14.12.2017, 19:54 
Евгений Машеров
Вы всё правильно понимаете, но человек же хотел конечную сумму (хотя учитывая что в этом случае там присутствует логарифм, это смысла имеет мало).
P.S. $\[\frac{{dJ}}{{dr}} = \frac{{{r^{n + 1}}}}{{1 - r}}\]$

 
 
 
 Re: Преобразовать ряд к конечной сумме
Сообщение14.12.2017, 21:36 
Аватара пользователя
Спасибо. Поправил.

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


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