MATLAB 多项式拟合代码解析 - 3次、4次、2次拟合曲线绘制 - 常规

MATLAB 多项式拟合代码解析 - 3次、4次、2次拟合曲线绘制

本文将详细解析一段 MATLAB 代码，展示如何使用 polyfit 和 polyval 函数进行多项式拟合，并绘制 3次、4次、2次拟合曲线。同时解释 poly2str 和 poly2sym 函数的用法，帮助您理解多项式拟合的基本原理和实现方法。

代码：

X = [0.0 0.1 0.2 0.3 0.5 0.8 1.0];
Y = [1.0 0.41 0.50 0.61 0.91 2.02 2.46];

p1 = polyfit(X, Y, 3);
p2 = polyfit(X, Y, 4);
Y1 = polyval(p1, X);
Y2 = polyval(p2, X);

poly2str(p1, 'x')  % 将 p1 多项式系数转换为字符串
poly2str(p2, 'x')  % 将 p2 多项式系数转换为字符串

plot(X, Y, 'r*', X, Y1, 'b-.', X, Y2, 'g--')

p3 = polyfit(X, Y, 2);
Y3 = polyval(p3, X);

poly2str(p3, 'x')  % 将 p3 多项式系数转换为字符串

f1 = poly2sym(p1);
f2 = poly2sym(p2);
f3 = poly2sym(p3);

plot(X, Y, 'r*', X, Y1, 'b-.', X, Y2, 'g--', X, Y3, 'm--')
legend('数据点', '3次多项式拟合', '4次多项式拟合', '2次多项式拟合')
xlabel('X轴'), ylabel('Y轴')

代码逐句解析：

X = [0.0 0.1 0.2 0.3 0.5 0.8 1.0]; Y = [1.0 0.41 0.50 0.61 0.91 2.02 2.46];
- 定义 X 和 Y 数组，分别代表自变量和因变量。
p1 = polyfit(X, Y, 3); p2 = polyfit(X, Y, 4);
- 使用 polyfit 函数拟合 X 和 Y，其中参数 3 和 4 分别表示拟合多项式的次数，p1 为 3 次多项式系数，p2 为 4 次多项式系数。
Y1 = polyval(p1, X); Y2 = polyval(p2, X);
- 使用 polyval 函数计算出 p1 和 p2 对应 X 值的 Y 值。
poly2str(p1, 'x') poly2str(p2, 'x')
- 将 p1 和 p2 多项式系数转换为字符串，方便查看多项式表达式。
plot(X, Y, 'r*', X, Y1, 'b-.', X, Y2, 'g--')
- 绘制 X 和 Y 的散点图，以及 p1 和 p2 对应 X 值的拟合曲线。
- 'r*' 表示红色星形，'b-.' 表示蓝色虚线，'g--' 表示绿色虚线。
p3 = polyfit(X, Y, 2); Y3 = polyval(p3, X);
- 使用 polyfit 函数拟合 X 和 Y，p3 为 2 次多项式系数，使用 polyval 函数计算出 p3 对应 X 值的 Y 值。
poly2str(p3, 'x')
- 将 p3 多项式系数转换为字符串。
f1 = poly2sym(p1); f2 = poly2sym(p2); f3 = poly2sym(p3);
- 将 p1、p2、p3 多项式系数转换为符号函数，方便进一步的数学运算和分析。
plot(X, Y, 'r*', X, Y1, 'b-.', X, Y2, 'g--', X, Y3, 'm--')
- 绘制 X 和 Y 的散点图，以及 p1、p2、p3 对应 X 值的拟合曲线。
- 'm--' 表示洋红色虚线。
legend('数据点', '3次多项式拟合', '4次多项式拟合', '2次多项式拟合') xlabel('X轴'), ylabel('Y轴')
- 添加图例和坐标轴标签，使图形更加清晰易懂。

总结：

这段代码展示了 MATLAB 中使用 polyfit 和 polyval 函数进行多项式拟合，以及使用 poly2str 和 poly2sym 函数对拟合结果进行处理和展示的完整过程。通过该代码示例，您可以学习到多项式拟合的基本原理和实现方法，并了解如何使用 MATLAB 工具箱进行数据分析和可视化。