curvefit07.m

数值分析与计算的PPT幻灯片与MAT源码

%curvefit07.m

x0=[19 25 31 38 44];
y0=[19.0 32.3 49.0 73.3 97.8];
x1=x0.^2;
x2=[ones(5,1),x1'];
ab=x2\y0'
a=ab(1)
b=ab(2)

x=[19.0:0.2:44];
y=a+b*x.^2;
plot(x0,y0,'.r',x,y,'-b','markersize',30)
xlabel('x')
ylabel('y')
legend('original data','fitted curve',2)
title('The fitted curve of polynomial')