《计算方法》实验课程指导书.doc

《计算方法》实验课程指导书 徐中宇 朱洪秀 2011年3月 目 录 第一章 插值法 1 §1拉格朗日插值 1 §2牛顿插值 5 第二章 线性方程组的解法 8 §1 高斯消去法 8 §２　列主元消去法 12 §3　线性方程组的迭代解法 16 第三章 方程求根 21 §1 二分法 21 §2 牛顿迭代法 24 §3 埃特金（Aitken）迭代法 26 第一章 插值法 目的与要求： 掌握不同的输入、输出语句，注意节约内存方法。 熟悉拉格朗日插值和牛顿插值公式，并体会它们不同的特点。 §1拉格朗日插值 1．方法概要： 拉格朗日n次插值多项式 ，其中 可用双重循环来实现；内循环为连乘；外循环为连加。用数组表下标变量，给出。 2．程序流程图 图1-1 拉格朗日插值法程序流程图 3．程序及例 例1．已知函数表： xi 0.56160 0.56280 0.56401 0.56521 yi 0.82741 0.82659 0.82577 0.82495 用三次拉格朗日插值多项式求x=0.5635时的函数值。 程序为 #include stdafx.h #include iostream.h const int n=3; const double xx=0.5635; void main() { double x[n+1]={0.56160,0.56280,0.56401,0.56521}; double y[n+1]={0.82741,0.82659,0.82577,0.82495}; double L=0,Li=1; for (int i=0;i=n;i++) { Li=1; for (int j=0;j=n;j++) { if (j!=i) { Li=Li*(xx-x[j])/(x[i]-x[j]); } } L+=Li*y[i]; } coutThe result is:Lendl; } 例2． 已知函数表，用拉格朗日插值多项式求0.5，0.7，0.85三点处的函数值。 程序的结构：在程序的说明部分将已知数据分别存入一维数组x和y。在主函数main中提示人工键入插值点的数目m，继而在m次循环计算中分别提示键入插值点，计算插值，并输出结果。 程序中的主要常量，变量及函数：n为插值基点的最大下标。x[n+1],y[n+1]分别存放插值基点及其函数值。M为插值点个数，X存放插值点的值。 程序清单： #include stdafx.h #include stdio.h #define n 4 /*插值基点的最大下标*/ void main() { double x[n+1]={0.4,0.55,0.8,0.9,1}; /*插值基点的值*/ double y[n+1]={0.41075,0.57815,0.88811,1.02652,1.17520}; int m,k,i,j; float X,L,P; printf(\n Please enter m=); scanf(%d,m); /*键入插值点的个数*/ for(k=1;k=m;k++) { printf(\n Please enter X%d=,k); scanf(%f,X); P=0.0; for (i=0;i=n;i++) { L=1.0; for (j=0;j=n;j++) { if (j!=i) { L=L*(X-x[j])/(x[i]-x[j]); } } P=P+y[i]*L; } printf(P(%f)=%f\n,X,P); } } 计算结果：（当提示“Please enter m=”时，键入插值点个数3，再分别根据提示X1=或X2=，X3=依次键入插值点0.5,0.7,0.85），输出的结果分别为 P（0.500000）=0.521090 P(0.700000)=0.758589 P(0.850000)=0.