A=[216.85 232.49 245.96 259.87 267.9 287.62 310.47 332.21 354.99 384.52];%人口数
B3=[727.96 889.86 1000.84 1265.27 1590.56 1565.71 1690.83 2166.14 3521.67 4238.79];%交通投入
p1=polyfit(B3,A,1);%数据拟合,一次
p2=polyfit(B3,A,2);%数据拟合,此处为2次
p3=polyfit(B3,A,3);%数据拟合,此处为3次
p4=polyfit(B3,A,4);%数据拟合,此处为4次
p5=polyfit(B3,A,5);%数据拟合,此处为5次
p6=polyfit(B3,A,6);%数据拟合,此处为6次
z1=poly2sym([p1(1),p1(2)]);%得出1次函数
z2=poly2sym([p2(1),p2(2),p2(3)]);%得出二次函数
z3=poly2sym([p3(1),p3(2),p3(3),p3(4)]);%得出三次函数
z4=poly2sym([p4(1),p4(2),p4(3),p4(4),p4(5)]);%得出四次函数
z5=poly2sym([p5(1),p5(2),p5(3),p5(4),p5(5),p5(6)]);%得出五次函数
z6=poly2sym([p6(1),p6(2),p6(3),p6(4),p6(5),p6(6),p6(7)]);%得出六次函数
y1=polyval(p1,B3);%将1次拟合参数带入求得函数值
y2=polyval(p2,B3);%将二次拟合参数带入求得函数值
y3=polyval(p3,B3);%将三次拟合参数带入求得函数值
y4=polyval(p4,B3);%将四次拟合参数带入求得函数值
y5=polyval(p5,B3);%将五次拟合参数带入求得函数值
y6=polyval(p6,B3);%将六次拟合参数带入求得函数值
figure(1)
plot(B3,A,'r',B3,y2,'b')
title('二次拟合与原曲线比较')
xlabel('交通投入')
ylabel('人口数')
hold on;
figure(2)
plot(B3,A,'r',B3,y3,'b')
title('三次拟合与原曲线比较')
xlabel('交通投入')
ylabel('人口数')
hold on;
figure(3)
plot(B3,A,'r',B3,y4,'b')
title('四次拟合与原曲线比较')
xlabel('交通投入')
ylabel('人口数')
hold on;
figure(4)
plot(B3,A,'r',B3,y5,'b')
title('五次拟合与原曲线比较')
xlabel('交通投入')
ylabel('人口数')
hold on;
figure(5)
plot(B3,A,'r',B3,y6,'b')
title('六次拟合与原曲线比较')
xlabel('交通投入')
ylabel('人口数')
hold on;
figure(6)
plot(B3,A,'r',B3,y1,'b')
title('一次拟合与原曲线比较')
xlabel('交通投入')
ylabel('人口数')
hold on