一元回归.txt

来自「一元回归预测法在很多地方被采用,本代码主要对其参数的确定」· 文本 代码 · 共 44 行

#include<stdio.h>
#include<math.h>

main()
{
	int i,n;
	char ch;
	double x[30],y[30],xy=0,x2=0,y2=0,x_=0,y_=0,r,a,b;
	printf("n:\n");
	scanf("%d",&n);
	printf("x:\n");
	for(i=0;i<n;i++)
	{
		scanf("%lf",&x[i]);
	}
	printf("y:\n");
	for(i=0;i<n;i++)
	{
		scanf("%lf",&y[i]);
	}
	
	for(i=0;i<n;i++)
	{
		xy=xy+x[i]*y[i];
		x_=x_+x[i]/n;
		y_=y_+y[i]/n;
		x2=x2+x[i]*x[i];
		y2=y2+y[i]*y[i];
		
	}
	
	b=(xy-n*x_*y_)/(x2-n*x_*x_);
	a=y_-b*x_;
	r=(n*xy-n*x_*n*y_)/sqrt((n*x2-n*x_*x_*n)*(n*y2-n*y_*n*y_));
	printf("a=%lf,b=%lf,r=%lf,xy=%lf,x_=%lf\ny_=%lf,x2=%lf,y2=%lf\n",a,b,r,xy,x_,y_,x2,y2);
	for(i=0;i<n;i++)
	{
		y[i]=b*x[i]+a;
		printf("%lf  ",y[i]);
	}
	system("pause");
	
}