stm32使用dsp库，结合Matlab进行FIR滤波器设计

首先我们打开Matlab。在命令串口输入fdatool，按回车。

Response Type ：这里可以设置滤波器类型 lowpass(低通)，highpass(高通)，banpass(带通)，bandstop(带阻)。
Design Method（设计模式）：有IIR和FIR两大类。（这里我们用FIR）

Specify order 滤波器的阶数
Window 窗口 这里可以设置滤波器类型。

Units 这里可以设置频率的单位。
Fs 你的采样率。
Fc 截止频率。
点击Design Filter，就会显示出所设计的滤波器。这里我设计的是30阶的低通滤波器，截止频率为125Hz。我们按照下图点击。
按照下图配置，然后点击生成.h文件。

打开文件就可以看到一个数组，这个数组就是stm32dsp库里面要用到的滤波系数。
接下来是stm32的配置。
首先我们先引用dap库。
点击上方的四角形
勾选dsp库，并在main.c添加两个文件。

#include "arm_math.h"
#include "arm_const_structs.h"

在下面添加USE_HAL_DRIVER,STM32F429xx,ARM_MATH_CM4,ARM_MATH_MATRIX_CHECK,ARM_MATH_ROUNDING （我用的是f429）
接下来就是stm32的代码部分了，首先先导入滤波系数，将数组改成浮点型。

const int BL = 30;
 float B[30] = {
 0.0005415165797, 0.001727426774, 0.003113061655, 0.003883380443, 0.002264012117,
  -0.003352924949, -0.01260423567,  -0.0218139533, -0.02428940125, -0.01267496031,
    0.01740843058,  0.06416342407,     0.11868231,   0.1671814919,   0.1957704276,
     0.1957704276,   0.1671814919,     0.11868231,  0.06416342407,  0.01740843058,
   -0.01267496031, -0.02428940125,  -0.0218139533, -0.01260423567,-0.003352924949,
   0.002264012117, 0.003883380443, 0.003113061655, 0.001727426774,0.0005415165797
};

按照上面的步骤配置正确，我下面封装的函数是能直接用的。（input,output数组需定义在main函数以外，不让程序会卡死）

void FIR(int Length ,int  BLOCK_SIZE ,int BL ,float32_t *input,float32_t *output,float *B)// Length 代表采的点数，BLOCK_SIZE大于1小于你的采样点数即可，BL代表滤波阶数，inout,output,代表输入数据和存储输出数据的数组，B是滤波系数。
{
uint32_t i;
uint32_t blockSize = BLOCK_SIZE;
uint32_t numBlocks = Length /BLOCK_SIZE;
float32_t State[BLOCK_SIZE + BL  - 1];
float32_t  *inputF32, *outputF32;
arm_fir_instance_f32 S;
	inputF32 = &input[0];
 	outputF32 = &output[0];
	
  	arm_fir_init_f32(&S, BL , (float32_t *)&B[0], &State[0], blockSize);
	
		for(int i=0; i < numBlocks; i++)
		{
			arm_fir_f32(&S, inputF32 + (i * blockSize), outputF32 + (i * blockSize), blockSize);
			
		}
}

记得在main.h里添加下面的文件，不然会有一大堆报错

下面是我的测试程序，这是通过串口直接给Matlab发送数据并画图，让我们更能直观的看到其滤波效果。具体怎么配置，可以查看我的另一篇文章，超简单的。单片机与matlab之间的串口通信

for(int i=0;i<1000;i++)
{
	input[i]=sin(2*i*PI*5/1000)+sin(2*i*PI*200/1000);//5hz和200hz的正弦波信号
	
}
  while (1)
  {
		FIR(1000,100,30,input,output,B);
		HAL_UART_Receive_IT(&huart1,(uint8_t*)aRxBuffer,1);// 串口持续接受信号
  //  printf("%d\r\n",aRxBuffer[0]);
		if(aRxBuffer[0]==100)// 判断握手信号是否到来
		{
				for(int i=0;i<1000;i++)
	{
		D_data.data=output[i];//将数据赋给共同体
		HAL_UART_Transmit(&huart1,D_data.ComInfo.com, 4, 0xffff);//将4个字节发送出去
		HAL_Delay(5);
		
	}
  aRxBuffer[0]=1;// 清空标志位
		}
	}

将D_data.data=output[i];改成D_data.data=input[i],在运行一次Matlab的程序，就可以得到一张很直观的图。下面是我测试的效果，红色的是200Hz的和5hz叠加正弦信号，蓝色的是滤波后的信号，滤波效果还是不错的。
创造不易，感谢关注与点赞。