23.多项式与非多项式曲线拟合对比（matlab程序）

1.简述

      

拟合标准：

(1)原始数据向量与拟合向量之间的距离最小，该距离的度量一般使用误差平方和表示，即均方误差：R=||Q-Y||22

(2)当均方误差最小时，说明构造的拟合向量与原始向量最为接近，这种曲线拟合的方法称为最小二乘法

(3)计算均方误差最小时的拟合系数，可以通过微积分中求解极值的方法实现

多项式拟合

(1) 多项式拟合在matlib中可以使用ployfit函数求解多项式系数，通过ployval函数求解拟合多项式在某数据点处的值

(2) 比如：拟合正弦函数曲线，使用4次多项式拟合

2.代码

%% 多项式拟合正弦函数曲线
% 原始数据点
x0 = -pi:0.1:pi;
y0 = sin(x0);
% 4次多项式拟合
p0 = polyfit(x0, y0, 4);
y1 = polyval(p0, x0);
plot(x0, y0, x0, y1, 'r');

%% 多项式拟合效果
% 阶次越高，并不一定拟合效果越好
x1 = -1:0.2:1;
y1 = 1./(1 + 25*x1.^2);
% 分别使用3次、5次、8次多项式拟合
p3 = polyfit(x1, y1, 3);
p5 = polyfit(x1, y1, 5);
p8 = polyfit(x1, y1, 8);

x  = -1:0.01:1;
y  = 1./(1 + 25*x.^2);
y3 = polyval(p3, x);
y5 = polyval(p5, x);
y8 = polyval(p8, x);
figure
plot(x, y, x, y3, 'r-', x, y5, 'm:', x, y8, 'b--');
legend('原始', '3次', '5次', '8次');

%% 使用非多项式拟合方法
% 首先建立拟合选项结构体
options = fitoptions('Method', 'NonlinearLeastSquare');
options.Lower = [-Inf, -Inf, -Inf];
options.Upper = [Inf, Inf, Inf];
% 通过fittype建立非线性拟合模型
type    = fittype('a/(b + c*x^n)', 'problem', 'n', 'options', options);
% 拟合
[cfun gof] = fit(x1', y1', type, 'problem', 2);
% 拟合效果
ynp = feval(cfun, x);
figure
hl = plot(x, y, 'k');
set(hl, 'LineWidth', 5);
hold on
plot(x, ynp, 'r');

3.运行结果