数值分析实验报告’.docVIP

实验一 误差分析 实验（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。 对数值方法的研究而言， 所谓坏问题就是问题本身对扰动敏感者， 反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。 病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价 （如耗用更多的机器时间、 占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 20 p(x) ( x 1)( x 2) (x 20) (x k) (1.1) k 1 显然该多项式的全部根为 1,2, ,20 共计 20 个，且每个根都是单重的。现考虑该多项式的一个扰动 p( x) x19 0 (1.2) 其中 是一个非常小的数。这相当于是对（）中 x19 的系数作一个小的扰动。我 们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。 实验内容： 为了实现方便，我们先介绍两个 Matlab 函数：“roots ”和“ poly ”。 u roots(a) 其中若变量 a 存储 n+1 维的向量，则该函数的输出 u 为一个 n 维的向量。设 a 的元素依次为 a1, a2 , ,an 1 ，则输出 u 的各分量是多项式方程 a1 xn a2 xn 1 an x an