你的浏览器版本过低,可能导致网站不能正常访问!
为了你能正常使用网站功能,请使用这些浏览器。

【经验之谈】基于STM32在单片机上做插值算法——拉格朗日插值的经验分享

[复制链接]
STMCU小助手 发布时间:2022-12-12 14:00

第一步,学你所学,不懂也得懂,最枯燥的数学公式来了     

在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式。数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现,不久后(1783年)由莱昂哈德·欧拉再次发现。1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起。







第二步,构造拉格朗日插值算法函数,并使用C/C++ 语言实现


  1. //预先定义插值节点的个数为1000个,根据控制台输入的个数num从而确定插值节点的个数

  2. const int N=1000;

  3. // 拉格朗日插值算法

  4. float lglr(float x[], float y[],int n,float t)

  5. {

  6.      float yResult=0.0;

  7.      //LValue[N]存放的是每次求解的插值基函数的通项

  8.      float LValue[N];

  9.      //循环变量k,m

  10.      int k,m;

  11.      //插值基函数中的上下累乘temp1,temp2

  12.      float temp1,temp2;

  13.      for(k=0;k<n;k++)

  14.      {

  15.          temp1=1.0;

  16.          temp2=1.0;

  17.          for(m=0;m<n;m++)

  18.          {

  19.              if(m==k)

  20.              {

  21.                  continue;

  22.              }

  23.              temp1 *= (t-x[m]);

  24.              temp2 *= (x[k]-x[m]);

  25.          }

  26.          LValue[k]=temp1/temp2;

  27.      }

  28.      for(int i=0;i<n;i++)

  29.      {

  30.          yResult += y*LValue;

  31.      }

  32.      return yResult;

  33. }

复制代码


第三步,测试线性插值效果

1.给写入的算法送入斜线的3个点,来预测此区间的其它点,发现预测出的线性值很好,线性度很饱满。


2.同样是给写入的算法送入正弦曲线3个点,来预测此区间的其它点,发现预测出的正弦曲线不太理想,相关度很差。


3.增加给写入的算法送入正弦曲线的点数到10个点,来预测此区间的其它点,发现预测出的正弦曲线已达到要求,相关度很好。

4.下面是使用matlab在绘制的图像,分别是线性3点、正弦3、5、10的和原始值对比的图像。


第四步,移植到单片机使用串口输出测试插值效果,效果还可以哦!



第四步  进阶完善,有点不甘心,使用Qt 搭了一个界面,做了一下可视化,你别说效果还可以,Qt 界面代码 放到了 gitee
有兴趣的来一起完善  



插值题外话:
关于插值算法:可以继续完善  牛顿插值、三次样条插值等,都很好玩!
---------------------
作者:一路向北lm

收藏 1 评论0 发布时间:2022-12-12 14:00

举报

0个回答

所属标签

相似分享

官网相关资源

关于
我们是谁
投资者关系
意法半导体可持续发展举措
创新与技术
意法半导体官网
联系我们
联系ST分支机构
寻找销售人员和分销渠道
社区
媒体中心
活动与培训
隐私策略
隐私策略
Cookies管理
行使您的权利
官方最新发布
STM32N6 AI生态系统
STM32MCU,MPU高性能GUI
ST ACEPACK电源模块
意法半导体生物传感器
STM32Cube扩展软件包
关注我们
st-img 微信公众号
st-img 手机版