在数学与工程领域中，数值积分是一种重要的工具。其中，复化梯形公式因其简单易用而广受欢迎。它通过将区间分成若干小段，并用梯形面积近似每一段的积分值来提高精度。今天，我们用Matlab编写一个程序来实现这一算法！🎉
首先，我们需要定义函数和积分区间。假设我们要计算函数`f(x) = x^2`从`a=0`到`b=1`的积分。接着，我们将区间划分为`n`等分，构建梯形并累加面积。代码逻辑清晰，易于调试。💡
以下是核心代码片段：
```matlab
function I = trapz_comp(f, a, b, n)
h = (b - a) / n;
x = a:h:b;
y = f(x);
I = h (sum(y) - 0.5 (y(1) + y(end)));
end
```
运行后，你会发现结果非常接近真实值！🎯
如果你对代码有疑问或需要进一步优化，请随时留言~ 📩💬
Matlab编程 数值积分 学习分享