安卓手机扫描二维码安装App

奥斯特罗格拉茨基积分法


奥斯特罗格拉茨基积分法

微积分处理技巧之二

z张z


这篇文章是对上一篇文章的补充与拓展,参见「分项分式积分法」,同样用来处理多项式积分。


相关概念


(一)真分式


如果一个分式的分子多项式的次数小于分母多项式的次数,就称它为真分式。


(二)即约分式


如果分式p(x)/q(x) p(x)/q(x)  的分子和分母除了常数因子外,没有其它公因式,即 p(x)p(x)q(x)q(x) 互质,则此分式叫做即约分式。


正文


部分分式可以解决大量多项式问题,但是仍有一些积分处理后计算步骤亢长,如


1(x23x+3)2dx\int \dfrac{1}{\left(x^{2}-3 x+3\right)^{2}} d x


我们知道可以处理成


Px23x+3\int \dfrac{P}{x^{2}-3 x+3} +Q(x23x+3)2dx+\dfrac{Q}{\left(x^{2}-3 x+3\right)^{2}} d x


其中PPQQ 是一次式,这样通分后得到P=0P=0 ,Q=1Q=1 ,这是无法进一步简化处理的,只能逐步运算,或许,也可以换以种思路,这就用到了本文讲述的方法-奥斯特罗格拉茨基积分法。


内容简述



P(x)Q(x)\dfrac{P(x)}{Q(x)} 为【即约真分式】,那么存在【即约真分式】P1Q1\dfrac{P_{1}}{Q_{1}}P2Q2\dfrac{P_{2}}{Q_{2}} 满足


PQ=(P1Q1)+P2Q2\dfrac{P}{Q}=\left(\dfrac{P_{1}}{Q_{1}}\right)^{\prime}+\dfrac{P_{2}}{Q_{2}} 这里 Q=Q1Q2Q=Q_{1} Q_{2}


Q1Q_1 是导数 QQ^{\prime}QQ 的最大公因数,数学语言 Q1=gcd(Q,Q)Q_{1}=\text{gcd}\left(Q, Q^{\prime}\right)


P1P_1P2P_2 是次数比Q1Q_1Q2Q_2 低的多项式,写成积分形式就是


PQdx=P1Q1+P2Q2dx\int \dfrac{P}{Q} d x=\dfrac{P_{1}}{Q_{1}}+\int \dfrac{P_{2}}{Q_{2}} d x


它的形状有些类似分部积分 udv=uvvdu\int \mathrm{udv}=\mathrm{uv}-\int \mathrm{vdu} 很协调不是吗?我们只需要记住这个重要的积分形式就好,证明i 将在文末给出。



例1


现在来用它解决我们开头提到的微积分。


1(x23x+3)2dx\int \dfrac{1}{\left(x^{2}-3 x+3\right)^{2}} d x 用上述公式可以写成


1(x23x+3)2dx\int \dfrac{1}{\left(x^{2}-3 x+3\right)^{2}} d x =P1Q1+P2Q2dx=\dfrac{P_{1}}{Q_{1}}+\int \dfrac{P_{2}}{Q_{2}} d x


这里,


P=1P=1

Q=(x23x+3)2Q=\left(x^{2}-3 x+3\right)^{2}

Q=2(2x3)(x23x+3)Q^{\prime}=2(2 x-3)\left(x^{2}-3 x+3\right)


Q1=deg(Q,Q)Q_{1}=\text{deg}\left(Q, Q^{\prime}\right) =(x23x+3)=\left(x^{2}-3 x+3\right)


Q2=QQ1Q_{2}=\dfrac{Q}{Q_{1}} =(x23x+3)=\left(x^{2}-3 x+3\right)


这样就有


1(x23x+3)2dx\int \dfrac{1}{\left(x^{2}-3 x+3\right)^{2}} d x


=P1x23x+3=\dfrac{P_{1}}{x^{2}-3 x+3} +P2x23x+3dx+\int \dfrac{P_{2}}{x^{2}-3 x+3} d x



P1P_1P2P_2 是次数小于分母的多项式,不妨设 P1=Ax+BP_1=Ax+BP2=Cx+DP_2=Cx+D ,依旧通过通分求出,即


1(x23x+3)2dx\int \dfrac{1}{\left(x^{2}-3 x+3\right)^{2}} d x


=Ax+Bx23x+3=\dfrac{A x+B}{x^{2}-3 x+3} +Cx+Dx23x+3dx+\int \dfrac{C x+D}{x^{2}-3 x+3} d x (1式)


两边求导


