Nice work, bit for those not familiar with this theory (like myself) could you provide links to some background info, or provide a playground article?
coaxed into adding a graph showing the frequency response?
What I am actually interested in is not a graph of the theoretical response (you can find in the net) but a simulation of the calculation, so I can see if the integer calculation optimizations are ok for the kind of signals I got.
long tmp = ((((data * 662828L) >> 4) //= ( 7.901529699e-2 * x) btw, should it be 662828 / 2^23 = 0.079015255? + ((_v * -540791L) >> 1) //+( -0.5157387562*v) + (_v * 628977L) //+( 1.1996775682*v) )+262144) >> 19; // round and downshift fixed point /524288 _v= (short)tmp;
int _order; short _v;
int _order; int _v; or int16_t _V;