您的位置:

c语言二元插值,二维线性插值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;

}