Метод Остроградского

Метод Остроградского — метод интегрирования рациональных функций с кратными неприводимыми множителями в знаменателе. Метод позволяет одними лишь алгебраическими операциями свести задачу интегрирования произвольной рациональной функции к задаче интегрирования рациональной функции без кратных корней в знаменателе.

ИсторияПравить

Метод Остроградского назван по имени М. В. Остроградского, впервые предложившего его 22 ноября 1844 года на заседании физико-математического отделения Академии наук, опубликован в следующем году на французском языке, статья переведена на русский в 1958 г.

Описание методаПравить

Любой интеграл от рациональной функции можно представить в виде

$\int <!--\; \int \; -->\frac{P(x)}{Q(x)}dx=\frac{P_1(x)}{Q_1(x)}+\int <!--\; \int \; -->\frac{P_2(x)}{Q_2(x)}dx$$\int {\frac {P(x)}{Q(x)}}dx={\frac {P_{1}(x)}{Q_{1}(x)}}+\int {\frac {P_{2}(x)}{Q_{2}(x)}}dx$ .

Здесь $Q_2(x)$$Q_{2}(x)$  представляет собой произведение всех неприводимых множителей многочлена $Q(x)$$Q(x)$  без учёта кратности (то есть каждый неприводимый множитель многочлена $Q(x)$$Q(x)$  встречается в разложении многочлена $Q_2(x)$$Q_{2}(x)$  один раз), $Q_1(x)$$Q_{1}(x)$  — произведение всех неприводимых множителей многочлена $Q(x)$$Q(x)$  с пониженной на 1 кратностью (каждый неприводимый множитель многочлена $Q(x)$$Q(x)$  кратности $n$$n$  встречается в разложении многочлена $Q_1(x)$$Q_{1}(x)$  $n-<!--\; -\; -->1$$n-1$  раз). Дробь $\frac{P_2(x)}{Q_2(x)}$${\frac {P_{2}(x)}{Q_{2}(x)}}$  является правильной. Эта формула называется формулой Остроградского. $\frac{P_1(x)}{Q_1(x)}$${\frac {P_{1}(x)}{Q_{1}(x)}}$  здесь есть алгебраическая (рациональная) часть интеграла от рациональной функции $\int <!--\; \int \; -->\frac{P(x)}{Q(x)}dx$$\int {\frac {P(x)}{Q(x)}}dx$ , а $\int <!--\; \int \; -->\frac{P_2(x)}{Q_2(x)}dx$$\int {\frac {P_{2}(x)}{Q_{2}(x)}}dx$  — трансцендентная.

Суть метода заключается в следующем. Запишем многочлены $P_1(x)$$P_{1}(x)$  и $P_2(x)$$P_{2}(x)$  с неопределёнными коэффициентами:

$P_1(x)=A_s{x}^s+\dots <!--\; \dots \; -->+A_{1}x+A_0$$P_{1}(x)=A_{s}x^{s}+\ldots +A_{1}x+A_{0}$ 

$P_2(x)=B_r{x}^r+\dots <!--\; \dots \; -->+B_{1}x+B_0$$P_{2}(x)=B_{r}x^{r}+\ldots +B_{1}x+B_{0}$ .

Степени многочленов можно выяснить позже, а можно заранее взять наверняка. Пусть далее $\mathrm{deg}<!--\; \; -->P=n,\mathrm{deg}<!--\; \; -->Q=m,\mathrm{deg}<!--\; \; -->Q_2=p,\mathrm{deg}<!--\; \; -->Q_1=m-<!--\; -\; -->p$$\operatorname {deg} {P}=n,\ \operatorname {deg} {Q}=m,\ \operatorname {deg} {Q_{2}}=p,\ \operatorname {deg} {Q_{1}}=m-p$ . Дробь под интегралом должна получиться правильной, поэтому степень $P_2$$P_2$  можно взять за $p-<!--\; -\; -->1$$p-1$ . Если первоначальная дробь была правильной, то и $\frac{P_1}{Q_1}$${\frac {P_{1}}{Q_{1}}}$  правильная и можно взять степень $P_1$$P_1$  как $m-<!--\; -\; -->p-<!--\; -\; -->1$$m-p-1$ . Если же она неправильная, то выделить целую часть и свести дробь к правильной (или же взять степень такую, чтобы степени целых частей слева и справа совпали).

