分部求和法简介

分部求和法

首先我们知道
d d x ( u ( x ) v ( x ) ) = u ( x ) v ( x ) + u ( x ) v ( x ) \cfrac{\mathrm{d}}{\mathrm{d}x}(u(x)v(x))=u'(x)v(x)+u(x)v'(x)

d ( u ( x ) v ( x ) ) = u ( x ) v ( x ) d x + u ( x ) v ( x ) d x \mathrm d(u(x)v(x))=u'(x)v(x)\mathrm{d}x+u(x)v'(x)\mathrm{d}x

d ( u ( x ) v ( x ) ) = u ( x ) v ( x ) d x + u ( x ) v ( x ) d x \int\mathrm d(u(x)v(x))=\int u'(x)v(x)\mathrm dx+\int u(x)v'(x)\mathrm dx

移项并注意到 d ( u ( x ) v ( x ) ) = u ( x ) v ( x ) \int \mathrm d(u(x)v(x))=u(x)v(x)
u ( x ) v ( x ) d x = u ( x ) v ( x ) u ( x ) v ( x ) d x \int u(x)v'(x)\mathrm dx=u(x)v(x)-\int u'(x)v(x)\mathrm dx
这启示我们:在求解积分 f ( x ) d x \int f(x)\mathrm dx 时,如果能找到合适的 u ( x ) , v ( x ) u(x),v(x) ,使得被积函数 f ( x ) f(x) 可以表示为 u ( x ) v ( x ) u(x)v'(x) 的形式,并且 u ( x ) v ( x ) u(x)v(x) u ( x ) v ( x ) d x \int u'(x)v(x)\mathrm dx 比较好求,就可以用后两项相减求得 f ( x ) f(x) 的积分。

这个式子可以简单地写作
u d v = u v v d u \int u\mathrm dv=uv-\int v\mathrm du
它称为分部积分法

举个栗子:求 x cos x d x \int x\cos x\mathrm dx

解:设 u = x , d v = cos x d x u=x,\mathrm dv=\cos x\mathrm dx ,那么 v = sin x , d u = d x v=\sin x,\mathrm du=\mathrm dx

x cos x d x = u d v = u v v d u = x sin x sin x d x = x sin x + cos x + C \int x\cos x\mathrm dx=\int u\mathrm dv=uv-\int v\mathrm du=x\sin x-\int\sin x\mathrm dx=x\sin x+\cos x+C

这种方法同样可以运用到求和式上,我们有对应的分部求和法

先介绍若干记号。

Δ f ( x ) = f ( x + 1 ) f ( x ) \Delta f(x)=f(x+1)-f(x) ,为函数 f ( x ) f(x) 的差分函数。 Δ \Delta 称为差分算子。

扫描二维码关注公众号,回复: 6669954 查看本文章

f ( x ) δ x \sum f(x)\delta x 不定和式,它的含义如下: F ( x ) = f ( x ) δ x F(x)=\sum f(x)\delta x 当且仅当 Δ F ( x ) = f ( x ) \Delta F(x)=f(x)

f ( x ) a b = f ( b ) f ( a ) f(x)|_a^b=f(b)-f(a)

a b f ( x ) δ x \sum_a^b f(x)\delta x 定和式,它的含义如下:若 F ( x ) = f ( x ) δ x F(x)=\sum f(x)\delta x ,则 a b f ( x ) δ x = F ( b ) F ( a ) \sum_a^b f(x)\delta x=F(b)-F(a)

利用传统的 \sum 表示方法,我们有 a b f ( x ) δ x = k = a b 1 f ( k ) \sum_a^b f(x)\delta x=\sum\limits_{k=a}^{b-1}f(k) 。注意这里求和的上界是 b 1 b-1 而不是 b b

在有些参考书上 a b \sum_a^b 也写作 a b \sum\limits_a^b

为了推导分部求和法,我们先推导两个函数的积的差分。
Δ ( u ( x ) v ( x ) ) = u ( x + 1 ) v ( x + 1 ) u ( x ) v ( x ) \Delta(u(x)v(x))=u(x+1)v(x+1)-u(x)v(x)

= ( u ( x ) + Δ u ( x ) ) v ( x + 1 ) u ( x ) ( v ( x + 1 ) Δ v ( x ) ) =(u(x)+\Delta u(x))v(x+1)-u(x)(v(x+1)-\Delta v(x))

= Δ u ( x ) v ( x + 1 ) + u ( x ) Δ v ( x ) =\Delta u(x)v(x+1)+u(x)\Delta v(x)

