C++实现最小二乘法对y=exp(x)的二次拟合
生活随笔
收集整理的這篇文章主要介紹了
C++实现最小二乘法对y=exp(x)的二次拟合
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
#include<iostream>
#include<math.h>
const int N=20;
const int n=3;
double A[n][n];double B[n];double start=0, end=1.0;//積分范圍
const double PI=3.1415926;
using namespace std;
void HibertMatrix();
double func1(double x);
double func2(double x);
void projectAarry();
void jacobi();
int main()
{double A[n][n];double B[n];double start=0, end=1.0;//積分范圍HibertMatrix(); //生成希爾伯特矩陣projectAarry();//生成投影矩陣jacobi();//求解正規方程
return 0;
}
void HibertMatrix()//希爾伯特矩陣的生成
{for(int i=0;i<n;i++)for(int j=0;j<n;j++)A[i][j]=1.0/(i+j+1);
}
double func1(double x) //Expansion function
{ return exp(x);
}
double func2(double x) //Base function
{return x;}
void projectAarry()
{double h=(end-start)/N;for(int i=0;i<n;i++){B[i]=(pow(func2(start),i)*func1(start)+pow(func2(end),i)*func1(end))*h/2;for(int k=1;k<N;k++) B[i]+=pow(func2(start+k*h),i)*func1(start+k*h)*h;cout<<B[i]<<'\t';
}
cout<<endl;
}
void jacobi()
{// double A[N][N]={{10,-1.0,-2.0},{-1.0,10.0,-2.0},{-1.0,-1.0,5.0}}; // double B[N]={72,83,42}; double X[N]={0,0,0,};
for (int k=0;k<100;k++)
{
for(int i=0;i<n;i++)
{
double sum=0;
for(int j=0;j<n;j++)
{ if(j==i) continue; sum+=A[i][j]*X[j];
}
X[i]=(B[i]-sum)/A[i][i];
}
}
for (int i=0;i<n;i++) cout<<X[i]<<'\t'; cout<<endl;
}
總結
以上是生活随笔為你收集整理的C++实现最小二乘法对y=exp(x)的二次拟合的全部內容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： C++理解全局变量在调用函数中的应用
- 下一篇： C++ 类的保护继承与构造函数