c语言 线性回归

来源:百度知道 编辑:UC知道 时间:2024/06/25 11:37:52
通过测试实验获得以下数据,要求利用所学的线性回归问题方法求解该直线方程。
X 0 1 2 3 4 5
Y 1.1 1.9 3.1 3.9 4.4 4.9
答案Y=1.266666x+0.7800002

根据:y=a+bx 可知

∑y=na+b∑x

∑x∑y=na∑x+b(∑x)^2 (1)

xy=ax+bx^2

∑xy=a∑x+b∑x^2

n∑xy=na∑x+nb∑x^2 (2)

联立(1)(2)即可求出a和b的表达式

main()
{
float x[6]={0.0,1.0,2.0,3.0,4.0,5.0};
float y[6]={1.1,1.9,3.1,3.9,4.4,4.9};
float a,b,mxy,xx,yy,x2,x22;
int i;
a=b=mxy=xx=yy=x2=x22=0.0;
for(i=0;i<6;i++)
{
mxy=6.0*x[i]*y[i]+mxy;
xx=1.0*x[i]+xx;
yy=1.0*y[i]+yy;
x2=1.0*x[i]*x[i]*6.0+x2;
x22=1.0*x[i]+x22;
}
b=1.0*(mxy-xx*yy)/(x2-x22*x22);
a=1.0*yy/6.0-b*xx/6.0;
printf("Y=%fx+%f\n",a,b);
}

这个是不是得用matlab啊