c语言插值法编程,Lagrange插值的C语言编程.pdf

Lagrange插值的C语言编程

线性插值

给出两点(x , y ),i 0,1 ,则有线性插值函数

i i

x −x x −x

1 0

L (x ) l (x )y =+l (x )y y =+ y

1 0 0 1 1 0 1

x −x x −x

0 1 1 0

抛物插值

给出三点(x , y ),i 0,1,2 ,则有三点抛物线插值函数

i i

(x −x1)(x −x 2 ) (x −x0 )(x −x2 ) (x −x0 )(x −x1 )

L2 (x ) y0 =+ y1 + y2

(x0 −x1 )(x0 −x2 ) (x1 −x0 )(x1 −x2 ) (x2 −x0 )(x2 −x1 )

Lagrange 插值函数的一般形式

N

两点线性插值函数和三点抛物线插值函数属于 Lagrange 插值中比较常用的低次插值,对于 个点

(x , y ),i 0,1,2, …,N - 2,N -1 的Lagrange 插值函数的一般形式为:

i i

(x −x )(x −x )…(x −x )(x −x ) (x −x )(x −x )…(x −x )(x −x )

Ln−1(x ) 1 2 n−2 n−1 y 0 =+ 0 2 n−2 n−1 y1 +…

(x −x )(x −x )…(x −x )(x −x ) (x −x )(x −x )…(x −x )(x −x )

0 1 0 2 0 n−2 0 n−1 1 0 1 2 1 n−2 1 n−1

(x −x )(x −x )…(x −x )(x −x )

1 2 n−3 n−2

+ y n−1

(x −x )(x −x )…( − )( − )

x x x x

n−1 1 n−1 2 n−1

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
#include #include #include typedef strUCt data { float x; float y; }Data;//变量x和函数值y的结构 Data d[20];//最多二十组数据 float f(int s,int t)//牛顿插值法,用以返回插商 { if(t==s+1) return (d[t].y-d[s].y)/(d[t].x-d[s].x); else return (f(s+1,t)-f(s,t-1))/(d[t].x-d[s].x); } float Newton(float x,int count) { int n; while(1) { cout<>n; if(n<=count-1)// 插值次数不得大于count-1次 break; else system("cls"); } //初始化t,y,yt。 float t=1.0; float y=d[0].y; float yt=0.0; //计算y值 for(int j=1;j<=n;j++) { t=(x-d[j-1].x)*t; yt=f(0,j)*t; //cout<<f(0,j)<<endl; y=y+yt; } return y; } float lagrange(float x,int count) { float y=0.0; for(int k=0;k<count;k++)//这儿默认为count-1次插值 { float p=1.0;//初始化p for(int j=0;j<count;j++) {//计算p的值 if(k==j)continue;//判定是否为同一个数 p=p*(x-d[j].x)/(d[k].x-d[j].x); } y=y+p*d[k].y;//求和 } return y;//返回y的值 } void main() { float x,y; int count; while(1) { cout<>count; if(count<=20) break;//检查输入的是否合法 system("cls"); } //获得各组数据 for(int i=0;i<count;i++) { cout<<"请输入第"<<i+1<>d[i].x; cout<<"请输入第"<<i+1<>d[i].y; system("cls"); } cout<>x; while(1) { int choice=3; cout<<"请您选择使用哪种插值法计算:"<<endl; cout<<" (0):退出"<<endl; cout<<" (1):Lagrange"<<endl; cout<<" (2):Newton"<<endl; cout<>choice;//取得用户的选择项 if(choice==2) { cout<<"你选择了牛顿插值计算方法,其结果为:"; y=Newton(x,count);break;//调用相应的处理函数 } if(choice==1) { cout<<"你选择了拉格朗日插值计算方法,其结果为:"; y=lagrange(x,count);break;//调用相应的处理函数 } if(choice==0) break; system("cls"); cout<<"输入错误!!!!"<<endl; } cout<<x<<" , "<<y<<endl;//输出最终结果 }

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值