数值分析实验指导书

1、数值分析实验指导书编者：陈晓红 殷明实验一11 水手、猴子和椰子问题：五个水手带了一只猴子来到南太平洋的一个荒岛上，发现那里有一大堆椰子。由于旅途的颠簸，大家都很疲惫，很快就入睡了。第一个水手醒来后，把椰子平分成五堆，将多余的一只给了猴子，他私藏了一堆后便又去睡了。第二、第三、第四、第五个水手也陆续起来，和第一个水手一样，把椰子分成五堆，恰多一只猴子，私藏一堆，再去入睡，天亮以后，大家把余下的椰子重新等分成五堆，每人分一堆，正好余一只再给猴子，试问原先共有几只椰子？试分析椰子数目的变化规律，利用逆向递推的方法求解这一问题（15621）。12 设，（1）从尽可能精确的近似值出发，利用递推公式：计算机从到的近似值；（2）从较粗糙的估计值出发，用递推公式：计算从到的近似值；（3）分析所得结果的可靠性以及出现这种现象的原因。13 绘制Koch分形曲线问题描述：从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的另两条边代替，形成具有5个结点的新的图形（图1-4）；在新的图形中，又将图中每一直线段中间的三分之一部分都用一个等边三角形的另两条边代替，再次形成新的图形（图1-5），这时，图形

2、中共有17个结点。这种迭代继续进行下去可以形成Koch分形曲线。在迭代过程中，图形中的结点将越来越多，而曲线最终显示细节的多少取决于所进行的迭代次数和显示系统的分辨率。Koch分形曲线的绘制与算法设计和计算机实现相关。 图14 第一次迭代 图15 第二次迭代问题分析：考虑由直线段（2个点）产生第一个图形（5个点）的过程，设和分别为原始直线段的两个端点。现在需要在直线段的中间依次插入三个点产生第一次迭代的图形（图1-4）。显然，位于点右端直线段的三分之一处，点绕旋转60度（逆时针方向）而得到的，故可以处理为向量经正交变换而得到向量，形成算法如下：（1）；（2）；（3）；在算法的第三步中，A为正交矩阵。；这一算法将根据初始数据（和点的坐标），产生图1-4中5个结点的坐标。这5个结点的坐标数组，组成一个52矩阵。这一矩阵的第一行为为的坐标，第二行为的坐标，第二行为的坐标第五行为的坐标。矩阵的第一列元素分别为5个结点的x坐标 ，第二列元素分别为5个结点的y坐标。问题思考与实验：（1）考虑在Koch分形曲线的形成过程中结点数目的变化规律。设第k次迭代产生结点数为，第迭代产生结点数为，试写出和之间

3、的递推关系式；（2）参考问题分析中的算法，考虑图1-4到图1-5的过程，即由第一次迭代的5个结点的结点坐标数组，产生第二次迭代的17个结点的结点坐标数组的算法；（3）考虑由第k次迭代的个结点的结点坐标数组，产生第次迭代的个结点的结点坐标数组的算法；（4）设计算法用计算机绘制出如下的Koch分形曲线（图1-6）。实验二21 用高斯消元法的消元过程作矩阵分解。设消元过程可将矩阵A化为上三角矩阵U，试求出消元过程所用的乘数、并以如下格式构造下三角矩阵L和上三角矩阵U验证：矩阵A可以分解为L和U的乘积，即A=LU。22 用矩阵分解方法求上题中A的逆矩阵。记分别求解方程组由于三个方程组系数矩阵相同，可以将分解后的矩阵重复使用。对第一个方程组，由于A=LU，所以先求解下三角方程组,再求解上三角方程组，则可得逆矩阵的第一列列向量；类似可解第二、第三方程组，得逆矩阵的第二列列向量的第三列列向量。由三个列向量拼装可得逆矩阵。2.3 验证希尔伯特矩阵的病态性：对于三阶矩阵取右端向量，验证：（1）向量是方程组的准确解；（2）取右端向量b的三位有效数字得，求方程组的准确解，并与X的数据作比较 。说明矩阵的病态

