用MATLAB建立一个M文件。
clear
clc
time=[0:0.5:4 5 6 6.1 7]
temperature=[300 281 261 244 228 214 202 191 181 164 151 149 141]
p=polyfit(time,temperature,3);
x=[0:0.5:4 5 6 6.1 7];
a=p(1);b=p(2);c=p(3);d=p(4);
y=a*x.^3+b*x.^2+c*x+d;
plot(time,temperature,'<',x,y,'g')
w=a*x.^3+b*x.^2+c*x+d
N=13;
MEAN=sum(temperature)/N
mean(temperature)
S=sum((temperature-MEAN).^2)
A=sum((w-temperature).^2)
r2=1-A/S
r2 =
0.9999
结果r2非常接近1,因此是很完美的拟合。