在各类积分运算中,我们常见到这样一类函数,它们有貌离神合的计算方法——有理函数。

例如
∫1xdx=ln⁡∣x∣+C,  x≠0 \int \frac{1}{x}\mathrm{d}x=\ln|x|+C, \ \ x \ne 0 x1dx=lnx+C,  x=0
又如
∫xndx=xn+1n+1+C,  n≠−1 \int x^n\mathrm{d}x=\frac{x^{n+1}}{n+1}+C, \ \ n\ne -1 xndx=n+1xn+1+C,  n=1
再如
∫11+x2dx=arctan⁡x+C \int \frac{1}{1+x^2}\mathrm{d}x=\arctan x+C 1+x21dx=arctanx+C
上述三例倒也不是随意举出来的,利用上述三个最简积分公式加上一些技巧,我们可以计算出任一有理函数的(不定)积分,下面我们就来介绍具体的范式。

一些定积分的技巧我们在这里不做讨论,下面所言积分都是不定积分

有理函数

首先,我们明确一下何谓有理函数,形如:
f(x)=Pm(x)Qn(x) f(x)=\frac{P_m(x)}{Q_n(x)} f(x)=Qn(x)Pm(x)
其中
Pm(x)=∑i=0maixi,Qn(x)=∑j=0nbjxj,Qn(x)≠0 P_m(x)=\sum_{i=0}^ma_ix^i,Q_n(x)=\sum_{j=0}^nb_jx^j,Q_n(x)\ne0 Pm(x)=i=0maixi,Qn(x)=j=0nbjxj,Qn(x)=0
也就是两个多项式的商称为有理函数,这个概念显然是有理数某种意义上的推广:
Q={f=pq∣p,q∈Z,q≠0} \mathbb{Q}=\left\{f=\frac{p}{q}\mid p,q\in \mathbb{Z}, q\ne 0\right\} Q={f=qpp,qZ,q=0}
并且类似的,我们也定义真分式:
proper fraction={f(x)=Pm(x)Qn(x)∣m<n} \text{proper fraction}=\left\{f(x)=\frac{P_m(x)}{Q_n(x)}\mid m<n\right\} proper fraction={f(x)=Qn(x)Pm(x)m<n}
和假分式:
improper fraction={f(x)=Pm(x)Qn(x)∣m≥n} \text{improper fraction}=\left\{f(x)=\frac{P_m(x)}{Q_n(x)}\mid m\ge n\right\} improper fraction={f(x)=Qn(x)Pm(x)mn}

上面写成了集合形式,有些显而易见的条件就略去了,并且在我们讨论有理函数积分的场景下,不妨假设PmP_mPmQnQ_nQn最高次项的系数都为1

分解技术

从上面的定义看得出来,有理函数和多项式的关系非常密切,实际上为了求解一般有理函数的积分我们需要下面两个分解:

因式分解

因式分解是对于一个多项式而言的,根据代数基本定理,Qn(x)Q_n(x)Qn(x)复数域内有如下分解形式:
Qn(x)=∏i=1n(x−xi) Q_n(x)=\prod_{i=1}^n(x-x_i) Qn(x)=i=1n(xxi)
其中xi∈Cx_i\in \mathbb{C}xiCQn(x)=0Q_n(x)=0Qn(x)=0的根。
更进一步,我们有:
Qn(x)=∏i=1n1(x−xi)∏j=1(n−n1)/2(x−xj)(x−xj‾) Q_n(x)=\prod_{i=1}^{n_1}(x-x_i)\prod_{j=1}^{(n-n_1)/2}(x-x_j)(x-\overline{x_j}) Qn(x)=i=1n1(xxi)j=1(nn1)/2(xxj)(xxj)
其中xj‾\overline{x_j}xj代表复数xj,j=1,⋯ ,n−n12x_j,j=1,\cdots,\frac{n-n_1}{2}xj,j=1,,2nn1的共轭,xi,i=1,⋯ ,n1x_i,i=1,\cdots,n_1xi,i=1,,n1都是实数。

能做出第二个形式是因为,实系数多项式的零点若不是实数,那么它的共轭必然也是零点,也就是说非实数零点是成对出现的

