f# - 低通滤波器在 MathNet.Filter、F# 中不起作用

我正在尝试使用 MathNet.Filter 在 F# 中使用低通滤波器，当我绘制结果时，结果不起作用。
我在C#中发现了一些类似的问题，但是当我在F#中重现代码时结果是错误的。

这是我使用 MathNet.Filter 的方式:

let fs ,fcut, order = 5000., 2.5, 5  // set the filter parameters
let lowPass2 = MathNet.Filtering.OnlineFilter.CreateLowpass(MathNet.Filtering.ImpulseResponse.Finite,fs,fcut,order)
let disp_beam_f = disp_beam |> lowPass2.ProcessSamples // apply the filter

// I print max/min to do a quick check
(disp_beam_f |> Array.min, disp_beam |> Array.min )

然后，当我绘制一些数据时，我的结果是:

当我只绘制过滤后的数据时，该图看起来像原始数据，但具有不同的比例:

我已经尝试过改变过滤器的顺序，并没有使用过滤器的顺序，结果更好，但它不是过滤后的数据。

我做错了什么？我知道 fs 是频率采样(以赫兹为单位)， fcut 是高频被切断的频率(以赫兹为单位)， 顺序 表示 fcut 以上频率的信号衰减的程度(或速度)。

有关更多信息，这里有一个脚本:
https://gist.github.com/josesoyo/f45b59a035c3e2ec2cdf0fe1f95cc84b
数据样本在:https://drive.google.com/open?id=1rA4B21i-S5g2wrB4eK7xk3_A6Rs0tLHf
结果应该是这个文件:https://drive.google.com/open?id=1ukmWtzFm_I4ou3zuQcLWCrFIOJs2eV41

编辑基于 Gene Belitski answer 我相信它澄清了这个案例:

运行下一个脚本可以更好地显示我遇到的问题:

#I __SOURCE_DIRECTORY__
#r @"..\packages\MathNet.Numerics.3.8.0\lib\net40\MathNet.Numerics.dll"
#r @"..\packages\MathNet.Filtering.0.4.0\lib\net40\MathNet.Filtering.dll"
#load @"..\packages\FSharp.Charting.0.91.1\lib\net45\FSharp.Charting.fsx"

let length, samplingRate,  amplitude = 750, 44100.,20.
let flow, fhigh = 200., 15000.
let lo = Generate.Sinusoidal(length,samplingRate,flow,amplitude)
let hi = Generate.Sinusoidal(length,samplingRate,fhigh,amplitude/4.)
let sumLoHi = lo |> Array.map2 (+) hi
let fs ,fcut, order = float(samplingRate), 5000., 10

// first low pass filter
let lowPass2 = MathNet.Filtering.OnlineFilter.CreateLowpass(MathNet.Filtering.ImpulseResponse.Finite,fs,fcut,order)
let filtered = sumLoHi |> lowPass2.ProcessSamples

// second low pass filter
let lowPass3 = MathNet.Filtering.OnlineFilter.CreateLowpass(MathNet.Filtering.ImpulseResponse.Finite,fs,fcut/5.,order)
let filtered_2 = sumLoHi |> lowPass3.ProcessSamples

// plot
Chart.Combine([Chart.Line(lo, Name="lo", Color=Color.DarkBlue);
               Chart.Line(sumLoHi, Name="sumLoHi", Color=Color.LightSkyBlue);
               Chart.Line(filtered, Name="filtered ratio fcut/SamplingRatio = 0.1", Color=Color.Purple);
               Chart.Line(filtered_2, Name="filtered ratio fcut/SamplingRatio = 0.02", Color=Color.Fuchsia)])
|> Chart.WithLegend(Title="Combined",Docking=ChartTypes.Docking.Bottom)
|> Chart.WithTitle("OnlineFilter")

结果是:


我已经指定了samplingRate和过滤频率之间的比率，因为我认为问题问题与这个参数有关。是我不应该以这种方式使用 OnlineFilter 还是其他什么？

我也尝试使用 IIR 和 FIR 滤波器，结果是一样的。相反，python/scipy 的类似情况效果很好。

最佳答案

在不了解样本数据的情况下，很难确定观察到的行为是对还是错。如果过滤器是否有效，它可能有助于快速收集样本；然后检查你的代码。

例如，我们可以取一个较低频率的波数据，将另一个较高频率的波加入其中，然后将滤波器应用于过滤掉较高频率的总和，并期望过滤后的数据将类似于较低的波。

下面的脚本实现了上述场景:

#I __SOURCE_DIRECTORY__
#r @"..\packages\MathNet.Numerics.3.8.0\lib\net40\MathNet.Numerics.dll"
#r @"..\packages\MathNet.Filtering.0.4.0\lib\net40\MathNet.Filtering.dll"
#load @"..\packages\FSharp.Charting.0.91.1\lib\net45\FSharp.Charting.fsx"

open FSharp.Charting
open MathNet.Numerics
open MathNet.Filtering
open System.Drawing

let length, samplingRate,  amplitude = 250, 44100.,20.
let lo = Generate.Sinusoidal(length,samplingRate,1000.,amplitude)
let hi = Generate.Sinusoidal(length,samplingRate,15000.,amplitude/3.)
let sumLoHi = lo |> Array.map2 (+) hi
let fs ,fcut, order = 44100., 5000., 10
let lowPass2 = MathNet.Filtering.OnlineFilter.CreateLowpass(MathNet.Filtering.ImpulseResponse.Finite,fs,fcut,order)
let filtered = sumLoHi |> lowPass2.ProcessSamples
Chart.Combine([Chart.Line(lo, Name="lo", Color=Color.DarkBlue);
               Chart.Line(sumLoHi, Name="sumLoHi", Color=Color.LightSkyBlue);
               Chart.Line(filtered, Name="filtered", Color=Color.Purple)])
|> Chart.WithLegend(Title="Combined",Docking=ChartTypes.Docking.Bottom)

sine 为 1kHz，lo 为 15kHz，采样率、数据长度和频率相同的两个 hi 波组合成 sumLoHi 。然后将滤波器应用于后者，切断 5kHz 以上的频率，获取 filtered 数据。得到的 lo 、 sumLoHi 和 filtered 组合图表显示过滤器工作正常:

关于f# - 低通滤波器在 MathNet.Filter、F# 中不起作用，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/48888244/