中北大学数值14实验汇总.docxVIP

实验类别：数值分析专 业：信息与计算科学班 级 号：1308024120姓 名：杨燕中北大学理学院实验一 函数插值方法【实验内容】给定一元函数的个节点值，数据如下：0.40.550.650.800.951.050.410750.578150.696750.901.001.25382求五次Lagrange多项式或分段三次插值多项式或Newton插值多项式，并计,的值。（提示：结果为≈0.625732,≈1.05423）【实验方法与步骤】利用Lagrange插值公式 ，用C语言编写出插值多项式程序如下：#include stdio.h#define N 5 float x[]={0.4,0.55,0.65,0.80,0.95,1.05}; float y[]={0.41075,0.57815,0.69675,0.90,1.00,1.25382}; float p(float xx) { int i,k; float pp=0,m1,m2; for(i=0;i=N;i++) { m1=1;m2=1; for(k=0;k=N;k++) if(k!=i) { m1*=xx-x[k]; m2*=x[i]-x[k]; } pp+=y[i]*m1/m2; } return pp; } main() { printf(f(0.596)=%lf\n,p(0.596)); printf(f(0.99)=%lf\n,p(0.99)); } 【实验结果】【思考】1、给出的程序求行不行，精度高不高？2、五次Lagrange多项式与Newton插值多项式是同一个多项式吗？五次Lagrange多项式与Newton插值多项式是同一个多项式。3、为什么高次插值不能令人满意？一般来说，节点个数越多，插值函数和被插值函数就有越多的地方相等。但是随着插值节点个数的增加，两个插值节点之间插值函数并不一定能够很好地逼近被插值函数。再次，从舍入误差看，高次插值由于计算量大，可能会产生更严重的误差积累，所以，稳定性得不到保证。此时就会出现龙格现象。实验类别：数值分析专 业：信息与计算科学班 级 号：1308024120姓 名：杨燕中北大学理学院实验二 函数逼近与曲线拟合【实验内容】1、编写出Legendre、Chebyshev多项式的程序；2、从随机的数据中找出其规律性，给出其近似表达式，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。例如在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量与时间的拟合曲线。(分)051015202530354045505501.72.162.863.443.874.154.374.514.584.024.64111111111111【实验方法与步骤】1、用C编写语言出Legendre多项式的程序如下：#include stdio.hdouble p(int n, double x){ if(n==0) return 1; else if(n==1) return x; else return ((2*n-1)*x*p(n-1,x)-(n-1)*p(n-2,x))/n;}int main(){ int n; double x; double y; printf(input n, x:\n); scanf(%d%lf,n,x); y=p(n,x); printf(%lf\n,y); return 0;}2、给出表格数据的近似解析表达式为；，选取基函数，，则得到，，，，，于是得方程组，用C语言编写曲线拟合的程序如下：#includestdio.h#define Max_N 25main(){int i,n;double x[Max_N],y[Max_N];double a11,a12,a21,a22,d1,d2;double a,b;printf(\n Input n value:);do{scanf(%d,n);if(nMax_N)printf(\n please re_input n value:); }while(nMax_N||n=0);printf(input x[i],i=0,...%d:\n,n-1);for(i=0;in;i++)scanf(%lf,x[i]);printf(input y[i],i=0,...%d:\n,n-1);for(i=0;in;i++)scanf(%lf,y[i]);for(i=0;in;i++){a21+=x[i];a22+=x[i]*x