计算数值方法实验指导书.docVIP

PAGE 第 PAGE \* Arabic 8 页 《计算机数值方法B》 实验指导书 2012-02-20 目录 TOC \o 1-3 \h \z 目录 1 一、实验目的 2 二、实验内容 2 三、实验要求 4 四、示例程序 4 五、实验报告要求 8 六、参考资料 8 一、实验目的 《数值计算方法》是一门实践性和实用性都很强的课程，学习的目的在于应用已有的理论基础和技术方法，实现对实际计算问题的计算机求解。本实验配合《计算数值方法》课堂教学，通过上机实践实验环节，巩固课堂和书本上所学的知识、加强实践能力、提高解决实际计算问题的水平、启发创新思想。其基本目的是： (1) 培养理论联系实际的思路，训练综合运用已经学过的理论和实际知识去分析和解决实际问题的能力。 (2) 帮助学生全面消化已学的相关课程内容，深刻理解计算数值方法课程的内涵，培养使用电子计算机进行科学计算和解决问题的能力。 (3) 进行基本技能训练和巩固。使学生得到选择算法、编写程序、分析数值结果、写数值试验报告、课堂讨论等环节的综合训练。 二、实验内容 本实验指导书共包含6个基本实验，其中实验1——实验4为必做实验，要求学生在8个实验课时内完成，实验5和实验6为选作实验，供学有余力的同学完成。为使实验更为有成效，需要写出实验报告，以此可作为《计算机数值方法B》课程成绩评定的参考。 实验1 方程求根 熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度|x*-xn|0.5×10-5 实验2 线性方程组的直接解法 合理利用Gauss消元法、LU分解法或追赶法求解下列方程组： ① ② ③ ④（n=5,10,100,…） 实验3 线性方程组的迭代解法 使用雅可比迭代法或高斯-赛德尔迭代法对下列方程组进行求解。 实验4矩阵特征值与特征向量问题 使用幂法求A模为最大的特征值及其相应的特征向量。 实验5代数插值 使用拉格朗日插值法或牛顿插值法求解：已知f(x)在6个点的函数值如下表所示，运用插值方法，求f(0.596)的近似值。 x 0.40 0.55 0.65 0.80 0.90 1.05 f(x) 0.41075 0.57815 0.69675 0.88811 1.02652 1.25386 实验6 最小二乘法拟合多项式 给定数据点（xi ,yi），用最小二乘法拟合数据的多项式，并求平方误差。 xi 0 0.5 0.6 0.7 0.8 0.9 1.0 yi 1 1.75 1.96 2.19 2.44 2.71 3.00 三、实验要求 (1) 应用结构化程序设计编出通用程序，源程序要有详细的注释和说明； (2) 比较计算结果，分析数值解误差的原因； (3) 实验完成，要求提交实验结果并写出报告，分析计算结果是否符合问题的要求，找出计算成功的原因或计算失败的教训。 四、示例程序 （1）Gauss消元法 #include stdafx.h #include stdio.h double a[15][15],a0[15][15]; double b[15],b0[15],l[15]; int n; void displayA() { printf(\n); for(int j=1;j=n;j++) { for(int i=1;i=n;i++) printf(a[%d][%d]=%f ,j,i,a[j][i]); printf(b[%d]=%f\n,j,b[j]); } for(j=1;j=n;j++) printf(l[%d]=%f ,j,l[j]); printf(\n); } void main() { FILE *f; int i,j,k; f=fopen(Gaussdata.txt,r); fscanf(f,%d,n); for(i=1;i=n;i++) { for(j=1;j=n;j++) { fscanf(f,%lf,a[i][j]); a0[i][j]=a[i][j]; } fscanf(f,%lf,b[i]); b0[i]=b[i]; } displayA(); fclose(f); k=1; do { for(i=1;i=n;i++) { if(i==k) continue; l[i]=a0[i][k]/a0[k][k]; } for (j=k+1;j=n;j++) a[k][j]=a0[k][j]/a0[k][k]; b[k]=b0[k]/a0[k][k];