热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

c语言二元插值,二维线性插值c语言

本文目录一览:1、用C语言实现拉格朗日插值、牛顿插值、等距结点插值算法2、用C语言编写一个线性插值程序3、c语言中插值排序法怎么打程序4、用拉格朗日插值法(二次插值)来

本文目录一览:


  • 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;

}


推荐阅读
  • 题目描述:给定n个半开区间[a, b),要求使用两个互不重叠的记录器,求最多可以记录多少个区间。解决方案采用贪心算法,通过排序和遍历实现最优解。 ... [详细]
  • golang常用库:配置文件解析库/管理工具viper使用
    golang常用库:配置文件解析库管理工具-viper使用-一、viper简介viper配置管理解析库,是由大神SteveFrancia开发,他在google领导着golang的 ... [详细]
  • 本文详细探讨了KMP算法中next数组的构建及其应用,重点分析了未改良和改良后的next数组在字符串匹配中的作用。通过具体实例和代码实现,帮助读者更好地理解KMP算法的核心原理。 ... [详细]
  • 本文将介绍如何编写一些有趣的VBScript脚本,这些脚本可以在朋友之间进行无害的恶作剧。通过简单的代码示例,帮助您了解VBScript的基本语法和功能。 ... [详细]
  • 1:有如下一段程序:packagea.b.c;publicclassTest{privatestaticinti0;publicintgetNext(){return ... [详细]
  • 本文介绍了Java并发库中的阻塞队列(BlockingQueue)及其典型应用场景。通过具体实例,展示了如何利用LinkedBlockingQueue实现线程间高效、安全的数据传递,并结合线程池和原子类优化性能。 ... [详细]
  • 本文详细介绍了Java编程语言中的核心概念和常见面试问题,包括集合类、数据结构、线程处理、Java虚拟机(JVM)、HTTP协议以及Git操作等方面的内容。通过深入分析每个主题,帮助读者更好地理解Java的关键特性和最佳实践。 ... [详细]
  • 本文探讨了如何在给定整数N的情况下,找到两个不同的整数a和b,使得它们的和最大,并且满足特定的数学条件。 ... [详细]
  • 从 .NET 转 Java 的自学之路:IO 流基础篇
    本文详细介绍了 Java 中的 IO 流,包括字节流和字符流的基本概念及其操作方式。探讨了如何处理不同类型的文件数据,并结合编码机制确保字符数据的正确读写。同时,文中还涵盖了装饰设计模式的应用,以及多种常见的 IO 操作实例。 ... [详细]
  • 本实验主要探讨了二叉排序树(BST)的基本操作,包括创建、查找和删除节点。通过具体实例和代码实现,详细介绍了如何使用递归和非递归方法进行关键字查找,并展示了删除特定节点后的树结构变化。 ... [详细]
  • 本文详细介绍了C语言中链表的两种动态创建方法——头插法和尾插法,包括具体的实现代码和运行示例。通过这些内容,读者可以更好地理解和掌握链表的基本操作。 ... [详细]
  • 深入理解Java泛型:JDK 5的新特性
    本文详细介绍了Java泛型的概念及其在JDK 5中的应用,通过具体代码示例解释了泛型的引入、作用和优势。同时,探讨了泛型类、泛型方法和泛型接口的实现,并深入讲解了通配符的使用。 ... [详细]
  • 本文详细介绍了Linux系统中init进程的作用及其启动过程,解释了运行级别的概念,并提供了调整服务启动顺序的具体步骤和实例。通过了解这些内容,用户可以更好地管理系统的启动流程和服务配置。 ... [详细]
  • 本文介绍了几种不同的编程方法来计算从1到n的自然数之和,包括循环、递归、面向对象以及模板元编程等技术。每种方法都有其特点和适用场景。 ... [详细]
  • 并发编程:深入理解设计原理与优化
    本文探讨了并发编程中的关键设计原则,特别是Java内存模型(JMM)的happens-before规则及其对多线程编程的影响。文章详细介绍了DCL双重检查锁定模式的问题及解决方案,并总结了不同处理器和内存模型之间的关系,旨在为程序员提供更深入的理解和最佳实践。 ... [详细]
author-avatar
_快樂Smile_903
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有