先输入x=[1 1.5 .......];
y=[-1.4 2.7 ......];
p=polyfit(x,y,2);%二次多项式
p=polyfit(x,y,3);%三次多项式
......
p=polyfit(x,y,10);%十次多项式
其中p是一个行向量,p(1)储存最高次幂的多项式系数,知道最后一个储存常数项
例如 p=polyfit(x,y,2);就返回长度为3的p,代表多项式
p(x)=p(1)*x^2+p(2)*x+p(3);
用函数polyval(p,x)就可以计算多项式的值
输入时间t序列
输入温服度T序列
interp1(t,T,[10 12],'spline')