Logiciel de filtrage passe-bas?
-
19-08-2019 - |
Question
Je recherche un code / bibliothèque / classe de filtre passe-bas numérique pour un projet de formulaires Windows .net, de préférence écrit en c, c ++ ou c #. Je dois probablement définir le nombre de pôles, de coefficients, de fenêtres, etc. Je ne peux utiliser aucun des codes gpl disponibles, et je ne sais pas ce qu'il y a d'autre. Toutes les suggestions appréciées.
La solution
Voici un filtre Butterworth Low Pass que j'ai écrit pour un projet récent.
Il y a des nombres magiques comme constantes qui m’ont été donnés. Si vous savez comment créer les nombres magiques avec vos pôles, vos coefficients, etc., cela pourrait être utile.
using System;
using System.Collections.Generic;
using System.Text;
namespace Filter
{
public class ButterworthLowPassFilter
{
//filter fc = 2hz, fs = 10hz
private const int LowPassOrder = 4;
private double[] inputValueModifier;
private double[] outputValueModifier;
private double[] inputValue;
private double[] outputValue;
private int valuePosition;
public ButterworthLowPassFilter()
{
inputValueModifier = new double[LowPassOrder];
inputValueModifier[0] = 0.098531160923927;
inputValueModifier[1] = 0.295593482771781;
inputValueModifier[2] = 0.295593482771781;
inputValueModifier[3] = 0.098531160923927;
outputValueModifier = new double[LowPassOrder];
outputValueModifier[0] = 1.0;
outputValueModifier[1] = -0.577240524806303;
outputValueModifier[2] = 0.421787048689562;
outputValueModifier[3] = -0.0562972364918427;
}
public double Filter(double inputValue)
{
if (this.inputValue == null && this.outputValue == null)
{
this.inputValue = new double[LowPassOrder];
this.outputValue = new double[LowPassOrder];
valuePosition = -1;
for (int i=0; i < LowPassOrder; i++)
{
this.inputValue[i] = inputValue;
this.outputValue[i] = inputValue;
}
return inputValue;
}
else if (this.inputValue != null && this.outputValue != null)
{
valuePosition = IncrementLowOrderPosition(valuePosition);
this.inputValue[valuePosition] = inputValue;
this.outputValue[valuePosition] = 0;
int j = valuePosition;
for (int i = 0; i < LowPassOrder; i++)
{
this.outputValue[valuePosition] += inputValueModifier[i] * this.inputValue[j] -
outputValueModifier[i] * this.outputValue[j];
j = DecrementLowOrderPosition(j);
}
return this.outputValue[valuePosition];
}
else
{
throw new Exception("Both inputValue and outputValue should either be null or not null. This should never be thrown.");
}
}
private int DecrementLowOrderPosition(int j)
{
if (--j < 0)
{
j += LowPassOrder;
}
return j;
}
private int IncrementLowOrderPosition(int position)
{
return ((position + 1) % LowPassOrder);
}
}
}
Keith
Autres conseils
Bien, j'ai découvert comment obtenir les coefficients que vous avez utilisés. J'ai téléchargé Octave pour Windows et ai exécuté la commande butter (comme dans MatLab) comme suit:
[b, a] = beurre (3, .4, 'faible')
Je peux maintenant utiliser ce code avec d'autres paramètres fs et fc.