微积分表

微分表

积分表

`sinh x = ("e"^x - "e"^-x)/2` `quad cosh x = ("e"^x + "e"^-x)/2`, `quad tanh x = ("e"^x - "e"^-x)/("e"^x + "e"^-x)`;
`"arsh "x = ln (x+sqrt(x^2+1))`, `quad x in RR`,
`"arch "x = ln (x+sqrt(x^2-1))`, `quad x ge 1`,
`"arth "x = 1/2 ln ((1+x)/(1-x))`, `quad -1 lt x lt 1`.

级数与无穷乘积

幂级数

设 `|x| lt 1`, 则
1. `sum_(n ge 0) x^n = 1/(1-x)` (几何级数);
2. `sum_(n ge 1) x^n/n = -ln(1-x)`;
3. `sum_(n ge 1, "odd") x^n/n = "arth "x = 1/2 ln|(1+x)/(1-x)|`;

1. 由等式 `(1-x)(1 + x + cdots + x^n) = 1 - x^(n+1)` 得 `1 + x + cdots + x^n = (1-x^(n+1))/(1-x)`, 两边令 `n to oo` 即得结论.
2. 由几何级数两边从 `0` 到 `x` 积分即可.
3. 由 `sum_(n ge 0) t^(2n) = 1/(1-t^2)` 两边从 `0` 到 `x` 积分即可.

1. 在 `1/sqrt(1-x^2) = sum_(n ge 0) (-1/2; n) (-x^2)^n` `= sum_(n ge 0) ((-1/2)(-3/2) cdots (-(2n-1)/2))/(n!) (-x^2)^n` `= sum_(n ge 0) ((2n-1)!!)/((2n)!!) x^(2n)` 两边逐项积分, 就得到 `arcsin x` 的展式.
2. 将 `"erf"(x) = 2/sqrt pi int_0^x "e"^(-t^2) dt` 的被积函数展开并逐项积分, 就得到 `"erf"(x)` 的展式.
3. 用 sympy 可以得到更多级数表达式, 如 series((1+x)**(1/x), x) 输出:

   E - E*x/2 + 11*E*x**2/24 - 7*E*x**3/16 + 2447*E*x**4/5760 - 959*E*x**5/2304 + O(x**6)
   

1. 双阶乘. `n` 为正整数时, `(2n)!! = 2 * 4 * 6 * cdots * 2n`, `quad (2n-1)!! = 1 * 3 * 5 * cdots * (2n-1)`. 特别约定 `0!! = (-1)!! = 1`.
2. `B_(2n)` 是 Bernoulli 数. 满足 `B_0 = 1, B_1 = -1/2, B_3 = B_5 = B_7 = cdots = 0`. 事实上由生成函数知道 `x = ("e"^x-1) sum_(k ge 0) B_k x^k/(k!)` `= (sum_(j ge 1) x^j/(j!)) (sum_(k ge 0) B_k x^k/(k!))` `= sum_(n ge 1) x^n/(n!) sum_(k=0)^(n-1) B_k (n; k)`. 等号两边比较系数得 `sum_(k=0)^n (n+1; k) B_k = 0`, `quad n ge 1`.
3. `E_(2n)` 是 Euler 数, 它是正负号交错的整数: `E_0 = 1`, `E_2 = -1`, `E_4 = 5`, `E_6 = -61`, `E_8 = 1385`, `E_10 = -50521`. 由 `E_(2n)` 的生成函数知道, `1 = cosh x sum_(k ge 0) E_(2k)/((2k)!) x^(2k)` `= (sum_(j ge 0) x^(2j)/((2j)!)) (sum_(k ge 0) E_(2k)/((2k)!) x^(2k))` `= sum_(n ge 0) x^(2n)/((2n)!) sum_(k=0)^n (2n;2k) E_(2k)` 比较等式两边的各次项系数, 得到递推公式: `sum_(k=0)^n (2n;2k) E_(2k) = 0`, `quad n ge 1`.
4. 三角函数与其对应的双曲三角函数有简单联系 `cos x = cosh("i"x)`, `quad cosh x = cos("i"x)`,
   `sin x = -"i" sinh("i"x)`, `quad sinh x = -"i" sin("i"x)`,
   `tan x = -"i" tanh("i"x)`, `quad tanh x = -"i" tan("i"x)`,
   `cot x = "i" coth("i"x)`, `quad coth x = "i" cot("i"x)`,
   `arccos x = "arch "x("i"x)`, `quad "arch "x = arccos("i"x)`,
   `arcsin x = -"i arsh "x("i"x)`, `quad "arsh "x = -"i"arcsin("i"x)`,
   `arctan x = -"i arth "x("i"x)`, `quad "arth "x = -"i"arctan("i"x)`, 因此只需将一个三角函数的级数从第二项起交错改变符号, 就得到对应的双曲三角函数的级数, 反之亦然. 如 `cos x = sum_(n ge 0) (-1)^n x^(2n)/((2n)!)`,
   `"arsh "x = sum_(n ge 0) (-1)^n ((2n-1)!!)/((2n)!!) x^(2n+1)/(2n+1)`,
   `cot x = 1/x sum_(n ge 0) (-1)^n B_(2n)/((2n)!) (2x)^(2n)`.

