[英]Calculating a derivative with Math.NET C#
我正在寻找一个简单的函数,它将接受一个 double 值数组并返回数学导数。
Math.NET 似乎有这样的函数,但它要求以下语法:
double FirstDerivative(Func<double, double> f, double x)
我不确定为什么我需要指定一个函数。 我只想要一个可以将数据传递给的预先存在的函数。
获取您的数据点并创建一个 Math.NET Numerics Cubic Spline对象。 然后使用.Differentiate()
方法获取您想要的每个点的斜率。
试试下面的代码:
static class Program
{
const int column_width = 12;
/// <summary>
/// The main entry point for the application.
/// </summary>
[STAThread]
static void Main(string[] args)
{
var xvec = new DenseVector(new double[] { 0.0, 1.0, 2.0, 3.0, 4.0 });
var yvec = new DenseVector(new double[] { 3.0, 2.7, 2.3, 1.6, 0.2 });
Debug.WriteLine("Input Data Table");
Debug.WriteLine($"{"x",column_width} {"y",column_width}");
for(int i = 0; i < xvec.Count; i++)
{
Debug.WriteLine($"{xvec[i],column_width:G5} {yvec[i],column_width:G5}");
}
Debug.WriteLine(" ");
var cs = CubicSpline.InterpolateNatural(xvec, yvec);
var x = new DenseVector(15);
var y = new DenseVector(x.Count);
var dydx = new DenseVector(x.Count);
Debug.WriteLine("Interpoaltion Results Table");
Debug.WriteLine($"{"x",column_width} {"y",column_width} {"dy/dx",column_width}");
for(int i = 0; i < x.Count; i++)
{
x[i] = (4.0*i)/(x.Count-1);
y[i] = cs.Interpolate(x[i]);
dydx[i] = cs.Differentiate(x[i]);
Debug.WriteLine($"{x[i],column_width:G5} {y[i],column_width:G5} {dydx[i],column_width:G5}");
}
}
}
并查看调试输出:
Input Data Table
x y
0 3
1 2.7
2 2.3
3 1.6
4 0.2
Interpoaltion Results Table
x y dy/dx
0 3 -0.28214
0.28571 2.919 -0.28652
0.57143 2.8354 -0.29964
0.85714 2.7469 -0.3215
1.1429 2.6509 -0.35168
1.4286 2.5454 -0.38754
1.7143 2.429 -0.42864
2 2.3 -0.475
2.2857 2.154 -0.55809
2.5714 1.9746 -0.7094
2.8571 1.7422 -0.92894
3.1429 1.4382 -1.1979
3.4286 1.0646 -1.4034
3.7143 0.64404 -1.5267
4 0.2 -1.5679
如果您不反对 Math.Net 以外的库,则可以尝试AlgLib及其spline1ddiff函数
构建样条曲线很容易,而且 Akima 样条曲线通过点看起来很好且平滑。 如果你想要一个方法来接收一组数据并返回导数,这里有一个使用 AlgLib 数学库的例子:
public static void CalculateDerivatives(this Dictionary<double, double> inputPoints, out Dictionary<double, double> firstDerivatives, out Dictionary<double, double> secondDerivatives)
{
var inputPointsXArray = inputPoints.Keys.ToArray();
var inputPointsYArray = inputPoints.Values.ToArray();
spline1dinterpolant akimaSplineToDifferentiate;
alglib.spline1dbuildakima(inputPointsXArray, inputPointsYArray, out akimaSplineToDifferentiate);
firstDerivatives = new Dictionary<double, double>();
secondDerivatives = new Dictionary<double, double>();
foreach (var pair in inputPoints)
{
var xPoint = pair.Key;
double functionVal, firstDeriv, secondDeriv;
alglib.spline1ddiff(akimaSplineToDifferentiate, xPoint, out functionVal, out firstDeriv, out secondDeriv);
firstDerivatives.Add(point, firstDeriv);
secondDerivatives.Add(point, secondDeriv);
}
}
警告:Akima Spline 在数据集范围之外具有不可预测的行为。
感谢您的回复。
我相信我需要遍历数组,从 Math.NET 调用 Differentiate 函数,或者可能只是编写我自己的(运行时上升)计算。
在这里你有函数的衍生扩展
public static Func<double[], double> Derivative(this Func<double[], double> func, int derivativeIndex)
{
double step = 0.001;
return income =>
{
double[] increasedIncome = (double[])income.Clone();
increasedIncome[derivativeIndex] += step;
double[] decreasedIncome = (double[])income.Clone();
decreasedIncome[derivativeIndex] -= step;
return (func(increasedIncome) - func(decreasedIncome)) / (2 * step);
};
}
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.