数值分析实一.docVIP

实验一 误差分析 实验1.1（病态问题） 实验目的： 算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 实验内容： 考虑一个高次的代数多项式 （1.1） 显然该多项式的全部根为1,2,,,,,,20共计20个，且每个根都是单重的。考虑多项式的一个扰动 （1.2） 其中是一个非常小的数，这相当于是对（1.1）中的系数作一个小的扰动。我们希望比较（1.2）和（1.2）根的差别，从而分析方程（1.1）的解对扰动的敏感性。 实验要求： 选择充分小的ess，反复进行上述实验，记录结果的变化并分析； 将方程中的扰动项改成或其他形式，试验中又有怎样的现象出现？ 请从理论上分析产生这一问题的根源。注意我们可以将方程（1.2）写成展开的形式 （1.3） 同时将方程的解x看成是系数的函数，考察方程的某个解关于的扰动是否敏感，与研究它关于的倒数的大小有何关系？为什么？你发现了什么现象，那些跟关于的变化更敏感 实验步骤： 选择充分小的ess，反复进行上述实验，分析结果变化； 对不同扰动项进行试验，比较所得结果； 理论分析，并与实验结果进行对比。 程序代码： function charpt1_1 result=inputdlg({请输入正确的扰动项：[0 20]之间的整数}, charpt1_1,1,{19}); Numb=str2num(char(result)); if((Numb20)|(Numb0))errordlg(请输入正确的扰动项：[0 20]之间的整数！); return; end result=inputdlg({请输入（0 1）之间的扰动常数：}, charpt1_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)); 实验结果： 扰动项为时，加不同的扰动所得的根如下表所示 表1 扰动项为时，加不同的扰动所得的根 ess 1e-012 1e-010 1e-006 0.001 20 19.9952+0i 21.3025+1.56717i 27.0817+5.03812i 19.0001 19.033+0i 21.3025-1.56717i 27.0817-5.03812i 17.999 17.8638+0i 18.5028+3.6004i 19.5337+9.1664i 17.0053 17.2299+0i 18.5028-3.6004i 19.5337-9.1664i 15.984 15.5044+0.155574i 15.1651+3.76125i 13.8235+7.77167i 15.0303 15.5044-0.155574i 15.1651-3.76125i 13.8235-7.77167i 13.9577 13.6984+0i 12.4866+2.88279i 10.7211+5.4609i 13.0419 13.2332+0i 12.4866-2.88279i 10.7211-5.4609i 11.9701 11.9006+0i 10.5225+1.71964i 8.91282+3.47317i 11.0166 11.05+0i 10.5225-1.71964i 8.91282-3.47317i 9.99354 9.98287+0i 9.04482+0.59469i 7.69268+1.89884i 9.00187 9.00502+0i 9.04482-0.59469i 7.69268-1.89884i 7.99962 7.99893+0i 7.94886+0i 6.75761+0.65471i 7.00005 7.00017+0i 7.00247+0i 6.75761-0.65471i 6 5.99998+0i 5.99995+0i 5.95208+0i 5 5+0i 5+0i 5.00061+0i 4 4+0i 4+0i 4+0i 3 3+0i 3+0i 3+0i 2 2+0i 2+0i 2+0i 1 1+0i 1+0i 1+0i 根在复平面上的位置如下图所示：