这是显然成立的,Q(x)=0Q(x)=0Q(x)=0两边取共轭有,Q(x)‾=0\overline{Q(x)}=0Q(x)=0,根据复数的性质(加、减、乘运算和共轭运算可交换),所以有Q(x‾)=0Q(\overline{x})=0Q(x)=0,这样就说明了非实数零点成对出现。

如果取xj=αj+βjix_j=\alpha_j + \beta_j \mathrm{i}xj=αj+βji 其中αj,βj∈R\alpha_j,\beta_j \in \mathbb{R}αj,βjR,那么很容易得到下面的分解形式:
Qn(x)=∏i=1n1(x−xi)∏j=1(n−n1)/2(x2−2αjx+αj2+βj2) Q_n(x)=\prod_{i=1}^{n_1}(x-x_i)\prod_{j=1}^{(n-n_1)/2}(x^2-2\alpha_j x+\alpha_j^2+\beta_j^2) Qn(x)=i=1n1(xxi)j=1(nn1)/2(x22αjx+αj2+βj2)
这也就是Qn(x)Q_n(x)Qn(x)实数域内的因式分解。

如果我们把重根合并在一起,就有
f(x)=∏i=1s(x−ai)di∏j=1t(x2−2αjx+αj2+βj2)ej f(x)=\prod_{i=1}^{s}(x-a_i)^{d_i}\prod_{j=1}^{t}(x^2-2\alpha_j x+\alpha_j^2+\beta_j^2)^{e_j} f(x)=i=1s(xai)dij=1t(x22αjx+αj2+βj2)ej
其中did_idi为实根aia_iai的重数,eje_jej为共轭复根对αj+βji\alpha_j + \beta_j \mathrm{i}αj+βji的重数。

部分分式分解

有了上述Qn(x)Q_n(x)Qn(x)的实数域内的分解形式,我们就可以得到下面有理真分式的部分分式分解形式:
P(x)Q(x)=∑i=1s∑k=1diAik(x−ai)k+∑j=1t∑k=1ejBjkx+Cjk(x2−2αjx+αj2+βj2)k \begin{aligned} &\frac{P(x)}{Q(x)}\\ =&\sum_{i=1}^{s}\sum_{k=1}^{d_i}\frac{A_{ik}}{(x-a_i)^k}+\sum_{j=1}^{t}\sum_{k=1}^{e_j}\frac{B_{jk}x+C_{jk}}{(x^2-2\alpha_j x+\alpha_j^2+\beta_j^2)^k}\\ \end{aligned} =Q(x)P(x)i=1sk=1di(xai)kAik+j=1tk=1ej(x22αjx+αj2+βj2)kBjkx+Cjk
其中Ai,Bj,Cj∈RA_i,B_j,C_j \in \mathbb{R}Ai,Bj,CjR

可以证明,上述待定系数一定有解。

我见过两种证明方法,一种是使用多项式的理论,详见

傅莺莺.有理真分式部分分式分解的证明及系数公式[J].大学数学,2014,30(02):82-87.

另外一种是使用线性空间的理论,详见

高翔宇,杨洪福,王世鹏.有理分式函数分解定理的新证明[J].大学教育,2020(12):103-104+154.

我们这里不多赘述,下面直接利用部分分式分解的形式来计算有理函数的积分。

积分范式

我们用一个例题来呈现有理函数积分的一般过程。

例一

F(x)=∫x4x3+1dx F(x)=\int \frac{x^4}{x^3+1}\mathrm{d}x F(x)=x3+1x4dx

分离整式

如果f(x)f(x)f(x)是一个有理假分式, 我们使用带余除法可以把f(x)f(x)f(x)分解为一个整式和有理分式, 积分也就分为两块, 整式部分的积分非常简单, 真分式的部分进入下一流程. 另外, 如果分母不是首一多项式, 直接把系数提出来即可.

