《数值分析》课程设计实验指导书.docVIP

数值分析实验指导书 考核标准： 及格：独立完成12—15题，其中八组实验中每组至少做1题； 中： 独立完成16—23题，其中八组实验中每组至少做1题； 良： 独立完成24—31题，其中八组实验中每组至少做2题； 优： 独立完成32—40题，其中八组实验中每组至少做3题。 结束课程时，抽查上机考核。 实验一 1.1 水手、猴子和椰子问题：五个水手带了一只猴子来到南太平洋的一个荒岛上，发现那里有一大堆椰子。由于旅途的颠簸，大家都很疲惫，很快就入睡了。第一个水手醒来后，把椰子平分成五堆，将多余的一只给了猴子，他私藏了一堆后便又去睡了。第二、第三、第四、第五个水手也陆续起来，和第一个水手一样，把椰子分成五堆，恰多一只猴子，私藏一堆，再去入睡，天亮以后，大家把余下的椰子重新等分成五堆，每人分一堆，正好余一只再给猴子，试问原先共有几只椰子？ 试分析椰子数目的变化规律，利用逆向递推的方法求解这一问题。 1.2 当时，计算积分. 1.3 绘制静态和动态的Koch分形曲线 问题描述：从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的另两条边代替，形成具有5个结点的新的图形；在新的图形中，又将图中每一直线段中间的三分之一部分都用一个等边三角形的另两条边代替，再次形成新的图形，这时，图形中共有17个结点。这种迭代继续进行下去可以形成Koch分形曲线。在迭代过程中，图形中的结点将越来越多，而曲线最终显示细节的多少取决于所进行的迭代次数和显示系统的分辨率。Koch分形曲线的绘制与算法设计和计算机实现相关。 图1.1 Koch曲线的形成过程 2.1 小行星轨道问题：一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，在五个不同的对小行星作了五次观察，测得轨道上五个点的坐标数据（单位：万公里）如下表所示： P1 P2 P3 P4 P5 X坐标 53605 58460 62859 66662 68894 Y坐标 6026 11179 16954 23492 68894 由开普勒第一定律知，小行星轨道为一椭圆，椭圆的一般方程可表示为： 现需要建立椭圆的方程以供研究。 (1) 分别将五个点的数据代入椭圆一般方程中，写出五个待定系数满足的等式，整理后写出线性方程组 以及方程组的系数矩阵和右端项b； (2) 用MARLAB求低阶方程的指令A＼b求出待定系数; (3) 分别用直接法、Jacobi迭代法、Gauss-Seidel迭代法求出待定系数. 2.2 (1) 用Gauss列主元消去法、Gauss按比例列主元消去法、Cholesky分解求解线性方程组，并验证(2) 判断用Jacobi 迭代法Gauss-Seidel迭代法SOR法(分别取)解下列线性方程组. 若收敛，再用Jacobi 迭代法Gauss-Seidel迭代法SOR法(分别取)分别解线性方程组并. (3) 用Cholesky分解求解线性方程组 （方程组的精确解是） 2.3 设 (1) 将矩阵A进行LU分解A=LU，其中U是上三角矩阵，L是主对角线上的元素都是1的下三角矩阵。 (2) 利用上述分解分别求解方程组 并由此求出逆矩阵。 (3) 用LU分解求下列线性方程组的解 （方程组的精确解是） 2.4 (1) 用追赶法求解方程组 (a) ， (b) ， (2) 设计实验验证Hilbert矩阵的病态性. 实验三 3.1 用Taylor级数的第项多项式，分别在区间和上逼近正弦函数，并用计算机绘出上面个函数的图形。 .3 追赶曲线的计算机模拟 问题描述：欧洲文艺复兴时期的著名人物达·芬奇曾经提出一个有趣的“狼追兔子”问题，当一只兔子正在它的洞穴南面60码处觅食时，一只饿狼出现在兔子正东的100码处。兔子急忙奔向自己的洞穴，狼立即以快于兔子一倍的速度紧追兔子不放。兔子一旦回到洞穴便逃脱厄，问狼是否会追赶上兔子？ 这一问题的研究方法可以推广到如鱼雷追击潜艇、地对空导弹击飞机等问题上去。 在对真实系统做实验时，可能时间太长、费用太高、危险太大、甚至很难进行。计算机模拟是用计算机模仿实物系统，对系统的结构和行为进行动态演示，以评价或预测系统的行为效果。根据模拟对象的不同特点，分为确定性模拟和随机性模拟两大类。模拟通常所用的是时间步长法，即按照时间流逝的顺序一步一步对所研究的系统进行动态演示，以提取所需要的数据。 问题分析：首先计算狼的初始位置到兔子洞穴的直线距离： 由于狼奔跑的速度是兔子速度的两倍，兔子跑60码的时间狼可以跑120码。如果狼沿直线奔向兔窝，应该是可以追上兔子的。但是，有人推导出狼在追赶兔子过程中的运动曲线为 根据曲线方程，当时，。也就是说，在没有兔窝的情况下兔子一直往北跑，在跑到大约66码处将被狼追上。由此可知，在有兔窝时狼是追