计算方法实验报告格式.docVIP

计算方法上机实验指导书(06-07上) PAGE PAGE 1 计算方法上机实验指导 一、数值实验报告格式及要求： 1、实验目的： 首先要求每一个做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出。 2、实验题目： 在下面分若干个实验详细给出，实验者可根据报告形式需要适当改写或重述。 3、实验原理与基本理论： 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及的理论基础，算法原理详尽列出。 4、实验内容： 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法流程图等。 5、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格形式实现，较为简单的结果可以与实验结果分析合并出现。 6、实验结果分析： 实验结果分析是数值实验的重要环节，只有对实验结果认真分析，才能对实验目的、实验方法进一步理解对实验的重要性充分认识，明确数值计算方法的使用范围及其优缺点。 7、实验体会 要求：每个实验都应在计算机上实现或演示，由实验者独立用Matlab语言编程实现 (使之尽量具有通用性)，程序清单以附录形式给出，程序中至少1/3行要加注释，特别要对程序中的主要变量给出说明。 二、数值实验类型：（可以根据课时安排在讲对应章节之后完成实验） 实验一 误差传播与算法稳定性 1.1 实验目的 体会稳定性在选择算法中的地位．误差扩张的算法是不稳定的，是我们所不期望的；误差衰竭的算法是稳定的．是我们努力寻求的，这是贯穿本课程的目标．通过上机计算，了解舍入误差所引起的数值不稳定性。 1.2 算法描述 概要：舍人误差在计算方法中是—个很重要的概念。在实际计算中，如果选用了不同的算法，由于舍入误差的影响．将会得到截然不同的结果。因此，选取稳定的算法，在实际计算中是十分重要的。 考虑一个简单的由积分定义的序列 ， (E.1) 利用分部积分易得 , 得递推公式 ， (E.2) 注意到 取 由利用(E.2)变形得到 (E.3) 计算方法：先估计一个,再反推要求的 ( )。 可取 （）。 1.3 实验内容 由递推关系(E.2)，可以得到计算(E.1)积分序列的两种算法．其一为(E.2)的直接应用，即 ， ， (E.4) 另外一种算法则是利用(E.2)变形得到的 ， ， (E.5) 1.4 实验要求 （1）分别用算法(E.4)、(E.5)并在计算中分别采用5位、6位和7位有效数字，请判断哪种算法能给出更精确的结果． （2）两种算法的优劣，与你的第一感觉是否吻合．请从理论上证明你实验得出的结果，解释实验的结果．设(E.4)中的计算误差为，由递推计算的误差为；算法(E.5)中的误差为，由向前递推计算的误差为．如果在上述两算法中都假定后面的计算不再引入其他误差，试给出与的关系和与的关系． （3）算法(E.4)中通常会很小，当增大时，的变化趋势如何?算法(E.5)中通常相对比较大，当减小时，误差又是如何传播的?也就是说比较一下上述两个算法，当某一步产生误差后．该误差对后面的影响是衰减还是扩张的． （4）通过理论分析与计算实验，针对算法(E.4)和(E.5)的稳定性，给出你的结论。 实验二 列主元高斯消去法 2.1 实验目的 掌握列主元高斯消去法的基本思路和迭代步骤； 培养编程与上机调试能力。 2.2 算法描述 2.2.1 高斯（列主元）消去法基本思路 设有方程组，设是可逆矩阵。高斯消去法的基本思想就是将矩阵的初等行变换作用于方程组的增广矩阵，将其中的变换成一个上三角矩阵，然后求解这个三角形方程组。列主元高斯消去法的基本思想是每次在系数矩阵中依次按列在主对角线及以下的元素中，挑选绝对值最大的元素作为主元，将它调至主对角线上，然后用它消去主对角线以下的元素，最后变为同解的上三角形方程组去求解。 2.2.2 列主元高斯消去法计算步骤 将方程组用增广矩阵表示。 步骤1：消元过程，对 选主元，找使得 如果，则矩阵奇异，程序结束；否则执行（3）。 如果，则交换第行与第行对应元素位置，，。 消元，对，