河北工业大学数值分析实验要求.docVIP

实验一 舍入误差与数值稳定性 （2学时） 1.实验名称2.实验目的3.算法描述4.源程序 5.运行结果6.对算法的理解与分析（包括改进与建议） 程序与实例 例 1 对 n = 0,1,2,…,20 计算定积分 = 算法 1 利用递推公式 = - 5 n = 1,2,…,20 取 ln6- ln5 0.182 322 算法 2 利用递推公式 n = 20,19,…,1 注意到 取 实验二 拉格朗日插值与牛顿插值 （4学时） 目的与要求： 熟悉拉格朗日插值多项式和牛顿插值多项式，注意其不同特点； 实验内容： 通过拉格朗日插值和牛顿插值多项式的两个实例的计算，了解两种求解方法，分析他们的优缺点。 程序与实例 算法 输入x,y(i=0,1,2,(,n),令L(x)=0; 对=0,1,2,(,n计算 ll(x)= L( L+l(x)y 程序与实例 已知函数表 x 0.56160 0.56280 0.56401 0.56521 y 0.82741 0.82659 0.82577 0.82495 用三次拉格朗日多项式求x=0.5635的函数近似值。 牛顿插值多项式 算法 输入n,x,y (i=0,1,2(,n); 对k=1,2,3(,n, i=1,2, (,k计算各阶差商f(x,x (,x); 计算函数值 N(x)=f(x)+f[x, x](x- x)+(+f[x, x,(,x](x- x)(x- x)((x-x) 程序与实例 已知函数表 x 0.4 0.55 0.65 0.8 0.9 y 0.41075 0.57815 0.88811 1.02652 1.02652 用牛顿插值多项式求N(0.596)和N(0.895)。 实验三 复化辛卜生法，龙贝格法 （4学时） 目的与要求： 通过实际计算体会各种方法的精确度； 会编写用复化辛卜生、龙贝格算法求定积分的程序。 实验内容： 通过实际计算体会各种方法的精确度并且会编写用复化辛卜生、龙贝格算法求定积分的程序 程序与实例 复化辛卜生公式 算法：复化辛卜生公式为Sn=h/6,计算过程为： 1．令 2．对 计算 3． 。 程序与实例 例 用复化辛卜生法计算积分 运行结果为 s(2)=0.785392 s(4)=0.785398 s(8)=0.785398 说明：本例运行了三次，当时，就与时有6位数字相同，若用复化梯形法计算，当n=512时有此结果。 龙贝格算法计算 算法 用事后估计法控制精度 。 实验四 改进欧拉法，二分法，牛顿法 （4学时） 目的与要求： 熟悉求解常微分方程初值问题的有关方法和理论，主要是改进欧拉法 会编制上述方法的计算程序 针对实习题编制程序，并上机计算其所需要的结果； 实验内容： 熟悉求解常微分方程初值问题的有关方法和理论，主要是改进欧拉法,体会其解法的功能。程序与实例 改进欧拉方法 算法概要 解一阶常微分方程初值问题 将区间[a,b]作n等分，取步长。 欧拉公式为 梯形公式为 改进欧拉法，采用公式 或表为 实验题：  一、目的与要求： 通过对二分法和牛顿迭代法作编程练习和上机运算，进一步体会它们在方程求根中的不同特点； 比较二者的计算速度和计算精度。 实验内容： 通过对二分法和牛顿迭代法作编程练习和上机运算，进一步体会它们在方程求根中的不同特点 程序与实例 二分法 算法：给定区间[a,b]，并设与符号相反，取为根的容许误差，为的容许误差。 令c=(a+b)/2 如果(c-a)或，则输出，结束；否则执行（3）， 如果，则令；否则则令，重复（1），（2），（3）。 书上课后习题1。 牛顿迭代法 算法：给定初值，为根的容许误差，为的容许误差，N为迭代次数的容许值。 如果=0或迭代次数大于N，则算法失败，结束；否则执行（2）。 计算=- 若或，则输出，程序结束；否则执行（4）。 令=，转向（1）。 书上课后习题7的（1）。 实验五 矩阵的LU分解法，雅可比迭代 目的与要求： 熟悉求解线性方程组的有关理论和方法； 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序； 通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。 实验内容： 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序，进一步了解各种方法的优缺点。 程序与实例 列主元高斯消去法 算法：将方程用增广矩阵[A∣b]=(表示 消元过程 对k=1,2,…,n-1 ①选主元，找使得 = ②如果，则