関数型プログラミング言語 F# (fsharp)の数値計算ライブラリ Math.NET Numericsを使う(2) 高速フーリエ変換（FFT）を実行（NumPyの例題）してPlotly.NETで表示する。

さて今回はFFT。こちらの例題をF#でやってみる。

Visual C# (C_sharp)の数学ライブラリ Math.NET Numericsを使う(3) 高速フーリエ変換（FFT）を実行する。FourierOptionsにMatlabとNumerical Recipesがあるのが意外。 

まずはNuGetで

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

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

をしたのち、

open MathNet.Numerics

open MathNet.Numerics.IntegralTransforms

open Plotly.NET

open System.Numerics

let t = [| for i in 0..255 -> (float)i|]

let y = t |> Array.map(fun t -> Complex(Math.Sin(t), 0.0))

Fourier.Forward(y, FourierOptions.Matlab)

let yreal = y |> Array.map(fun y -> y.Real)

let yimag = y |> Array.map(fun y -> y.Imaginary)

let freq = Fourier.FrequencyScale(256, 1.0)

[Chart.Line(freq, yreal, ShowLegend = true, Name = "Real");

Chart.Line(freq, yimag, ShowLegend = true, Name = "Imaginary")]

|> Chart.combine |> Chart.withXAxisStyle("Frequency") |> Chart.withTitle("FFT")

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

とすればいいだけ。めちゃくちゃ簡単。関数型プログラミング言語面白くなってきた（初歩の機能しかつかってないけど…）。

結果はNumPyとも一緒。