如何在matlab中集成这个函数

时间:2017-01-08 00:01:11

标签: matlab integration

我是matlab的新手。

如何整合这行代码? ?

p2= polyfit(x,y,length(x));
from= x(1);
to= x(length(x));

我需要整合p2

我尝试使用积分功能:

 value = integral(p2,from,to);

但我得到了

  

使用积分时出错(第82行)第一个输入参数必须是函数   处理

Error in poly_integral (line 5)
 value = integral(p2,from,to);

2 个答案:

答案 0 :(得分:0)

这是因为代码中的p2不是函数。它只是系数的向量。 integral的第一个参数需要处理要集成的函数。

从您的代码判断,您似乎想要定义一个评估多项式p2的函数。如果是这样,您可以执行以下示例:

% take an example set of x and y 
x = linspace(0, pi, 1000); % uniform samples between 0 to pi
y = sin(x); % assume, for sake of example, output is sine function of input

% polynomial fit
p2 = polyfit(x,y,4); % 4th order polynomial
% Note that, in general, the order should be much smaller than length(x).
% So you probably should review this part of your code as well.

% define a function to evaluate the polynomial
fn = @(x) polyval(p2, x);
% this means: fn(x0) is same as polyval(p2, x0)

% compute integral
value = integral(fn,x(1),x(end));

答案 1 :(得分:0)

您可以使用polyint函数获取多项式系数,以便精确积分多项式:

p2  = polyfit(x,y,length(x));
int = diff(polyval(polyint(p2),x([1 end])));