[题源 brilliant.org] 求 `(sum_(n=1)^oo n/((2n+1)!))/(sum_(n=1)^oo n/((2n-1)!))`.

将求和写成 `(2/(3!) + 4/(5!) + 6/(7!) + cdots)/ (2/(1!) + 4/(3!) + 6/(5!) + cdots)` 分子等于 `(3-1)/(3!) + (5-1)/(5!) + (7-1)/(7!) + cdots` `= 1/(2!) - 1/(3!) + 1/(4!) - 1/(5!) + 1/(6!) - 1/(7!) + cdots` `= "e"^-1 - 1 + 1/(1!)` `= "e"^-1`. 类似可得分母等于 `"e"`, 于是结果为 `"e"^-2`.

由 digamma 函数的展式得 `sum_(n ge 1)(1/n - 1/(n+s)) = psi(s+1) + gamma`.

从 Euler 的正弦无穷乘积 `prod_(n ge 1) (1-x^2/n^2) = (sin pi x)/(pi x)` (参见 Gamma 函数一节) 出发, 可以得到一些推论:
1. 在正弦乘积公式中令 `x = "i" y`, 就得到双曲正弦的公式 `prod_(n ge 1) (1 + y^2/n^2) = (sinh pi y)/(pi y)`.
2. 在正弦乘积公式两边取对数再求导, 整理得到余切函数的展开: `sum_(n ge 1) 1/(n^2 - x^2) = 1/(2x^2) (1-pi x cot pi x)`. 同样令 `x = "i" y` 得到 `sum_(n ge 1) 1/(n^2+y^2) = 1/(2y^2)(pi y coth pi y - 1)`.
3. 一个与上式有几分神似的求和 (可以通过裂项得到): `sum_(n ge 2) 1/(n^2-1) = 3/4`. 事实上, 这个结果也可以通过求极限得到: `lim_(x to 1) ((1-pi x cot pi x)/(2x^2) - 1/(1-x^2)) = 3/4`.
4. 在余切函数的展式中令 `z = pi x`, 稍加整理可以得到更漂亮的公式 `sum_(n=-oo)^oo 1/(z-n pi) = cot z`, 这里的求和取主值, 即 `1/z + sum_(n ge 1) (1/(z+n pi) + 1/(z-n pi))`. 两边求导立得 Jumping 恒等式: `sum_(n=-oo)^oo 1/(z-pi n)^2 = 1/(sin^2 z)`.

求 `f(x) = sum_(n ge 0) x^n/(n!!)` (`0!! = 1`).

`f(0) = 1`, 且 `f` 适合微分方程 `y' = x y + 1`. 解得 `f(x) = sqrt(2/pi) int_(-oo)^x "e"^((x^2-t^2)/2) dt` `= "e"^(x^2/2) ("erf"(sqrt 2 x)+1)`. 利用 Wallis 公式知道, 收敛半径为 `oo`.

Fourier 级数