两边求和
u ( x ) v ( x ) = Δ u ( x ) v ( x + 1 ) d x + u ( x ) Δ v ( x ) d x u(x)v(x)=\sum\Delta u(x)v(x+1)\mathrm dx+\sum u(x)\Delta v(x)\mathrm dx

u ( x ) Δ v ( x ) d x = u ( x ) v ( x ) Δ u ( x ) v ( x + 1 ) d x \sum u(x)\Delta v(x)\mathrm dx=u(x)v(x)-\sum\Delta u(x)v(x+1)\mathrm dx

这个式子可以简写成
u Δ v δ x = u v Δ u E v δ x \sum u\Delta v\delta x=uv-\sum\Delta u\mathrm Ev\delta x
其中 E \mathrm E 移位算子,表示 E f ( x ) = f ( x + 1 ) \mathrm Ef(x)=f(x+1)

很多时候我们要求的都是定和式,那么定和式的分部求和法形式就是
a b u Δ v δ x = ( u v ) a b a b Δ u E v δ x \sum_a^bu\Delta v\delta x=(uv)|_a^b-\sum_a^b\Delta u\mathrm Ev\delta x
举个栗子:化简 S = k = 0 n 1 C k m H k S=\sum\limits_{k=0}^{n-1}C_k^mH_k 。其中 H n H_n 调和级数,定义为 H 0 = 0 , H n = i = 1 n 1 i H_0=0,H_n=\sum\limits_{i=1}^n\cfrac{1}{i}

这里我们使用广义组合数,意味着当 k < m k<m C k m = 0 C_k^m=0

解: S = 0 n C x m H x δ x S=\sum_0^nC_x^mH_x\delta x

u = H x , v = C x m + 1 u=H_x,v=C_x^{m+1} ,那么 Δ v = C x + 1 m + 1 C x m + 1 = C x m , Δ u = 1 x + 1 \Delta v=C_{x+1}^{m+1}-C_x^{m+1}=C_x^m,\Delta u=\cfrac{1}{x+1}

S = 0 n u Δ v δ x = ( u v ) 0 n 0 n 1 x + 1 C x + 1 m + 1 S=\sum_0^nu\Delta v\delta x=(uv)|_0^n-\sum_0^n\frac{1}{x+1}C_{x+1}^{m+1}

= C n m + 1 H n k = 0 n 1 1 k + 1 C k + 1 m + 1 =C_n^{m+1}H_n-\sum\limits_{k=0}^{n-1}\frac{1}{k+1}C_{k+1}^{m+1}

右边这个和式再熟悉不过了,利用吸收恒等式 1 k + 1 C k + 1 m + 1 = 1 m + 1 C k m \frac{1}{k+1}C_{k+1}^{m+1}=\frac{1}{m+1}C_k^m ,那么

S = C n m + 1 H n 1 m + 1 k = 0 n 1 C k m = C n m + 1 H n 1 m + 1 C n m + 1 = C n m + 1 ( H n 1 m + 1 ) S=C_n^{m+1}H_n-\frac{1}{m+1}\sum\limits_{k=0}^{n-1}C_k^m=C_n^{m+1}H_n-\frac{1}{m+1}C_n^{m+1}=C_n^{m+1}(H_n-\frac{1}{m+1})

还有一个我们很熟悉的栗子:化简 T = k = 1 n k C m + k m + 1 T=\sum\limits_{k=1}^nkC_{m+k}^{m+1} 。下面我们用分部求和法来试一试。

T = 1 n + 1 x C m + x m + 1 δ x T=\sum_1^{n+1}xC_{m+x}^{m+1}\delta x 。设 u = x , v = C m + x m + 2 u=x,v=C_{m+x}^{m+2} ,那么 Δ u = 1 , Δ v = C m + x m + 1 \Delta u=1,\Delta v=C_{m+x}^{m+1}

于是 T = 1 n + 1 u Δ v δ x = ( u v ) 1 n + 1 1 n + 1 Δ u E v δ x T=\sum_1^{n+1}u\Delta v\delta x=(uv)|_1^{n+1}-\sum_1^{n+1}\Delta u\mathrm Ev\delta x

= ( n + 1 ) C m + n + 1 m + 2 k = 1 n C m + k + 1 m + 2 =(n+1)C_{m+n+1}^{m+2}-\sum\limits_{k=1}^nC_{m+k+1}^{m+2}

= ( n + 1 ) C m + n + 1 m + 2 C m + n + 2 m + 3 =(n+1)C_{m+n+1}^{m+2}-C_{m+n+2}^{m+3}

比传统做法更简便。

猜你喜欢

转载自blog.csdn.net/qq_39677783/article/details/93842036