matlab 计算方法 实验指导 误差分析

《matlab 计算方法 实验指导 误差分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验一误差分析实验1（病态问题）实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。问题提出：考虑一个高次的代数多项式显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动其中是一个非常小的数。这相当于是对（1.

2、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的帮助。思考题二：（二进制产生的误差）用MATLAB计算。结果居然有误差！因为从十进制数角度分析，这一计算应该是准确的。实验反映了计算机内部的二进制本质。思考题三：（一个简单公式中产生巨大舍入误差的例子）可以用下列式子计算自然对数的底数这个极限表明随着n的增加，计算e值的精度是不确定的。现编程计算与exp(1)值的差。n大到什么程度的时候误差最大？你能解释其中的原因吗？相关MATLAB函数提示：poly(a)求

5、给定的根向量a生成其对应的多项式系数（降序）向量roots(p)求解以向量p为系数的多项式（降序）的所有根poly2sym(p)将多项式向量p表示成为符号多项式（降序）sym(arg)将数字、字符串或表达式arg转换为符号对象symsarg1arg2argk将字符arg1,arg2,argk定义为基本符号对象solve('eq1')求符号多项式方程eq1的符号解实验二非线性方程求根实验2（迭代法、初始值与收敛性）实验目的：初步认识非线性问题的迭代法与线性问题迭代法的差别，探讨迭代法及初始值与迭代收敛性的关系。问题提出：迭代法是求解非线性方程的基本思想方法，与线性方程的情况一

6、样，其构造方法可以有多种多样，但关键是怎样才能使迭代收敛且有较快的收敛速度。实验内容：考虑一个简单的代数方程针对上述方程，可以构造多种迭代法，如在实轴上取初始值x0，请分别用迭代（7.1）-（7.3）作实验，记录各算法的迭代过程。实验要求：（1）取定某个初始值，分别计算（7.1）-（7.3）迭代结果，它们的收敛性如何？重复选取不同的初始值，反复实验。请自选设计一种比较形象的记录方式（如利用MATLAB的图形功能），分析三种迭代法的收敛性与初值选取的关系。（2）对三个迭代法中的某个，取不同的初始值进行迭代，结果如何？试分析迭代法对不同的初值是否有差异？（3）线性方程组迭代法的

7、收敛性是不依赖初始值选取的。比较线性与非线性问题迭代的差异，有何结论和问题。思考题一：用Newton法求方程在区间[-3,3]上误差不大于的根，分别取初值进行计算，比较它们的迭代次数。相关MATLAB函数提示：x=fzero(fun,x0)返回一元函数fun的一个零点，其中fun为函数句柄，x0为标量时，返回在x0附近的零点；x0为向量[a,b]时，返回函数在[a,b]中的零点[x,f,h]=fsolve(fun,x0)返回一元或多元函数x0附近fun的一个零点，其中fun为函数句柄，x0为迭代初值；f返回fun在