数值计算(或计算方法)试验教学讲义.docVIP

第一章 实验的目的和要求 1.1 实验目的 为了掌握计算方法的基本思想、原理和方法，要注意计算方法的处理技巧与计算机实现的结合，需将各种数值方法设计成算法，并编制好程序，拿到计算机上实现，最后得到可行性的验证。 1.2 实验要求 ⑴ 用C或C++、Java、FORTRAN、Matlab等计算机程序设计语言编写程序。 ⑵ 上机前充分准备，复习相关知识，选用合适的数据结构并详细设计算法，尽量写出具有通用性的程序，反复检查程序。 ⑶ 上机时快速输入程序；首先排除语法错误；然后采用多组数据，详细测试，排除逻辑错误；最后将程序调试成功，运行程序得到准确结果。 ⑷ 完成计算后，反复体会和分析，试着改善计算复杂性，使程序或算法更加完美。 1.3 实验环境 1.3.1 硬件环境 CPU : Pentium 4以上 内存：256MB以上 1.3.2 软件环境 (1)操作系统：Microsoft Windows XP 和 2000 (2)编译器 ：C或C++、Java、FORTRAN、Matlab 1.4 本实验课程与其它课程的关系 本课程的前导课程有高等数学、线性代数（或高等代数）、C语言或FORTRAN语言等，最好事先开设数据结构；后续课程有计算机图形学、图像处理、模式识别等。 第二章 实验的计划和内容 2.1 实验计划 计算方法实验课共安排20学时。计算方法实验计划如下 ⑴Lagrange插值多项式 ⑵Newton插值多项式 ⑶Hermite插值多项式 ⑷最小二乘法 ⑸复化求积公式 ⑹Romberg求积公式 ⑺数值微分的外推算法 ⑻Gauss消元法 ⑼直接三角分解法 ⑽解方程组的迭代法 2.2 实验内容 共十个实验题目，每次课(2学时)一个题目。 2.2.1实验一 实验题目：Lagrange插值多项式 相关知识：通过n+1个节点的次数不超过n的Lagrange插值多项式为： 其中，Lagrange插值基函数，k=0,1,…,n。 另外，补充C语言绘制图形方面的内容如下 屏幕坐标系 坐标原点在屏幕的左上角，x轴水平向右，y轴垂直向下。 常用的绘图函数(绘图库函数所在的头文件 graphics.h) 初始化图形系统的函数 void initgraph(int *graphdriver,int *graphmode, char *pathtodriver); 画点函数 void putpixel(int x,int y,int pixelcolor); 移“画笔”函数 void moveto(int x,int y); 画直线函数 void line(int x1,int y1,int x2,inty2); void lineto(int x,int y); 设置前景颜色函数 void setcolor(int color); 设置背景颜色函数 void setbkcolor(int color); 设置画线宽度和类型函数 void setlinestyle(int linestyle, unsigned upattern,int thickness); 关闭图形系统函数 void closegraph(void); 绘图程序的设计模式 #include graphics.h main() {int graphdriver=DETECT,graphmode; initgraph(graphdriver,graphmode, ); 调用绘图函数进行绘图 closegraph(); } 数据结构：两个一维数组或一个二维数组 算法设计：（略） 编写代码：（略） 实验用例： 已知函数y=f(x)的一张表： x 0 10 20 30 40 50 60 70 80 90 100 110 120 y 5 1 7.5 3 4.5 8.8 15.5 6.5 -5 -10 -2 4.5 7 试验要求：利用Lagrange插值多项式求被插值函数f(x)在点x=65处的近似值。建议：画出Lagrange插值多项式的曲线。 2.2.2 实验二 实验题目：Newton插值多项式 相关知识：通过n+1个节点的次数不超过n的Newton插值多项式为： 数据结构：两个一维数组或一个二维数组 算法设计：（略） 编写代码：（略） 实验用例： 已知函数y=f(x)的一张表(同上一个试验) 试验要求：利用Newton插值多项式求被插值函数f(x)在点x=65处的近似值。建议：画出Newton插值多项式的曲线。 2