计算方法B上机题目.docVIP

2016年《计算方法B》上机题目 一.计算机语言要求 使用语言不做限制，可以使用C、C++、FORTRAN、VC、VB、C#、Matlab、PHP、JavaScript等语言。 二.实习报告内容 上机题目完成后，必须交一份上机报告。上机报告中应对每道题目包括： （1）上机题目内容; （2）详细说明实现的思想、算法依据、算法实现的结构; （3）详细完整的源程序，并附相关的注释说明; （4）给出必要的计算结果（若数据量大，则只需列出主要的数据内容），并对结果进行分析; （5）对上机中出现的问题进行分析总结； 三.实习报告要求 1.提供一份完整的上机报告的电子文档；然后再提供一份与电子文档对应的纸质上机报告。 2.电子文档中提供上机过程中的所有源程序、输出数据，以及可以运行的文件。如果程序的运行环境特殊，请附上运行程序所需要的软件环境。 3.上机报告严禁抄袭，如发现有抄袭现象，所有涉及抄袭的上机报告将被以作弊处理，并按零分处理,不再另行通知。 4.上机报告电子版一、二、三班发送到邮箱:xjtujsff@126.com,上机作业纸面作业请送到:理科楼338。 上机实习题目 1.某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所示： 分点 0 1 2 3 4 5 6 深度 9.01 8.96 7.96 7.97 8.02 9.05 10.13 分点 7 8 9 10 11 12 13 深度 11.18 12.26 13.28 13.32 12.61 11.29 10.22 分点 14 15 16 17 18 19 20 深度 9.15 7.90 7.95 8.86 9.81 10.80 10.93 (1)请用合适的曲线拟合所测数据点； (2)估算所需光缆长度的近似值，并作出铺设河底光缆的曲线图； 2.假定某天的气温变化记录如下表所示，试用数据拟合的方法找出这一天的气温变化的规律;试计算这一天的平均气温,并试估计误差。 时刻 0 1 2 3 4 5 6 7 8 9 10 11 12 平均气温 15 14 14 14 14 15 16 18 20 20 23 25 28 时刻 13 14 15 16 17 18 19 20 21 22 23 24 平均气温 31 34 31 29 27 25 24 22 20 18 17 16 3.线性方程组求解。 （1）编写程序实现大规模方程组的高斯消去法程序，并对所附的方程组进行求解。所附方程组的类型为对角占优的带状方程组。 （2）针对本专业中所碰到的实际问题，提炼一个使用方程组进行求解的例子，并对求解过程进行分析、求解。 附：数据文件格式说明 1.数据文件的文件名后缀为.dat,形式为：文件名.dat； 2.数据文件中的数据均为二进制记录结构，因此必须使用二进制方式进行读取； 3.数据文件的结构，分为以下四个部分： (1)文件标示部分,该部分用于存放数据文件的描述信息 结构如下(用C语言格式进行描述)： typedef struct FileInfo { long int id; // 数据文件标示 long int ver; // 数据文件版本号 long int id1; // 备用标志 } FILEINFO; 其中： ①id:为该数据文件的标识，值为0xF1E1D1A0即为：十六进制的F1E1D1A0 ②ver:为数据文件的版本号，值为16进制数据, 版本号 说明 0x101 系数矩阵为非压缩格式稀疏矩阵 0x102 系数矩阵为非压缩格式带状矩阵 0x201 系数矩阵为压缩格式稀疏矩阵 0x202 系数矩阵为压缩格式带状矩阵 ③id1:为备用标志字段，暂时未用; (2)矩阵描述部分:此部分中包括矩阵的阶数和上下带宽,如果是稀疏矩阵，则上下带宽值为0。 结构如下： typedef struct HeadInfo { long int n; // 方程组的阶数 long int q; // 带状矩阵上带宽 long int p; // 带状矩阵下带宽 } HEADINFO; (3) 系数矩阵数据部分：该部分存放方程组系数矩阵中的所有元素 ①如存贮格式为非压缩格式,则按行方式顺序存贮系数矩阵中的每一个元素