数值分析实验报告大学论文.docVIP

数值分析实验报告 学 院： 学生姓名： 学号： 2016年11月 目录 实验一 插值法 1 一 实验内容 1 二 实验步骤 1 三 所编程序清单 2 四 例题验证 3 实验二 数据拟合 5 一 实验内容 5 二 实验步骤 5 三 难点分析 6 四 所编算法清单 6 五 例题验证 6 实验三 数值积分 8 一 实验内容 8 二 实验步骤 8 三 难点分析 8 四 所编算法清单 9 五 例题验证 10 实验四 线性方程组的直接求法之高斯消去法 12 一 实验内容 12 二 实验步骤 12 三 难点分析 13 四 所编算法清单 13 五 例题验证 13 实验五 线性方程组的迭代解法之雅可比迭代法 15 一 实验内容 15 二 实验步骤 15 三 难点分析 16 四 所编算法序清单 16 四 例题验证 17 实验六 非线性方程的数值解法之Newton迭代法 19 一 实验内容 19 二 实验步骤 19 三 实验难点 19 四 所编程序详细清单 19 五 例题验证 20 实验七 矩阵特征值问题的数值计算之幂法 21 一 实验内容 21 二 实验步骤 21 三 难点分析 21 四 具体所编算法清单 21 五 例题验证 22 实验八 常微分方程的数值解法之改进的Euler法 24 一 实验内容 24 二 实验步骤 24 三 实验难点 24 四 具体所编算法 24 五 例题验证 25 实验一 插值法 本人选择Newton插值法进行编程。选择的例题为：教材（P25，例2.7）进行验证： 题目：设，并已知的数据如表1-1试用2次Newton插值多项式计算计算的近似值，并讨论其误差。 表1-1 x 2.0 2.1 2.2 1.414214 1.449138 1.483240 一 实验内容 根据Newton插值多项式 (*) 使用Matlab编程实现Newton插值算法。 2.使用教材上的例题验证所编算法的正确性。 二 实验步骤 Newton插值法的关键是求出的k阶均差，然后将k阶均差代入上式(*)，求出插值多项式，最后将所选例题的条件代入Newton插值多项式，与正确答案进行比对，进而得出结论。 1 难点分析 在编写程序时需要考虑的难点之一是如何确定所选用的Newton插值法的次数，本人拟选用直接从Matlab中m文件中直接人工输入的方法，设一个未知量n作为选用Newton插值法的次数确定选用几次的Newton插值多项式。 在编写程序时的另一个难点是经过Newton插值公式计算之后，输出是什么？是直接列写多项式，还是以系数矩阵的形式输出？本人考虑到编程的简便性，拟直接以均差表形式输出。 2 具体步骤 详细步骤如下： (1)已知k阶均差的求解公式如下： k=1,2,....; (1) 通过上式分别求取一阶均差、二阶均差、...、k阶均差，并保存在与函数自变量x同阶的同一个大矩阵中。在大矩阵中，第一列用来存放自变量x，第二列用来存放因变量，第三列用来存放一阶均差，...，第（k+2）列用来存放k阶均差。具体矩阵表示如下表1-2所示 表1-2 一阶均差 二阶均差 三阶均差 ... ... ... ... ... 注：表1-2未标示出的元素都以0填充。 (2)根据公式(*)求解n次(n为人为设定的Newton插值多项式系数)Newton插值多项式，并且在Matlab中Command Window输出该插值多项式。 (3)在教材中随机选取某一例题来验证所编程序的正确性。 三 所编程序清单 所编程序对应的文件为Newton.m %Newton插值法 function [y]=Newton(Xi,Yi,x,n) %%------------------------------------------------- %计算均差表 %根据教材，均差表一般列举到三阶均差即可 %一阶均差 f1=zeros(1,n);%一阶均差数值保存在f1中 for i=1:n f1(1,i)=(Yi(1,i+1)-Yi(1,i))/(Xi(1,i+1)-Xi(1,i)); end disp(一阶均差数值为：); disp(f1); %二阶均差 f2=zeros(1,n-1);%二阶均差数值保存在f2中 for i=1:(n-1) f2(1,i)=(f1(1,i+1)-f1(1,i))/(Xi(1,i+2)-Xi(1,i)); end disp(二阶均差数值为：); disp(f2); %三阶均差 if n2 f3=zeros(1,n-2);%三阶均差数值保存