例1 步骤1, 分离整式:
F(x)=∫x4x3+1dx=∫[x−xx3+1]dx=12x2−∫xx3+1dx \begin{aligned} &F(x)\\ =&\int \frac{x^4}{x^3+1}\mathrm{d}x\\ =&\int \left[x-\frac{x}{x^3+1}\right]\mathrm{d}x\\ =&\frac{1}{2}x^2-\int \frac{x}{x^3+1}\mathrm{d}x \end{aligned} ===F(x)x3+1x4dx[xx3+1x]dx21x2x3+1xdx

部分分式分解

前文提到任何一个有理真分式都可以做部分分式分解, 其中待定系数的求法较为灵活. 主要包括三种, 其一是极限法: 对分解式取不同极限可以得到系数的值, 其二是列出待定系数的线性方程组求解(通用), 其三是利用导数求解.

例1 步骤2, 有理真分式的部分分解
x3+1=(x+1)(x2−x+1) x^3+1 = (x+1)(x^2-x+1) x3+1=(x+1)(x2x+1)
于是有部分分解
xx3+1=x(x+1)(x2−x+1)=Ax+1+Bx+Cx2−x+1 \begin{aligned} &\frac{x}{x^3+1}\\ =&\frac{x}{(x+1)(x^2-x+1)}\\ =&\frac{A}{x+1}+\frac{Bx+C}{x^2-x+1} \end{aligned} ==x3+1x(x+1)(x2x+1)xx+1A+x2x+1Bx+C
通分得到
x=A(x2−x+1)+(Bx+C)(x+1)(*) \begin{aligned} x=A(x^2-x+1)+(Bx+C)(x+1) \tag{*} \end{aligned} x=A(x2x+1)+(Bx+C)(x+1)(*)
下面求待定系数.

法一(极限法)

(*)式取x→−1x\to -1x1立得A=−13A=-\frac{1}{3}A=31;

再取x→0x \to 0x0立得C=−A=13C=-A=\frac{1}{3}C=A=31;

最后取x→1x \to 1x1可得B=13B=\frac{1}{3}B=31.

法二(线性方程组)

比较(*)式中各次项系数即可得到线性方程组:
x2 ∣ A+B=0x1 ∣ −A+B+C=1x0 ∣ A+C=0 \begin{aligned} x^2& \ | \ A+B=0\\ x^1& \ | \ -A+B+C=1\\ x^0& \ | \ A+C=0 \end{aligned} x2x1x0  A+B=0  A+B+C=1  A+C=0
解得−A=B=C=13-A=B=C=\frac{1}{3}A=B=C=31

法三(导数)

对(*)求0,1,2阶导数可得:
0 ∣ x=A(x2−x+1)+(Bx+C)(x+1)1 ∣ 1=2Ax−A+B(x+1)+Bx+C2 ∣ 0=2A+B+B \begin{aligned} 0& \ | \ x=A(x^2-x+1)+(Bx+C)(x+1)\\ 1& \ | \ 1=2Ax-A+B(x+1)+Bx+C\\ 2& \ | \ 0=2A+B+B \end{aligned} 012  x=A(x2x+1)+(Bx+C)(x+1)  1=2AxA+B(x+1)+Bx+C  0=2A+B+B
同样比较各次项系数可得方程组, 进而解得−A=B=C=13-A=B=C=\frac{1}{3}A=B=C=31.

得到部分分式分解形式
F(x)=12x2−∫xx3+1dx=12x2−∫[−13x+1+13x+13x2−x+1]dx \begin{aligned} &F(x)\\ =&\frac{1}{2}x^2-\int \frac{x}{x^3+1}\mathrm{d}x\\ =&\frac{1}{2}x^2-\int \left[ \frac{-\frac{1}{3}}{x+1}+\frac{\frac{1}{3}x+\frac{1}{3}}{x^2-x+1}\right]\mathrm{d}x \end{aligned} ==F(x)21x2x3+1xdx21x2[x+131+x2x+131x+31]dx

两类积分

做完部分分式分解之后,我们需要处理的只有两种积分

分母为一次项

也就是
∫Ax−adx=Aln⁡∣x−a∣+C \int \frac{A}{x-a}\mathrm{d}x=A\ln|x-a|+C xaAdx=Alnxa+C

分母为二次项

