python代码:
import numpy as nparrX=np.random.randint(0,100,30)#print(arrX)arrY=arrX*5+3#print(arrY)for i in range(30):rand=2*np.random.random()-1arrY[i]+=randsum=0for i in arrX:sum+=iexp_x=sum/30;#x的平均值sum=0for i in arrY:sum+=iexp_y=sum/30#y的平均值sum=0for i in range(0,len(arrX)):sum+=arrX[i]*arrY[i]exp_xy=sum/30;#xy的平均值sum=0for i in range(0,len(arrX)):sum+=arrX[i]*arrX[i]exp_xx=sum/30;#xx的平均值a=(exp_xy-exp_x*exp_y)/(exp_xx-exp_x*exp_x)b=exp_y-a*exp_xprint(a,b)