计算方法实验指导书.docVIP

《计算方法》 实验指导书 数学实验室 2013年3月 实验一 病态问题 2 实验二 误差传播与算法稳定性 6 实验三 线方程组的直接解法 9 实验四 解线性方程组的迭代法 12 实验五 非线性方程求根 13 实验六 函数插值方法 14 提高：〖多项式插值的振荡现象〗 15 扩展：〖样条插值的收敛〗 19 实验七 数值积分与数值微分 23 实验八 常微分方程初值问题数值解法 24 实验九 函数逼近与曲线拟合 26 附实验一：〖多项式最小二乘拟合〗 27 附实验二：〖正交化多项式的最小二乘拟合〗 29 附录一：实验报告内容要求 32 附录二：部分程序示例 33 1、〖线方程组的直接解法系列程序〗 33 2、〖Jacobi迭代法解线性方程组程序〗 38 3、〖牛顿法非线性方程求根程序〗 39 4、〖多项式系列程序〗 40 5、〖Romberg算法计算数值积分程序〗 42 6、〖四阶经典的龙格-库塔方法解常微分方程程序〗 42 7、〖最小二乘法曲线拟合程序〗 43 8、〖幂法求矩阵特征值程序〗 46 前 言 结合课程教学，配备适当的上机实验(16个学时)以便加深课堂教学的实践性，同时通过实验可以加强对数学模型的总体分析，算法选取，程序结构，上机调试和结果分析等环节的训练。本实验指导书共包含8个实验，要求学生在16个实验课时内完成。为使实验更为有成效，需要写出实验报告（格式要求见附录）,以此可作为《计算方法》课程成绩评定的参考。 实验 病态问题 一、实验目的： 算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题是指问题本身对扰动敏感，反之属于好问题。本实验通过对一个高次多项式方程的求解，初步认识病态问题。 二、实验设备和实验环境 安装有C、C++或MATLAB的计算机。 三、实验内容及要求： 考虑一个高次的代数多项式 (1.1) 显然该多项式的全部根为1，2，…，20，共计20个，且每个根都是单重的(也称为简单的)。现考虑该多项式的一个扰动 (1.2) 其中，ε是一个非常小的数。这相当于是对方程(1)中x19的系数作一个小的扰动。比较方程(1.1)和方程(1.2)根的差别，从而分析方程(1.1)的解对扰动的敏感性。 四、实验步骤与结果分析 (一) 实验源程序 function t_charpt1_1 % 数值实验1病态问题 % 输入：[0 20]之间的扰动项及小的扰动常数 % 输出：加扰动后得到的全部根 clc result=inputdlg({请输入扰动项:在[0 20]之间的整数:},charpt 1_1,1,{19}); Numb=str2num(char(result)); if((Numb20)|(Numb0))errordlg(请输入正确的扰动项:[0 20]之间的整数!);return;end result=inputdlg({请输入(0 1)之间的扰动常数:},charpt 1_1,1,{0.00001}); ess=str2num(char(result)); ve=zeros(1,21); ve(21-Numb)=ess; root=roots(poly(1:20)+ve); x0=real(root); y0=imag(root); plot(x0,y0, *); disp([对扰动项 ,num2str(Numb),加扰动,num2str(ess),得到的全部根为:]); disp(num2str(root)); (二) 实验结果分析 对于x19项的扰动ess，不同的取值对应的结果如下所示。 对扰动项 19加扰动1e-010得到的全部根为: 19.9961，19.0257，17.9085，17.1508，15.7982，15.181，13.8995，13.0571，11.9753，11.0109，9.99608，9.00111，7.99978，7.00003，6，5，4，3，2，1。 对扰动项 19加扰动1e-009得到的全部根为: 19.952，19.2293，17.6573+0.692896i，17.6573-0.692896i，15.4524+0.875524i，15.4524-0.875524i，13.3527+0.486992i，13.3527-0.486992i，11.8578，11.0427，9.9916，9.00201，7.99952，7.00009，5.99999，5，4，3，2，1。 对扰动项 19加扰动1e-007得到的全部根为: 20.422+0.999203i，20.422-0.999203i，18.1572+2.4702i，18.1