AKima 插值实现

最新推荐文章于 2024-08-11 14:33:36 发布
最新推荐文章于 2024-08-11 14:33:36 发布
阅读量1.1w 收藏 13
点赞数 3
分类专栏： Algorithm
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/bluels01/article/details/43561131
版权

测量数据的内插已有各种方法,如线性内插、多项式内插、样条函数插值等,但这里的Akima插值法具有独特的优点。线性内插只顾及其附近两点的影响。多项式内插时,低阶多项式由于参数较少,内插精度很低,而使用高阶多项式又会使解不稳定,出现“龙格”现象,即内插函数在插值点与实际数据符合得很好,而在插值点外出现较大的偏差。因此人们又在多项式的基础上发展了分片多项式,即样条函数。样条函数既保持了多项式运算简单的特点,又避免了多项式阶数较高时数值不稳定的缺点,因而得到了广泛的应用。但在样条函数插值中,确定任何一个小区间上的多项式,都要考虑所有数据点对它的影响。这不仅扩大了误差传播的范围,还增加了不少工作量。有时只用内插点附近的几个数据点作为控制点来内插。Aikma插值法和三次样条函数一样考虑了要素导数值的效应,因而得到的整个插值曲线是光滑的。三次样条函数插值法具有最小模、最佳最优逼近和收敛的特性,而Aikma插值法所得曲线比样条函数插值曲线更光顺,更自然。[1]

Akima相关实现的代码：

Akima's original paper:``A new method of interpolation and smooth curve fitting based on local procedures'', Journal of ACM 17, 4 (1970), 589-602

http://student.ndhu.edu.tw/~u9111023/akima.pdf

C implementation

http://packages.debian.org/source/squeeze/spline

C++ implementation

http://www.koders.com/cpp/fid1393B9D668316C1700966643DE0609660B9CB13A.aspx?s=%22Brian+Smith%22

C# implementation

https://github.com/mathnet/mathnet-numerics/blob/master/src/Numerics/Interpolation/Algorithms/AkimaSplineInterpolation.cs

Delphi implementation (see procedure BuildAkimaSpline in delphi/src/spline3.pas)

http://www.alglib.net/translator/re/alglib-2.6.0.delphi.zip

Akima's Fortran 66 implementation

http://cran.r-project.org/web/packages/akima/

Fortran 90 implementation

http://miyoshi.googlecode.com/svn-history/r72/trunk/common/common.f90

Java implementation

http://code.google.com/p/msnlab4/source/browse/trunk/src/br/edu/ufcg/msn/interpolacao/spline/AkimaCubicSplineInterpolator.java?spec=svn83&r=83

Lisp implementation "for AutoCAD 2d-Polylines"

http://autocad.xarch.at/code/candido/akima.lsp

Matlab implementation

http://www.mathworks.se/matlabcentral/fileexchange/1814-akima-interpolation

Pascal implementation (program description)

http://jean-pierre.moreau.pagesperso-orange.fr/Pascal/akima_pas.txt

Python implementation

http://www.lfd.uci.edu/~gohlke/code/akima.py.html

VB6 implementation (see subroutine BuildAkimaSpline in vb6/src/spline3.bas)

http://www.alglib.net/translator/re/alglib-2.6.0.vb6.zip  [2]

参考文献：

[1]聂桂根, and 万剑华. "Akima 插值法在测量中的应用." 测绘科技动态 3 (1998): 31-34.

[2]http://www.mathworks.com/matlabcentral/fileexchange/1814-akima-interpolation

