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)的二次拟合相关推荐

MATLAB | 绘制复指数函数 y = exp(j*w*n)的三维图像
MATLAB | 绘制复指数函数 y = exp(j×w×n)的三维图像 1 思路 2 代码 3 结果 4 说明概要: 本文内容是大二下学期<信号与系统>课上做的一道小题,目的是让我们直 ...
python输入y继续运行_Python二三事 - 接触Python(x,y)
注释:本文原来转自博客大巴的一篇文章,向原作者致谢!我也做了相应的修改,更加完善对于Pythonx,y的介绍 Python 二三事面向初学者介绍Python相关的一些工具,以及可能遇到的常见问题. ...
最小二乘法曲线拟合y=c1+c2x+c3x^2——java版
import java.util.Scanner; public class Nllsf {/*** 拟合目标P(x)=c1+c2x+c3x^2* c1+c2x+c3x^2=f(x)*/static ...
opencv图像处理学习（五十二）——拟合（最小二乘法）
怀着沉痛的心情,拖着疲惫的身心,为了拟合好圆,我实在不得不上最小二乘法了(我上班写的代码不要想了,不可能发在blog里的),现在进入正题. (1)基本原理名称自变量:x - ...
基于最小二乘法与matlab软件的动态称重数据拟合,结合面参数最小相对误差二乘法拟合...
论文写作指导:请加QQ229366758 摘要:结合文献资料,利用最小相对误差二乘法对原模型进行求解,得到新的回归方程组.结果证明:用最小相对误差二乘法解出的回归方程组相对误差在3.5%以内,相对于 ...
计算机对化工的影响PPT,化工计算机应用ppt课件.ppt
化工计算机应用 1 教材及主要参考书方利国陈砺计算机在化学化工中的应用化学工业出版社 2006AspenPlus10 0用户指南各种网络资源 2 教学安排第1章拟合方法 2学时第2章非线性 ...
matlab多项式拟合要求系数项大于零,matlab多项式系数
要求一高阶多项式的根往往须借助数值方法,所幸MATLAB已将这些数值方法写成一函数 roots(p),我们只要输入多项式的各阶系数 (以 p 代表)即可求解到对应的根 >...... 2. ...
matlab拟合直线交点,多条直线如何求交点
clc clear all close all % imread('4001.bmp') x=[2.7058976e+002 9.1422067e+002 9.8335400e+002]; y=[1. ...
Matlab 多项式拟合
文章目录 1. 按 1.1. 语法 2. 一次多项式 3. 二次多项式 1. 按 polyfit函数用于多项式的曲线拟合. 1.1. 语法 p = polyfit(x,y,n) [p,S] = pol ...

C++实现最小二乘法对y=exp(x)的二次拟合

C++实现最小二乘法对y=exp(x)的二次拟合相关推荐

最新文章

热门文章