以编程方式拟合

适用于多项式模型的 MATLAB 函数

有两个 MATLAB^® 函数可通过多项式对您的数据建模。

如果您尝试对物理情况建模，务必考虑特定阶次的模型是否对您的情况有意义。

带非多项式项的线性模型

此示例说明如何使用含有非多项式项的线性模型拟合数据。

若多项式函数并未得出适合您数据的满意模型，您可尝试使用带非多项式项的线性模型。以如下函数为例，它在参数 $$a_0$$、$$a_1$$ 和 $$a_2$$ 中为线性，而在 $$t$$ 数据中为非线性：

您可通过构建及求解一组联立方程并为参数求解，计算未知系数 $$a_0$$、$$a_1$$ 和 $$a_2$$。以下语法通过构建一个设计矩阵实现此目的，该矩阵中的每一列代表用于预测响应（模型中的项）的变量，每一行对应于这些变量的一个观测值。

输入 t 和 y 作为列向量。

t = [0 0.3 0.8 1.1 1.6 2.3]';
y = [0.6 0.67 1.01 1.35 1.47 1.25]';

构建设计矩阵。

X = [ones(size(t))  exp(-t)  t.*exp(-t)];

计算模型系数。

a = X\y

a = 3×1

    1.3983
   -0.8860
    0.3085

因此，该数据的模型由以下公式提供：

现在以固定间隔点评估模型，并以原始数据绘制模型。

T = (0:0.1:2.5)';
Y = [ones(size(T))  exp(-T)  T.*exp(-T)]*a;
plot(T,Y,'-',t,y,'o'), grid on
title('Plot of Model and Original Data')

多重回归

此示例说明如何使用多重回归对具有多个预测变量的函数进行数据建模。

若 y 是具有多个预测变量的函数，则必须对表示各变量之间关系的矩阵方程进行扩展，以容纳额外的数据。这称为多重回归。

为多个 $$x_1$$ 和 $$x_2$$ 值测量对应的量 $$y$$。将这些值分别存储在向量 x1、x2 和 y 中。

x1 = [.2 .5 .6 .8 1.0 1.1]';
x2 = [.1 .3 .4 .9 1.1 1.4]';
y  = [.17 .26 .28 .23 .27 .24]';

此数据的模型采用以下形式：

多重回归可通过最大限度地减小数据与模型偏差的平方和（最小二乘拟合），对未知系数 $$a_0$$、$$a_1$$ 和 $$a_2$$ 求解。

通过构建设计矩阵 X，构建和求解一组联立方程。

X = [ones(size(x1))  x1  x2];

使用反斜杠运算符对参数求解。

a = X\y

a = 3×1

    0.1018
    0.4844
   -0.2847

数据的最小二乘拟合模型为

为了验证该模型，请求出数据与模型偏差绝对值的最大值。

Y = X*a;
MaxErr = max(abs(Y - y))

MaxErr = 
0.0038

该值远小于任何数据值，表明该模型能够准确贴合数据。

以编程方式拟合

此示例说明如何使用 MATLAB 函数执行以下操作：

• 计算相关系数

• 对数据进行多项式拟合

• 绘制和计算置信边界

从 census.mat 加载样本人口普查数据，它含有 1790 年至 1990 年的美国人口数据。

load census

这会向 MATLAB 工作区添加以下两个变量。

绘制数据图。

plot(cdate,pop,'ro')
title('U.S. Population from 1790 to 1990')

该图呈现出一个明显的模式，表明各变量之间存在高相关性。

corrcoef(cdate,pop)

ans = 2×2

    1.0000    0.9597
    0.9597    1.0000

[p,ErrorEst] = polyfit(cdate,pop,2);

pop_fit = polyval(p,cdate,ErrorEst);

plot(cdate,pop_fit,'-',cdate,pop,'+');
title('U.S. Population from 1790 to 1990')
legend('Polynomial Model','Data','Location','NorthWest');
xlabel('Census Year');
ylabel('Population (millions)');

res = pop - pop_fit;
figure, plot(cdate,res,'+')
title('Residuals for the Quadratic Polynomial Model')

[pop_fit,delta] = polyval(p,cdate,ErrorEst);

plot(cdate,pop,'+',...
     cdate,pop_fit,'g-',...
     cdate,pop_fit+2*delta,'r:',...
     cdate,pop_fit-2*delta,'r:'); 
xlabel('Census Year');
ylabel('Population (millions)');
title('Quadratic Polynomial Fit with Confidence Bounds')
grid on