(由于不方便上传公式,有部分直接截图)
2. Polyval函数
对应多项式拟合,多项式求值函数也相当重要,可以通过设定向量,快速定义如下多项式
据此,可以大概猜测到polyval调用的形式了
y = polyval(p,x)
如上所述
[y,delta] = polyval(p,x,S)
S即为polyfit函数返回的误差估计,具体delta的数学意义本人暂时不清楚
y = polyval(p,x,[],mu)
[y,delta] = polyval(p,x,S,mu)
结合polyfit的例子,可以尝试将拟合曲线画出来
x = -10 : 0.2 : 10;
y = 0.59*x.*x + 1.4*x+2+2*(rand(1,length(x))-0.5);
plot(x,y,'.')
[p, S] = polyfit(x,y,2)
yfit = polyval(p,x);
hold on
plot(x,yfit, 'r');
结果如下图