4、性。实验三31 用泰勒级数的有限项逼近正弦函数用计算机绘出上面四个函数的图形。32 绘制飞机的降落曲线一架飞机飞临北京国际机场上空时，其水平速度为540km/h，飞行高度为1 000m。飞机从距机场指挥塔的横向距离12 000m处开始降落。根据经验，一架水平飞行的飞机其降落曲线是一条三次曲线。建立直角坐标系，设飞机着陆点为原点O，降落的飞机为动点，则表示飞机距指挥塔的距离，表示飞机的飞行高度，降落曲线为该函数满足条件：（1）试利用满足的条件确定三次多项式中的四个系数；（2）用所求出的三次多项式函数绘制出飞机降落曲线。33 追赶曲线的计算机模拟问题描述：欧洲文艺复兴时期的著名人物达芬奇曾经提出一个有趣的“狼追兔子”问题，当一只兔子正在它的洞穴南面60码处觅食时，一只饿狼出现在兔子正东的100码处。兔子急忙奔向自己的洞穴，狼立即以快于兔子一倍的速度紧追兔子不放。兔子一旦回到洞穴便逃脱厄，问狼是否会追赶上兔子？这一问题的研究方法可以推广到如鱼雷追击潜艇、地对空导弹击飞机等问题上去。在对真实系统做实验时，可能时间太长、费用太高、危险太大、甚至很难进行。计算机模拟是用计算机模仿实物系统，对系统的

5、结构和行为进行动态演示，以评价或预测系统的行为效果。根据模拟对象的不同特点，分为确定性模拟和随机性模拟两大类。模拟通常所用的是时间步长法，即按照时间流逝的顺序一步一步对所研究的系统进行动态演示，以提取所需要的数据。问题分析：首先计算狼的初始位置到兔子洞穴的直线距离：由于狼奔跑的速度是兔子速度的两倍，兔子跑60码的时间狼可以跑120码。如果狼沿直线奔向兔窝，应该是可以追上兔子的。但是，有人推导出狼在追赶兔子过程中的运动曲线为根据曲线方程，当时，。也就是说，在没有兔窝的情况下兔子一直往北跑，在跑到大约66码处将被狼追上。由此可知，在有兔窝时狼是追赶不上兔子的。用计算机模拟的方法也可以得到同样的结论。取时间步长为1s，随时间步长的增加，考虑这一系统中的各个元素（狼和兔子）所处的位置变化规律，用计算机作出模拟。最后，根据第60s时狼所在的位置的坐标，判断狼是否能追上兔子。问题思考与实验：（1）设兔子奔跑的速度为，则狼运动的速度为。建立平面直角坐标系，若当时刻，兔子位于点处，狼位于点处。试根据，的坐标确定一个单位向量描述狼在时段内的运动方向。（2）根据狼的运动方向和速度推导到的坐标的具体表达式；

6、（3）用计算机绘制追赶曲线的图形（包括静态和动态的图形）。实验四41 曾任英特尔公司董事长的摩尔先生早在1965年时，就观察到一件很有趣的现象：集成电路上可容纳的零件数量，每隔一年半左右就会增长一倍，性能也提升一倍。因而发表论文，提出了大大有名的摩尔定律（Moores Law），并预测未来这种增长仍会延续下去。下面数据中，第二行数据为晶片上晶体数目在不同年代与1959年时数目比较的倍数。这些数据是推出摩尔定律的依据：年代19591962196319641965增加倍数13456试从表中数据出发，推导线性拟合的函数表达式。42 参考算法4.2设计绘制Bezier曲线的程序，选取四个点的坐标数据作为控制点绘制飞机机翼剖面图草图的下半部分图形；结合例4.4中上半部分图形绘出完整的机翼草图。最后写出机翼剖面图曲线上20个点处的坐标数据。43 神经元模型用于蠓的分类识别（MCM1989A题）问题描述：生物学字试图对两类蠓虫（Af与Apf）进行鉴别，依据的资料是蠓虫的触角和翅膀的长度，已经测得9只Af和6只Apf的数据（触角长度用x表示，翅膀长度用y表示）Af数据X1241361381381381