Теперь мы можем найти коэффициенты этих многочленов методом неопределённых коэффициентов. Продифференцируем это равенство.

$\frac{P(x)}{Q(x)}=\frac{P_1^{\prime }(x)Q_1(x)-<!--\; -\; -->P_1(x)Q_1^{\prime }(x)}{Q_1^2(x)}+\frac{P_2(x)}{Q_2(x)}$${\frac {P(x)}{Q(x)}}={\frac {P_{1}'(x)Q_{1}(x)-P_{1}(x)Q_{1}'(x)}{Q_{1}^{2}(x)}}+{\frac {P_{2}(x)}{Q_{2}(x)}}$ 

Умножим обе части на $Q(x)$$Q(x)$ .

$P(x)=\frac{P_1^{\prime }(x)Q(x)}{Q_1(x)}-<!--\; -\; -->\frac{P_1(x)Q_1^{\prime }(x)Q(x)}{Q_1^2(x)}+\frac{P_2(x)Q(x)}{Q_2(x)}=P_1^{\prime }(x)Q_2(x)-<!--\; -\; -->P_1(x)\frac{Q_1^{\prime }(x)Q_2(x)}{Q_1(x)}+P_2(x)Q_1(x)$$P(x)={\frac {P_{1}'(x)Q(x)}{Q_{1}(x)}}-{\frac {P_{1}(x)Q_{1}'(x)Q(x)}{Q_{1}^{2}(x)}}+{\frac {P_{2}(x)Q(x)}{Q_{2}(x)}}=P_{1}'(x)Q_{2}(x)-P_{1}(x){\frac {Q_{1}'(x)Q_{2}(x)}{Q_{1}(x)}}+P_{2}(x)Q_{1}(x)$ 

В обоих частях равенства стоят многочлены. $P_1(x)\frac{Q_1^{\prime }(x)Q_2(x)}{Q_1(x)}$$P_{1}(x){\frac {Q_{1}'(x)Q_{2}(x)}{Q_{1}(x)}}$  здесь тоже многочлен, так как $Q_1^{\prime }(x)Q_2(x)$$Q_{1}'(x)Q_{2}(x)$  делится на $Q_1(x)$$Q_{1}(x)$ . Приравниваем коэффициенты при равных степенях и получаем систему линейных алгебраических уравнений. Решая её, получаем в итоге коэффициенты многочленов $P_1(x)$$P_{1}(x)$  и $P_2(x)$$P_{2}(x)$ .

В итоге мы представили первоначальный интеграл в виде $\frac{P_1(x)}{Q_1(x)}+\int <!--\; \int \; -->\frac{P_2(x)}{Q_2(x)}dx$${\frac {P_{1}(x)}{Q_{1}(x)}}+\int {\frac {P_{2}(x)}{Q_{2}(x)}}dx$ . Задача свелась к интегрированию дроби без кратных неприводимых множителей в знаменателе.

Формула $P(x)=P_1^{\prime }(x)Q_2(x)-<!--\; -\; -->P_1(x)\frac{Q_1^{\prime }(x)Q_2(x)}{Q_1(x)}+P_2(x)Q_1(x)$$P(x)=P_{1}'(x)Q_{2}(x)-P_{1}(x){\frac {Q_{1}'(x)Q_{2}(x)}{Q_{1}(x)}}+P_{2}(x)Q_{1}(x)$  позволяет более точно подобрать степени для многочленов $P_1(x)$$P_{1}(x)$  и $P_2(x)$$P_{2}(x)$ . Если приравнять степени всех слагаемых, то получим $s=n-<!--\; -\; -->p+1$$s=n-p+1$  и $r=n+p-<!--\; -\; -->m$$r=n+p-m$ .

