最近在弄STM32自带的DSP库里的滤波,记录一下:
arm_fir_instance_q15 instance_q15_S;
#define NUM_TAPS 16 //滤波系数的个数
#define BLOCK_SIZE 32
q15_t firStateF32[BLOCK_SIZE + NUM_TAPS];
q15_t Fir_Coeff[NUM_TAPS] = {-79, -136, 312, 654, -1244, -2280, 4501, 14655, 14655, 4501, -2280, -1244, 654, 312, -136, -79};
q15_t firStateF32[BLOCK_SIZE + NUM_TAPS];//1、初始化S
arm_fir_init_q15(&instance_q15_S,NUM_TAPS,(q15_t*)&Fir_Coeff[0],&firStateF32[0],BLOCK_SIZE);//2、用三角函数生成一个数组,
for (i = 0; i < SAMPLE_NUMBER; i++)
{f = 0.5 * arm_cos_f32(2 * 3.1415926f * 50 * i / SAMPLE_NUMBER) + 0.5 * arm_cos_f32(2 * 3.1415926f * 6000 * i / SAMPLE_NUMBER);arm_float_to_q15(&f, &testInput[i], 1); //f * 32768
}//3、滤波
for (i = 0; i < SAMPLE_NUMBER / BLOCK_SIZE; i++)
{arm_fir_q15(&instance_q15_S, testInput + (i * BLOCK_SIZE), testInput + (i * BLOCK_SIZE), BLOCK_SIZE);
}
整个滤波的过程比较简单,先是调用arm_fir_init_q15函数来初始化instance_q15_S,然后用三角函数生成一个50Hz+6000Hz的正弦数组,放在testInput数组中,最后调用arm_fir_q15函数来滤波。
这个里面需要注意的是,滤波系数的个数,也就是NUM_TAPS必须是偶数。这个翻了些文档,都没讲,若是奇数的话一运行就进HARDFAULT。打开arm_fir_init_q15这个函数的原型,里面讲到了参数的设置。
下面是函数的路径:
用软件打开可以看到
还有就是firStateF32这个数组,它的大小numTaps+blockSize,不要减一,因为我用的STM32F407的芯片,减一的话也会进hardfault。
这里可以看到,内核是M4的话,数组大小是numTaps+blockSize;内核是M3的话,数组大小是numTaps+blockSize-1。