2014 dxdy logo

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

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


Правила форума


Посмотреть правила форума



Начать новую тему Ответить на тему
 
 Определенный интеграл
Сообщение30.04.2022, 20:10 
Аватара пользователя


08/10/09
835
Херсон
Необходимо вычислить определенный интеграл: $$\int_{0}^{a}\,\frac{dx}{\sqrt{1-\frac{1}{\left(1-\frac{1}{a}+\frac{1}{x}\right)^2}}}$$ при произвольном $a>0$. Maple через два раза на третий выдает горбатое выражение в элементарных функциях. Однако оно содержит корень $\sqrt{2a-1}$, который перестает быть действительным числом при $a<1/2$. Wolfram вообще выдает выражение с мнимой единицей. Однако, результат очевидно должен быть положительным действительным числом. Прошу помочь разобраться с этой несуразицей

 Профиль  
                  
 
 Re: Определенный интеграл
Сообщение30.04.2022, 20:13 
Заслуженный участник


20/12/10
8858
Сделайте замену $1-1/a+1/x=y$. И обозначьте $1-1/a=b$, так что $b<1$.

 Профиль  
                  
 
 Re: Определенный интеграл
Сообщение01.05.2022, 04:29 
Аватара пользователя


07/01/16
1426
Аязьма
Трудно этот способ рекомендовать, но в качестве борьбы с вольфрамом за результат он по-моему забавен.
Введя $k=1/a-1,\dfrac1{1/x-k}=\sin t$ можно свести интеграл к $$I(k)=\int\limits_0^{\pi/2}\dfrac{dt}{(1+k\sin t)^2}$$Вольфрам чихает при попытке вычислить его в параметрическом виде, но бодро рассчитывает точные значения для конкретных $k$; посчитав его для $k\in\{2,3,4,5,7,11\}$ можно угадать, что для $k>1$$$I(k)=\frac1{(k^2-1)^2}\left\{k(k^2-1)+\sqrt{k^2-1}\ln\left(k-\sqrt{k^2-1}\right)\right\}$$Не возьмусь утверждать, но похоже это работает и для $-1<k\leqslant1$; при этом, для $k\! =\!1\, (a\! =\!1/2)$ - надо брать предел, $\lim\limits_{k\rightarrow1}{I(k)}=2/3$

 Профиль  
                  
 
 Re: Определенный интеграл
Сообщение01.05.2022, 05:38 
Аватара пользователя


07/01/16
1426
Аязьма
(тут видимо имеет смысл дописать, что при $|k|<1$ выражение для $I(k)$ все равно даст действительное число; в этом можно убедиться подставив $k=\ch{w}$ и рассмотрев, что будет при мнимых значениях $w$, соответствующих данному диапазону $k$)

 Профиль  
                  
 
 Re: Определенный интеграл
Сообщение01.05.2022, 13:53 
Заслуженный участник


12/07/07
4438
Если рассматривать неопределённый интеграл, то указанный в начальном сообщении интеграл — это интеграл от квадратичной иррациональности. И он может быть найден в элементарных функциях. (См. в учебниках по началам анализа. Например, в первом томе книги Ильин, Позняк «Начала математического анализа».)

Если вопрос по Maple. Если интеграл содержит параметр, то в разных версиях присутствуют в результатах разные ошибки при нахождении интеграла. Чтобы избежать комплексных или отрицательных результатов (в рассматриваемом в начальном сообщении несобственном интеграле) можно вычислять интеграл отдельно для разных значений диапазонов параметра.
Maple 7:
Код:
> assume(a<1/2); additionally(a>0);
> F1:= int(1/sqrt(1-1/(1-1/a+1/x)^2), x=0..a);
      F1 := a*(a^2*ln(a)-2*a^3*ln(a)+a*(-2*a+1)^(3/2)-a^2*ln(1-sqrt(-2*a+1)-a)+2*a^3*ln(1-sqrt(-2*a+1)-a)-(-2*a+1)^(3/2))/(-1+2*a)/(-2*a+1)^(3/2)
> P1:= plot(F1, a=0..1/2):
> assume(a>1/2); additionally(a< 1);
> F2:= int(1/sqrt(1-1/(1-1/a+1/x)^2), x=0..a);
      F2 := 1/2*a*(-a^2*Pi+2*a^3*Pi+2*a*(-1+2*a)^(3/2)-2*a^2*arcsin(1/a*(a-1))+4*a^3*arcsin(1/a*(a-1))-2*(-1+2*a)^(3/2))/(-1+2*a)^(5/2)
> P2:= plot(F2, a=1/2..1):
> assume(a>=1);
> F3:= int(1/sqrt(1-1/(1-1/a+1/x)^2), x=0..a);
      F3 := 1/2*a*(-a^2*Pi+2*a^3*Pi+2*a*(-1+2*a)^(3/2)-2*a^2*arcsin(1/a*(a-1))+4*a^3*arcsin(1/a*(a-1))-2*(-1+2*a)^(3/2))/(-1+2*a)^(5/2)
