东南大学计算方法实验报告.doc

计算方法与实习 实验报告 学 院：电气工程学院 指导老师：李翠平 班 级：160093 姓 名：黄芃菲 学 号实习题一 实验1? 拉格朗日插值法 一、方法原理 n次拉格朗日插值多项式为：Ln(x)=y0l0(x)+y1l1(x)+y2l2(x)+…+ynln(x) n=1时，称为线性插值，L1(x)=y0(x-x1)/(x0-x1)+ y1(x-x0)/(x1-x0)=y0+(y1-x0)(x-x0)/(x1-x0) n=2时，称为二次插值或抛物线插值，精度相对高些 L2(x)=y0(x-x1)(x-x2)/(x0-x1)/(x0-x2)+y1(x-x0)(x-x2)/(x1-x0)/(x1-x2)+y2(x-x0)(x-x1)/(x2-x0)/(x2-x1) 二、主要思路 使用线性方程组求系数构造插值公式相对复杂，可改用构造方法来插值。 对节点xi(i=0,1,…,n)中任一点xk(0=k=n)作一n 次多项式lk(xk)，使它在该点上取值为1，而在其余点xi(i=0,1,…,k-1,k+1,…,n)上为0，则插值多项式为Ln(x)=y0l0(x)+y1l1(x)+y2l2(x)+…+ynln(x) 上式表明：n 个点xi(i=0,1,…,k-1,k+1,…,n)都是lk(x)的零点。可求得lk 三．计算方法及过程：1.输入节点的个数n ?????????????????????2.输入各个节点的横纵坐标 ??????????????????? ?3.输入插值点 ???????????????????? 4.调用函数，返回z 函数语句与形参说明 程序源代码如下： 形参与函数类型 参数意义 int n 节点的个数 double x[n]（double *x） 存放n个节点的值 double y[n]（double *y） 存放n个节点相对应的函数值 double p 指定插值点的值 double fun() 函数返回一个双精度实型函数值，即插值点p处的近似函数值 ? ? #includeiostream #includemath.h using namespace std; #define N 100??? double fun(double *x,double *y, int n,double p); void main() {int i,n; cout输入节点的个数n：; ?cinn; double x[N], y[N],p; coutplease input xiangliang x= endl; for(i=0;in;i++)cinx[i]; coutplease input xiangliang y= endl; for(i=0;in;i++)ciny[i]; coutplease input LagelangrichazhiJieDian p= endl; ?cinp; coutThe Answer= fun(x,y,n,p)endl; system(pause) ;} double fun(double x[],double y[], int n,double p) {double z=0,s=1.0; int k=0,i=0; double L[N]; while(kn) { if(k==0) { for(i=1;in;i++)s=s*(p-x[i])/(x[0]-x[i]); ?? ??L[0]=s*y[0]; ? ???k=k+1;} ?? else ??? {s=1.0; ? ??for(i=0;i=k-1;i++)s=s*((p-x[i])/(x[k]-x[i])); ? ??for(i=k+1;in;i++) s=s*((p-x[i])/(x[k]-x[i])); ? ??L[k]=s*y[k]; ?k++;} } for(i=0;in;i++)z=z+L[i]; return z; } 五．实验分析 n=2时，为一次插值，即线性插值 n=3时，为二次插值，即抛物线插值 n=1，此时只有一个节点，插值点的值就是该节点的函数值 n1时，结果都是返回0的；这里做了n=0和n=-7两种情况 3n100时，也都有相应的答案 常用的是线性插值和抛物线插值，显然，抛物线精度相对高些 n次插值多项式Ln（x）通常是次数为n的多项式，特殊情况可能次数小于n.例如：通过三点的二次插值多项式L2(x),如果三点共线，则y=L2(x)就是一条直线，而不是抛物线，这时L2(x)是一次式。 拟合曲线光顺性差 ? ? 实验2? 牛顿插值法 ? 一、方法原理及基本思路 在拉格朗日插值方法中，若增加一个节点数据，其插值的多项式需重新计算。现构造一个插值多项式Nn(x)，只需对