也就是
∫Bx+Cx2−2αx+βdx=∫B(x−α)+C+αB(x−α)2+β−α2d(x−α)=∫Bt+C+αBt2+β−α2dt=∫12Bt2+β−α2dt2+∫C+αBt2+β−α2dt=B2ln⁡∣t2+β−α2∣+∫C+αBt2+β−α2dt \begin{aligned} &\int \frac{Bx+C}{x^2-2\alpha x+\beta}\mathrm{d}x\\ =&\int \frac{B(x-\alpha)+C+\alpha B}{(x-\alpha)^2+\beta-\alpha^2}\mathrm{d}(x-\alpha)\\ =&\int \frac{Bt+C+\alpha B}{t^2+\beta-\alpha^2}\mathrm{d}t\\ =&\int \frac{1}{2}\frac{B}{t^2+\beta-\alpha^2}\mathrm{d}t^2+\int \frac{C+\alpha B}{t^2+\beta-\alpha^2}\mathrm{d}t\\ =&\frac{B}{2}\ln|t^2+\beta-\alpha^2|+\int \frac{C+\alpha B}{t^2+\beta-\alpha^2}\mathrm{d}t \end{aligned} ====x22αx+βBx+Cdx(xα)2+βα2B(xα)+C+αBd(xα)t2+βα2Bt+C+αBdt21t2+βα2Bdt2+t2+βα2C+αBdt2Blnt2+βα2+t2+βα2C+αBdt
β−α2=0\beta-\alpha^2=0βα2=0
∫C+αBt2+β−α2dt=−C+αBt+C \begin{aligned} &\int \frac{C+\alpha B}{t^2+\beta-\alpha^2}\mathrm{d}t\\ =&-\frac{C+\alpha B}{t}+C \end{aligned} =t2+βα2C+αBdttC+αB+C
β−α2<0\beta-\alpha^2<0βα2<0
∫C+αBt2+β−α2dt=C+αB2−β+α2∫[1t+−β+α2−1t−−β+α2]dt=C+αB2−β+α2[ln⁡∣t+−β+α2∣+ln⁡∣t−−β+α2∣]+C \begin{aligned} &\int \frac{C+\alpha B}{t^2+\beta-\alpha^2}\mathrm{d}t\\ =&\frac{C+\alpha B}{2\sqrt{-\beta+\alpha^2}}\int \left[\frac{1}{t+\sqrt{-\beta+\alpha^2}}-\frac{1}{t-\sqrt{-\beta+\alpha^2}}\right]\mathrm{d}t\\ =&\frac{C+\alpha B}{2\sqrt{-\beta+\alpha^2}} \left[\ln|t+\sqrt{-\beta+\alpha^2}|+\ln|t-\sqrt{-\beta+\alpha^2}|\right]+C \end{aligned} ==t2+βα2C+αBdt2β+α2 C+αB[t+β+α2 1tβ+α2 1]dt2β+α2 C+αB[lnt+β+α2 +lntβ+α2 ]+C
β−α2>0\beta-\alpha^2>0βα2>0
∫C+αBt2+β−α2dt=C+αBβ−α2arctan⁡(tβ−α2)+C \begin{aligned} &\int \frac{C+\alpha B}{t^2+\beta-\alpha^2}\mathrm{d}t\\ =&\frac{C+\alpha B}{\sqrt{\beta-\alpha^2}}\arctan (\frac{t}{\sqrt{\beta-\alpha^2}})+C \end{aligned} =t2+βα2C+αBdtβα2 C+αBarctan(βα2 t)+C
以上三种情况注意回代t=x−αt=x-\alphat=xα

利用以上两种积分公式, 我们就可以完全解决部分分解之后的积分问题了.

