第二章值法与数值微分.docVIP

第二章 插值法与数值微分 1. 设,在三处的值是很容易求得的,试以这三个点建立的二次插值多项式,并用此多项式计算的近似值,且给出误差估计.用其中的任意两点,构造线性插值函数,用得到的三个线性插值函数,计算的近似值,并分析其结果不同的原因. 解: 已知, 建立二次Lagrange插值函数可得: 所以. 误差,所以 利用前两个节点建立线性插值函数可得: 所以. 利用后两个节点建立线性插值可得: 所以. 利用前后两个节点建立线性插值可得: 所以. 与的真实值比较,二次插值比线性插值效果好,利用前两个节点的线性插值比其他两个线性插值效果好.此说明,二次插值比线性插值效果好,内插比外插效果好. 2. 利用(2.9)式证明 证明: 由(2.9)式 当时, , 所以 3. 若为互异节点,且有 证明 证明: 由于 且和都为k次多项式,而且在k+1个不同的节点处的函数值都相同 , 所以马上有 . 4. 设给出在上的数值表,用二次插值进行计算,若希望截断误差小于,问函数表的步长最大能取多少? 解: 记插值函数为p(x),则 所以 ;令,设,得 又的最大值为,所以有 所以 . 5. 用拉格朗日插值和牛顿插值找经过点的三次插值公式. 解: Lagrange插值函数: 牛顿插值: 首先计算差商 也可以利用等距节点构造,首先计算差分 可得前插公式 和后插公式 6. 确定一次数不高于4的多项式,使. 解: 利用重节点计算差商 则可构造Hermite插值函数满足题设条件: 7. 寻找过个点的次多项式,满足条件: 解: 和Lagrange插值函数的构造类似,可将插值函数写成 其中,基函数满足条件 (1); (2) 则可由已知条件,可得 ; . 所以可得 8. 过0,1两点构造一个三次Hermite插值多项式,满足条件: 解: 计算重节点的差商 马上可得 9. 过给定数组 75 76 77 78 79 80 2.768 2.833 2.903 2.979 3.062 3.153 作一分段线性插值函数. 取第二类边界条件,作三次样条插值多项式. 用两种插值函数分别计算的函数值. 解: (1)做分段线性插值函数可得: 其中, (2)把已知节点值带入M关系式可得: 由边界条件可得,所以上面方程组变为可求解方程组 解得. 所以可得在每个区间上的三次样条函数的表达式: (3)当时, ; 当时, ; 10. 若给出的函数表: 1.567 1.568 1.569 1.570 0.999 992 8 0.999 996 1 0.999 998 4 0.999 999 7 0.003 796 3 0.002 796 3 0.001 796 3 0.000 796 3 263.411 25 357.611 06 556.690 98 1 255.765 59 用表上的数据和任一插值公式求: 用表格直接计算. 用和来计算.并讨论这两个结果中误差变化的原因. 解: 利用Lagrange插值直接用tan表计算得 ; 利用Lagrange插值计算sin得 ; 利用Lagrange插值计算cos得 ; 最后利用sin/cos计算tan得 . 出现小除数,误差被放大. 11. 求三次样条函数,已知 0.25 0.30 0.39 0.45 0.53 0.500 0 0.547 7 0.624 5 0.670 8 0.728 0 和边界条件 解: 把表中数据带入M关系式可得 由边界条件还可得到两个方程: 联立两个方程组可解得: 带入M表达式便可得所求三次样条函数. 12. 称n阶方阵具有严格对角优势,若 试证明:具有严格对角优势的方阵必可逆. 证明:方程组(2.62)解存在唯一. 证明: (1)设矩阵A按行严格对角占优,如果A奇异,则存在非零向量x使得Ax=0,写成分量形式为 令指标使得,则 因