I want to design a simple IIR filter with mathematica. For example, consider the filter deigned with this online IIR calculator:
y[n] = x[n- 2] + 2*x[n- 1] + x[n- 0] - 0.8752145483*y[n- 2] + 1.8668922797*y[n- 1]
Butterworth, Lowpass, filter order: 2, sample rate: 1000, corner frequency 1: 15
I can certainly envision just using entering this equation into Mathematica with a
For loop, but that doesn't sound satisfactory to me. Especially for a higher order IIR filter.
Needless to say the mathematica page about this is way above my pay grade. I just want to implement a simple filer. It is easy in Matlab, so it should be easier in Mathematica, right?
My question is:
How can I implement a simple low pass IIR filter with Mathematica?
The correct answer will allow all the order to increase to a higher number than 2. It will can use the in-built Mathematica functions, or some other custom means. The correct answer will not have a procedural loop. I assume the equation above can somehow be entered using a Mathematica style list?