以函数y=(x+1)的平方为例 (未知)。假设测得5组数据,x=[-2,-1,0,1,2];y=[1,0,1,4,9]; polyfit(x,y,2);便可得到 ans=1.000 2.000 1.000 由此可知y=1*x平方 + 2*x + 1 ,
再根据这个函数关系,定义y=x平方+2*x+1;x=(-2:0.1:2);plot(x,y);便可得到所测数据拟和为函数后的图象。
x=0:0.1:1;
y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; A=polyfit(x,y,2); z=polyval(A,x);
plot(x,y,'r*',x,z,'b') x=VarName5'; y=VarName6'; A=polyfit(x,y,3); a=0:0.01:100; z=polyval(A,a);
plot(x,y,'r*',a,z,'b')(%小二乘法拟合数据曲线) a’为a的转置,’为英文输入法下单引号, import numpy as np import pylab as pl
# Use numpy to load the data contained in the file # ’fakedata.txt’ into a 2-D array called data data = np.loadtxt(’fakedata.txt’)
# plot the first column as x, and second column as y pl.plot(data[:,0], data[:,1], ’ro’)
pl.xlabel(’x’) pl.ylabel(’y’) pl.xlim(0.0, 10.) pl.show()
因篇幅问题不能全部显示,请点此查看更多更全内容