関数型プログラミング言語 F# (fsharp)の数値計算ライブラリ Math.NET Numericsを使う(3) 補間を行う（Interpolate) リニア、３次スプライン、有理関数などいろいろ使える。Plotly.NETもF#から簡単に使える。

今回はこのC#の例をF#でやってみる。使うのはPolyglot Notebooks。

Visual C# (C_sharp)の数学ライブラリ Math.NET Numericsを使う(2) 補間を行う（Interpolate) リニア、３次スプライン、有理関数などいろいろ使える。 

グラフにもしたいのでPlotly.NETも導入する。F#から簡単に使えるのにびっくり。

https://plotly.net/

まずはMath.NET numericsとPlotly.NETをNuGetで導入するためにセルにこれを書いてSHIFT＋リターン。

#r "nuget: MathNet.Numerics.FSharp, 5.0.0"

#r "nuget: Plotly.NET.Interactive, 5.0.0"

そしてコードは関数型らしくパイプラインやmapがあるのでめちゃくちゃ簡単にかける。C#よりはるかに簡単。

open MathNet.Numerics

open Plotly.NET

let x = [| for i in 0..9 -> (float)i * 2.0 * Math.PI / 9.0 |]

let y = x |> Array.map(fun x -> Math.Sin(x))

let interpRational = Interpolate.RationalWithoutPoles(x, y)

let interpLinear = Interpolate.Linear(x, y)

let xval =  [| for i in 0..49 -> (float)i * 2.0 * Math.PI / 49.0 |]

let yinterpRational = xval |> Array.map(fun x -> interpRational.Interpolate(x))

let yinterpLinear = xval |> Array.map(fun x -> interpLinear.Interpolate(x))

[Chart.Point(x, y, ShowLegend = true, Name = "Original");

Chart.Line(xval, yinterpRational, ShowLegend = true, Name = "Rational");

Chart.Line(xval, yinterpLinear, ShowLegend = true, Name = "Linear") ]

|> Chart.combine |> Chart.withLayoutStyle(Width=800, Height=500)

|> Chart.withTitle("Interpolation")

|> Chart.withXAxisStyle("x")

|> Chart.withYAxisStyle("y")

結果はこちら。

ちゃんとC#の結果を再現している。

しかしPolyglot Notebooksを使うとPythonより簡単にいろいろな計算ができるような気がしてきた。