设 `x in (0, 2pi)`, 则
1. `sum_(n ge 1) (sin n x)/n = (pi-x)/2`;
2. `sum_(n ge 1) (cos n x)/n^2 = x^2/4 - (pi x)/2 + pi^2/6`.
3. `sum_(n ge 1) (sin n x)/n^3 = x^3/12 - (pi x^2)/4 + (pi^2 x)/6`.
在 1. 中替换 `x |-> pi-x`, 就得到 `x` 在 `|x| lt pi` 上的 Fourier 展开. 类似地, 由 2, 3 分别得 `x^2, x^3` 的 Fourier 展开.

1. 直接计算 `(pi-x)//2` 的 Fourier 级数即可.
2. 上题结果在 `[0, x]` 上逐项积分得 `sum_(n ge 1) (1-cos n x)/n^2 = pi/2 x - x^2/4`, 再利用 `sum_(n ge 1) n^-2 = pi^2//6` 即得结论.
3. 由上题结果在 `[0, x]` 上积分即得.

类似上题有
1. `sum_(n ge 1) (cos n x)/n = -ln(2 sin{:x/2:})`;
2. `sum_(n ge 1) (sin n x)/n^2 = -int_pi^x ln(2 sin{:t/2:}) dt`;
   
令 `x = pi-x`, 又得到
3. `sum_(n ge 1) (-1)^n (cos n x)/n = -ln(2 cos{:x/2:})`;
4. `sum_(n ge 1) (-1)^(n-1) (sin n x)/n^2 = int_0^x ln(2 cos{:t/2:}) dt`;
   `x = 0, pi` 时, 等式两边等于 `0`; `x = pi/2` 时, 等于 Calatan 常数 `K`.

一般地, `f` 在 `(-pi,pi)` 上的 Fourier 级数可以由以下公式计算: `f(x) = a_0/2 + sum_(n ge 1) (a_n cos n x + b_n sin n x)`,
`a_n = 1/pi int_(-pi)^pi f(x) cos n x dx`, `quad b_n = 1/pi int_(-pi)^pi f(x) sin n x dx`.

提示: 1. `x sin x` 和 `x cos x` 积分时需借助积化和差; 2. `"e"^(a x) = sinh a x + cosh a x`.

在 `x` 的展式中取 `x = pi/2`, 得到 `pi//4 = sum_(n ge 1) (-1)^(n-1)/(2n-1)`; 在 `x^2` 的展式中取 `x = pi`, 得到 `sum_(n ge 1) n^-2 = pi^2//6`.

积分变换

详细的讨论参见 积分变换.

Laplace 变换

用 Gamma 函数的定义验证表中的第一条公式; 后面几条都是它的推论.

1. `cc L[delta(t)]`. 严格说, `delta(t)` 是广义函数, 我们暂不追究详情; 直观上它是下面函数族当 `epsi to 0` 的极限: `f_epsi(t) = { 1//epsi, if |t| lt epsi; 0, otherwise :}` 于是 `cc L[delta(t)]` `= lim_(epsi to 0) int_0^epsi 1/epsi "e"^(-s t) dt` `= lim_(epsi to 0) 1/(epsi s) (1-"e"^(-epsi s)) = 1`.
2. `cc L[t^(a-1) ln^n t]`. 由 `cc L[t^(a-1)] = (Gamma(a))/s^a` 两边对 `a` 求 `n` 阶导即得结论. 特别 `a = n = 1` 时 `cc L[ln t] = -(gamma + ln s)/s`

Mellin 变换

Fourier 变换

连分数

Euler 连分数公式 `sum_(i=0)^n prod_(j=0)^i a_j` `= a_0/(1-) a_1/(1+a_1-) a_2/(1+a_2-) cdots` `a_(n-1)/(1+a_(n-1)-) a_n/(1+a_n)`.

