找回密码
 快速注册
搜索
查看: 19|回复: 3

[函数] 用二次多项式近似$(1 - x^2)^{1/2}$使误差的最大值最小化

[复制链接]

3149

主题

8387

回帖

6万

积分

$\style{scale:11;fill:#eff}꩜$

积分
65396
QQ

显示全部楼层

hbghlyj 发表于 2024-10-29 22:00 |阅读模式
https://web.archive.org/web/2022 ... n_chap4.pdf#page=12
在区间 $[-1,1]$ 上,用二次多项式近似函数 $f(x) = (1 - x^2)^{1/2}$,使误差的最大值最小化。

首先,我们注意到如果 $x=\cos \theta$,那么 $f(\cos \theta)=\left(1-\cos ^{2} \theta\right)^{1 / 2}=\sin \theta$,并且区间 $x \in[-1,1]$ 变为 $\theta \in[0, \pi]$。
$\sin \theta$ 在 $[0, \pi]$ 上的傅里叶余弦级数为
$$\sin \theta=\frac{2}{\pi}-\frac{4}{\pi}\left[\frac{\cos 2 \theta}{3}+\frac{\cos 4 \theta}{15}+\frac{\cos 6 \theta}{35}+\ldots\right]$$
代入 $x =\cos θ$ 得$$\left(1-x^{2}\right)^{1 / 2}=\frac{2}{\pi}-\frac{4}{\pi}\left[\frac{T_{2}(x)}{3}+\frac{T_{4}(x)}{15}+\frac{T_{6}(x)}{35}+\ldots\right],\quad\text{其中 }-1\le x\le1$$
(因为要求二次多项式)我们取前两项,
\begin{align*}
p_{2}(x) &=\frac{2}{\pi}-\frac{4}{\pi} \frac{T_{2}(x)}{3}=\frac{2}{\pi}-\frac{4}{3 \pi}\left(2 x^{2}-1\right) \\
&=\frac{2}{3 \pi}\left(3-2\left(2 x^{2}-1\right)\right)=\frac{2}{3 \pi}\left(5-4 x^{2}\right).
\end{align*}
误差
$$f(x)-p_{2}(x) \approx-\frac{4}{\pi} \frac{T_{4}(x)}{15},$$
在 $[-1,1]$ 上振荡 4+1=5 次。为了使误差的最大值最小化,至少需要 4 个等振幅点。
$T_{4}(x)$ 在 $8 x^{4}-8 x^{2}+1=\pm 1$ 时有极值,即在
$$x=0, x=1, x=-1, x=1 / \sqrt{2} \text { 和 } x=-1 / \sqrt{2} .$$
\begin{array}{l|c|c|c}
& \left(1-x^{2}\right)^{1 / 2} & p_{2}(x) & \text { 误差 } \\
\hline
x=0 & 1 & 10 / 3 \pi & -0.0610 \\
x=\pm 1 / \sqrt{2} & 1 / \sqrt{2} & 2 / \pi & 0.0705 \\
x=\pm 1 & 0 & 2 / 3 \pi & -0.2122
\end{array}
所以误差振荡但不等振幅。因此,$p_2(x)$ 不是 $f(x) = (1 − x^2)^{1/2}$ 的使误差的最大值最小化的近似二次多项式,但它是一个好的初步近似。

$\left(1-x^{2}\right)^{1 / 2}$ 的真正的使误差的最大值最小化的近似二次多项式实际上是 $\left(\frac{9}{8}-x^{2}\right)=\left(1.125-x^{2}\right)$,因此我们估计的 $\left(1.061-0.8488 x^{2}\right)$ 还不错。

3149

主题

8387

回帖

6万

积分

$\style{scale:11;fill:#eff}꩜$

积分
65396
QQ

显示全部楼层

 楼主| hbghlyj 发表于 2024-10-29 22:06
hbghlyj 发表于 2024-10-29 14:00
至少需要 4 个等振幅点。

为什么需要使上述振荡点的振幅相等

61

主题

155

回帖

1241

积分

积分
1241

显示全部楼层

血狼王 发表于 2024-11-5 00:12
hbghlyj 发表于 2024-10-29 22:06
为什么需要使上述振荡点的振幅相等

这是“切比雪夫交错点组定理”(又名“一致最优逼近多项式判定定理”)

血狼王者,格罗特克斯(Grotex)是也。
AOPS的id:Grotex

3149

主题

8387

回帖

6万

积分

$\style{scale:11;fill:#eff}꩜$

积分
65396
QQ

显示全部楼层

 楼主| hbghlyj 发表于 2024-11-5 01:18
血狼王 发表于 2024-11-4 16:12
这是“切比雪夫交错点组定理”

谢谢回复,补充一下,NA_2023
定理 9.17 (Chebyshev 交错定理) 设函数 $f \in C[a, b]$ 且 $f \notin \mathbb{P}_n[x]$, 则 $p^*$ 是 $f$ 的 $n$次最佳一致逼近多项式的充分必要条件是 $f-p^*$ 在 $[a, b]$ 上存在有 $n+2$ 个点组成的交错点组,即有 $a \leqslant x_0<x_1<\cdots<x_{n+1} \leqslant b$ 使得
$$
f\left(x_i\right)-p^*\left(x_i\right)=(-1)^i \sigma\left\|f-p^*\right\|_{\infty}, \quad i=0,1, \cdots, n+1 .
$$
证明 (充分性) 设 $f-p^*$ 在 $[a, b]$ 上存在有 $n+2$ 个点组成的交错点组 $\left\{x_i\right\}_{i=0}^k \subset \mathcal{L}_M$.用反证法. 如果存在多项式 $p \in \mathbb{P}_n[x]$ 使得条件 (9.28) 成立, 那么 $p$ 显然至少有 $n+1$个零点. 另一方面, 利用代数基本定理知, 次数不超过 $n$ 次的非零多项式 $p$ 至多有 $n$ 零点,矛盾!因此,由定理(9.16)知, $p^*$ 是 $f$ 的最佳一致逼近多项式。
(必要性) 设 $p^*$ 是 $f$ 的最佳一致逼近多项式, 记 $E_n(f)=\left\|f-p^*\right\|_{\infty}$, 因 $f \notin \mathbb{P}_n[x]$,故 $E_n(f)>0$. 用反证法. 设 $f-p^*$ 在 $[a, b]$ 上的任意交错点组为 $\left\{x_i\right\}_{i=0}^k$, 均有 $k \leqslant n$.因此, 存在区间 $[a, b]$ 的一个分割 $\left\{t_i\right\}_{i=0}^{k+1}$, 其中 $t_0=a, t_{k+1}=b$, 使得
$$
S_i=\left\{x \in\left[t_{i-1}, t_i\right]:\left|\left(f-p^*\right)(x)\right|=E_n(f)\right\}, \quad i=1,2, \cdots, k+1,
$$
非空,且函数 $f-p^*$ 在 $S_i$ 上保持符号不变。由于 $f-p^*$ 在相邻 $S_i$ 和 $S_{i+1}$ 上交替取 $\pm E_n(f)$, 故 $\left[t_{i-1}, t_{i+1}\right]$ 上至少存在一个 $f-p^*$ 的零点, 不失一般性, 可以调整 $t_i$ 使得
$$
\left(f-p^*\right)\left(t_i\right)=0, \quad i=1,2, \cdots, k .
$$
由于 $k \leqslant n$, 可构造如下 $k$ 次多项式
$$
p(x)=\sigma \prod_{i=1}^k\left(x-t_i\right) \in \mathbb{P}_n[x]
$$
使得 $p$ 满足条件 (9.28). 而这与定理(9.16)矛盾!因此,假设错误,原命题成立.

手机版|悠闲数学娱乐论坛(第3版)

GMT+8, 2025-3-4 16:07

Powered by Discuz!

× 快速回复 返回顶部 返回列表