[工学]计算方法复习题.ppt

近似解的误差估计 设　　 的近似解为　 ，则一般有 cond (A) 残向量 注意：对于病态方程组，不能用残向量的大小衡量近似解的精度． 定义1. 因为一般情况下,求解公式的每一步都存在误差,因此有 定义2. 定义3. 一．２． 相对误差 相对误差限 绝对误差限 设 为精确值， 为近似值， * 总复习 一.填空题 1．对于较大的正数 ,计算 时为使计算 结果更准确应将其变形为 . 2．设 是某个数四舍五入的近似值，其有效 相对误差限为 数字的位数为 , . . 3．已知 , 则 , 4．已知函数 的数据表为 1.6 1.5 1.2 y 3 2 1 x 用复化梯形公式计算 用Simpson公式计算 . , 5．求解线性代数方程组 的迭代格式 充要条件是 收敛的 . 条件是其代数精度至少为 . 10 数值求积公式 是插值型 的充要 二． 确定常数 使得迭代法 局部收敛到 并有尽可能高的收敛阶。 解： 要 为 的根，则 要 ，则有 要 ，则有 将方程联立得 可验证 所以是三阶收敛。 三． 四． 用LU分解法求下列方程组 解 (1)求A的LU分解 五．方程组的系数矩阵为 问J迭代法当a为何值时收敛？ 下面方程组J迭代法与GS迭代法是否收敛？ J迭代法： Jacobi迭代法收敛! GS迭代法收敛？ Gauss-Seidel迭代法发散! 六 解: 且 七. 试确定下面积分公式中的参数使其代数精确度尽量高. 解: 因此 所以该积分公式具有3次代数精确度 八. 用改进的Euler法求解下面的初值问题 解 改进的Euler法 计算如下 九．已知求解 的公式如下： 试确定常数 ，使该公式具有尽可能高的精度。 解 求局部截断误差 即 此时公式为 方法阶至少为2。 是2阶精度的。 证明求微分方程 数值解的梯形公式 十. 证明: 越小，误差就越小。对吗？ 设 则 1. 不相容方程组 的最小二乘解为 2. 3.对于线性方程组 的近似解 , 十一. 设求积公式 为插值型求积公式 （1）试确定其求积系数 （2）试确定该公式的代数精度； （3）利用该公式计算 的近似值 十二. 解 Newton迭代格式为 十三.用Newton迭代法求方程 在0.5附近的根, 求迭代公式的局部收敛速度. 十四.用Seidel迭代法求解方程组 解： Seidel迭代格式为 因为系数矩阵严格对角占优，故Seidel　方 法对任意 收敛。 2范数为最小即 线性最小二乘拟合问题的提法 给定m+1对数据 (xi,yi),i=0,1,2,……,m 在选定的有限维 维函数空间Ф中求p(x)使偏差 如果函数 p*(x)能使偏差达到最小，则称函数 p*(x)为 最小二乘解！ 由此确定线性组合的系数，就可以求出最小二乘解！ 显然Q是关于线性组合的系数的多元函数： 它是关于 a0,a1,a2,……,an的一个二次齐次多项式函数，故 任意阶可导，又由于Q是正定的，从而必有最小值， 最小值点可由驻点求出： 它是关于 a0,a1,a2,……,an的一个二次齐次多项式函数，故 任意阶可导，又由于Q是正定的，从而必有最小值， 最小值点可由驻点求出： 即 这其实是一个关于 ak 方程组!它共有n+1个方程，n+1个未知量！ 我们称它为最小二乘问题的法方程组！ 引入向量 则它们都是向量空间Rm的向量，由内积的概念可知: 显然内积满足交换律 将其表示成矩阵形式 这是法方程组的内积形式！ 构造一个(m+1) ×(n+1) 矩阵： 则有： G的第j列就是向量ψj! 从而法方程组又可以写成： 这是法方程组的矩阵形