数值分析实验解析.docVIP

实验课程：数值分析 专 业： 班 级： 学 号： 姓 名： 中北大学理学院 实验一 函数插值方法 【实验类型】验证性 【实验课时】2 【】的的节点值。试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。 数据如下： 0.4 0.55 0.65 0.80 0.95 1.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次Lagrange多项式或分段三次插值多项式，并计算，的值。（提示：结果为，） 【实验方法及步骤】 利用Lagrange插值公式编写出插值多项式程序； 【实验结果】 #define N 5 float x[]={0.4,0.55,0.65,0.80,0.95,1.05}; double y[]={0.41075,0.57815,0.69675,0.90,1.00,1.25382}; float L(float xx) { int i,k; float f=0,m1,m2; for(k=0;k=N;k++) {m1=1;m2=1; for(i=0;i=N;i++) if(i!=k) {m1*=xx-x[i]; m2*=x[k]-x[i];} f+=y[k]*m1/m2; } return f; } main() { printf(f(0.596)=%lf\n,L(0.596)); printf(f(0.99)=%lf\n,L(0.99)); } 实验二 正交多项式与曲线拟合 【实验类型】检验性 【实验学时】2 【实验内容】 编写出legendre、Chebyshev多项式的程序； 从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。 在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量与时间的拟合曲线。 0 5 10 15 20 25 30 35 40 45 50 55 0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.02 4.64 1 1 1 1 1 1 1 1 1 1 1 1 【实验前的预备知识】 1、掌握曲线拟合的最小二乘法； 2、最小二乘法亦可用于解超定线代数方程组； 3、探索拟合曲线的选择与拟合精度间的关系。 【实验方法与步骤】 1、用最小二乘法进行曲线拟合； 2、近似解析表达式为： 选取基函数，，则得到 ，，， ，，于是得方程组 另外选取一个近似表达式，尝试拟合效果的比较。 【实验程序与结果】 legendre多项式的程序: #includestdio.h double L(int n,double x) { if(n==0) return 1; if(n==1) return x; else return ((2*n-1)*x*L(n-1,x)-(n-1)*L(n-2,x))/n; } int main() { int n; double x; double y; printf(Input n,x:\n); scanf(%d%lf,n,x); y=L(n,x); printf(%lf\n,y); //输出第n个勒让德正交多项式 return 0; //输入错误返回0 } 利用数据的最小二乘法求得拟合曲线的程序： #includestdio.h #define N 25 main() { int i,n; double x[N],y[N]; double a11,a12=0,a21=0,a22=0,d1=0,d2=0; double a,b; printf(\n input n value:); do {scanf(%d,n); if(nN) printf(\n please re-input n value!); }while(nN||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[i]; d1+=y[i]; d2+=x[i]*y[i];} a12