记 `s_k = sum_(i=k)^n prod_(j=k)^i a_j`, 则 `s_k = sum_(i=k)^n a_k prod_(j=k+1)^i a_j` `= a_k(1 + sum_(i=k+1)^n prod_(j=k+1)^i a_j)` `= a_k(1+s_(k+1))`, `k = 0, 1, cdots, n-1`. 因此 `s_0 = a_0 (1+s_1)` `= a_0/(1-s_1//(1+s_1))`,
`s_k/(1+s_k)` `= a_k/(a_k + a_k//s_k)` `= a_k/(a_k + 1//(1+s_(k+1)))` `= a_k/(a_k + 1 - s_(k+1)//(1+s_(k+1)))`. 连续运用上式即得原式成立.

级数化为连分数 用 `u_i` 替换 Euler 连分数中的乘积 `prod_(j=0)^i a_j`, 得到 `sum_(i=0)^n u_i` `= a_0/(1-) a_1/(b_1-) a_2/(b_2-) cdots a_n/b_n`, 其中 `a_0 = u_0`, `b_0 = 1`, `a_n = u_n//u_(n-1)`, `quad b_n = 1 + a_n`, `quad n ge 1`. 故一个各项不为零的级数, 形式上可以表示为连分数, 使得其第 `n` 个部分和等于连分数的第 `n` 项近似.

无穷乘积化为连分数 利用下式将无穷乘积写为级数: `prod_(n ge 1) (1 + a_n/b_n)` `= 1 + a_1/b_1 + a_2/b_2(1+a_1/b_1) + a_3/b_3(1+a_1/b_1)(1+a_2/b_2) + cdots` 再将级数化为连分数, 得 `1/(1-) (a_1//b_1)/(1+a_1//b_1-)` `(a_2//b_2 (1+b_1//a_1))/(1+a_2//b_2(1+b_1//a_1)-)` `(a_3//b_3 (1+b_2//a_2))/(1+a_3//b_3(1+b_2//a_2)-) cdots`.

利用 Taylor 展开,
1. 指数函数 `"e"^x = 1 + 1 * x + 1 * x * x^2/2 + cdots`
   `= 1/(1-) x/(1+x-) (x//2)/(1+x//2-) (x//3)/(1+x//3-) cdots`
   `= 1/(1-) x/(1+x-) x/(2+x-) (2x)/(3+x-) (3x)/(4+x-) cdots`, 特别 `"e" = 1/(1-) 1/(2-) 1/(3-) 2/(4-) 3/(5-) cdots`. 上面的连分数只是幂级数的另一种写法, 其各项截断等于幂级数的各个部分和, 因此其收敛性由幂级数的收敛性保证. 和 `"e"` 的无穷级数一样, 该连分数收敛很快, 取 n = 12 就得到 `"e"` 的近似值 2.718281828703704.
2. 正弦函数 `sin x = x + x * (-x^2)/3 + x * (-x^2)/3 * (-x^2)/5 + cdots`
   `= x/(1-) ((-x^2)//(2*3))/(1-x^2//(2*3)-) ((-x^2)//(4*5))/(1-x^2//(4*5)-) cdots`
   `= x/(1+) x^2/(2*3-x^2+) (2*3 x^2)/(4*5-x^2+) cdots`
3. 反正切函数 `arctan x = x + x * (-x^2)/3 + x * (-x^2)/3 * (-3x^2)/5 + cdots`
   `= x/(1-) (-x^2//3)/(1-x^2//3-) (-3x^2//5)/(1-3x^2//5-) cdots`
   `= x/(1+) x^2/(3-x^2+) (3x)^2/(5-3x^2+) (5x)^2/(7-5x^2+) cdots`, 特别 `pi/4 = 1/(1+) 1^2/(2+) 3^2/(2+) 5^2/(2+) cdots`. 和 `arctan` 的展式一样, 上面的连分数收敛很慢. 取 n = 999, 只得到 `pi` 的近似值 3.143584671557847.

