matlab编写拟合

可运行的拟合matlab程序，附详细解释

%%拟合
%n阶多项式拟合
%参数：数据、阶数

clc
clear
close all

p = [17,38.7;16,41.4;21,51.5;20,51.8;25,52.3;25,53.9;26,54.7;27,56.6;29,57.1;28,55.2;29,57.5;30,56.5;31,58.5;32,59.4;34,60.2;35,61.5;35,62.4;36,63.8;37,65.8;38,66.4]; %数据
px = p(:,1); %数据x轴坐标
py = p(:,2); %数据y轴坐标
m = length(p); %数据个数

figure %绘制数据点图
plot(px,py,'ob')
hold on
grid on
axis equal
axis tight
title('数据点图')
xlabel('X')
ylabel('Y')

n = 1; %阶数
qy1 = (0); %某阶数，拟合多项式单项式值表1，初始化
qy2 = (0); %某数据，拟合多项式值单项式值表2，初始化
qy = (0); %拟合多项式值表，初始化
q = polyfit(px,py,n); %拟合参数
for i = 0:n %阶数
    qy3 = q(i+1).*px.^(n-i); %当前阶数，拟合多项式单项式值
    for j = 1:m %数据个数
        qy1(j,i+1) = qy3(j,1); %当前阶数，拟合多项式单项式值，记录在拟合多项式单项式值表1中
    end
end
for k = 1:m %数据个数
    qy2 = qy1(k,:); %当前数据，拟合多项式值单项式值，记录在拟合多项式单项式值表2中
    qy(k,1) = sum(qy2); %当前数据，拟合多项式值，即拟合多项式值单项式值总和，记录在拟合多项式值表中
end

figure %绘制拟合图
plot(px,py,'ob')
hold on
plot(px,qy,'-r')
hold on
grid on
axis equal
axis tight
title('拟合图')
xlabel('X')
ylabel('Y')

运行结果