本文目录一览:
- 1、用C语言实现拉格朗日插值、牛顿插值、等距结点插值算法
- 2、用C语言编写一个线性插值程序
- 3、c语言中插值排序法怎么打程序
- 4、用拉格朗日插值法(二次插值)来计算某个f(x)的值,我用了C语言编,但是编译的时候很多错误,求修改啊!
- 5、两元lagrange插值如何用C语言表示
- 6、求双线性插值法的C语言程序!帮帮忙!拜托各位了!
用C语言实现拉格朗日插值、牛顿插值、等距结点插值算法
#includestdio.h
#includestdlib.h
#includeiostream.h
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值(即n次插值):";//获得插值次数
cinn;
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;
//coutf(0,j)endl;
y=y+yt;
}
return y;
}
float lagrange(float x,int count)
{
float y=0.0;
for(int k=0;kcount;k++)//这儿默认为count-1次插值
{
float p=1.0;//初始化p
for(int j=0;jcount;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"请输入x[i],y[i]的组数,不得超过20组:";//要求用户输入数据组数
cincount;
if(count=20)
break;//检查输入的是否合法
system("cls");
}
//获得各组数据
for(int i=0;icount;i++)
{
cout"请输入第"i+1"组x的值:";
cind[i].x;
cout"请输入第"i+1"组y的值:";
cind[i].y;
system("cls");
}
cout"请输入x的值:";//获得变量x的值
cinx;
while(1)
{
int choice=3;
cout"请您选择使用哪种插值法计算:"endl;
cout" (0):退出"endl;
cout" (1):Lagrange"endl;
cout" (2):Newton"endl;
cout"输入你的选择:";
cinchoice;//取得用户的选择项
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;
}
coutx" , "yendl;//输出最终结果
}
用C语言编写一个线性插值程序
#include stdio.h
double Lerp(double x0,double y0,double x1,double y1,double x)
{
double dy = y1 - y0;
if(dy == 0){
printf("除0错误!\n");
return 0;
}
return x * (x1 - x0) / dy;
}
int main()
{
double x0,x1,y1,y0,x,y;
printf("Inptu x0 y0 x1 y1 x:");
scanf("%lf %lf %lf %lf %lf",x0,y0,x1,y1,x);
y = Lerp(x0,y0,x1,y1,x);
printf("y = %lf\n",y);
return 0;
}
c语言中插值排序法怎么打程序
(1)“冒泡法”
冒泡法大家都较熟悉。其原理为从a[0]开始,依次将其和后面的元素比较,若a[0]a[i],则交换它们,一直比较到a[n]。同理对a[1],a[2],...a[n-1]处理,即完成排序。下面列出其代码:
void bubble(int *a,int n) /*定义两个参数:数组首地址与数组大小*/
{
int i,j,temp;
for(i=0;in-1;i++)
for(j=i+1;jn;j++) /*注意循环的上下限*/
if(a[i]a[j]) {
temp=a[i];
a[i]=a[j];
a[j]=temp;
}
}
冒泡法原理简单,但其缺点是交换次数多,效率低。
下面介绍一种源自冒泡法但更有效率的方法“选择法”。
(2)“选择法”
选择法循环过程与冒泡法一致,它还定义了记号k=i,然后依次把a[k]同后面元素比较,若a[k]a[j],则使k=j.最后看看k=i是否还成立,不成立则交换a[k],a[i],这样就比冒泡法省下许多无用的交换,提高了效率。
void choise(int *a,int n)
{
int i,j,k,temp;
for(i=0;in-1;i++) {
k=i; /*给记号赋值*/
for(j=i+1;jn;j++)
if(a[k]a[j]) k=j; /*是k总是指向最小元素*/
if(i!=k) { /*当k!=i是才交换,否则a[i]即为最小*/
temp=a[i];
a[i]=a[k];
a[k]=temp;
}
}
}
选择法比冒泡法效率更高,但说到高效率,非“快速法”莫属,现在就让我们来了解它。
(3)“快速法”
快速法定义了三个参数,(数组首地址*a,要排序数组起始元素下标i,要排序数组结束元素下标j). 它首先选一个数组元素(一般为a[(i+j)/2],即中间元素)作为参照,把比它小的元素放到它的左边,比它大的放在右边。然后运用递归,在将它左,右两个子数组排序,最后完成整个数组的排序。下面分析其代码:
void quick(int *a,int i,int j)
{
int m,n,temp;
int k;
m=i;
n=j;
k=a[(i+j)/2]; /*选取的参照*/
do {
while(a[m]kmj) m++; /* 从左到右找比k大的元素*/
while(a[n]kni) n--; /* 从右到左找比k小的元素*/
if(m=n) { /*若找到且满足条件,则交换*/
temp=a[m];
a[m]=a[n];
a[n]=temp;
m++;
n--;
}
}while(m=n);
if(mj) quick(a,m,j); /*运用递归*/
if(ni) quick(a,i,n);
}
(4)“插入法”
插入法是一种比较直观的排序方法。它首先把数组头两个元素排好序,再依次把后面的元素插入适当的位置。把数组元素插完也就完成了排序。
void insert(int *a,int n)
{
int i,j,temp;
for(i=1;in;i++) {
temp=a[i]; /*temp为要插入的元素*/
j=i-1;
while(j=0tempa[j]) { /*从a[i-1]开始找比a[i]小的数,同时把数组元素向后移*/
a[j+1]=a[j];
j--;
}
a[j+1]=temp; /*插入*/
}
}
(5)“shell法”
shell法是一个叫 shell 的美国人与1969年发明的。它首先把相距k(k=1)的那几个元素排好序,再缩小k值(一般取其一半),再排序,直到k=1时完成排序。下面让我们来分析其代码:
void shell(int *a,int n)
{
int i,j,k,x;
k=n/2; /*间距值*/
while(k=1) {
for(i=k;in;i++) {
x=a[i];
j=i-k;
while(j=0xa[j]) {
a[j+k]=a[j];
j-=k;
}
a[j+k]=x;
}
k/=2; /*缩小间距值*/
}
}
上面我们已经对几种排序法作了介绍,现在让我们写个主函数检验一下。
#includestdio.h
/*别偷懒,下面的"..."代表函数体,自己加上去哦!*/
void bubble(int *a,int n)
{
...
}
void choise(int *a,int n)
{
...
}
void quick(int *a,int i,int j)
{
...
}
void insert(int *a,int n)
{
...
}
void shell(int *a,int n)
{
...
}
/*为了打印方便,我们写一个print吧。*/[code]
void print(int *a,int n)
{
int i;
for(i=0;in;i++)
printf("%5d",a[i]);
printf("\n");
}
main()
{ /*为了公平,我们给每个函数定义一个相同数组*/
int a1[]={13,0,5,8,1,7,21,50,9,2};
int a2[]={13,0,5,8,1,7,21,50,9,2};
int a3[]={13,0,5,8,1,7,21,50,9,2};
int a4[]={13,0,5,8,1,7,21,50,9,2};
int a5[]={13,0,5,8,1,7,21,50,9,2};
printf("the original list:");
print(a1,10);
printf("according to bubble:");
bubble(a1,10);
print(a1,10);
printf("according to choise:");
choise(a2,10);
print(a2,10);
printf("according to quick:");
quick(a3,0,9);
print(a3,10);
printf("according to insert:");
insert(a4,10);
print(a4,10);
printf("according to shell:");
shell(a5,10);
print(a5,10);
}
用拉格朗日插值法(二次插值)来计算某个f(x)的值,我用了C语言编,但是编译的时候很多错误,求修改啊!
#include stdio.h
float lagrange(float x[],float y[],float xx) // 不写 i,不要 int n
{
int i,j;
float a[3],yy=0.0; // 定态数组 a[3]
for(i=0;i=2;i++) { a[i]=y[i];
for(j=0;j=2;j++)
if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]);
yy+=a[i]; }
return yy; }
void main() {
int i;
// 去掉函数原型声明
float x[3],y[3],xx,yy; // 定态数组 x[3],y[3]
for(i=0;i=2;i++) {
printf("x[%d]:",i); scanf("%f",x[i]); }
printf("\n");
for(i=0;i=2;i++)
{
printf("y[%d]:",i);scanf("%f",y[i]);
}
printf("\n"); printf("Input xx:");
scanf("%f",xx);
yy=lagrange(x,y,xx);
printf("x=%f,y=%f\n",xx,yy);
}
两元lagrange插值如何用C语言表示
//Lagrange插值多项式
//算法描述:
// 1、输入:插值节点控制数n,插值点序列(x_i,y_i),i=0,1,...n,要计算的函数点x。
// 2、for(i=0,1,2,...,n)
// {
// temp=1;
// for(j=0,1,...i,i+1,...n)
// !x要事先给定
// temp=temp*(x-x_j)/(x_i-x_j);
// }
// fx=fx+temp*y_i;
// }
#includestdio.h
#includestring.h
#define MAX_n 20
typedef struct tagPOINT
{
double x;
double y;
}POINT;
double Lagrange()
{
int n,i,j;
double x,temp,fx=0;
POINT points[MAX_n];
printf("Now,please input the n value: \n");
scanf("%d",n);
if(n=1||nMAX_n)
{
printf("The value of n should be between 2 and %d\n",MAX_n);
return 1;
}
printf("Now,please input the (x_i,y_i),i=0,...,%d\n",n-1);
for(i=0;in;i++)
scanf("%lf%lf",points[i].x,points[i].y);
printf("Now,please input the x value:\n");
scanf("%lf",x);
for(i=0;in;i++)
{
temp=1;
for(j=0;jn;j++)
if(i==j)continue;
else temp=temp*(x-points[j].x)/(points[i].x-points[j].x);
fx=fx+temp*points[i].y;
}
printf("So,when x=%lf,the Lagrange(%lf)=%lf\n",x,x,fx);
}
int main()
{
char s[10];
Lagrange();
gets(s);
while(strcmp(s,"exit"))
{
if(strcmp(s,"con")==0)
{
Lagrange();
getchar();
}
printf("继续(输入con),退出(输入exit)!\n");
gets(s);
}
return 0;
}
求双线性插值法的C语言程序!帮帮忙!拜托各位了!
a b
t
c d
就是两次线性插值,先在x方向插出t上下方的_t1、_t2,然后再用它们插出t来
float test(float x,float y)
{
float _t1,_t2,t;
_t1 = a+(b-a)*(x-ax)/(bx-ax);
_t2 = c+(d-c)*(x-cx)/(dx-cx);
t = _t1 +(_t2-_t1)*(y - ay);
return t;
}