我是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);
答案 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])));