数值分析课程设计题目与要求(10级应数及创新班).docx

数值分析课程设计题目与要求 （10级应数及创新班） ［设计题一］ 编写顺序Gauss消去法和列主元Gauss消去法的函数，再分别调用这两个函数求解下面的84阶方 程组： 6 10 0 8 6 10 0 8 6 1 0 0 8 6 然后考虑将方程组的阶数取为10至100之间多个值进行求解。将你的计算结果与方程组的精确解进行 比较。从44快气“准二“省円三个方面分析以上两个算法，试提出改进的算法并加以实现和验证。 ［设计题二］ 编写平方根法和改进的平方根法（参见教材《计算方法》P54的例题2.5）的函数，然后分别调用 这两个函数求解对称正定方程组Ax二b,其中A和b分别为： （1）系数矩阵A为矩阵（阶数取为10至100之间多个值）： 10 1 0 0 1 10 1 0 ? ? ? 0 0 1 10 1 0 0 ■ ? 1 1 ? ? ? ? ■ ■ ? ■ ■ ? ? ? 10 ? 1 0 0 0 ? ? ? 1 10 1 0 0 1 10 1 - 0 0 1 10 向量b随机地选取； （2）系数矩阵A为Hilbert矩阵（阶数取为5至40之间多个值），即A的第i行第j列元素a,= 占，向量b的第i个分量取为S 将你的计算结果与方程组的精确解进 9 计 j—1 ” 计 j—1 行比较。 若出现问题，分析其原因，提出改进的设想并尝试实现之。 ［设计题三］ 对于迭代法 兀利=0?99心-球伙=0丄2,??????),它显然有不动点“ =0。试设计2个数值实验 得到收敛阶数的大概数值(不利用判定收敛阶的判据定理): 直接用收敛阶的定义； 用最小二乘拟合的方法。 ［设计题四］ 湖水在夏天会出现分层现象，接近湖面温度较高，越往下温度变低。这种上热下冷的现象影响了水 的对流和混合过程，使得下层水域缺氧.导致水生鱼类的死亡。如果把水温T看成深度x的函数T(x),有 某个湖的观测数据如下： T(°C) 22.8 22.8 22.8 20.6 13.9 11 7 11.1 11.1 x(m) 0 2.3 4.9 9.1 13.7 18.3 22.9 27.2 环境工程师希望： 用三次样条插值求出T(x)。 rrp 1 2 rp 求在什么深度处》的绝对值达到最大(即於 ［设计题五］ 某飞机头部的光滑外形曲线的型值点坐标由下表给出: ? I 0 1 2 3 4 5 6 7 8 9 10 xi 0 70 130 210 337 578 776 1012 1142 1462 1841 yi 0 57 78 103 135 182 214 244 256 272 275 试建立其合适的模拟曲线(未必是用拟合方法)，并求在点X = 1OO, 250, 400, 500, 800处的函数 ? ? ? 值y及一阶、二阶导数值y\ y\绘出模拟曲线的图形。 ［设计题六］ 给定初值问题 给定初值问题 / = -y + x2e\ 1 x 2 x W) = o 其精确解为卩=/（「一 °）,分别按下列方案求它在节点% = 1+O k （比二12，…,10）处的数值解及 误差。比较各方法的优缺点，并将计算结果与精确解做比较（列表、画图）。 （方案I）欧拉法，步长/?二0.025,力二0.1 ； （方案II）改进的欧拉法，步长/?二0.05,/?二0.1 ； （方案川）四阶经典龙格一库塔法，步长/7=0.1o ［设计题七］ 生态环境学家在研究自然界中两个生物种群数目变化时得到一组常微分方程。 假设有两种生物（例如一种是蓝鲸，另一种是南极磷虾），前者在时刻r时的数量为“（/）,后者在 时刻/时的数目为x2（r）,并假设它们都是f的连续可微函数。蓝鲸是以磷虾为主要食物的。当没有食物 来源时蓝鲸数目会减少，其减少速度与当时蓝鲸的数目成线性关系，即 牛= . （1） at 当有食物来源时，蓝鲸的数目会增加。增加的速度和它捕食的数目有关，即 ⑵⑶⑷⑸dx{ ⑵ ⑶ ⑷ ⑸ ■- d“（/） x2［t）. at 合并（1）和（2）,得到蓝鲸变化速度满足的微分方程 dx{ ,、 = -CX,（0 + d%,（r）X2（r）. 同样，在没有蓝鲸时，磷虾的增加速度满足 dXy 訂吒⑴. 考虑到被捕食情况，则磷虾的数目满足 JXy 万二“2⑴一 b坷⑴心⑴. 合并（3）和（5）,得到著名的Lotka-Volterra方程 —=-ex + dx.x^di - = ax2 -bxlx1 .dt 其中u,b,c,d均为正常数。 （6）是一个非线性常微分方程组，不可能有解析解。 假设 a = 1.2,/? = 0.6,c = 0.〃 = 0.3,而且初始值为.勺(0)二2, x2 (0)=1. 1） 分别用欧拉法、改进的欧拉法和四阶经典龙格一库塔法，取多种步长求解（6）o把“（Q和心 （/）画在