课程实验及其解决方案.docVIP

?数值计算?实验指导书 要求：　　 用Matlab语言或你熟悉的其它算法语言语言编程序，使之尽量具有通用性． 上机前兖分准备，复习有关算法，写出计算步骤，反复查对程序． 完成计算后写出实验报告报告，内容包括：计算机型号及所用算法语言，CPU时间，结果分析和小结等． 熟练掌握Marlab的使用． 实验－ 实验目的：掌握Matlab在数值计算求解中的基本方法． 实验要求：会使用Matlab工具．并能用它求解以下问题问题． 分别用库函数quad8和quad求以下定积分： 库函数eig求下列矩阵的特征值和特征向量： 作下列函数在指定区间上的图形，并用库函数roots求它的根： 实验二 实验目的：掌握三次样条插值的三弯矩法． 实验要求：运用三弯矩法编制程序求下列问题的解． 实验题： 设已知数据如下： 0.2 4 0.6 0.8 1.0 0.9798652 0.9177710 0.8080348 0.6386093 0.3843735 求的三次样条插值函数 ，满足： 自然边界条件　 第一种边界条件　要求输出用追赶法解出的弯矩向量和的值，并画出的图形． 　　　　　　　　　　　　实验三 实验目的：复化求积公式及高斯－勒让德公式计算定积分. 实验题： （1） (2) 实验要求： 分别用复化梯形公式，复化Simpson公式求解，要求绝对误差限为． 用三点的五点的Gauss-Legendre求积公式求解上述积分，并输出I的近似值． 分析比较各种计算结果． 实验四 实验目的：迭代法求解非线性方程的根． 实验题： 　　　　　　　 实验要求： 用你自己设计的一种线性收敛迭代法求方程的根，然后用斯蒂芬森加速迭代重新计算． 用牛顿法求方程的根，输出迭代初值．各次迭代值及迭代次数，并与（1）比较． 实验五 实验目的：观察的理解高斯消元过程中出现小主元即很小时，引起方程组解数值不稳定性． 实验题：求解线性方程组． 　　　　， 实验要求： 计算矩阵的条件数，判断系数矩阵是良态的还是病态的． 用高斯列主元消去法求得L和U及解向量． 用不选主元的高斯消去法求得L和U及解向量． 观察小主元并分析对计算结果的影响． 实验六 实验目的：认识迭代法收敛的含义以及迭代初值和方程组系数矩阵性质对收敛速度的影响． 实验题：用迭代法求解方程组，它的每条对角线元素是常数． 实验要求： 选取不同的初始向量和不同的方程组右端向量，给定迭代误差要求．用Jacobi和Gauss-Seidel迭代法计算，观测得到的迭代向量序列是否收敛，若收敛，记录迭代次数，分析计算结果并得出你的结论． 取定右端向量和初始向量，将的主对角线元素成倍增长若干次，非主对角线元素不变，每次用Jacobi迭代法计算，要求迭代误差满足，比较收敛速度，分析现象并得出你的结论． 实验七 实验目的：求初值问题的常微分方程的数值解 实验题： 给定初值问题为 实验要求： 用改进欧拉法（取h=0.05）求它的数值解，并输出的数值解 用四阶R-K方法（取h=0.1）求它的数值解，并输出的数值解 分析结果，对两种算法进行比较． 综合实验 实验目的：由实际问题建立数学模型，然后会用相应的数值算法求解． 实验题：生态环境学家在研究自然界中两个生物种群数目变化时得到一组常微分方程。 假设有两种生物（例如一种是蓝鲸，另一种是南极磷虾），前者在时刻时的数量为()，后者在时刻时的数目为()，并假设它们都是的连续可微函数。蓝鲸是以磷虾为主要食物的。当没有食物来源时蓝鲸数目会减少，其减少速度与当时蓝鲸的数目成线性关系，即 （1） . 当有食物来源时，蓝鲸的数目会增加。增加的速度和它捕食的数目有关，即 = d() () （2） . 合并（1）和（2），得到蓝鲸变化速度满足的微分方程 d() (). (3) 同样，在没有蓝鲸时，磷虾的增加速度满足 =(). (4) 考虑到被捕食情况，则磷虾的数目满足 =()－b() (). (5) 合并（3