用Python学《微积分B》（换元法与分部积分）

一、知识点

1，换元法
  换元法主要是根据“积分形式不变性”，通过“替换被积变量”的方式将“被积函数”转换为“积分表”中的函数及其组合，待求出原函数之后，再将变量替换回来的求积分的方法。
1）复合函数换元法（第一类换元法）
  复合函数换元法主要是利用复合函数微分来进行“凑微分”，即将“微分符”外面的东西变到“微分符”里面。其关键点是要找到合适的中间变量$u = u{\left (x \right )}$ ，使得 ：

∫f(x)dx=∫g(u)du=G(u)+C=G[u(x)]+C

$$\int f{\left (x \right )}\, dx = \int g{\left (u \right )}\, du = G{\left (u \right )} + C = G{\left [u{\left (x \right )} \right ]} + C$$
例：
$$\int tan{\left (x \right )}\, dx = \int \frac{sin(x)}{cos(x)}dx=-\int \frac{d[cos(x)]}{cos(x)}=-ln|cos(x)|+C$$
2）反函数换元法（第二类换元法）
  反函数换元法与复合函数换元法刚好相反，它需要做的是将“微分符”里面的东西变到“微分符”外面。其关键点是找到新变量t，使得$x=\varphi (t)$，积分过程如下：
$$\int f(x)dx=\int f[\varphi (t)]d[\varphi (t)]=\int f[\varphi (t)]*\varphi '(t)dt=G[\varphi ^{-1}(x)]+C$$
  最后又要用反函数将x替换回来
例：
$$\int \frac{1}{\sqrt{x^{2}-a^{2}}}dx=?\; when\; a>0$$
解：令$sec(t)=\frac{x}{a}$，则$tan(t)=\sqrt{sec^{2}(t)-1}=\frac{1}{a}\sqrt{x^{2}-a^{2}}$，代入上式得
$$\int \frac{1}{\sqrt{x^{2}-a^{2}}}dx=ln|\frac{x}{a}+\frac{\sqrt{x^{2}-a^{2}}}{a}|+C$$

1，分部积分法
  分部积分法是从“函数乘积的微分”，即$d(u*v)=udv+vdu$，推导出来的一种求积分的方法。如果说“换元法”是将微分符外面的东西变到微分符里面，或反过来；那么“分部积分法”则是直接交换微分符里面和外面的东西。

∫udv=uv−∫vdu

$$\int udv=uv-\int vdu$$ 或 $$\int u(x)v'(x)dx=u(x)v(x)-\int v(x)u'(x)dx$$
  分部积分法主要是将“求函数的积分转变为求导函数的积分”（又是“降维”的方法），所以主要应用于这两种情况：一是函数的比较复杂，而它的导数比较简单，比如：ln x、arctan(x)、arcsin(x)等；二是函数与导数差不多的情况，比如：$e^{x}$、sin(x)、cos(x)等。
例：
$$\int x^{2}ln(x)dx= \frac{x^{3}}{3}*ln(x)-\int \frac{x^{3}}{3}d(ln(x))=\frac{x^{3}}{3}*ln(x)-\int \frac{x^{3}}{3}*\frac{1}{x}dx$$
$$\int x^{2}e^{-2x}dx=\int x^{2}d[\frac{e^{-2x}}{-2}]=-\frac{1}{2}x^{2}e^{-2x}+\int xe^{-2x}dx$$
  从上面这两个例题的推导过程来看，分部积分法实际上是一种“递归”（recursion）算法。
$$\int e^{x}sin(x)dx = e^{x}sin(x) - \int e^{x}cos(x)dx= e^{x}sin(x) - e^{x}cos(x) - \int e^{x}sin(x)dx$$
  上式是一种特殊的递归，通过两次分部积分，又得到了自身，故，它的递归结束条件是重回自身。

二、课后习题

换元法：
3，求积分$\int e^{- \left|{x}\right|}\, dx$
解：包含绝对值的积分式，不能用python直接解。可以将它分段表示，再用“凑微分法”求解，最后取x=0点连续，求出共同的常数C。

∫e−|x|dx={−e−x+Cex−2+Cx≥0x<0

$$\int e^{- \left|{x}\right|}\, dx=\left\{\begin{matrix} -e^{-x}+C & x\geq 0\\ e^{x}-2+C&x<0 \end{matrix}\right.$$
下面再给出5-10题的python求解过程（注，对于“换元”，sympy.subs函数给与了很好的支持）：

from sympy import *
init_printing()
#Exercise 3
x = Symbol('x')
integrate(exp(-abs(x)), x), integrate(exp(-x), x), integrate(exp(x), x)

$$\left ( \int e^{- \left|{x}\right|}\, dx, \quad - e^{- x}, \quad e^{x}\right )$$

#Exercise 5
x = Symbol('x')
f = log(x) / x ** 2
f, integrate(f, x)

$$\left ( \frac{1}{x^{2}} \log{\left (x \right )}, \quad - \frac{1}{x} \log{\left (x \right )} - \frac{1}{x}\right )$$

#Exercise 6
x = Symbol('x')
Df = exp(-x).diff(x)
g = Df.subs(x, log(x))
g, integrate(g / x, x)

$$\left ( - \frac{1}{x}, \quad \frac{1}{x}\right )$$

#Exercise 7
x = Symbol('x')
F = x ** 2
f = F.diff(x)
g = f.subs(x, 1 - x ** 2)
f, g, integrate(x * g, x)

$$\left ( 2 x, \quad - 2 x^{2} + 2, \quad - \frac{x^{4}}{2} + x^{2}\right )$$

#Exercise 8
x = Symbol('x')
f = x / (1 - x) ** 3
f, integrate(f, x)

$$\left ( \frac{x}{\left(- x + 1\right)^{3}}, \quad \frac{2 x - 1}{2 x^{2} - 4 x + 2}\right )$$

#Exercise 9
x, u = symbols('x u')
f = u / (1 + u ** 4)
Gu = integrate(f, u)
f, Gu, Gu.subs(u, sin(x))

$$\left ( \frac{u}{u^{4} + 1}, \quad \frac{1}{2} \operatorname{atan}{\left (u^{2} \right )}, \quad \frac{1}{2} \operatorname{atan}{\left (\sin^{2}{\left (x \right )} \right )}\right )$$

#Exercise 10
x = Symbol('x')
f = 2 ** x + x ** 2
f, integrate(f, x)

$$\left ( 2^{x} + x^{2}, \quad \frac{2^{x}}{\log{\left (2 \right )}} + \frac{x^{3}}{3}\right )$$