例1 步骤3, 部分分式分解完之后应用积分公式求解
F(x)=12x2−∫[−13x+1+13x+13x2−x+1]dx=12x2+13ln⁡∣x+1∣−∫13(x−0.5)+12(x−0.5)2+0.75d(x−0.5)=12x2+13ln⁡∣x+1∣−∫13t+12t2+0.75dt=12x2+13ln⁡∣x+1∣−16ln⁡∣(x−12)2+34∣+33arctan⁡(233(x−12))+C=16{3x2+2ln⁡∣x+1∣−ln⁡(x2−x+1)−23arctan⁡(2x−13)}+C \begin{aligned} &F(x)\\ =&\frac{1}{2}x^2-\int \left[ \frac{-\frac{1}{3}}{x+1}+\frac{\frac{1}{3}x+\frac{1}{3}}{x^2-x+1}\right]\mathrm{d}x\\ =&\frac{1}{2}x^2+\frac{1}{3}\ln|x+1|-\int \frac{\frac{1}{3}(x-0.5)+\frac{1}{2}}{(x-0.5)^2+0.75}\mathrm{d}(x-0.5)\\ =&\frac{1}{2}x^2+\frac{1}{3}\ln|x+1|-\int \frac{\frac{1}{3}t+\frac{1}{2}}{t^2+0.75}\mathrm{d}t\\ =&\frac{1}{2}x^2+\frac{1}{3}\ln|x+1|-\frac{1}{6}\ln|(x-\frac{1}{2})^2+\frac{3}{4}|+\frac{\sqrt{3}}{3}\arctan(\frac{2\sqrt{3}}{3}(x-\frac{1}{2}))+C\\ =&\frac{1}{6}\left\{3x^2+2\ln|x+1|-\ln(x^2-x+1)-2\sqrt{3}\arctan(\frac{2x-1}{\sqrt{3}})\right\}+C \end{aligned} =====F(x)21x2[x+131+x2x+131x+31]dx21x2+31lnx+1(x0.5)2+0.7531(x0.5)+21d(x0.5)21x2+31lnx+1t2+0.7531t+21dt21x2+31lnx+161ln(x21)2+43+33 arctan(323 (x21))+C61{3x2+2lnx+1ln(x2x+1)23 arctan(3 2x1)}+C

例题

我们再给出几个例题

例二(分母无重因子)

∫dx(x2+1)(x2+x) \int \frac{\mathrm{d}x}{(x^2+1)(x^2+x)} (x2+1)(x2+x)dx
解:
∫dx(x2+1)(x2+x)=∫[1x−12(x+1)−1+x2(x2+1)]dx=ln⁡∣x∣−12ln⁡∣x+1∣−12arctan⁡x−14∫d(x2+1)x2+1=ln⁡∣x∣−12ln⁡∣x+1∣−12arctan⁡x−14ln⁡(x2+1)+C \begin{aligned} &\int \frac{\mathrm{d}x}{(x^2+1)(x^2+x)}\\ &=\int\left[\frac{1}{x}-\frac{1}{2(x+1)}-\frac{1+x}{2\left(x^{2}+1\right)}\right] \mathrm{d} x \\ &=\ln |x|-\frac{1}{2} \ln |x+1|-\frac{1}{2} \arctan x-\frac{1}{4} \int \frac{\mathrm{d}\left(x^{2}+1\right)}{x^{2}+1} \\ &=\ln |x|-\frac{1}{2} \ln |x+1|-\frac{1}{2} \arctan x-\frac{1}{4} \ln \left(x^{2}+1\right)+C \end{aligned} (x2+1)(x2+x)dx=[x12(x+1)12(x2+1)1+x]dx=lnx21lnx+121arctanx41x2+1d(x2+1)=lnx21lnx+121arctanx41ln(x2+1)+C

例三(分母有重因子)

∫x2+1(x+1)2(x−1)dx \int \frac{x^{2}+1}{(x+1)^{2}(x-1)} \mathrm{d} x (x+1)2(x1)x2+1dx
解:
∫x2+1(x+1)2(x−1)dx=∫[12(x+1)−1(x+1)2+12(x−1)dx]=12ln⁡∣x2−1∣+1x+1+C \begin{aligned} &\int \frac{x^{2}+1}{(x+1)^{2}(x-1)} \mathrm{d} x\\ =&\int\left[\frac{1}{2(x+1)}-\frac{1}{(x+1)^{2}}+\frac{1}{2(x-1)} \mathrm{d}x\right]\\ =&\frac{1}{2} \ln \left|x^{2}-1\right|+\frac{1}{x+1}+C \end{aligned} ==(x+1)2(x1)x2+1dx[2(x+1)1(x+1)21+2(x1)1dx]21lnx21+x+11+C

