clear all;
close all;
a=[-0.2142,-0.3087,-0.3938,-0.8153,-0.9323,-1.5011];
N=[13945,19395,26250,33967,36914,51328];
plot(a,N,'*','markersize',3)
grid on %画图并观察离散数据的特性
p=polyfit(a,N,1); %用1次多项式进行拟合
f = polyval(p,a);
hold on
plot(a,f,'r');
xlabel('a')
ylabel('N')
title('数据拟合图')
p %显示拟合多项式系数
p =
1.0e+004 *
-2.7068 1.1508
得到拟合直线的方程为:
y=-27068*x+11508