莱布尼茨公式 (Leibniz formula) 指的是几种不同的公式,最常见的是指计算 π 的莱布尼茨公式,也称为莱布尼茨级数或莱布尼茨-格雷戈里级数 (Leibniz–Gregory series)。 它是一个无限级数,其和等于 π/4。

1. 计算 π 的莱布尼茨公式:

这个公式表示为:

π4=∑n=0∞(−1)n2n+1=1−13+15−17+19−⋯ \frac{\pi}{4} = \sum_{n=0}^{\infty} \frac{(-1)^n}{2n+1} = 1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \frac{1}{9} - \cdots 4π​=n=0∑∞​2n+1(−1)n​=1−31​+51​−71​+91​−⋯

推导过程 (简略): 这个公式的推导通常依赖于反正切函数的泰勒展开式。 反正切函数的泰勒级数展开式为:

arctan⁡(x)=∑n=0∞(−1)nx2n+12n+1=x−x33+x55−x77+⋯ \arctan(x) = \sum_{n=0}^{\infty} \frac{(-1)^n x^{2n+1}}{2n+1} = x - \frac{x^3}{3} + \frac{x^5}{5} - \frac{x^7}{7} + \cdots arctan(x)=n=0∑∞​2n+1(−1)nx2n+1​=x−3x3​+5x5​−7x7​+⋯

该级数在 |x| ≤ 1 时收敛。 当 x = 1 时,我们得到:

arctan⁡(1)=∑n=0∞(−1)n2n+1=1−13+15−17+19−⋯ \arctan(1) = \sum_{n=0}^{\infty} \frac{(-1)^n}{2n+1} = 1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \frac{1}{9} - \cdots arctan(1)=n=0∑∞​2n+1(−1)n​=1−31​+51​−71​+91​−⋯

由于 arctan(1) = π/4,因此我们就得到了莱布尼茨公式。

2. 莱布尼茨公式的其他形式:

除了计算 π 的公式,莱布尼茨还提出了其他一些公式,例如:

  • 微分法则 (莱布尼茨法则): 这可能是莱布尼茨最著名的贡献之一,它描述了两个函数乘积的导数:

ddx(f(x)g(x))=f′(x)g(x)+f(x)g′(x) \frac{d}{dx}(f(x)g(x)) = f'(x)g(x) + f(x)g'(x) dxd​(f(x)g(x))=f′(x)g(x)+f(x)g′(x)

这个公式可以推广到 n 个函数的乘积,以及高阶导数。

  • 莱布尼茨积分法则: 这个公式描述了参数积分的导数:

ddx∫a(x)b(x)f(x,t)dt=f(x,b(x))dbdx−f(x,a(x))dadx+∫a(x)b(x)∂∂xf(x,t)dt \frac{d}{dx} \int_{a(x)}^{b(x)} f(x,t) dt = f(x, b(x)) \frac{db}{dx} - f(x, a(x)) \frac{da}{dx} + \int_{a(x)}^{b(x)} \frac{\partial}{\partial x} f(x,t) dt dxd​∫a(x)b(x)​f(x,t)dt=f(x,b(x))dxdb​−f(x,a(x))dxda​+∫a(x)b(x)​∂x∂​f(x,t)dt

这个公式在微积分学中非常重要,它允许我们对积分进行微分运算。

3. 莱布尼茨公式的收敛速度:

需要注意的是,莱布尼茨公式计算 π 的收敛速度非常慢。 要得到 π 的一个比较精确的近似值,需要计算大量的项。 这是因为该级数是一个交错级数,其收敛速度受制于调和级数的收敛速度。 因此,在实际应用中,通常使用其他更快速的算法来计算 π。

让我们分别用几个例子来说明莱布尼茨的不同公式:

例1:使用莱布尼茨级数近似计算 π

用莱布尼茨公式计算 π 的近似值,取前五项:

