[:es]Objetivos:

• Implementar un filtro Butterworh en Arduino

Diseño de filtro:

1. Segundo Orden
2. Butterworth
3. Filtro pasa bajas
4. Frecuencia de corte 400 hz
5. Frecuencia de muestreo 100kHz
6. Periodo de muestreo 10 us

 

Usamos la pagina http://www.schwietering.com/jayduino/filtuino/

para generar la ecuación en diferencias

 

http://www.schwietering.com/jayduino/filtuino/index.php?characteristic=bu&passmode=lp&order=2&usesr=usesr&sr=10000&frequencyLow=400&noteLow=&noteHigh=&pw=pw&calctype=double&run=Send

 

Código Arduino:

#include <SPI.h> // Remember this line!
#include <DAC_MCP49xx.h>
#define SS_PIN 10

DAC_MCP49xx dac(DAC_MCP49xx::MCP4921, SS_PIN);

#include “TimerOne.h” // libreria sirve para generar interrupciones por tiempo

int row; // señal de entrada
float yn; // señal de salida

//Low pass butterworth filter order=2 alpha1=0.04
class FilterBuLp2
{
public:
FilterBuLp2()
{
v[0]=0.0;
v[1]=0.0;
}
private:
double v[3];
public:
double step(double x) //class II
{
v[0] = v[1];
v[1] = v[2];
v[2] = (1.335920002785651040e-2 * x)
+ (-0.70089678118840259557 * v[0])
+ (1.64745998107697655399 * v[1]);
return
(v[0] + v[2])
+2 * v[1];
}
};

FilterBuLp2 f;

void setup() {

Timer1.initialize(10); // Dispara cada 100 us
Timer1.attachInterrupt(muestrea); // Activa la interrupcion y la asocia a muestrea

dac.setSPIDivider(SPI_CLOCK_DIV16);

dac.setPortWrite(true);

}

void muestrea () { // se ejecuta cada 100 us

int raw = analogRead(A0);

float filtered = f.step(raw);

yn=map(filtered,0,1024,0,4096); // se mapa xn de 0 1024 a 0 4096 ya que el DAC es de 12 bits

dac.output(yn); // se manda yn al DAC

}

void loop() {

}

 

Resultados:

 

Para una frecuencia de entrada de 100 Hz

 

 

Para una frecuencia de salida de 800 Hz

 

[:en]Objetivos:

• Implementar un filtro Butterworh en Arduino

Diseño de filtro:

1. Segundo Orden
2. Butterworth
3. Filtro pasa bajas
4. Frecuencia de corte 400 hz
5. Frecuencia de muestreo 100kHz
6. Periodo de muestreo 10 us

 

Usamos la pagina http://www.schwietering.com/jayduino/filtuino/

para generar la ecuación en diferencias

 

http://www.schwietering.com/jayduino/filtuino/index.php?characteristic=bu&passmode=lp&order=2&usesr=usesr&sr=10000&frequencyLow=400&noteLow=&noteHigh=&pw=pw&calctype=double&run=Send

 

Código Arduino:

 

#include // Remember this line!
#include
#define SS_PIN 10

DAC_MCP49xx dac(DAC_MCP49xx::MCP4921, SS_PIN);

#include “TimerOne.h” // libreria sirve para generar interrupciones por tiempo

int row; // señal de entrada
float yn; // señal de salida

//Low pass butterworth filter order=2 alpha1=0.04
class FilterBuLp2
{
public:
FilterBuLp2()
{
v[0]=0.0;
v[1]=0.0;
}
private:
double v[3];
public:
double step(double x) //class II
{
v[0] = v[1];
v[1] = v[2];
v[2] = (1.335920002785651040e-2 * x)
+ (-0.70089678118840259557 * v[0])
+ (1.64745998107697655399 * v[1]);
return
(v[0] + v[2])
+2 * v[1];
}
};

FilterBuLp2 f;

void setup() {

Timer1.initialize(10); // Dispara cada 100 us
Timer1.attachInterrupt(muestrea); // Activa la interrupcion y la asocia a muestrea

dac.setSPIDivider(SPI_CLOCK_DIV16);

dac.setPortWrite(true);

}

void muestrea () { // se ejecuta cada 100 us

int raw = analogRead(A0);

float filtered = f.step(raw);

yn=map(filtered,0,1024,0,4096); // se mapa xn de 0 1024 a 0 4096 ya que el DAC es de 12 bits

dac.output(yn); // se manda yn al DAC

}

void loop() {

}

 

Resultados:

 

Para una frecuencia de entrada de 100 Hz

 

 

Para una frecuencia de salida de 800 Hz

 

[:]