leastsquares.cpp

数值分析的简单算法VC++实现

// Leastsquares.cpp: implementation of the Leastsquares class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "Leastsquares.h"

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

double Leastsquares::cal(double DATA[100][2],double var[100],double x,double t)
{double a11=0,a12=0,a21=0,a22=0,result;
 int i,j;
 double y1=0,y2=0;
 double str[100][100];
 for(i=0;i<100;i++)
	 for(j=0;j<100;j++)
	 {str[i][j]=0;}
 for(i=0;DATA[i][0]>-999999&&DATA[i][0]<999999;i++)
 {str[0][0]=str[0][0]+var[i];
  str[1][0]=str[1][0]+var[i]*DATA[i][0];
  str[0][1]=str[1][0];
  str[1][1]=str[1][1]+var[i]*DATA[i][0]*DATA[i][0];
  y1=y1+var[i]*DATA[i][1];
  y2=y2+var[i]*DATA[i][1]*DATA[i][1];
 }
 a=(y1*str[1][1]-str[1][0]*y2)/(str[0][0]*str[1][1]-str[1][0]*str[0][1]);
 b=(y1-str[0][0]*a)/str[1][0];
 result=a+b*t;
 return(a+b*t);
}

/////////////////////////////////////////////////////////////////////
/*

double Leastsquares::line(double var)
{double result;
 result=a+var*b;
// result=1+1*var;  
return result;
}*/