连分数的第 `n` 项近似 记 `b_1/(a_1+) b_2/(a_2+) cdots b_n/a_n = y_n/x_n` (不约分), 则分子和分母均满足二阶线性递推关系: `x_0 = 1`, `quad x_1 = a_1`, `quad x_n = a_n x_(n-1) + b_n x_(n-2)`,
`y_0 = 0`, `quad y_1 = b_1`, `quad y_n = a_n y_(n-1) + b_n y_(n-2)`. 上式给出连分数第 `n` 项近似的算法. 反之, 任意两个满足上面递推关系的数列, 补充定义 `x_0, y_0` 后, 它们商可以展开为连分数. 例如 Fibonacci 数列满足 `F_0 = 0`, `F_1 = 1`, `F_n = F_(n-1) + F_(n-2)`, 因此 `F_n/F_(n+1)` `= F_n/(F_n + F_(n-1))` `= 1/(1+F_(n-1)//F_n)` `= cdots` `= 1/(1+) 1/(1+) cdots 1/1`.

容易验证边界条件成立. 假设成立 `y_n/x_n = (a_n y_(n-1) + b_n y_(n-2))/(a_n x_(n-1) + b_n x_(n-2))`, 用 `a_n + b_(n+1)//a_(n+1)` 替换公式中的 `a_n` 得 `y_(n+1)/x_(n+1)` `= ((a_n + b_(n+1)//a_(n+1)) y_(n-1) + y_n - a_n y_(n-1)) / ((a_n + b_(n+1)//a_(n+1)) x_(n-1) + x_n - a_n x_(n-1))` `= (a_(n+1) y_n + b_(n+1) y_(n-1))/(a_(n+1) x_n + b_(n+1) x_(n-1))`.

考虑超几何函数 `{::}_2 F_1(a,b";"c";"z)` `= 1 + (a b)/c z + (a(a+1)b(b+1))/(c(c+1)) z^2/(2!) + cdots`, 取 `f_(2n) = {::}_2 F_1(a+n,b+n";"c+2n-1";"z)`,
`f_(2n+1) = {::}_2 F_1(a+n+1,b+n";"c+2n";"z)`,
`n = 0, 1, 2, cdots`, 则 `f_n = f_(n+1) + k_n z f_(n+2)`, 其中 `k_(2n) = ((a-c-n+1)(b+n))/((c+2n-1)(c+2n))`, `quad k_(2n+1) = ((b-c-n-1)(a+n+1))/((c+2n)(c+2n+1))`. 于是 `f_0/f_1` 可以展开为 `f_0/(f_1)` `= (f_1 + k_0 z f_2)/(f_1)` `= 1 + (k_0 z)/(f_1 // f_2)` `= cdots` `= 1 + (k_0 z)/(1+) (k_1 z)/(1+) cdots`. 特别取 `a = 0`, 则 `f_1 = 1`, 我们就得到 `f_0` 的连分数. 鉴于许多初等函数都能写成超几何函数, 如 `arctan z = z {::}_2 F_1(1,1/2";"3/2";"-z^2)`. 我们得到一些初等函数的连分数:

`tanh x` 的连分数的另一证明: 记 `a_(n,k) = 1//(2n)!(2n+1)(2n+3)cdots(2n+2k-1)`, 则 `a_(n,k-1) - (2k-1) a_(n,k)` `= { a_(n-1,k-1), if n ge 1; 0, if n = 0 :}` 进而有 `s_k - (2k-1)` `:= (sum_(n ge 0) a_(n,k-1) x^(2n))/(sum_(n ge 0) a_(n,k) x^(2n)) - (2k-1)` `= (sum_(n ge 1) a_(n-1,k+1)x^(2n))/(sum_(n ge 0) a_(n,k) x^(2n))` `= x^2(sum_(n ge 0) a_(n,k+1)x^(2n))/(sum_(n ge 0) a_(n,k) x^(2n))` `= x^2 // s_(k+1)`, 即 `s_k = (2k-1) + x^2//s_(k+1)`. 将 `sinh x` 和 `cosh x` 的 Taylor 展式相除, `tanh x = x (1 + x^2/(3!) + x^4/(5!) + cdots)/(1 + x^2/(2!) + x^4/(4!) + cdots)` `= x (sum_(n ge 0) a_(n,1) x^(2n))/(sum_(n ge 0) a_(n,0) x^(2n))` `= x // s_1`, 递推得 `tanh x = x/(1+) x^2/(3+) x^2/(5+) x^2/(7+) cdots`.