1(x23x+3)2\dfrac{1}{\left(x^{2}-3 x+3\right)^{2}}


=A(x23x+3)(Ax+B)(2x3)(x23x+3)2=\dfrac{A\left(x^{2}-3 x+3\right)-(A \mathrm{x}+\mathrm{B})(2 x-3)}{\left(x^{2}-3 x+3\right)^{2}} +Cx+Dx23x+3+\dfrac{C x+D}{x^{2}-3 x+3}


通分解得


A=23A=\dfrac{2}{3} ,B=1B=-1 ,C=0C=0 ,D=23D=\dfrac{2}{3}


代入(1式)


1(x23x+3)2dx\int \dfrac{1}{\left(x^{2}-3 x+3\right)^{2}} d x


== 23x1x23x+3\dfrac{\dfrac{2}{3} x-1}{x^{2}-3 x+3} +23x23x+3dx+\int \dfrac{\dfrac{2}{3}}{x^{2}-3 x+3} d x


== 2x33(x23x+3)\dfrac{2 x-3}{3\left(x^{2}-3 x+3\right)} +43d(2x3)(2x3)2+3+\dfrac{4}{3} \int \dfrac{d(2 x-3)}{(2 x-3)^{2}+3}


== 2x33(x23x+3)\dfrac{2 x-3}{3\left(x^{2}-3 x+3\right)} +433arctan(2x33)+\dfrac{4}{3 \sqrt{3}} \arctan \left(\dfrac{2 x-3}{\sqrt{3}}\right) +C+C


实际上过程只有三步,只需要一点准备工作,就能简化这样的积分,参考下图。




例2



再来举一个例子,


12x+16x2+4x3+4x4+8(x+1)2(x2+1)2dx\int \dfrac{12 x+16 x^{2}+4 x^{3}+4 x^{4}+8}{(x+1)^{2}\left(x^{2}+1\right)^{2}} d x


这样的积分无论有哪种方法计算步骤都会很多,因为复杂的微积分本身计算量偏多,而奥斯特罗格拉茨基积分法ii我认为是一种清晰明了的方法。


过程如下



12x+16x2+4x3+4x4+8(x+1)2(x2+1)2dx\int \dfrac{12 x+16 x^{2}+4 x^{3}+4 x^{4}+8}{(x+1)^{2}\left(x^{2}+1\right)^{2}} d x


== Ax2+Bx+C(x+1)(x2+1)\dfrac{A x^{2}+B x+C}{(x+1)\left(x^{2}+1\right)} +Dx2+Ex+F(x+1)(x2+1)dx+\int \dfrac{D x^{2}+E x+F}{(x+1)\left(x^{2}+1\right)} d x


== xx24(x+1)(x2+1)\dfrac{ x-x^{2}-4}{(x+1)\left(x^{2}+1\right)} ++ 3x+1(x+1)(x2+1)dx3 \int \dfrac{x+1}{(x+1)\left(x^{2}+1\right)} d x


== xx24(x+1)(x2+1)\dfrac{ x-x^{2}-4}{(x+1)\left(x^{2}+1\right)} +3arctanx+3 \arctan x +C+C


只用了四步,再来看一下计算机给出的标准解法,对比下


12x+16x2+4x3+4x4+8(x+1)2(x2+1)2dx\int \dfrac{12 x+16 x^{2}+4 x^{3}+4 x^{4}+8}{(x+1)^{2}\left(x^{2}+1\right)^{2}} d x



=43x+4x2+x3+x4+2(x+1)2(x2+1)2dx=4 \int \dfrac{3 x+4 x^{2}+x^{3}+x^{4}+2}{(x+1)^{2}\left(x^{2}+1\right)^{2}} d x


=4x+22(x2+1)2=4 \int \dfrac{x+2}{2\left(x^{2}+1\right)^{2}} +14(x2+1)+\dfrac{1}{4\left(x^{2}+1\right)} +34(x+1)2dx+\dfrac{3}{4(x+1)^{2}} d x


=arctanx=\arctan x +2x+2(x2+1)2dx+2 \int \dfrac{x+2}{\left(x^{2}+1\right)^{2}} d x +31(x+1)2dx+3 \int \dfrac{1}{(x+1)^{2}} d x


x+2(x2+1)2dx\int \dfrac{x+2}{\left(x^{2}+1\right)^{2}} d x 做代换 x=tanux=\tan u


=arctanx+=\arctan x+ 2cos2u(tanu+2)du2 \int \cos ^{2} u(\tan u+2) d u +31(x+1)2dx+3 \int \dfrac{1}{(x+1)^{2}} d x


