数值分析实验应用数学研究生..docVIP

数值分析实验题 理学院 应用数学研究生 实验1 病态问题与数值稳定性 实验1.1 病态问题 实验目的：算法有“优”与“劣”之分，问题也有“好”和“坏”之别.对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题.希望读者通过本实验对此有一个初步的体会. 数值分析的大部分研究课题中，如线性代数方程组，矩阵特征值问题，非线性方程及方程组等都存在病态的问题.病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间，占用更多的存储空间等） 问题提出：考虑一个高次的代数多项式 ， （1） 显然该多项式的全部根为1,2，，20共计20个，且每个根都是单重的（也称为简单的）.现考虑该多项式的一个扰动 , （2） 式中，是一个非常小的数.这相当于是对（1）中的系数作一个小的扰动.希望比较式（1）和（2）根的差别，从而分析方程（1）的解对扰动的敏感性. 实验内容：为了实现方便，介绍两个Matlab函数：roots和poly. u=roots(a)，其中若变量a存储n+1维的向量，则该函数的输出u为一个n维的向量.设a的元素依次为，则输出u的各分量是多项式方程 的全部根；而函数b=poly(v)的输出b是一个维变量，它是以维变量的各分量为根的多项式系数.可见roots和poly是两个互逆的运算函数. ve=zeros(1,21); ve(2)=ess; roots(poly(1:20)+ve) 上述简单的Matlab程序便得到（2）的全部根，程序中的“ess”即是（2）中的. 实验要求： （1）选择充分小的ess，反复进行上述实验，记录结果的变化并分析它们.如果扰动项的系数很小，我们自然感觉式（1）和式（2）的解应当相差很小.计算中你有什么出乎意料的发现？表明有些解关于如此的扰动敏感性如何？ （2）将方程（2）中的扰动项改成或其他形式，实验中又有怎样的现象出现？ （3）请从理论上分析产生这一问题的根源.注意我们可以将方程（2）写成展开的形式 . （3） 同时将方程的解看成是系数的函数，考察方程的某个解关于的扰动是否敏感，与研究它关于的导数的大小有何关系？为什么？你发现了什么现象，哪些根关于的变化更敏感？ 实验1.2 误差传播与算法稳定性 实验目的：体会稳定性在选择算法中的地位，误差扩张的算法是不稳定的，是我们所不期望的；误差衰竭的算法是稳定的，是我们努力寻求的，这是贯穿本课程的目标. 问题提出：考虑一个简单的由积分定义的序列 . 显然利用分部积分易得, . （1） 又有 . 实验内容：由递推关系（1），可以得到计算积分序列{}的两种算法 算法1： (2) 算法2： (3) 实验要求： （1）分别用算法1、算法2并在计算中分别采用5位、6位和7位有效数字，请判断哪种算法能给出更精确的结果. （2）两种算法的优劣，与你的第一感觉是否吻合.请从理论上证明你实验得出的结果，解释实验的结果.设算法1中的计算误差为,由递推计算的误差为；算法2中的误差为，由向前递推计算误差为.如果在上述两算法中都假定后面的计算不再引入其他误差，试给出与的关系和的关系. （3）算法1中通常会很小，当n增大时，的变化趋势如何？算法2中通常相对比较大，当n减小时，误差又是如何传播的？也就是说比较一下上述两个算法，当某一步产生误差后，该误差对后面的影响是衰减还是扩张的. （4）通过理论分析与计算实验，针对算法1和算法2的稳定性，给出你的结论. 实验2 非线性方程（组） 实验2.1 迭代法、初始值与收敛性 实验目的：初步认识非线性问题的迭代法与线性问题迭代法的差别，探讨迭代法及初始值与迭代收敛性的关系. 问题提出：迭代法是求解非线性方程及方程组的基本思想方法，与线性方程的情况一样，其构造方法可以有多种多样，但关键是怎样才能使迭代收敛且有较快的收敛速度. 实验内容：考虑一个简单的代数方程可以构造多种迭代法，如 (1) (2) (3) 在实轴上取初始值,请分别用迭代（1）～（3）做实验，记录各算法的迭代过程. 实验要求： （1）取定某个初始值，分别计算（1）～（3）迭代结果，收敛性如何？重复选取不同的初始值，反复实验，请读者自行设计一种比较