sample6_12.cs

来自「C#的数值计算的几乎所有算法」· CS 代码 · 共 34 行

/*
 * 示例程序Sample6_12: Interpolation类的光滑等距插值
 */

using System;
using CSharpAlgorithm.Algorithm;

namespace CSharpAlgorithm.Sample
{
	class Class1
	{
		[STAThread]
		static void Main(string[] args)
		{
            // 数据
            int n = 10;
            double x0 = -1;
            double xstep = 0.2;
            double[] y = {0.0384615,0.0588236,0.1,0.2,0.5,1.0,0.5,0.2,0.1,0.0588236,0.0384615};
            double[] t = {-0.65,0.25};
		
            // 插值运算
            for (int i=0; i<t.Length; ++i)
            {
                double[] coef = new double[5];
                double yt = Interpolation.GetValueAkima(n, x0, xstep, y, t[i], coef, -1);
                Console.WriteLine("f(" + t[i] + ") = " + yt);
                for (int j=0; j<4; ++j)
                    Console.WriteLine("s" + j + " = " + coef[j]);
            }		
        }
	}
}