=arctanx+=\arctan x+ 2(1sin2u)(tanu+2)du2 \int\left(1-\sin ^{2} u\right)(\tan u+2) d u +31(x+1)2dx+3 \int \dfrac{1}{(x+1)^{2}} d x


=arctanx+=\arctan x+ 22sin2u+tanu2 \int-2 \sin ^{2} u+\tan u tanusin2u+2du-\tan u \sin ^{2} u+2 d u +31(x+1)2dx+3 \int \dfrac{1}{(x+1)^{2}} d x


=arctanx+=\arctan x+ 2sin2utanudu2 \int \sin ^{2} u \tan u d u +2tanudu+2 \int \tan u d u 4sin2udu-4 \int \sin ^{2} u d u +4du+4 \int d u +31(x+1)2dx+3 \int \dfrac{1}{(x+1)^{2}} d x


== \cdots \cdots


== x2+3(x3+x2+x+1)arctanx+x4(x+1)(x2+1)\dfrac{-x^{2}+3\left(x^{3}+x^{2}+x+1\right) \arctan x+x-4}{(x+1)\left(x^{2}+1\right)} +C+C


中间还有很长的步骤可以明显的看出差距。



i 证明


PQdx=P1Q1+P2Q2dx\int \dfrac{P}{Q} d x=\dfrac{P_{1}}{Q_{1}}+\int \dfrac{P_{2}}{Q_{2}} d x ,求导


PQ=P1Q1P1Q1Q2\dfrac{P}{Q}=\dfrac{P_{1}^{\prime} Q_{1}-P_{1} Q_{1}^{\prime}}{Q^{2}} +P2Q2+\dfrac{P_{2}}{Q_{2}} ,由于 Q=Q1Q1Q=Q_{1} Q_{1}


P=P1Q2P=P_{1}^{\prime} Q_{2} Q1Q1Q2P1+Q1P2-\dfrac{Q_{1}}{Q_{1}} Q_{2} P_{1}+Q_{1} P_{2}


deg(P)deg(Q)1\text{deg}(P) \leq \text{deg}(Q)-1


deg(P1Q2)(deg(Q1)1)\text{deg}\left(P_{1}^{\prime} Q_{2}\right) \leq\left(\text{deg}\left(Q_{1}\right)-1\right) +(deg(Q)deg(Q1))+\left(\text{deg}(Q)-\text{deg}\left(Q_{1}\right)\right)


=deg(Q)1=\text{deg}(Q)-1



deg(Q1Q1Q2P)(deg(Q1)1)\text{deg}\left(\frac{Q_{1}^{\prime}}{Q_{1}} Q_{2} P\right) \leq\left(\text{deg}\left(Q_{1}\right)-1\right) +(deg(Q)deg(Q1))+\left(\text{deg}(Q)-\text{deg}\left(Q_{1}\right)\right) +(deg(Q1)1)+\left(\text{deg}\left(Q_{1}\right)-1\right) deg(Q1)-\text{deg}\left(Q_{1}\right)


=deg(Q)2=\text{deg}(Q)-2


deg(Q1P2)(deg(Q)deg(Q1)1)\text{deg}\left(Q_{1} P_{2}\right) \leq\left(\text{deg}(Q)-\text{deg}\left(Q_{1}\right)-1\right) +deg(Q1)+\text{deg}\left(Q_{1}\right)


=deg(Q)1=\text{deg}(Q)-1


P1,P2P_1,P_2 有解,证毕。



ii 奥斯特罗格拉茨基


米哈伊尔·瓦西里耶维奇·奥斯特罗格拉德斯基(也译作奥斯特罗格拉茨基,1801年9月24日-1862年1月1日)生平介绍


俄国数学家,他的工作包括变分学、代数函数的积分、数学物理、古典力学,是当时俄国数学界的重要人物。


他出生于现时乌克兰境内的Пашенная,在1816年入读哈尔科夫大学,1820年他参加毕业考试,但宗教兼国民教育主任要求他重考。官方原因是他没有上哲学和神学课,实际上的原因是他的老师奥西波夫斯基之前因宗教问题被停职,而校方认为这事也和奥西波夫斯基的学生有关。结果,他拒绝重考,不取学位便离开俄罗斯到巴黎求学。他上著名数学家的课,包括勒让德、泊松、柯西。他在法国科学院发表论文,内容以物理和积分学为主。


1828年他回到俄罗斯,到了圣彼得堡。曾在海军学院、通讯研究所、教育学院讲课。



苹果手机扫描二维码安装App