数值计算方法I实验报告.docVIP

实 验 报 告 实验课程名称 数值计算方法I 开课实验室 数学实验室 学 院 理学院 年级 2012 专业班 信息与计算科学 2班 学 生 姓 名 学 号 开 课 时 间 2012 至 2013 学年第 2 学期 评分细则 评分 报告表述的清晰程度和完整性（20分） 程序设计的正确性（40分） 实验结果的分析（30分） 实验方法的创新性（10分） 总成绩 教师签名 实验一 误差分析 试验1.1(病态问题） 问题提出：考虑一个高次的代数多项式 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 其中是一个非常小的数。这相当于是对（1.1）中的系数作一个小的扰动。我们希望比较（1.1）和（1.2）根的差别，从而分析方程（1.1）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个MATLAB函数：“roots”和“poly”。 其中若变量a存储n+1维的向量，则该函数的输出u为一个n维的向量。设a的元素依次为，则输出u的各分量是多项式方程 的全部根；而函数 的输出b是一个n+1维向量，它是以n维向量v的各分量为根的多项式的系数。可见“roots”和“poly”是两个互逆的运算函数。 上述简单的MATLAB程序便得到（1.2）的全部根，程序中的“ess”即是（1.2）中的。 实验要求：（1）选择充分小的ess，反复进行上述实验，记录结果的变化并分析它们。如果扰动项的系数很小，我们自然感觉（1.1）和（1.2）的解应当相差很小。计算中你有什么出乎意料的发现？表明有些解关于如此的扰动敏感性如何？ （2）将方程（1.2）中的扰动项改成或其它形式，实验中又有怎样的现象？ （3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将方程（1.2）写成展开的形式， 同时将方程的解x看成是系数的函数，考察方程的某个解关于的扰动是否敏感，与研究它关于的导数的大小有何关系？为什么？你发现了什么现象，哪些根关于的变化更敏感？ 实验过程：打开matlab程序，在command window输入依次以下程序，按回车键即可。 程序： 以i=2为例 ess=0.000001; ve=zeros(1,21); ve(2)=ess; roots(poly(1:20)+ve) ess=0.0000001; ve=zeros(1,21); ve(2)=ess; roots(poly(1:20)+ve) ess=0 ve=zeros(1,21); ve(2)=ess; roots(poly(1:20)+ve) ess=0.000000001; ve=zeros(1,21); ve(2)=ess; roots(poly(1:20)+ve) ess=0.0000000001; ve=zeros(1,21); ve(2)=ess; roots(poly(1:20)+ve) 实验结果: 0.000001 0.0000001 00.000000001 0.0000000001 2 21.3025 + 1.5672i 21.3025 - 1.5672i 18.5028 + 3.6004i 18.5028 - 3.6004i 15.1651 + 3.7612i 15.1651 - 3.7612i 12.4866 + 2.8828i 12.4866 - 2.8828i 10.5225 + 1.7196i 10.5225 - 1.7196i 9.0448 + 0.5947i 9.0448 - 0.5947i 7.9489 7.0025 6.0000 5.0000 4.0000 3.0000 2.0000 1.0000 20.4220 + 0.9992i 20.4220 - 0.9992i 18.1572 + 2.4702i 18.1572 - 2.4702i 15.3147 + 2.6987i 15.3147 - 2.6987i