7、40148154156Y127174164182190170182182208Apf数据x1.141.181.201.261.281.30Y1.781.961.862.002.001.96现需要解决三个问题：（1）如何凭借原始资料（15对数据，被称之为学习样本）制定一种方法，正确区分两类蠓虫；（2）依据确立的方法，对题目提供的三个样本：（1.24,1.80）,(1.28,1.84),(1.40,2.04)加以识别；（3）设Af是宝贵的传粉益虫，Apf是某种疾病的载体，是否应该修改分类方法。问题分析：首先画出15对数据的散点图，其中，Af用*标记，Apf用标记。观察图形，可以发现，Af的点集中在图中右下角，而Apf的点集中在图中左上角。应该存在一条直线L位于两类点之间，作为Af和Apf分界线，这条直线L的确定应依据问题所给的数据，即学习样本。设这条直线的方程为对于平面上任意一点P（x,y），如果该点在直线上，将其坐标代入直线方程则使方程成为恒等式，即使方程左端恒为零；如果点不在直线上，将其坐标代入直线方程，则方程左端不为零。由于f和pf的散点都不在所求的直线上，故将问题所提供的数据代入直

8、线方程左端应该得到表达式的值大于零或者小于零两种不同的结果。这需要建立一个判别系统，引入判别函数,当属于f类时，否则。为了对判别系统引入学习机制，在学习过程中将两种不同的状态，以“”和“”表示。当属于f类时，否则。取于是由所给数据形成约束条件，这是关于判别函数中的三个待定系数的线性方程组：这是包括三个未知数共15个方程的超定方程组，可以求方程组的最小二乘解。（1）根据上面分析写出对应的正规方程组并求解。（2）确定分类边界直线的方程。由所给数据用判别函数判别三个新蠓虫的类属，即当时，判为Af类：当时，判为Apf类。实验五51 用几种不同的方法求积分的值。（1）牛顿-莱布尼茨公式；（2）梯形公式；（3）辛卜生公式；（4）复合梯形公式。52 设X为标准正态随机变量，即XN（0，1）。现分别取，试设计算法计算30个不同的概率值；，并将计算结果与概率论教科书中的标准正态分布函数表作比较。（提示： 53 设某城市男子的身高XN（170，36）（单位：cm），应如何选择公共汽车门的高度H使男子与车门碰头的机会小于1%。问题分析：由题设男子身高XN（170，36）（单位：cm），应如何选择公共汽车门的

9、高度H使男子与车门碰头的机会小于1%。问题分析：由题设男子身高数据服从平均值为170（cm），方差为6（cm）的正态分布，其分布密度函数为按正态分布的分布规律（原则），这个城市的男子身高超过188（cm）的人数极少。故可以对H=188，187，186，求出概率的值，观察使概率不超过1%的H，以确定公共汽车门应该取的高度。概念值的计算实际上是求定积分（1）选用一种数值求积公式或用数学软件分别计算出H=180、181、188时定积分近似值。（2）根据上面计算的积分值，按题目要求确定公共汽车门的高度取值（答案184cm）。如果将汽车门的高度取180cm，是否满足大多数市民的利益？（3）用计算机模拟的方法来检验你的结论，计算机产生10 000个正态随机数（它们服从均值为170，方差为6的正态分布）来模拟这个城市中10 000个男子的身高，然后统计出这10 000人中身高超过180（cm）的男子数量所占的百分比。54 用蒙特卡罗方法计算二重积分，其中，55 求空间曲线L：弧长公式为实验六61 用欧拉公式和四阶龙格-库塔法分别求解下列初值问题；62 用求解常微分方程初值问题的方法计算积分上限函数的值，实际上是将上面表达式两端求导化为常微分方程形式，并用初值条件。

《数值分析实验指导书》由会员新**分享，可在线阅读，更多相关《数值分析实验指导书》请在金锄头文库上搜索。