bluels01
关注
  • 3
    点赞
  • 踩
  • 13
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
Python Akima 插值法教程 | 关键原理和实践
叶庭云 成为自己的光
06-30 2796
Python Akima 插值法教程 | 关键原理和实践
commons-math3-3.6.1源码学习
Poo_Chai的博客
01-15 3848
commons-math3-3.6.1是一个很实用的工具包，JAVA中矩阵运算，云计算CloudSim平台都能用到。 学习commons-math3-3.6.1可到网站去。 以下为包中所含函数： build.xmlcheckstyle.xmlclirr-ignored.xmlfindbugs-exclude-filter.xmlpmd-ruleset.xmlpom.xmlsr
各种插值方法的实现方式
twinkle的博客
12-03 1805
每种插值方法都有其适用场景和局限性。选择哪种方法取决于数据的特性、问题的复杂性以及性能要求。对于不同的应用，可能需要尝试多种插值方法，调整插值参数，以找到最适合特定数据集和需求的方法。对于实际应用，您可能还需要进行数据的预处理和后处理，以确保插值结果的准确性和可靠性。
二维Akima插值Matlab程序
12-27
自己实现的二维Akima插值，包含akima_interp2和makima_interp2两个函数，分别为原始版本（由Akima提出）的插值算法和Cleve修正的算法，同时提供了测试脚本（示例）
akima 插值拟合算法 Python/C++/C版本
silent_dusbin的博客
06-20 2260
鉴于C站上Akima算法文章大部分要VIP观看或者下载，即使是付费也有质量不佳，浪费Money也浪费时间。笔者更具查到的资料分享给大家。
样条之Akima光滑插值函数
weixin_34413357的博客
10-18 708
 核心代码： 1 ////////////////////////////////////////////////////////////////////// 2 // Akima光滑插值 3 // t - 存放指定的插值点的值 4 // s[] - 一维数组，长度为5，其中s(0)，s(1)，s(2)，s(3)返回三次多项式的系数， 5 // ...
