最小二乘法 方程组.cpp

来自「大家请多指教！！！挺不错的哦！」· C++ 代码 · 共 55 行

#include<iostream.h>
#include<math.h>
#include<stdio.h>
void main()
{   int n,i;
    double xx=0,xy=0,yy=0,xz=0,yz=0,x1,x2;
	double x[100],y[100],z[100];

	cout<<"输入数据组数：";
	cin>>n;
	for( i=0;i<n;i++)
    {   cout<<"x["<<i<<"]=";
		cin>>x[i];
	}
	for( i=0;i<n;i++)
	{	cout<<"y["<<i<<"]=";
		cin>>y[i];
    }
	for( i=0;i<n;i++)
	{	cout<<"z["<<i<<"]=";
		cin>>z[i];
	}
	  

	for( i=0;i<n;i++)
	{   xx=xx+x[i]*x[i];
	    xy=xy+x[i]*y[i];
		yy=yy+y[i]*y[i];
		xz=xz+x[i]*z[i];
		yz=yz+y[i]*z[i];
	}


x1=(xz*yy-xy*yz)/(xx*yy-xy*xy);
x2=(xx*yz-xz*xy)/(xx*yy-xy*xy);

cout<<"x1= "<<x1<<endl;
cout<<"x2= "<<x2;
}