π4≈1−13+15−17+19 \frac{\pi}{4} \approx 1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \frac{1}{9} 4π​≈1−31​+51​−71​+91​

计算结果:

π4≈1−0.3333+0.2−0.1429+0.1111=0.8349 \frac{\pi}{4} \approx 1 - 0.3333 + 0.2 - 0.1429 + 0.1111 = 0.8349 4π​≈1−0.3333+0.2−0.1429+0.1111=0.8349

π≈4×0.8349=3.3396 \pi \approx 4 \times 0.8349 = 3.3396 π≈4×0.8349=3.3396

这个结果非常粗略,因为我们只取了前五项。 要得到更高的精度,需要计算更多的项,但这会非常耗时,这也是莱布尼茨公式的一个缺点。

例2:应用莱布尼茨微分法则

求函数 f(x)=x2sin⁡(x)f(x) = x^2 \sin(x)f(x)=x2sin(x) 的导数。

我们可以使用莱布尼茨微分法则:

ddx(f(x)g(x))=f′(x)g(x)+f(x)g′(x) \frac{d}{dx}(f(x)g(x)) = f'(x)g(x) + f(x)g'(x) dxd​(f(x)g(x))=f′(x)g(x)+f(x)g′(x)

在这个例子中, f(x)=x2f(x) = x^2f(x)=x2 且 g(x)=sin⁡(x)g(x) = \sin(x)g(x)=sin(x)。 因此:

f′(x)=2x f'(x) = 2x f′(x)=2x

g′(x)=cos⁡(x) g'(x) = \cos(x) g′(x)=cos(x)

应用莱布尼茨法则:

ddx(x2sin⁡(x))=(2x)sin⁡(x)+x2cos⁡(x)=2xsin⁡(x)+x2cos⁡(x) \frac{d}{dx}(x^2 \sin(x)) = (2x)\sin(x) + x^2\cos(x) = 2x\sin(x) + x^2\cos(x) dxd​(x2sin(x))=(2x)sin(x)+x2cos(x)=2xsin(x)+x2cos(x)

例3:应用莱布尼茨积分法则 (一个简化例子)

求 ddx∫0xt2dt \frac{d}{dx} \int_0^x t^2 dt dxd​∫0x​t2dt

这里, a(x)=0a(x) = 0a(x)=0, b(x)=xb(x) = xb(x)=x, f(x,t)=t2f(x, t) = t^2f(x,t)=t2。 根据莱布尼茨积分法则:

ddx∫0xt2dt=f(x,x)ddx(x)−f(x,0)ddx(0)+∫0x∂∂xt2dt \frac{d}{dx} \int_0^x t^2 dt = f(x, x) \frac{d}{dx}(x) - f(x, 0) \frac{d}{dx}(0) + \int_0^x \frac{\partial}{\partial x} t^2 dt dxd​∫0x​t2dt=f(x,x)dxd​(x)−f(x,0)dxd​(0)+∫0x​∂x∂​t2dt

由于 ∂∂xt2=0\frac{\partial}{\partial x} t^2 = 0∂x∂​t2=0 (因为 t² 不依赖于 x),并且 ddx(0)=0\frac{d}{dx}(0) = 0dxd​(0)=0,所以:

ddx∫0xt2dt=x2(1)−0+0=x2 \frac{d}{dx} \int_0^x t^2 dt = x^2(1) - 0 + 0 = x^2 dxd​∫0x​t2dt=x2(1)−0+0=x2

这个结果也可以通过先计算积分,再求导来验证:

∫0xt2dt=x33 \int_0^x t^2 dt = \frac{x^3}{3} ∫0x​t2dt=3x3​

ddx(x33)=x2 \frac{d}{dx} (\frac{x^3}{3}) = x^2 dxd​(3x3​)=x2

例4:莱布尼茨级数的应用(误差分析)

我们想用莱布尼茨公式计算π的近似值,并估计误差。 假设我们计算到第N项,那么误差的绝对值小于等于第N+1项的绝对值。

