拉格朗日差值法求解.doc

#includestdio.h #includemath.h #define n 3 float x[3],y[3],c[3]={1,1,1}; //////////////////////////////////////////////////////////////////// void main() { float input(),output(),l(),f(float); float a,b; input(); output(); l(); printf(请输入变量的值 a:\n); scanf(%f,a); b=f(a); printf(输入所求的值 y:\n); printf(%f ,b); printf(\n); } //////////////////////////////////////////////////////////////////// float f(float a) { int i,j; float t,y=0.0; for(i=0;in;i++){ t=1.0; for(j=0;jn;j++){ if(x[i]-x[j]){ t=t*(a-x[j]); } } t=t*c[i]; y=y+t; } return y; } //////////////////////////////////////////////////////////////////// float l() { int i,j; for(i=0;in;i++){ for(j=0;jn;j++){ if(x[i]-x[j]){ c[i]=c[i]/(x[i]-x[j]);} } c[i]=c[i]*y[i]; } return 1; } //////////////////////////////////////////////////////////////////// float input() { int i; printf(输入x(n):\n); for(i=0;in;i++){ scanf(%f,x[i]); } printf(输入 y(n):\n); for(i=0;in;i++){ scanf(%f,y[i]); } return 1; } //////////////////////////////////////////////////////////////////// float output() { int i; printf(输出 x(n):\n); for(i=0;in;i++){ printf(%f ,x[i]); } printf(\n); printf(输出 y(n):\n); for(i=0;in;i++){ printf(%f ,y[i]); } printf(\n); return 1; } //////////////////////////////////////////////////////////////////// 作业5（1） 课本p70 用拉格朗日插值法求值 主函数：void main() 子函数：input(),output(),l(), f(). 输入所要求的变量的参数值；主函数的主要作用是调用各个子函数，计算各个值。 子函数f(float a);的主要作用是传递a的值然后计算初相应的结果。 利用for循环来计算函数的叠加和每个分量。 拉格朗日迭代法德每个步骤含有x的量。 计算拉格朗日迭代法的每个系数（权）的大小 利用全局变量待会计算值 子函数input(); 本函数的作用是输入各个变量的值， 输入x（n） n=3 输入y(n) n=3 利用子函数的形式可以减少书写，方便 子函数output() 输出 x(n) n=3 输出y(n) n=3 可以随时随地的输出，只需要写一个函数名即可。