数值分析思考题（钟尔杰版）参考解答——第五章
题目：
1. 用代数插值方法描述平面两点定一直线的数学问题。
假定给定区间
[
x
k
,
x
k
+
1
]
\left[x_k, x_{k+1}\right]
[xk,xk+1] 及端点函数值
y
k
=
f
(
x
k
)
,
y
k
+
1
=
f
(
x
k
+
1
)
y_k=f\left(x_k\right), y_{k+1}=f\left(x_{k+1}\right)
yk=f(xk),yk+1=f(xk+1),
要求线 性揷值多项式
L
1
(
x
k
)
=
y
k
,
L
1
(
x
k
+
1
)
=
y
k
+
1
y
=
L
1
(
x
)
L_1\left(x_k\right)=y_k, L_1\left(x_{k+1}\right)=y_{k+1}\quad y=L_1(x)
L1(xk)=yk,L1(xk+1)=yk+1y=L1(x) 的几何意义就是通过两点
(
x
k
,
y
k
)
\left(x_k, y_k\right)
(xk,yk) 与
(
x
k
+
1
,
y
k
+
1
)
\left(x_{k+1}, y_{k+1}\right)
(xk+1,yk+1) 的直线,
L
1
(
x
)
L_1(x)
L1(x) 的表达式可由几何意义直接给出
L
1
(
x
)
=
y
k
+
x
k
+
1
−
x
x
k
+
1
−
x
k
(
x
−
x
k
)
L_1(x)=y_k+\frac{x_{k+1}-x}{x_{k+1}-x_k}\left(x-x_k\right) \quad
L1(x)=yk+xk+1−xkxk+1−x(x−xk) 由上可得给定平面两点可确定一直线。
2. 如何用插值方法求得区间[–1，1]到[a，b]的变换公式？
给定区间t ϵ [-1,1],得函数值 x ϵ [a,b] 端点处函数值如下表所示
| t | -1 | 1 |
|---|---|---|
| x | a | b |
由两点式方程可以构造出变换公式为：
x
=
1
2
[
(
1
−
t
)
a
+
(
t
+
1
)
b
]
x=\frac{1}{2}[(1-t) a+(t+1) b]
x=21[(1−t)a+(t+1)b]
3. 插值问题的存在唯一性定理对插值条件是如何叙述的？
插值存在唯一性定理叙述插值条件如下所示:
设在区间
[
a
,
b
]
[\mathrm{a}, \mathrm{b}]
[a,b] 上给定
n
+
1
\mathrm{n}+1
n+1 个点，
a
≤
x
0
<
x
1
<
⋯
<
x
n
≤
b
a \leq x_0<x_1<\cdots<x_n \leq b
a≤x0<x1<⋯<xn≤b 上的函数值
y
i
=
f
(
x
i
)
,
(
i
=
0
,
1
,
…
n
)
y_i=f\left(x_i\right),(i=0,1, \ldots n)
yi=f(xi),(i=0,1,…n), 求次数不超过
n
\mathrm{n}
n 的多顶式
P
(
x
)
=
a
0
+
a
1
x
+
⋯
+
a
n
x
n
,
P(x)=a_0+a_1 x+\cdots+a_n x^n, \quad
P(x)=a0+a1x+⋯+anxn, 使
P
(
x
i
)
=
y
i
,
i
=
0
,
1
,
…
,
n
P\left(x_i\right)=y_i, i=0,1, \ldots, n
P(xi)=yi,i=0,1,…,n
即当代数插值问题中所给出的插值结点是 (n+1) 个互异的结点时, 则必然存 在唯一的一个
n
\mathrm{n}
n 次揷值多项式满足所给的揷值条件。代数插值问题本质上与一 个线性方程组求解问题等价。
4. 对二次多项式计算三点函数值并构造二次插值。插值与被插值函数有无区别？
无区别，证明如下：
由插值余项定理可知
R
(
x
)
=
f
(
x
)
−
P
n
(
x
)
=
f
(
n
+
1
)
(
ξ
)
(
n
+
1
)
!
ω
n
+
1
(
x
)
R(x)=f(x)-P_n(x)=\frac{f^{(n+1)}(\xi)}{(n+1) !} \omega_{n+1}(x)
R(x)=f(x)−Pn(x)=(n+1)!f(n+1)(ξ)ωn+1(x)
当
f
(
x
)
f(x)
f(x) 是
m
(
m
≤
n
)
\mathrm{m}(\mathrm{m} \leq n)
m(m≤n) 次多项式时
f
(
n
+
1
)
(
x
)
=
0
f^{(n+1)}(x)=0
f(n+1)(x)=0, 所以
R
(
x
)
=
0
R(x)=0
R(x)=0,
从而
f
(
x
)
=
P
n
(
x
)
f(x)=P_n(x)
f(x)=Pn(x), 所以当
f
(
x
)
f(x)
f(x) 是二次多项式时显然满足。所以无区别, 证明 完毕。
5. 拉格朗日插值基函数与牛顿插值基函数有何不同特点？
一、性质不同
- 1、牛顿插值：代数插值方法的一种形式。牛顿差值引入了差商的概念，使其在差值节点增加时便于计算。
- 2、拉格朗日插值：满足插值条件的、次数不超过n的多项式是存在而且是唯一的。
二、公式意义不同
- 1、牛顿插值：牛顿差值作为一种常用的数值拟合方法，由于其计算简单、计算点多、逻辑清晰、编程方便等特点，在实验分析中得到了广泛的应用。
特别是在实验中，当只能测量离散数据点或用数值解表示相应的关系时，可以用牛顿插值公式拟合离散点，得到更精确的函数解析值。 - 2、拉格朗日插值：在许多实际问题中，函数被用来表示某些内部关系或规律，许多函数只能通过实验和观察来理解。如果实际观测到一个物理量，并在多个不同的地点得到相应的观测值，拉格朗日插值法可以找到一个多项式，它可以精确地提取每个观测点的观测值。
综上：拉格朗日插值法与牛顿插值法都是二种常用的简便的插值法。但牛顿法插值法则更为简便，与拉格朗日插值多项式相比较，它不仅克服了“增加一个节点时整个计算工作必须重新开始的缺点，而且可以节省乘、除法运算次数。同时，在牛顿插值多项式中用到的差分与差商等概念，又与数值计算的其他方面有着密切的关系。
6. Runge反例说明一个什么样的问题？
Runge反例说明了用高次代数插值可能引起震荡，使得插值函数不收敛于被插函数，因而通常不用高次插值，而用分段低次插值。
7. 何谓切比雪夫插值? 切比雪夫插值误差有何改变？
- 使用平均分布的点作为插值多项式的基点X_i很普遍。在很多情况下，用于插值的数据点仅以这种形式存在，例如当数据由相同时间间隔分布的一起读取的数据所组成时。在其他情况下，我们可以在认为合适的地方自由选择基点。事实证明，基点间距选取的方式对于插值误差有很大的影响，切比雪夫插值是一种特定最优的点间距选取方式。
- 选择切比雪夫的根作为插值的基点，在区间[-1,1]中尽可能均匀地分散了插值误差，我们将使用切比雪夫根作为基点的插值多项式叫做切比雪夫插值多项式。
- 切比雪夫插值对于误差的改变：
- 总的来说，切比雪夫插值是一种特定最优的点间距选取方式，它和拉格朗日插值、牛顿差商公式相比，更能有效地避免龙格现象地出现，而且，能够更好地贴合原函数，也就是说，它的结果预测相对会更准确。
8.如何将区间[–1，1]上的切比雪夫插值结点映射到区间[a，b]上去？
移动基点使得它们在区间
[
a
,
b
]
[a, b]
[a,b] 上的相对位置和在区间
[
−
1
,
1
]
[-1,1]
[−1,1] 上一致。这 可以通过如下两步实现:
（1）使用因子
(
b
−
a
)
2
\frac{(b-a)}{2}
2(b−a) 拉伸点（这是两个区间长度的比值）
（2）将点平移
b
+
a
2
\frac{b+a}{2}
2b+a ，使得中心从 0 移动到区间
[
a
,
b
]
[\mathrm{a}, \mathrm{b}]
[a,b] 的中心。
换言之, 从原始点
:
x
k
=
cos
(
2
k
+
1
)
π
2
n
,
k
=
0
,
1
,
2
…
n
−
1
: x_k=\cos \frac{(2 k+1) \pi}{2 n}, k=0,1,2 \ldots n-1
:xk=cos2n(2k+1)π,k=0,1,2…n−1 移动到
x
i
=
b
−
a
2
cos
(
2
i
+
1
)
π
2
n
+
b
+
a
2
,
i
=
0
,
1
,
2
…
n
−
1
x_i=\frac{b-a}{2} \cos \frac{(2 i+1) \pi}{2 n}+\frac{b+a}{2}, i=0,1,2 \ldots n-1
xi=2b−acos2n(2i+1)π+2b+a,i=0,1,2…n−1
∣
(
x
−
x
i
)
…
(
x
−
x
n
)
∣
≤
(
b
−
a
2
)
n
2
n
−
1
\left|\left(x-x_i\right) \ldots\left(x-x_n\right)\right| \leq \frac{\left.\frac{(b-a}{2}\right)^n}{2^{n-1}}
∣(x−xi)…(x−xn)∣≤2n−12(b−a)n, 在区间
[
a
,
b
]
[\mathrm{a}, \mathrm{b}]
[a,b] 上成立。
9．埃尔密特插值与通常的插值方法有何不同？
埃尔密特插值: 有时候, 我们不仅要求插值函数在给定节点上函数值重 合, 而且要求若干阶导数重合; 即: 要求插值函数
φ
(
x
)
\varphi(x)
φ(x) 满足:
φ
(
x
i
)
=
f
(
x
i
)
φ
′
(
x
i
)
=
f
′
(
x
i
)
φ
′
′
(
x
i
)
=
f
′
′
(
x
i
)
φ
(
n
)
(
x
i
)
=
f
(
n
)
(
x
i
)
\begin{aligned} &\varphi\left(x_i\right)=f\left(x_i\right) \\ &\varphi^{\prime}\left(x_i\right)=f^{\prime}\left(x_i\right) \\ &\varphi^{\prime \prime}\left(x_i\right)=f^{\prime \prime}\left(x_i\right) \\ &\varphi^{(n)}\left(x_i\right)=f^{(n)}\left(x_i\right) \end{aligned}
φ(xi)=f(xi)φ′(xi)=f′(xi)φ′′(xi)=f′′(xi)φ(n)(xi)=f(n)(xi)
- 通常的插值方法：
不需要节点上的若干阶导数值也相等，只需要在插值节点上的函数值相等。
10．样条插值有何特点？根据三次样条的数学概念描述二次样条函数。
-
什么是样条插值：简言之，就是依然对每一个小区间进行插值，但是我们不需要依赖于导数的已知；由于我们要做的就是使得端点处满足某种条件的光滑（一般来说，三次样条插值就是要满足二阶导数连续），根据这个要求，我们就可以在未知导数的情况下推导出样条函数。
-
因此，样条插值特点：可以在未知导数的情况下推导出样条函数
-
二次样条揷值:
若 Δ : a = x 0 < x 1 < x 2 … < x n = b \Delta: \mathrm{a}=x_0<x_1<x_2 \ldots<x_n=b Δ:a=x0<x1<x2…<xn=b 是区间 [ a , b ] [\mathrm{a}, \mathrm{b}] [a,b] 的一个分割, 函数 S ( x ) \mathrm{S}(\mathrm{x}) S(x) 满足下列 条件:
（1） 在每个小区间 [ x i , x i + 1 ] ( i = 0 , 1 , ⋯ , n ) \left[x_i, x_{i+1}\right](\mathrm{i}=0,1, \cdots, \mathrm{n}) [xi,xi+1](i=0,1,⋯,n) 上 S ( x ) \mathrm{S}(\mathrm{x}) S(x) 是二次多项式;
（2） S ( x ) \mathrm{S}(\mathrm{x}) S(x) 在区间 [ a , b ] [\mathrm{a}, \mathrm{b}] [a,b] 上连续, 则称 S ( x ) \mathrm{S}(\mathrm{x}) S(x) 为区间 [ a , b ] [\mathrm{a}, \mathrm{b}] [a,b] 上的二次样条函数。
如果 S ( x ) \mathrm{S}(\mathrm{x}) S(x) 还满足 S ( x i ) = y i , ( i = 0 , 1 , … n ) S\left(x_i\right)=y_i,(i=0,1, \ldots n) S(xi)=yi,(i=0,1,…n) 则称 S ( x ) \mathrm{S}(\mathrm{x}) S(x) 为区间 [ a , b ] [\mathrm{a}, \mathrm{b}] [a,b] 上的二次样条 插值函数。
扫一扫
2127
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
