插值算法C实现（一元全区间）

一元全区间插值

算法：

double lgr(x,y,n,t)int n;double t,x[],y[];{ int i,j,k,m;double z,s;z=0.0;if (n<1) return(z);if (n==1) { z=y[0];return(z);}if (n==2){ z=(y[0]*(t-x[1])-y[1]*(t-x[0]))/(x[0]-x[1]);return(z);}i=0;while ((x[i]<t)&&(i<n)) i=i+1;k=i-4;if (k<0) k=0;m=i+3;if (m>n-1) m=n-1;for (i=k;i<=m;i++){ s=1.0;for (j=k;j<=m;j++)if (j!=i) s=s*(t-x[j])/(x[i]-x[j]);z=z+s*y[i];}return(z);}

例子：

#include "8lgr.c"#include "stdio.h"main(){ double t,z;static double x[10]={0.10,0.15,0.25,0.40,0.50,0.57,0.70,0.85,0.93,1.00};static double y[10]={0.904837,0.860708,0.778801,0.670320,0.606531,0.565525,0.496585,0.427415,0.394554,0.367879};t=0.63; z=lgr(x,y,10,t);printf("\n");printf("t=%6.3f,   z=%e\n",t,z);printf("\n");}

结果：

插值算法C实现（一元全区间）相关推荐

插值算法C实现（二元全区间）
二元全区间插值程序: double slgq(x,y,z,n,m,u,v)int n,m;double x[],y[],z[],u,v;{ int ip,ipp,i,j,l,iq,iqq,k;dou ...
阿当姆斯校正程序代码MATLAB,全区间积分的阿当姆斯预报校正法(常微分方程组的求解)...
/* 代码作者:不详代码整理者:设计天下 MySDN网站算法天下工作室功能:全区间积分的阿当姆斯预报校正法(常微分方程组的求解) */ #include "stdio.h" ...
20世纪最好的十大算法、算法笔记（2008-11-15 22:16:57、2011-04-21 19:29:05）
Algorithm(算法)一词与9世纪的阿拉伯学者al-Khwarizmi有关,他写的书<al-jabr w'al muqabalah>(代数学)演变成为现在中学的代数教科书.Ad-Khw ...
java数值计算算法编程,Java数值计算算法编程
第1章 Java与数值计算. 1．1 数值计算中存在的问题 1．2 用Java实现数值计算算法的要点 1．3 实数类设计与实现第2章复数运算 2．1 复数类设计 2．2 复数乘法 2．3 复数除法 ...
计算机数值计算原理,C#数值计算算法编程
第1章 C#与数值计算. 1．1 数值计算中存在的问题 1．2 用C#实现数值计算算法的要点第2章复数运算 2．1 复数类设计 2．2 复数乘法 2．3 复数除法 2．4 复数的模 2．5 复数的 ...
基于c#的相关性分析_不同区间衰退路径下锂离子电池的性能相关性及温度适用性分析...
锂离子电池在不同容量区间和不同工况应力下使用衰退性能不同.北京交通大学国家能源主动配电网技术研发中心.北京电动车辆协同创新中心.北京新能源汽车股份有限公司的研究人员孙丙香.刘佳.韩智强.任鹏博.张维戈 ...
激活函数和全连接层——基于Pytorch
1.激活函数 1.1.什么是激活函数? 神经网络中的每个神经元接受上一层的输出值作为本神经元的输入值,并将处理结果传递给下一层(隐藏层或输出层).在多层神经网络中,上层的输出和下层的输入之间具有一个函 ...
多项式插值之Lagrange、PCHIP与Spline以及BD-Rate和BD-PSNR的计算
做视频编解码算法的话一般都会接触到一个目前比较通用的客观评价指标,即 Bjøntegaard Delta(BD),如果关注码率变化的话就是 BD-Rate,如果关注 PSNR 变化的话就是BD-PSN ...
matlab的插值方法
1. Lagrange插值插值是在已知数据之间寻找估计值的过程.在信号处理和图像处理中,插值极其常用. 类型很多:比如多项式插值,一.二.三维插值,样条插值等. 方法介绍: 对给定的n个插值点x1, ...

插值算法C实现（一元全区间）

一元全区间插值

插值算法C实现（一元全区间）相关推荐

最新文章

热门文章