clc;
clear;
% 定义 x 和 y
x = [-3 -1 1 3];
y = [15 5 1 5];
% 使用 polyfit 函数进行拟合
p = polyfit(x, y, 3);
% 生成新的数据点用于绘图
xp = linspace(min(x), max(x), 100);
yp = polyval(p, xp);
% 绘制原始数据和拟合曲线
plot(x, y, 'o', xp, yp, '-');
legend('原始数据', '拟合曲线');
xlabel('x');
ylabel('y');
title('三次多项式');