(完整word版)数值分析上机作业1-1解析.docx

数值计算方法上机题目 1 1、 1. 病 目的 ： 算法有“ ”与“劣”之分， 也有“好”和“坏”之 。所 坏 就是 本身 的解 数据 化的比 敏感， 反之属于好 。 希望 者通 本 此有一个初步的体会。 数 分析的大部分研究 中， 如 性代数方程 、 矩 特征 、 非 性方程及方 程 等都存在病 的 。 病 要通 研究和构造特殊的算法来解决， 当然一般要付出一些代价（如耗用更多的机器 、占用更多的存 空 等）。 提出 ： 考 一个高次的代数多 式 20 p( x) (x 1)( x 2)...( x 20) ( x k) （ E1-1） k 1 然 多 式的全部根 l， 2，?， 20，共 20 个，且每个根都是 重的（也称 的）。 考 多 式方程的一个 p( x) x19 0 (E1-2) 其中 是一个非常小的数。 相当于是 （ E1-1）中 x19 的系数作一个小的 。我 希望 比 （ E1-1）和（ E1-2）根的差 ，从而分析方程（ E1-1）的解 的敏感性。 内容 ： 了 方便，我 先介 两个 Matlab 函数：“ roots ”和“ poly ”， 入函数 u ＝ roots （ a ） 其中若 量 a 存 n 1 的向量， 函数的 出 u 一个 n 的向量。 a 的元素依次 a1 , a2 ,..., an 1 ， 出 u 的各分量是多 式方程 a x n a x n 1 ...a x a n 1 0 1 2 n 的全部根，而函数 b=poly(v) 的 出 b 是一个 n＋1 量，它是以 n 量 v 的各分量 根的多 式的系数。 可 “ roots” 和“ Poly ”是两个互逆的运算函数 . ve=zeros(1,21); ve(2)=ess; roots(poly(1:20))+ve) 上述 的 Matlab 程序便得到 （ E1-2）的全部根， 程序中的 “ ess”即是（ E1-2）中的 。 要求 ： （ 1） 充分小的 ess，反复 行上述 ， 果的 化并分析它 。如果 的系数 很小，我 自然感 （ E1-1）和 (E1-2) 的解 当相差很小。 算中你有什么出乎意料的 ？表明有些解关于如此的 敏感性如何？ （ 2）将方程（ E1-2）中的 改成 x18 或其他形式， 中又有怎 的 象出 ？ 1 实验步骤： （1）程序 function t_charpt1_1 clc result=inputdlg({ 请输入扰动项 : 在 [0 20] 之间的整数 : }, charpt 1_1 ,1,{ 19 }); Numb=str2num(char(result)); if ((Numb20)|(Numb0))errordlg( 请输入正确的扰动项 :[0 20] 之间的整 数! ); return ; end result=inputdlg({ 请输入 (0 1) 之间的扰动常数 : }, charpt 1_1 ,1,{ 0.00001 }); ess=str2num(char(result)); ve=zeros(1,21); ve(21-Numb)=ess; root=roots(poly(1:20)+ve); x0=real(root); y0=imag(root); plot(x0,y0, * ); disp([ 对扰动项 ,num2str(Numb), 加扰动 ,num2str(ess), 得到的全部根 为: ]); disp(num2str(root)); 二、实验结果分析 ess 分别为 1e-6,1e-8.1e-10,1e-12. 对扰动项 19 加扰动 1e-006 得到的全部根为 : 21.3025+1.56717i 21.3025-1.56717i 18.5028+3.6004i 18.5028-3.6004i 15.1651+3.76125i 15.1651-3.76125i 12.4866+2.88278i 12.4866-2.88278i 10.5225+1.71959i 10.5225-1.71959i 9.04485+0.594589i 9.04485-0.594589i 7.9489+0i 7.00247+0i 5.99995+0i 5+0i 4+0i 3+0i 2+0i 1+0i 对扰动项 19 加扰动 1e-010 得到的全部根为 : 19.9953+0i 19.0323+0i 17.8696+0i 17.2186+0i 15.4988+0.0211828i 15.