akima_akima插值函数_
10-02
`akima`插值函数是其中一种方法，它由Akima Hidehiko于1970年提出，特别适合处理有噪声的数据。本文将深入探讨`akima`插值函数的概念、工作原理以及在实际应用中的优势。 `akima`插值方法是一种四次多项式插值算法...
Akima.rar_akima_akima插值_插值
07-14
在MATLAB环境中，实现Akima插值是一种高效且灵活的方法。 Akima插值的核心在于它能够处理数据的突变，同时保持插值函数的连续性和一阶导数的连续性。这使得Akima插值在处理具有局部变化的复杂数据时表现优秀，如...
Akima 插值：通过平面上的给定点对平滑曲线进行插值。-matlab开发
06-01
在 MATLAB 中，你可以使用内置函数 `akima` 来实现 Akima 插值。以下是一般的使用步骤： 1. 加载数据：首先，你需要将数据存储在两个向量 x 和 y 中，然后调用 `akima` 函数创建一个插值对象。 ```matlab x = ...
akimai.zip_Akima 插值 算法_Akima插值算法_akima_akima插值_非线性插值
07-14
可以利用Akima插值算法对非线性问题进行插值计算
Akima算法
06-11
最近有人问Akima算法，这里是源代码，已经测试。
插值法VB源程序.docx
10-25
插值法VB源程序.docx
matlab 拉格朗日插值及akima插值
03-11
插值程序代码matlab版function，可轻松改为fortran或c
matlab_包含用于akima插值的源代码和Matlab示例
05-29
在提供的项目中，`akima-interpolation-in-matlab`可能包含了实现Akima插值的MATLAB脚本，以及展示如何使用这些脚本的示例。这些示例可能涵盖了不同的数据集和应用场景，帮助用户理解算法的工作原理，并学习如何将其...
数学建模-插值算法原理笔记
Claire_cf的博客
07-31 611
目的 比赛中常常需要根据已知的函数点进行数据、模型的处理和分析，而有时候现有的数据是极少的，不足以支撑分析的进行，这时就可以使用一些方法“模拟产生”一些新的但又比较靠谱的值，这就是插值的作用。 概念 设函数y=f(x)y=f(x)y=f(x)​在区间[a,b][a,b][a,b]​上有定义，且已知在点a≤x0<x1<…<xn≤ba≤x_0<x_1<…<x_n≤ba≤x0​<x1​<…<xn​≤b​上的值分别为y0,y1,…,yny_0,y_1,…,y_n
Akima 插值拟合算法VBA版
duet的博客
09-05 151
源代码下载：Akima样条插值.bas。改写为VBA版，共需要的同学使用。本文参考silent_dusbin的。
阿克玛插值(2D)
yliang的专栏
01-01 1998
#ifndef INTERPOLATION_HEADER#define INTERPOLATION_HEADER #include#include#include#include using namespace std; #define DELTA 1e-6#define VERY_CLOSE(x,y) fabs((x)-(y)) templat
Python数值计算（22）——akima插值
最新发布
cdinten的专栏
08-11 599
前面介绍到的很多插值算法，在多项式阶次较高时，都会存在龙格现象，最直观的表现就是导数的变化很急剧，所以一般采用的方式是使用分段的三次插值，例如前面介绍到的三次样条曲线，这里再介绍Akima插值算法，通过该方法生成的插值曲线，具有更加平滑的特性。
Qt实现akima 插值
06-02
在Qt中实现Akima插值，可以使用Qwt库提供的插值函数。Qwt是一个基于Qt的数据可视化库，提供了许多用于绘制曲线、图表等的控件和函数。以下是一个使用Qwt实现Akima插值的示例代码： ```c++ #include <QApplication> #include <qwt_plot.h> #include <qwt_plot_curve.h> #include <qwt_point_data.h> #include <qwt_plot_grid.h> #include <qwt_legend.h> #include <qwt_scale_draw.h> #include <qwt_scale_widget.h> #include <qwt_plot_renderer.h> #include <qwt_plot_zoomer.h> #include <qwt_plot_panner.h> #include <qwt_picker_machine.h> #include <qwt_picker.h> #include <qwt_plot_magnifier.h> #include <qwt_text.h> #include <qwt_math.h> #include <qwt_scale_engine.h> #include <qwt_interval.h> #include <qwt_plot_renderer.h> #include <qwt_plot_canvas.h> #include <qwt_plot_marker.h> #include <qwt_plot_akima_curve.h> #include <math.h> const int N = 20; const double x[N] = { 0.0, 0.1625, 0.325, 0.4875, 0.65, 0.8125, 0.975, 1.1375, 1.3, 1.4625, 1.625, 1.7875, 1.95, 2.1125, 2.275, 2.4375, 2.6, 2.7625, 2.925, 3.0875 }; const double y[N] = { 0.0, 0.2624, 0.5051, 0.7054, 0.8513, 0.9332, 0.9434, 0.8751, 0.7238, 0.489, 0.1681, -0.289, -0.9014, -1.6304, -2.4764, -3.44, -4.5207, -5.717, -7.0262, -8.4459 }; int main(int argc, char *argv[]) { QApplication a(argc, argv); QwtPlot plot; plot.setTitle("Akima Interpolation"); plot.setCanvasBackground(Qt::white); QwtPlotGrid *grid = new QwtPlotGrid; grid->enableXMin(true); grid->enableYMin(true); grid->setMajorPen(Qt::gray, 0, Qt::DotLine); grid->attach(&plot); QwtLegend *legend = new QwtLegend; legend->setItemMode(QwtLegend::CheckableItem); plot.insertLegend(legend, QwtPlot::RightLegend); QwtPlotAkimaCurve *curve = new QwtPlotAkimaCurve("Akima"); curve->setSamples(x, y, N); curve->attach(&plot); plot.resize(600,400); plot.show(); return a.exec(); } ``` 在上面的代码中，我们使用了QwtPlotAkimaCurve类来实现Akima插值，该类继承自QwtPlotCurve类，具有更高的插值精度和更少的振荡。在代码中，我们首先定义了一组数据点，然后创建了一个QwtPlot对象，设置了标题和背景色，并创建了一个QwtPlotGrid对象用于绘制网格线。接下来，我们创建了一个QwtPlotAkimaCurve对象，并通过setSamples函数将数据点传递给曲线对象。最后，将曲线对象附加到QwtPlot对象上，并显示出来。 需要注意的是，在使用Qwt库前，需要将其添加到Qt项目中。可以在.pro文件中添加以下代码来实现： ```c++ QT += qwt ``` 以上是一个简单的Qt实现Akima插值的示例代码，可以根据自己的需要进行修改和扩展。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值