Метод Остроградского позволяет сразу же получить алгебраическую часть интеграла рациональной функции. Более того, для этого даже не нужно вычислять разложение $Q(x)$$Q(x)$  на неприводимые. Действительно, $Q_1(x)=\text{НОД}(Q,Q^{\prime })$$Q_{1}(x)={\text{НОД}}(Q,Q')$ , $Q_2(x)=\frac{Q}{\text{НОД}(Q,Q^{\prime })}$$Q_{2}(x)={\frac {Q}{{\text{НОД}}(Q,Q')}}$ . НОД многочленов же можно вычислить при помощи алгоритма Евклида. Таким образом, алгебраическая часть интеграла от рациональной функции может быть найдена при помощи метода Остроградского с использованием только лишь алгебраических операций.

ДоказательствоПравить

Доказательство того, что для любой рациональной дроби можно записать формулу Остроградского, получается сразу же из общего вида интеграла.

Запишем общий вид интеграла от рациональной функции.

$\int <!--\; \int \; -->\frac{P(x)}{\underset{i=1}{\overset{l}{\prod <!--\; \prod \; -->}}}dx=\frac{T(x)}{\underset{i=1}{\overset{l}{\prod <!--\; \prod \; -->}}}+\underset{i=1}{\overset{l}{\sum <!--\; \sum \; -->}}{A}_i\ln <!--\; \; -->|x-<!--\; -\; -->x_i|+\underset{i=1}{\overset{n}{\sum <!--\; \sum \; -->}}\left(B_i\ln <!--\; \; -->(x^2+p_{i}x+q_i)+C_i\mathrm{arctg}<!--\; \; -->\frac{x+r_i}{h}\right)$$\int {\frac {P(x)}{\displaystyle \prod _{i=1}^{l}(x-x_{i})^{k_{i}}\displaystyle \prod _{i=1}^{n}(x^{2}+p_{i}x+q_{i})^{m_{i}}}}\,dx={\frac {T(x)}{\displaystyle \prod _{i=1}^{l}(x-x_{i})^{k_{i}-1}\displaystyle \prod _{i=1}^{n}(x^{2}+p_{i}x+q_{i})^{m_{i}-1}}}+\sum _{i=1}^{l}{A_{i}\ln {|x-x_{i}|}}+\sum _{i=1}^{n}\left({B_{i}\ln(x^{2}+p_{i}x+q_{i})+C_{i}\operatorname {arctg} {\frac {x+r_{i}}{h}}}\right)$ 

$x+r_i$$x+r_{i}$  здесь линейный двучлен, получаемый выделением полного квадрата из $x^2+p_{i}x+q_i$$x^{2}+p_{i}x+q_{i}$ , т. е. $x^2+p_{i}x+q_i=(x+r_i{)}^2+h^2$$x^{2}+p_{i}x+q_{i}=(x+r_{i})^{2}+h^{2}$ . Занесём логарифмы и арктангенсы под интеграл.

$\int <!--\; \int \; -->\frac{P(x)}{\underset{i=1}{\overset{l}{\prod <!--\; \prod \; -->}}}dx=\frac{T(x)}{\underset{i=1}{\overset{l}{\prod <!--\; \prod \; -->}}}+\int <!--\; \int \; -->\left(\underset{i=1}{\overset{l}{\sum <!--\; \sum \; -->}}\frac{A_i}{x-<!--\; -\; -->x_i}+\underset{i=1}{\overset{n}{\sum <!--\; \sum \; -->}}\left(\frac{B_i(2x+p_i)+C_{i}h}{x^2+p_{i}x+q_i}\right)\right)dx$$\int {\frac {P(x)}{\displaystyle \prod _{i=1}^{l}(x-x_{i})^{k_{i}}\displaystyle \prod _{i=1}^{n}(x^{2}+p_{i}x+q_{i})^{m_{i}}}}\,dx={\frac {T(x)}{\displaystyle \prod _{i=1}^{l}(x-x_{i})^{k_{i}-1}\displaystyle \prod _{i=1}^{n}(x^{2}+p_{i}x+q_{i})^{m_{i}-1}}}+\int \left(\sum _{i=1}^{l}{\frac {A_{i}}{x-x_{i}}}+\sum _{i=1}^{n}\left({\frac {B_{i}(2x+p_{i})+C_{i}h}{x^{2}+p_{i}x+q_{i}}}\right)\right)\,dx$ 

Полученная формула и есть формула Остроградского. Дробь под интегралом правильная, поскольку является суммой правильных дробей.