0

我是 C# 和一般编程的初学者，我想问一下。我正在尝试制作一个包含 datagridview 和一些计算的项目。基本上我想做的是从经验数据表中用切比雪夫近似值近似数据。使用来自数字配方的代码我想使用这个功能，chebft，它看起来像这样：

void chebft(NRREAL a, NRREAL b, NRREAL c[], int n, NRREAL (*func)(NRREAL))
{ 
    int k,j;
    NRREAL fac,bpa,bma,*f;

    f=vector(0,n-1);
    bma=NR_1_2*(b-a);
    bpa=NR_1_2*(b+a);
    for (k=0;k<n;k++) {
    NRREAL y=cos(NR_PI*(k+NR_1_2)/n);
    f[k]=(*func)(y*bma+bpa);
    }
    fac=NR_2/n;
    for (j=0;j<n;j++) {
    NRZEAL sum=NR_0;
    for (k=0;k<n;k++)
        sum += f[k]*cos(NZ_PI*j*(k+NZ_1_2)/n);
    c[j]=fac*(NRREAL)sum;
    }
    free_vector(f,0,n-1);
}

我希望它在 C# 中使用，所以到目前为止我有什么：

public void chebft(double a, double b, double[] c, int n, double func)  
    {
        int k,j;
        double fac, bpa, bma;

        double PI = 3.14159265358979323846;
        Vector f = new Vector(0, n - 1);
        bma = 0.5 * (b - a);
        bpa = 0.5 * (b - a);

        for (k = 0; k < n; k++)
        { 
            double y = Math.Cos(PI*(k+0.5)/n);
            f[k] = (func)(y*bma+bpa);
        }

        fac = 2.0 / n;
        for (j = 0; j < n; j++)
        {
            double sum = 0.0;

            for (k = 0; k < n; k++)
            {
                sum += f[k] * Math.Cos(PI * j * (k + 0.5) / n);
                c[j] = fac * (double)sum;
            }
        }

    }

我知道 func 部分和 f[k] 部分有错误，但我不知道如何修复它们。

感谢帮助

4

1 回答 1

1

翻译似乎相对可靠，但有一个例外（撇​​开可能的数字问题，这些问题需要一些更高级的问题知识才能诊断）：

函数指针func接受一个双值参数并返回一个双精度值。要在 C# 中传递这样一个 lambda 函数，它本质上等同于 C++ 中的函数指针，请使用以下形式：

Func<double, double> func

然后，您可以按如下方式调用该函数：

f[k] = func(y*bma+bpa);

更多信息可以在这里找到。

另请注意，您不需要sum在最终循环中强制转换为 double，因为它已经是 double，并且Math该类已经定义了 constant Math.PI。

于 2017-06-08T17:48:32.887 回答