←
PDF 719 / 1160 This integral cannot be evaluated analytically. Using the trapezium rule with step sizes $\pi/2$
→
English · PDF 719
Original PDF page 719
中文 · PDF 719

该积分无法解析求值。使用步长为 $\pi/2$ , $\pi/4$ 和 $\pi/8$ 的梯形法则,结合理查森外推,我们得到

$$\bar{y} = 2.9148$$

现在我们寻求一个形如 $y(x)$ 的近似,

$$y(x) \simeq \bar{y} + \cos x + a \cos 2x$$

使得平方误差在一个完整周期上的积分尽可能小。即, $a$ 的选取使得

$$\frac{d}{da} \left{ \int_{-\pi}^{\pi} [\bar{y} + \cos x + a \cos 2x - y(x)]^2 dx \right} = 0$$

展开被积函数,得到

$$\begin{aligned} \frac{d}{da} \left{ \int_{-\pi}^{\pi} [\bar{y}^2 + a^2 \cos^2 2x + 8 + \cos^2 x + 2\bar{y}a \cos 2x \right. \ \left. - 2\bar{y}\sqrt{(8 + \cos^2 x)} - 2a \cos 2x \sqrt{(8 + \cos^2 x)}] dx \right} = 0 \end{aligned}$$

整理后得

$$\begin{aligned} \frac{d}{da} \left{ \int_{-\pi}^{\pi} [\bar{y}^2 + 8 + \cos^2 x - 2\bar{y}\sqrt{(8 + \cos^2 x)}] dx + 2a\bar{y} \int_{-\pi}^{\pi} \cos 2x dx \right. \ \left. - 2a \int_{-\pi}^{\pi} \cos 2x \sqrt{(8 + \cos^2 x)} dx + a^2 \int_{-\pi}^{\pi} \cos^2 2x dx \right} = 0 \end{aligned}$$

大括号内的第一个积分与 $a$ , 因此求导后为零。第二个积分值为零,最后一个积分的值为 $\pi$ 。因此对 $a$ 求导得到

$$-2 \int_{-\pi}^{\pi} \cos 2x \sqrt{(8 + \cos^2 x)} dx + 2a\pi = 0$$

于是

$$a = \frac{1}{\pi} \int_{-\pi}^{\pi} \cos 2x \sqrt{(8 + \cos^2 x)} dx$$

对该积分进行数值计算得到 $a = 0.0858$ 。考察近似值与 $y$ 在一个完整周期内的差异可知,最大误差的大小小于 0.0007。