例四(分母无实根)

∫x2+1x4+1dx \int \frac{x^2+1}{x^4+1}\mathrm{d}x x4+1x2+1dx
解:
F(x)=12∫dxx2−2x+1+12∫dxx2+2x+1=12∫dx(x−22)2+12+12∫dx(x+22)2+12=12arctan⁡(2x−1)+12arctan⁡(2x+1)+C \begin{aligned} F(x)=& \frac{1}{2} \int \frac{d x}{x^{2}-\sqrt{2} x+1}+ \frac{1}{2} \int \frac{d x}{x^{2}+\sqrt{2} x+1} \\ =& \frac{1}{2} \int \frac{d x}{\left(x-\frac{\sqrt{2}}{2}\right)^{2}+\frac{1}{2}}+ \frac{1}{2} \int \frac{d x}{\left(x+\frac{\sqrt{2}}{2}\right)^{2}+\frac{1}{2}}\\ =&\frac{1}{\sqrt{2}} \arctan (\sqrt{2} x-1)+ \frac{1}{\sqrt{2}} \arctan (\sqrt{2} x+1)+C \end{aligned} F(x)===21x22 x+1dx+21x2+2 x+1dx21(x22 )2+21dx+21(x+22 )2+21dx2 1arctan(2 x1)+2 1arctan(2 x+1)+C

例五(经典题)

∫dxx3+1 \int \frac{\mathrm{d} x}{x^{3}+1} x3+1dx
解:
∫dxx3+1=∫[13(x+1)−x−23(x2−x+1)]dx=13∫dxx+1−16∫2x−1x2−x+1 dx+12∫d(x−12)(x−12)2+34=16ln⁡(x+1)2x2−x+1+13arctan⁡2x−13+C. \begin{aligned} & \int \frac{\mathrm{d} x}{x^{3}+1}\\ =&\int\left[\frac{1}{3(x+1)}-\frac{x-2}{3\left(x^{2}-x+1\right)}\right] \mathrm{d} x \\ =& \frac{1}{3} \int \frac{\mathrm{d} x}{x+1}-\frac{1}{6} \int \frac{2 x-1}{x^{2}-x+1} \mathrm{~d} x+\frac{1}{2} \int \frac{\mathrm{d}\left(x-\frac{1}{2}\right)}{\left(x-\frac{1}{2}\right)^{2}+\frac{3}{4}} \\ =& \frac{1}{6} \ln \frac{(x+1)^{2}}{x^{2}-x+1}+\frac{1}{\sqrt{3}} \arctan \frac{2 x-1}{\sqrt{3}}+C . \end{aligned} ===x3+1dx[3(x+1)13(x2x+1)x2]dx31x+1dx61x2x+12x1 dx+21(x21)2+43d(x21)61lnx2x+1(x+1)2+3 1arctan3 2x1+C.

例六(最好别用我们的范式)

∫dxx4−1 \int \frac{\mathrm{d} x}{x^{4}-1} x41dx
解:
∫dxx4−1=12∫[1x2−1−1x2+1]dx=14ln⁡∣x−1x+1∣−12arctan⁡x+C \begin{aligned} &\int \frac{\mathrm{d} x}{x^{4}-1}\\ =&\frac{1}{2} \int\left[\frac{1}{x^{2}-1}-\frac{1}{x^{2}+1}\right] \mathrm{d} x\\ =&\frac{1}{4} \ln \left|\frac{x-1}{x+1}\right|-\frac{1}{2} \arctan x+C \end{aligned} ==x41dx21[x211x2+11]dx41lnx+1x121arctanx+C

结语

至此, 有理函数积分的范式就介绍完了. 但是切忌盲目使用此方法, 老实说这种范式更加适合计算机——没有任何技巧性但是运算量较大. 更多时候积分的题目还是需要运用一些技巧灵活应对, 不过当你手足无措的时候, 不要忘了还有这样一套办法.

祝你好运!

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