用于嵌入式系统的连续定积分离散化方法

最新推荐文章于 2024-06-09 15:03:53 发布

ayuan0211

最新推荐文章于 2024-06-09 15:03:53 发布

阅读量523

点赞数 9

文章标签：算法

本文链接：https://blog.csdn.net/weixin_38342580/article/details/134988726

版权

矩形积分

梯形积分

辛普森积分

矩形积分

1.优点：计算量小；缺点：不够精确

2.公式：

$\int_{a}^{b}f(x)dx \approx (f(b)/2)*(b-a)$

若在嵌入式应用过程中，则自变量为时间t，则公式可写为：

$\int_{a}^{b}f(t)dt \approx f(b)*Ts$

式中Ts为嵌入式系统程序执行间隔Ts。

代码实现：

sum +=f(N)*Ts;

梯形积分

1.与矩形积分计算量、精度均无明显差别

2.公式：

$\int_{a}^{b}f(x)dx \approx (f(b)+f(a))*(b-a)/2$

若自变量为时间t,则公式可写为： $\int_{a}^{b}f(t)dt \approx (f(b)+f(a))*Ts/2$ ，Ts为程序采样间隔。

代码：

sum += （f(N)+f(N-1)）*Ts/2; f(N)N时刻的采样值，f(N-1)上一时刻的采样值。

辛普森积分

1.优点：更准确；缺点，计算量大，占用CPU资源，适用于对某一性能指标精度要求较高的场合。

2.公式：

$\int_{a}^{b}f(x)dx \approx \frac{h}{3}[f(a))+f(b)+4\sum_{i=1}^{n/2}f(x_{2i-1})+2\sum_{i=1}^{n/2-1}f(x_{2i})]$

其中h=(b-a)/n,n是子区间数目，且为偶数。

n的具象示意：嵌入式系统程序调用时间间隔为T，积分调用时间间隔为T/2，则n=T/(T/2)=2;

实际举例：当n=2时，

3.程序代码：

// function
double f(double x) {
  return 1/(1+x*x);
}

// Simpson's rule
double simpson(double a, double b, int n) {
  double h = (b-a)/n; // width of each subinterval
  double sum = f(a) + f(b); // first and last terms
  for (int i = 1; i < n; i += 2) {
    sum += 4*f(a+i*h); // odd terms
  }
  for (int i = 2; i < n-1; i += 2) {
    sum += 2*f(a+i*h); // even terms
  }
  return sum*h/3; // final result
}

参考资料：c辛普森积分法 - 知乎 (zhihu.com)