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
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啊