例如,我们计算到第1000项:

π4≈∑n=01000(−1)n2n+1 \frac{\pi}{4} \approx \sum_{n=0}^{1000} \frac{(-1)^n}{2n+1} 4π​≈n=0∑1000​2n+1(−1)n​

误差的绝对值小于:

∣(−1)10012(1001)+1∣=12003≈0.000499 \left| \frac{(-1)^{1001}}{2(1001)+1} \right| = \frac{1}{2003} \approx 0.000499 ​2(1001)+1(−1)1001​​=20031​≈0.000499

因此,使用1000项计算出的π/4的误差小于0.0005,π的误差小于0.002。 这说明即使计算了大量项,莱布尼茨公式的收敛速度仍然很慢。

例5:莱布尼茨微分法则的高阶导数

求函数 f(x)=x3exf(x) = x^3 e^xf(x)=x3ex 的三阶导数。 我们可以反复使用莱布尼茨法则:

首先,令 g(x)=x3g(x) = x^3g(x)=x3 和 h(x)=exh(x) = e^xh(x)=ex。 则:

f′(x)=g′(x)h(x)+g(x)h′(x)=3x2ex+x3exf'(x) = g'(x)h(x) + g(x)h'(x) = 3x^2 e^x + x^3 e^xf′(x)=g′(x)h(x)+g(x)h′(x)=3x2ex+x3ex

f′′(x)=(6xex+3x2ex)+(3x2ex+x3ex)=6xex+6x2ex+x3exf''(x) = (6x e^x + 3x^2 e^x) + (3x^2 e^x + x^3 e^x) = 6x e^x + 6x^2 e^x + x^3 e^xf′′(x)=(6xex+3x2ex)+(3x2ex+x3ex)=6xex+6x2ex+x3ex

f′′′(x)=(6ex+6xex)+(12xex+6x2ex)+(3x2ex+x3ex)=6ex+18xex+9x2ex+x3exf'''(x) = (6e^x + 6x e^x) + (12x e^x + 6x^2 e^x) + (3x^2 e^x + x^3 e^x) = 6e^x + 18x e^x + 9x^2 e^x + x^3 e^xf′′′(x)=(6ex+6xex)+(12xex+6x2ex)+(3x2ex+x3ex)=6ex+18xex+9x2ex+x3ex

这个例子展示了莱布尼茨法则在计算高阶导数时的便捷性。

例6:莱布尼茨积分法则的更复杂应用

计算 ddx∫sin⁡xx2cos⁡(t2)dt \frac{d}{dx} \int_{\sin x}^{x^2} \cos(t^2) dt dxd​∫sinxx2​cos(t2)dt

这里我们应用莱布尼茨积分法则的完整形式。 令 a(x)=sin⁡xa(x) = \sin xa(x)=sinx, b(x)=x2b(x) = x^2b(x)=x2, f(x,t)=cos⁡(t2)f(x, t) = \cos(t^2)f(x,t)=cos(t2)。 则:

ddx∫sin⁡xx2cos⁡(t2)dt=cos⁡((x2)2)(2x)−cos⁡((sin⁡x)2)(cos⁡x) \frac{d}{dx} \int_{\sin x}^{x^2} \cos(t^2) dt = \cos((x^2)^2) (2x) - \cos((\sin x)^2) (\cos x) dxd​∫sinxx2​cos(t2)dt=cos((x2)2)(2x)−cos((sinx)2)(cosx)

这个例子展示了莱布尼茨积分法则在处理积分限为函数的情况下的应用。

这些例题涵盖了莱布尼茨公式的不同方面,从简单的近似计算到更复杂的微积分运算。 希望这些例子能够帮助你更好地理解莱布尼茨公式的应用。 记住,选择哪种方法来解决问题取决于问题的具体情况。 有时,直接运用莱布尼茨公式很方便,有时则需要结合其他微积分技巧。

Logo

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

更多推荐