数学公式之求 log2(1+x)-x的积分

为了算神奇数字0x5F3759DF,有一步绕不开,那就是求 ∫ 0 1 [ l o g 2 ( 1 + x ) − x ] d x (0) \int_0^1[log_2(1+x)-x]dx \tag{0} 01[log2(1+x)x]dx(0) 为此,我们将详细探究这个一下这个函数 F ( x ) = l o g 2 ( 1 + x ) − x (1) F(x)=log_2(1+x)-x \tag{1} F(x)=log2(1+x)x(1)

0. 准备工作

首先, 准备两个数字: l n 2 = 0.693147180 ln2=0.693147180 ln2=0.693147180 1 l n 2 = 1.442695041 \cfrac{1}{ln2}=1.442695041 ln21=1.442695041

1.画出 F ( x ) F(x) F(x)图像

数学公式之求 log2(1+x)-x的积分_第1张图片

2. 画出其导数图像

F ′ ( x ) = 1 l n 2 × 1 1 + x − 1 (2) F'(x)= \cfrac{1}{ln2} \times \cfrac{1}{1+x} - 1 \tag{2} F(x)=ln21×1+x11(2)
数学公式之求 log2(1+x)-x的积分_第2张图片

3. 求 F ′ ( x ) F'(x) F(x)在[0, 1]上的零点以及对应F(x)的极值

由(2)式,令 F ′ ( x ) = 0 F'(x)=0 F(x)=0,求出 x 0 = 1 l n 2 − 1 = 0.442695 (3) x_0=\cfrac{1}{ln2} - 1=0.442695 \tag{3} x0=ln211=0.442695(3) 对应的把x0带入F(x),就求出其在[0,1]上的极大值就是 F ( x 0 ) = 1 − l n ( l n 2 ) + 1 l n 2 = 0.086071 (4) F(x_0)=1-\cfrac{ln(ln2)+1}{ln2}=0.086071 \tag{4} F(x0)=1ln2ln(ln2)+1=0.086071(4) x 0 x_0 x0, F ( x 0 ) F(x_0) F(x0)均可以在以上两张图得到验证

4. 求 F ( x ) F(x) F(x)在[0, 1]上的积分

(数学好的直接跳过,看这一部分最后的结果(7)式即可)
∫ 0 1 F ( x ) d x = ∫ 0 1 [ l o g 2 ( 1 + x ) − x ] d x \int_0^1 F(x) dx= \int_0^1[log_2(1+x)-x]dx 01F(x)dx=01[log2(1+x)x]dx ∫ 0 1 F ( x ) d x = ∫ 0 1 l o g 2 ( 1 + x ) d x − ∫ 0 1 x d x (5) \int_0^1 F(x) dx =\int_0^1log_2(1+x)dx-\int_0^1xdx \tag{5} 01F(x)dx=01log2(1+x)dx01xdx(5)
因为 ∫ 0 1 x d x = x 2 2 ∣ x = 0 x = 1 = 1 2 \int_0^1xdx=\cfrac{x^2}{2} \biggm\vert_{x=0}^{x=1}=\cfrac{1}{2} 01xdx=2x2x=0x=1=21, 因此我们只需要计算的是前半部分 ∫ 0 1 F 1 ( x ) d x = ∫ 0 1 l o g 2 ( 1 + x ) d x \int_0^1 F_1(x) dx= \int_0^1log_2(1+x)dx 01F1(x)dx=01log2(1+x)dx即可。
∫ 0 1 F 1 ( x ) d x = ∫ 0 1 l o g 2 ( 1 + x ) d x \int_0^1 F_1(x) dx= \int_0^1log_2(1+x)dx 01F1(x)dx=01log2(1+x)dx = ∫ 0 1 l n ( 1 + x ) l n 2 d x =\int_0^1\cfrac{ln(1+x)}{ln2}dx =01ln2ln(1+x)dx = 1 l n 2 ∫ 0 1 l n ( 1 + x ) d x =\cfrac{1}{ln2}\int_0^1ln(1+x)dx =ln2101ln(1+x)dx
注意接下来用了分部积分法

= 1 l n 2 × [ l n ( 1 + x ) × x − ∫ x d [ l n ( 1 + x ) ] ] =\cfrac{1}{ln2}\times[ln(1+x) \times x - \int x d[ln(1+x)]] =ln21×[ln(1+x)×xxd[ln(1+x)]] = 1 l n 2 × [ l n ( 1 + x ) × x − ∫ x 1 + x d x ] ∣ x = 0 x = 1 =\cfrac{1}{ln2}\times[ln(1+x) \times x - \int \cfrac{x}{1+x} dx] \biggm\vert_{x=0}^{x=1} =ln21×[ln(1+x)×x1+xxdx]x=0x=1 = 1 l n 2 × [ l n ( 1 + x ) × x − ∫ ( 1 + x 1 + x − 1 1 + x ) d x ] ∣ x = 0 x = 1 =\cfrac{1}{ln2}\times[ln(1+x) \times x - \int (\cfrac{1+x}{1+x} - \cfrac{1}{1+x}) dx] \biggm\vert_{x=0}^{x=1} =ln21×[ln(1+x)×x(1+x1+x1+x1)dx]x=0x=1 = 1 l n 2 × [ l n ( 1 + x ) × x − ∫ 1 d x + ∫ 1 1 + x d x ] ∣ x = 0 x = 1 =\cfrac{1}{ln2}\times[ln(1+x) \times x - \int 1dx + \int\cfrac{1}{1+x} dx] \biggm\vert_{x=0}^{x=1} =ln21×[ln(1+x)×x1dx+1+x1dx]x=0x=1 = 1 l n 2 × [ l n ( 1 + x ) × x − x + l n ( 1 + x ) ] ∣ x = 0 x = 1 =\cfrac{1}{ln2}\times[ln(1+x)\times x - x + ln(1+x)] \biggm\vert_{x=0}^{x=1} =ln21×[ln(1+x)×xx+ln(1+x)]x=0x=1 = 1 l n 2 × [ l n ( 1 + x ) × ( 1 + x ) − x ] ∣ x = 0 x = 1 =\cfrac{1}{ln2}\times[ln(1+x)\times (1+x) - x ] \biggm\vert_{x=0}^{x=1} =ln21×[ln(1+x)×(1+x)x]x=0x=1 ∫ 0 1 F 1 ( x ) d x = 1 l n 2 × ( 2 l n 2 − 1 ) (6) \int_0^1 F_1(x) dx=\cfrac{1}{ln2} \times (2ln2 - 1) \tag{6} 01F1(x)dx=ln21×(2ln21)(6) 综上,结合(3)(4)式, ∫ 0 1 F ( x ) d x = 1 l n 2 × ( 2 l n 2 − 1 ) − 1 2 \int_0^1 F(x) dx =\cfrac{1}{ln2} \times (2ln2 - 1) - \cfrac{1}{2} 01F(x)dx=ln21×(2ln21)21 ∫ 0 1 F ( x ) d x = 3 2 − 1 l n 2 = 0.057304959 (7) \int_0^1 F(x) dx = \cfrac{3}{2} - \cfrac{1}{ln2} = 0.057304959 \tag{7} 01F(x)dx=23ln21=0.057304959(7)

你可能感兴趣的:(计算机科学,算法)