《数值分析实验报告110560》由会员分享，可在线阅读，更多相关《数值分析实验报告110560（27页珍藏版）》请在金锄头文库上搜索。
1、数 值 分 析 实 验 报 告 实验一 误差分析实验1.1（病态问题）实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。问题提出：考虑一个高次的代数多项式显然该多项式的全部根为1,2,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动其中是一个非常小的数。这相
2、当于是对（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，反复
3、进行上述实验，记录结果的变化并分析它们。如果扰动项的系数很小，我们自然感觉（1.1）和（1.2）的解应当相差很小。计算中你有什么出乎意料的发现？表明有些解关于如此的扰动敏感性如何？（2） 将方程（1.2）中的扰动项改成或其它形式，实验中又有怎样的现象出现？（3） （选作部分）请从理论上分析产生这一问题的根源。注意我们可以将方程（1.2）写成展开的形式，同时将方程的解x看成是系数的函数，考察方程的某个解关于的扰动是否敏感，与研究它关于的导数的大小有何关系？为什么？你发现了什么现象，哪些根关于的变化更敏感？思考题一：（上述实验的改进）在上述实验中我们会发现用roots函数求解多项式方程的精度不高，
4、为此你可以考虑用符号函数solve来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab的帮助。实验过程：程序：a=poly(1:20);rr=roots(a);for n=2:21 n for m=1:9 ess=10(-6-m);ve=zeros(1,21); ve(n)=ess;r=roots(a+ve); -6-m s=max(abs(r-rr) endend利用符号函数：（思考题一）a=poly(1:20);y=poly2sym(a);rr=solve(y)for n=2:21 n for m=1:8 ess=10(-6-m)
5、; ve=zeros(1,21); ve(n)=ess; a=poly(1:20)+ve; y=poly2sym(a); r=solve(y); -6-m s=max(abs(r-rr) endend数值实验结果及分析：27format long-6-m n-7-8-9-1022.797226874783311.867536320091581.060527623807480.2527314421904731.693766997674240.923106667069640.084716145697410.4080402640941140.854013934155360.1994102202006
6、10.03972935295834050.110311005388710.042965323628440060000700008000090000100000110000120000130000140000150000160000170000180000190000200000210000-6-m n-11-12-13-1420.038776764393800.162565848682800.13322664013598030.021642583175460004000050000600007000080000900001000001100001200001300001400001500001
7、60000170000180000190000200000210000讨论：利用这种方法进行这类实验，可以很精确的扰动敏感性的一般规律。即当对扰动项的系数越来越小时，对其多项式扰动的结果也就越来越小，即扰动敏感性与扰动项的系数成正比，扰动项的系数越大，对其根的扰动敏感性就越明显，当扰动的系数一定时，扰动敏感性与扰动的项的幂数成正比，扰动的项的幂数越高，对其根的扰动敏感性就越明显。实验总结： 利用MATLAB来进行病态问题的实验，虽然其得出的结果是有误差的，但是可以很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动，对其多项式的根会有一定的扰动的，所以对于这类病态问题可以借助于MATL
8、AB来进行问题的分析。学号：06450210姓名：万轩实验二 插值法实验2.1（多项式插值的振荡现象）问题提出：考虑一个固定的区间上用插值逼近一个函数。显然拉格朗日插值中使用的节点越多，插值多项式的次数就越高。 我们自然关心插值多项式的次数增加时，L(x)是否也更加靠近被逼近的函数。龙格给出了一个极著名例子。设区间-1，1上函数 f(x)=1(1+25x2)实验内容：考虑区间-1，1的一个等距划分，分点为： x(i)=-1+2i/n,i=0,1,2,n泽拉格朗日插值多项式为： L(x)=l(i)(x)/(1+25x(j)2 ) i=0,1,n其中l(i)(x), i=0,1,n,n是n次拉格朗
9、日插值基函数。实验要求： 选择不断增大的分点数目n=2,3，画出f(x)及插值多项式函数L(x)在-1，1上的图象，比较分析实验结果。（2）选择其它的函数，例如定义在区间-5，5上的函数 h(x)=x/(1+x4) , g(x)=arctanx 重复上述的实验看其结果如何。 （3）区间a,b上切比雪夫点的定义为： xk=（b+a）/2+(b-a)/2)cos(2k-1)/(2(n+1),k=1,2,n+1 以x1,x2x(n+1)为插值节点构造上述各函数的拉格朗日插值多项式，比较其结果。实验过程：程序：多项式插值的震荡现象（实验2.1）for m=1:6 subplot(2,3,m) %把窗口
10、分割成2*3大小的窗口 largrang(6*m) %对largrang函数进行运行 if m=1 title(longn=6) elseif m=2 title(longn=12) elseif m=3 title(longn=18) elseif m=4 title(longn=24) elseif m=5 title(longn=30) elseif m=6 title(longn=36) end %对每个窗口分别写上标题为插值点的个数end保存为：chazhi.mfunction largrang(longn)mm=input(please input mm(运行第几个函数就输入mm为
11、几):mm=)if mm=1 %d表示定义域的边界值 d=1;elseif mm=2|mm=3 d=5;endx0=linspace(-d,d,longn); %x的节点if mm=1 y0=1./(1.+25.*x0.2);elseif mm=2 y0=x0./(1.+x0.4);elseif mm=3 y0=atan(x0);endx=sym(x);n=length(x0); s=0.0;for k=1:n p=1.0; for j=1:n if j=k p=p*(x-x0(j)/(x0(k)-x0(j); end end s=p*y0(k)+s;endy=s;if mm=1 ezplot
12、(1/(1+25*x2)elseif mm=2 ezplot(x/(1+x4)elseif mm=3 ezplot(atan(x)endhold onezplot(y,-d,d)hold off保存为：largrang.m数值实验结果及分析：对于第一个函数f(x)=1/(1+25x2)对于第二个函数h(x)=x/(1+x4)对于第三个函数g(x)=arctan(x)讨论：通过对三个函数得出的largrang插值多项式并在数学软件中的运行，得出函数图象，说明了对函数的支点不是越多越好，而是在函数的两端而言支点越多，而largrang插值多项式不是更加靠近被逼近的函数，反而更加远离函数，在函数两端的跳动性更加明显，argrang插值多项式对函数不收敛。实验总结：利用MATLAB来进行函数的largrang插值多项式问题的实验，虽然其得出的结果是有误差的，但是增加支点的个数进行多次实验，可以找出函数的largrang插值多项式的一般规律，当支点增加时，largrang
