load dat2.mat
dat2_max=max(dat2(:,2));
dat2_min=min(dat2(:,2));
dat21=(dat2(:,2)-dat2_min)/(dat2_max-dat2_min);%归一化处理
plot(dat21)
axis([0 1200 0 1.5])
P2=[7.2714e-010 -2.583e-006 0.0029667 -0.111];
yss2=polyval(P2,0);
for k=1:1200
ys2= polyval(P2,k);
yss2=[yss2;ys2];
end