昨天在研究一个最优化问题时,遇到了一个这样的积分:
$$\int \frac{1}{\cos^3 \theta} d\theta$$

然后就顺便研究了一下这种类型的函数的积分。一般来讲,这类积分可以写成$\int cos^n \theta d\theta$或$\int sin^n \theta d\theta$,其中n是一个整数。

首先我们来解决n=1的情况,我们很容易就有$\int cos\theta d\theta=sin\theta +C$或$\int sin\theta d\theta=-cos\theta +C$,这是一个基本的结果。

如果n是大于1的正整数,那么可以用递推的方法来搞定:
$$\begin{aligned}\int \cos^n \theta d\theta=\int \cos^{n-1} \theta d(\sin\theta)=\cos^{n-1} \theta \sin\theta-\int \sin\theta d(\cos^{n-1}\theta) \\ =\cos^{n-1} \theta \sin\theta+(n-1)\int \sin^2\theta \cos^{n-2}\theta d\theta \\ =\cos^{n-1} \theta \sin\theta+(n-1)\int (1-\cos^2\theta) \cos^{n-2}\theta d\theta \\ =\cos^{n-1} \theta \sin\theta+(n-1)\int \cos^{n-2}\theta d\theta-(n-1)\int \cos^n\theta d\theta\end{aligned}$$
(这里运用了分部积分法)

即$n\int cos^n \theta d\theta=cos^{n-1} \theta sin\theta+(n-1)\int cos^{n-2}\theta d\theta$
$$\int \cos^n \theta d\theta=\frac{1}{n}[\cos^{n-1} \theta \sin\theta+(n-1)\int \cos^{n-2}\theta d\theta]$$
这样就完成了递推。关于正弦也有类似公式,这里直接写出:
$$\int \sin^n \theta d\theta=\frac{1}{n}[-\sin^{n-1} \theta \cos\theta+(n-1)\int \sin^{n-2}\theta d\theta]$$

不过,我昨晚想到将其变为多项式积分来算:
$$\int \cos^n \theta d\theta=\int \cos^{n-1} \theta d(\sin\theta)$$
如果令$sin\theta=x$,则积分成了
$$\int (1-x^2)^{\frac{n-1}{2}} dx$$

要是n是奇数(正或负都行)的话,这个积分是相当好办的。要是它是偶数或其它有理数的话,就需要用到超几何函数(Hypergeometric Function)来表达。这里只对n是奇数进行讨论。

当n是正数时,直接多项式展开进行积分,不再讨论。当n是负数时,相当于讨论积分
$$\int \frac{1}{(1-x^2)^k} dx$$
其中k是正整数。当k=1时,有
$$\int \frac{1}{1-x^2} dx=\frac{1}{2} ln|\frac{1+x}{1-x}|$$

当k > 1时,我们即$\int \frac{1}{(1-x^2)^k} dx=I_k$,有
$$\begin{aligned}I_k=\int \frac{1-x^2+x^2}{(1-x^2)^k} dx=I_{k-1}+\int \frac{x^2}{(1-x^2)^k} dx \\ =I_{k-1}-\frac{1}{2}\int \frac{x}{(1-x^2)^k} d(1-x^2) \\ =I_{k-1}+\frac{1}{2(k-1)}\int x d[(1-x^2)^{1-k}] \\ =I_{k-1}+\frac{1}{2(k-1)}[x (1-x^2)^{1-k}-\int (1-x^2)^{1-k}dx] \\ =(\frac{2k-3}{2k-2})I_{k-1}+\frac{x}{2(k-1)(1-x^2)^{k-1}}\end{aligned}$$
递归完毕。

这里直接写出我可能用到的结果:
$$\int \frac{1}{\cos^3 \theta} d\theta=\frac{1}{4}(\frac{2\sin \theta}{\cos^2 \theta}+ln|\frac{1+\sin \theta}{1-\sin \theta}|)$$
$\int \frac{1}{sin^3 \theta} d\theta=-\frac{1}{4}(\frac{2cos \theta}{sin^2 \theta}+ln|\frac{1+cos \theta}{1-cos \theta}|)$

转载到请包括本文地址:https://kexue.fm/archives/1422

更详细的转载事宜请参考:《科学空间FAQ》

如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。

如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!

如果您需要引用本文,请参考:

苏剑林. (Jul. 22, 2011). 《三角函数幂的积分 》[Blog post]. Retrieved from https://kexue.fm/archives/1422

@online{kexuefm-1422,
        title={三角函数幂的积分},
        author={苏剑林},
        year={2011},
        month={Jul},
        url={\url{https://kexue.fm/archives/1422}},
}