计算方法实验指导书12621.docVIP

《计算方法》 实 验 指 导 书 汪荣贵 编 计算机学院 计算机应用教研室 2006年2月 前 言 《计算方法》所涉及的都是可以直接在计算机上求解的数值方法。这些方法既有数学的抽象性与严密的科学性，又有应用的广泛性与实验的高度技术性。故结合课堂教学，本课程配备适当的上机实验(12个学时)，通过实验可以加强学生对数学模型的总体分析，算法选取，上机调试和结果分析等环节的训练。 本实验指导书共包含5个实验，要求学生在12个实验课时内完成，并提交实验报告,以此可作为《计算方法》课程成绩评定的一部分。 目 录 前 言 2 实验一 插值与拟合 4 1 实验目的 4 2 实验内容 4 3 算法基本原理 5 4 算法描述 6 5 计算用例的参考输出 6 6 思考题 6 实验二 数值积分 7 1 实验目的 7 2 实验内容 7 3 算法基本原理 8 4 算法描述 9 5 计算用例的参考输出 10 6思考题 10 实验三 非线性方程求根迭代法 11 1 实验目的 11 2 实验内容 11 3 算法基本原理 11 4 算法描述 12 5 计算用例的参考输出 12 6 思考题 12 实验四 求解线性方程组 13 1 实验目的 13 2 实验内容 13 3 算法基本原理 14 4 算法描述 15 5计算用例的参考输出 15 6 思考题 16 实验五 数值微分 17 1 实验目的 17 2 实验内容 17 3 算法基本原理 17 4 算法描述 18 5 计算用例的参考输出 18 6 思考题 19 实验报告内容要求 20 实验一 插值与拟合 1 实验目的 明确插值多项式和分段插值多项式各自的优缺点； 编程实现拉格朗日插值算法，分析实验结果体会高次插值产生的龙格现象； 理解最小二乘拟合，并编程实现线性拟合，掌握非线性拟合转化为线性拟合的方法 运用常用的插值和拟合方法解决实际问题。 2 实验内容 对于 要求选取11个等距插值节点，分别采用拉格朗日插值和分段线性插值，计算x为0.5, 4.5处的函数值并将结果与精确值进行比较。 输入：区间长度，n(即n+1个节点)，预测点 输出：预测点的近似函数值，精确值，及误差 （2）对于给定的一元函数 的n+1个节点值 。 数据如下： 0.4 0.55 0.65 0.80 0.95 1.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次拉格朗日多项式和分段三次插值多项式，计算 的值。 输入：数据点集及个数，预测点。 输出：预测点的近似函数值 （3）在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量与时间的拟合曲线。 分） 0 5 10 15 20 25 30 35 40 45 50 55 0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.02 4.64 输入： 数据点集及个数。 输出：拟合函数，并打印出与的误差 提示：近似解析表达式为 （4）用线性拟合编程实现给药方案 背景知识：一种新药用于临床之前，必须设计给药方案. 临床上，每种药物有一个最小有效浓度c1和一个最大有效浓度c2。 设计给药方案时，要使血药浓度 保持在c1~c2之间。 本题设c1=10，c2=25(ug/ml). 在实验方面，对某人用快速静脉注射方式一次注入该药物300mg后,在一定时刻t(小时)采集血药,测得血药浓度c(ug/ml)如下表: t (h) 0.25 0.5 1 1.5 2 3 4 6 8 c ((g/ml) 19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01 问题： 在快速静脉注射的给药方式下，研究血药浓度（单位体积血液中的药物含量）的变化规律 设计给药方案：每次注射剂量多大；间隔时间多长 提示： 血药浓度的变化规律符合指数关系 将指数关系转化为线性关系 3 算法基本原理 当精确函数 y = f(x) 非常复杂或未知，在一系列节点 x0 … xn 处测得函数值 y0 = f(x0), … yn = f(xn)， 希望由此构造一个简单易算的近似函数 g(x) ( f(x)，满足条件g(xi) = f(xi) (i = 0, … n)。这里的 g(x) 称为f(x) 的插值函数，由插值函数可以去近似估计f(x)在一些未知点处的函数值。最常用的插值函数是多项式插值。 所谓多项式插值即求 n 次多项