C#，数值计算——数据建模Fitab的计算方法与源程序

news2025/1/10 16:25:30

1 文本格式

using System;

namespace Legalsoft.Truffer
{
/// <summary>
/// Fitting Data to a Straight Line
/// </summary>
public class Fitab
{
private int ndata { get; set; }
private double a { get; set; }
public double b { get; set; }
private double siga { get; set; }
private double sigb { get; set; }
private double chi2 { get; set; }
private double q { get; set; }
private double sigdat { get; set; }
private double[] x { get; set; }
private double[] y { get; set; }
private double[] sig { get; set; }

public Fitab(double[] xx, double[] yy, double[] ssig)
{
this.ndata = xx.Length;
this.x = xx;
this.y = yy;
this.sig = ssig;
this.chi2 = 0.0;
this.q = 1.0;
this.sigdat = 0.0;
Gamma gam = new Gamma();

double ss = 0.0;
double sx = 0.0;
double sy = 0.0;
double st2 = 0.0;
b = 0.0;
for (int i = 0; i < ndata; i++)
{
double wt = 1.0 / Globals.SQR(sig[i]);
ss += wt;
sx += x[i] * wt;
sy += y[i] * wt;
}
double sxoss = sx / ss;
for (int i = 0; i < ndata; i++)
{
double t = (x[i] - sxoss) / sig[i];
st2 += t * t;
b += t * y[i] / sig[i];
}
b /= st2;
a = (sy - sx * b) / ss;
siga = Math.Sqrt((1.0 + sx * sx / (ss * st2)) / ss);
sigb = Math.Sqrt(1.0 / st2);
for (int i = 0; i < ndata; i++)
{
chi2 += Globals.SQR((y[i] - a - b * x[i]) / sig[i]);
}
if (ndata > 2)
{
q = gam.gammq(0.5 * (ndata - 2), 0.5 * chi2);
}
}

public Fitab(double[] xx, double[] yy)
{
this.ndata = xx.Length;
this.x = xx;
this.y = yy;
this.sig = xx;
this.chi2 = 0.0;
this.q = 1.0;
this.sigdat = 0.0;

double sx = 0.0;
double sy = 0.0;
double st2 = 0.0;
b = 0.0;
for (int i = 0; i < ndata; i++)
{
sx += x[i];
sy += y[i];
}
double ss = ndata;
double sxoss = sx / ss;
for (int i = 0; i < ndata; i++)
{
double t = x[i] - sxoss;
st2 += t * t;
b += t * y[i];
}
b /= st2;
a = (sy - sx * b) / ss;
siga = Math.Sqrt((1.0 + sx * sx / (ss * st2)) / ss);
sigb = Math.Sqrt(1.0 / st2);
for (int i = 0; i < ndata; i++)
{
chi2 += Globals.SQR(y[i] - a - b * x[i]);
}
if (ndata > 2)
{
sigdat = Math.Sqrt(chi2 / (ndata - 2));
}
siga *= sigdat;
sigb *= sigdat;
}
}
}

2 代码格式

using System;

namespace Legalsoft.Truffer
{
    /// <summary>
    /// Fitting Data to a Straight Line
    /// </summary>
    public class Fitab
    {
        private int ndata { get; set; }
        private double a { get; set; }
        public double b { get; set; }
        private double siga { get; set; }
        private double sigb { get; set; }
        private double chi2 { get; set; }
        private double q { get; set; }
        private double sigdat { get; set; }
        private double[] x { get; set; }
        private double[] y { get; set; }
        private double[] sig { get; set; }

        public Fitab(double[] xx, double[] yy, double[] ssig)
        {
            this.ndata = xx.Length;
            this.x = xx;
            this.y = yy;
            this.sig = ssig;
            this.chi2 = 0.0;
            this.q = 1.0;
            this.sigdat = 0.0;
            Gamma gam = new Gamma();

            double ss = 0.0;
            double sx = 0.0;
            double sy = 0.0;
            double st2 = 0.0;
            b = 0.0;
            for (int i = 0; i < ndata; i++)
            {
                double wt = 1.0 / Globals.SQR(sig[i]);
                ss += wt;
                sx += x[i] * wt;
                sy += y[i] * wt;
            }
            double sxoss = sx / ss;
            for (int i = 0; i < ndata; i++)
            {
                double t = (x[i] - sxoss) / sig[i];
                st2 += t * t;
                b += t * y[i] / sig[i];
            }
            b /= st2;
            a = (sy - sx * b) / ss;
            siga = Math.Sqrt((1.0 + sx * sx / (ss * st2)) / ss);
            sigb = Math.Sqrt(1.0 / st2);
            for (int i = 0; i < ndata; i++)
            {
                chi2 += Globals.SQR((y[i] - a - b * x[i]) / sig[i]);
            }
            if (ndata > 2)
            {
                q = gam.gammq(0.5 * (ndata - 2), 0.5 * chi2);
            }
        }

        public Fitab(double[] xx, double[] yy)
        {
            this.ndata = xx.Length;
            this.x = xx;
            this.y = yy;
            this.sig = xx;
            this.chi2 = 0.0;
            this.q = 1.0;
            this.sigdat = 0.0;

            double sx = 0.0;
            double sy = 0.0;
            double st2 = 0.0;
            b = 0.0;
            for (int i = 0; i < ndata; i++)
            {
                sx += x[i];
                sy += y[i];
            }
            double ss = ndata;
            double sxoss = sx / ss;
            for (int i = 0; i < ndata; i++)
            {
                double t = x[i] - sxoss;
                st2 += t * t;
                b += t * y[i];
            }
            b /= st2;
            a = (sy - sx * b) / ss;
            siga = Math.Sqrt((1.0 + sx * sx / (ss * st2)) / ss);
            sigb = Math.Sqrt(1.0 / st2);
            for (int i = 0; i < ndata; i++)
            {
                chi2 += Globals.SQR(y[i] - a - b * x[i]);
            }
            if (ndata > 2)
            {
                sigdat = Math.Sqrt(chi2 / (ndata - 2));
            }
            siga *= sigdat;
            sigb *= sigdat;
        }
    }
}

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.coloradmin.cn/o/1068196.html

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！