第四-八章 习题解答.pptVIP

第四章 线性方程组的迭代解法习 题 解 答 第五章 数值插值方法习 题 解 答 第六章 数据拟合与函数逼近习 题 解 答 第七章 数值积分习 题 解 答 第八章 常微分方程的数值解法习 题 解 答 1. 用最小二乘法求解超定方程组 解： 超定方程组的矩阵形式为 将方程两端同时乘以系数矩阵的转置矩阵，得正规方程组 求解正规方程组，得 2. 观测一个做直线运动的物体，测得以下数据： 解： 在表中，时间单位为s（秒），距离单位为m（米）. 假若加速度为常数，求物体的加速度和初速度. 设物体运动的初速度和加速度分别为v0和a，初始时刻距离为0，则距离函数为 30 1.9 50 3.0 110 80 10 0 距离S 5.0 3.9 0.9 0 时间t 用后五个点的数据作曲线拟合，可得超定方程组 求解得 解： 令 30 2 15 20 60 y 4 3 1 x 则 数据转换如下 3. 用最小二乘法求一个形如y=AeBx的经验公式，使其与下列数据相拟合： 3.4012 2 2.7081 2.9957 4.0943 z=lny 4 3 1 x 对上述数据作曲线拟合，可得超定方程组 求解得 所以 A=84.8528，故所求经验公式为 4. 已知实验观测数据（xi, yi)(i=1,2,…,m). 令 取拟合函数为 试利用曲线拟合的最小二乘法确定组合系数a0,a1（推导计算公式）. 解： 记 将实验数据代入拟合函数，得超定方程组 即 将方程组两端同时左乘 得超定方程组 因为 故 所以 6. 求f(x)=ex在区间[-1,1]上的三次最佳逼近多项式. 解： 利用勒让德正交多项式作基函数，即设 其中 利用正交性，得系数为 又因为 故 所以 1. 对给定节点x0,x1，插值型求积公式为 证明： 将f(x)=1,x分别代入公式，得 试证明：求积系数为 解得 2. 确定下列求积公式中的待定参数，使其代数精度尽可能高，并指明所构造出的求积公式所具有的代数精度. 解： （1）分别用f(x)=1,x,x2代入公式，得 解得 故求积公式至少具有2次代数精度. 再用f(x)=x3代入公式有 故求积公式成立. 但若用f(x)=x4代入则 故所求求积公式具有3次代数精度. （2）分别用f(x)=1,x,x2代入公式，得 解得 故求积公式至少具有2次代数精度. 再用f(x)=x3代入公式有 故求积公式成立. 但若用f(x)=x4代入则 故所求求积公式具有3次代数精度. （3）分别用f(x)=1,x,x2代入公式，得 解得 故求积公式至少具有2次代数精度. 再用f(x)=x3代入公式有 或 故求积公式具有2次代数精度. （4）由于求积公式中只有一个参数，用f(x)=1代入公式有 用f(x)=x代入公式有 当f(x)=x3时代入公式有 为恒等式. 仍为恒等式. 再用f(x)=x2代入公式有 解得 而当f(x)=x4时 故所求求积公式具有3次代数精度. 4. 推导下列三种矩形求积公式： 解： （1）将f(x)在x=a处进行泰勒展开有 两边同时在[a,b]上关于x积分，得 这是左矩形公式. 两边同时在[a,b]上关于x积分，得 这是右矩形公式. （2）将f(x)在x=b处进行泰勒展开有 两边同时在[a,b]上关于x积分，得 （3）将f(x)在x=(a+b)/2处进行泰勒展开有 由于 故 这是中矩形公式. 8. 试证明求积公式 证明： 当f(x)=1时 对于不高于5次的多项式准确成立，并计算积分 当f(x)=x5时 同理可验证，当f(x)= x , x2 , x3 , x4时，左边=右边， 但当f(x)=x6时 故求积公式具有5次代数精度. * * * 1. 设方程组： 解： 写出用Jacobi迭代法、Gauss-Seidel迭代法解此方程组的迭代格式，并讨论其收敛性. Jacobi迭代格式： G-S迭代格式： 因原方程组系数矩阵是严格对角占优矩阵，故Jacobi迭代和G-S迭代均收敛. 2. 设有方程组： 解： 试考察求解上述方程组的Jacobi迭代法及Gauss-Seidel迭代法的收敛性. （1）Jacobi迭代法的迭代矩阵为 解得 即ρ(BJ)1，故Jacobi迭代法解此方程组不收敛. Gauss-Seidel迭代法的迭代矩阵为 故Gauss-Seidel迭代法解此方程组收敛. （2）Jacobi迭代法的迭代矩阵为 故 故Jacobi迭代法收敛. Gauss-Seidel迭代法的迭代矩阵为 故 所以 所以Gauss-Seidel迭代法不收敛. 3. 设A=(aij)为二阶矩阵，且a11a22≠0，试证明求解方程组Ax=