cernroot对数据进行拟合例子

假设我们有一组数据点(x, y)，我们想要拟合一个函数f(x)来描述这些数据。使用CERN根软件包中的TF1类，我们可以轻松地进行拟合。

首先，我们需要定义一个函数来拟合我们的数据。这个函数应该有一些参数，我们可以通过调整这些参数来优化我们的拟合。例如，我们可以定义一个简单的线性函数：

Double_t linear(Double_t *x, Double_t *par)
{
  return par[0] + par[1]*x[0];
}

这个函数的参数par[0]和par[1]是我们想要拟合的参数，x[0]是我们输入的变量x。现在，我们可以使用TF1类来创建一个实例，并将我们的函数和参数传递给它：

TF1 *f1 = new TF1("f1", linear, 0, 10, 2);
f1->SetParameters(0, 1);

这将创建一个名为"f1"的函数，它将使用我们定义的linear函数来进行拟合。我们还设置了两个参数，初始值为0和1。

现在，我们可以将我们的数据点添加到图形中，并使用TF1类中的Fit方法来拟合我们的函数：

TGraph *gr = new TGraph(n, x, y);
gr->Fit("f1");

这将使用最小二乘法来拟合我们的数据，并将拟合结果存储在我们的函数f1中。

最后，我们可以绘制我们的数据和拟合结果：

gr->Draw("AP");
f1->Draw("same");

这将绘制我们的数据点和拟合曲线。完整的示例代码如下：

#include <iostream>
#include "TCanvas.h"
#include "TF1.h"
#include "TGraph.h"

using namespace std;

Double_t linear(Double_t *x, Double_t *par)
{
  return par[0] + par[1]*x[0];
}

int main()
{
  const Int_t n = 10;
  Double_t x[n] = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9};
  Double_t y[n] = {1, 3, 5, 7, 9, 11, 13, 15, 17, 19};

  TCanvas *c1 = new TCanvas("c1", "c1", 800, 600);

  TF1 *f1 = new TF1("f1", linear, 0, 10, 2);
  f1->SetParameters(0, 1);

  TGraph *gr = new TGraph(n, x, y);
  gr->Fit("f1");

  gr->Draw("AP");
  f1->Draw("same");

  c1->Print("linear_fit.png");

  return 0;
}

这将生成一个名为"linear_fit.png"的图像文件，显示我们的数据点和拟合曲线。