> P3:= plot(F3, a=1..3):
> plots[display]([P1, P2, P3], thickness=3, labels=['a', 'F']);
Для проверки положительности интеграла построен график.
Вложение:
Комментарий к файлу: График зависимости интеграла от параметра. Maple 7
IntQRut.PNG
IntQRut.PNG [ 4.61 Кб | Просмотров: 842 ]

 Профиль  
                  
 
 Re: Определенный интеграл
Сообщение01.05.2022, 15:16 
Заслуженный участник


03/01/09
1677
москва
После замены $1-\frac 1a+\frac 1x=\ch t, k=\frac 1a-1, k\in (-1,\infty )$ получается:$$I(k)=\int \limits _0^{\infty }\dfrac {\ch tdt}{(\ch t+k)^2}=R+k\dfrac {dR}{dk}, R(k)=\int \limits _0^{\infty }\dfrac {dt}{\ch t+k}$$Делаем еще замену $\th \frac t2=u$, тогда $R(k)=\dfrac 2{k+1}\int \limits _0^1\dfrac {du}{1-\frac k{k+1}u^2}.$
Если $\frac k{k+1}\geq 0$, то $R(k)=\dfrac 2{\sqrt {k^2+k}}\ln (\sqrt {k+1}+\sqrt {k}).$
Если $\frac k{k+1}<0$, то $$R(k)=\dfrac 2{k+1}\sqrt {\frac {k+1}{-k}}\arctg \sqrt {\frac {k+1}{-k}}$$

-- Вс май 01, 2022 17:02:57 --

В первом сообщении была ошибка. Должно быть так:После замены $1-\frac 1a+\frac 1x=\ch t, k=\frac 1a-1, k\in (-1,\infty )$ получается:$$I(k)=\int \limits _0^{\infty }\dfrac {\ch tdt}{(\ch t+k)^2}=R+k\dfrac {dR}{dk}, R(k)=\int \limits _0^{\infty }\dfrac {dt}{\ch t+k}$$Делаем еще замену $\th \frac t2=u$, тогда $R(k)=\dfrac 2{k+1}\int \limits _0^1\dfrac {du}{1+\frac {k-1}{k+1}u^2}.$Если $\frac {k-1}{k+1}>0$, тогда $R(k)=\dfrac 2{\sqrt {k^2-1}}\arctg \sqrt {\frac {k-1}{k+1}}$.
Если $\frac {k-1}{k+1}<0$, то $$R(k)=\dfrac 2{\sqrt {1-k^2}}\ln \dfrac {\sqrt {1+k}+\sqrt {1-k}}{\sqrt {1+k}-\sqrt {1-k}}$$

 Профиль  
                  
 
 Re: Определенный интеграл
Сообщение02.05.2022, 10:11 
Заслуженный участник


12/07/07
4438
У меня после замены $\th \frac t 2 = u$ получился интеграл $R(k) = 2\int_0^1 \frac {du} {(1+k) + (1-k)u^2}$.
Если $-1 < k < 1$, то $R(k) = \frac 2 {\sqrt {1-k^2}} \arctg \sqrt{\frac {1-k}{1+k}} $.
Если $1 < k$, то $R(k) = \frac 1 {\sqrt {k^2-1}} \ln \frac{\sqrt{k+1}+\sqrt{k-1}} {\sqrt{k+1}+\sqrt{k-1}}$.
(Если вычислить по формулам из предыдущего сообщения, то $\lim\limits_{k \to 1} I(k) \ne 2/3$.)
Использование дифференцирования по параметру сильно упрощает нахождение интеграла. Замечательно Вами подмечено. Спасибо.

 Профиль  
                  
 
 Re: Определенный интеграл
Сообщение02.05.2022, 16:01 
Аватара пользователя


08/10/09
835
Херсон
GAA в сообщении #1553769 писал(а):
У меня после замены $\th \frac t 2 = u$ получился интеграл $R(k) = 2\int_0^1 \frac {du} {(1+k) + (1-k)u^2}$.
Если $-1 < k < 1$, то $R(k) = \frac 2 {\sqrt {1-k^2}} \arctg \sqrt{\frac {1-k}{1+k}} $.
Если $1 < k$, то $R(k) = \frac 1 {\sqrt {k^2-1}} \ln \frac{\sqrt{k+1}+\sqrt{k-1}} {\sqrt{k+1}+\sqrt{k-1}}$.
(Если вычислить по формулам из предыдущего сообщения, то $\lim\limits_{k \to 1} I(k) \ne 2/3$.)
Использование дифференцирования по параметру сильно упрощает нахождение интеграла. Замечательно Вами подмечено. Спасибо.


В знаменателе аргумента логарифма между корнями должен стоять знак "минус", а так все верно. Спасибо огромное всем участникам топика!!

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

Модераторы: Модераторы Математики, Супермодераторы



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

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


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

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