%% 多项式曲线拟合
%自变量向量
x=[1 2 3 4 5];
%应变量向量
y=[5.6 40 150 250 498.9];
x=[1 2 3 4 5];
y=[5.6 40 150 250 498.9];
p=polyfit(x,y,1)
x=[1 2 3 4 5];
y=[5.6 40 150 250 498.9];
p=polyfit(x,y,1)
%计算多项式的值 (polyvalm 计算矩阵多项式)
y2=polyval(p,x2);
x=[1 2 3 4 5];
y=[5.6 40 150 250 498.9];
% n=3 表示使用三次多项式进行拟合
% 三次拟合方程为 a*x^3+bx^2+cx+d
p=polyfit(x,y,3)
x2=1:0.1:5;
y2=polyval(p,x2);
plot(x,y,'*',x2,y2);
help CFTOOL
%年份从 1900 到 2000,间隔为 10
year=1900:10:2000;
%人口数量
number=100*sort(random('logn',0,1,1,length(year)));
%知道了 1900,1910,……,2000 年,每个 10 年的人口数量
%通过插值方法获取 1901 或 1999 年人口的数据
x=1900:1:2000;
%采用样条差值方法
y=interp1(year,number,x,'spline');