[问答]

如何利用DSP库去实现一种低通滤波的设计呢

问答对人有帮助，内容完整，我也想知道答案 0 如何利用DSP库去实现一种低通滤波的设计呢？有哪些设计步骤呢？ 0
2021-11-19 08:03:16　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × heks 该类别下有 50 个回答。邀请回答 hgimtk 该类别下有 44 个回答。邀请回答新星之火12138 该类别下有 42 个回答。邀请回答 wang21cj 该类别下有 41 个回答。邀请回答 chm5 该类别下有 40 个回答。邀请回答 hjfjsdgfjdsf 该类别下有 37 个回答。邀请回答 h1654155957.9520 该类别下有 35 个回答。邀请回答 fdjslkjd 该类别下有 35 个回答。邀请回答 uvysdfydad 该类别下有 35 个回答。邀请回答 werywer 该类别下有 35 个回答。邀请回答凤毛麟角该类别下有 34 个回答。邀请回答 jenny042 该类别下有 34 个回答。邀请回答 dfzvzs 该类别下有 34 个回答。邀请回答维生素B2 该类别下有 34 个回答。邀请回答刘洋该类别下有 33 个回答。邀请回答 meihuacg 该类别下有 33 个回答。邀请回答 hfgdzc 该类别下有 33 个回答。邀请回答 ggfvxv 该类别下有 33 个回答。邀请回答 wenminglang 该类别下有 32 个回答。邀请回答风来吴山该类别下有 32 个回答。邀请回答举报郝埃连相关推荐 • 如何利用单片机去实现一种数字滤波呢 941 • 怎样利用单片机去实现一种数字加权平均滤波算法呢 892 • 如何去实现一种基于MPU6050及卡尔曼滤波的平衡小车呢 1092 • 如何去实现一种PLL环路滤波器的设计？ 1097 • 如何利用XR806开发板去实现一种UDP通信呢 901 • 如何利用头文件去实现滤波呢 675 • 如何去实现一种DSP与PC机串行通讯的设计？ 1393 • 如何利用C#去实现一种HMACSHA256加密算法呢 1342 • 如何去实现一种STM32标准库函数呢 708 • 如何去实现一种巴特沃斯低通滤波器的设计呢 1736 1个回答

答案对人有帮助，有参考价值 0 　　一、需要工具　　1. 一方面下载一个STM32的内核板（这里我用的野火的stm32f429igt6核心板）　　2.matlab 　　二、测试步骤　　1.移植ST官方DSP库　　（1）首先在官网上DSP库（这里是f4的DSP库，f1的有点不一样，f1移植的时候可以参考富莱的DSP教程）　　下不下来就用我这个版本的（可能不是最新版的）　　（2）　　首先要带着库模板，然后将下面两个文件复制到工程的库（也放放库的文件下，实际上什么时候放都没有提到，只要编译器能找到就行）　　　　　　这两个文件在下载的ST官方DSP库文件里面：　　　　　　然后在中　　　　　　最终效果图：　　　　然后像工程一直保持蓝色一直使用Single Precision（这样有一点坏处就是所有的小数与小数之间）的智力编译器默认会变成双型，所以在精确要求不高的前提下，在小数加一个，这样就不会出现警告）　　　　然后添加头路径　　　　在定义添加：__FPU_PRESENT=1，__TARGET_FPU_VFP ，ARM_MATH_CM4，__CC_ARM 　　　　最后再包含#include “arm_math.h” 　　然后测试一下是否可以用　　　　可以清晰的看到当打下臂时，出现了DSP库里面的函数，然后通知移植成功　　2.用matlab设计低通软件　　（1）这里将演示如何使用DSP库，采样频率设计成45K，到最后频率为6K，并用matlab生成1K与15K的正弦波数据，将可以用数据1excel中并形成表　　　　可见K的正弦波中夹杂着15K的谐波　　（2）matlab的fdatool工具设计　　　　设计设计后：电机Targets，然后生成C头，　　　　然后将生成文件的数据复制到工程的数组里面（这里用的ST官方示例程序完成）　　　　注意：这里的工程是移植完成DSP库的工程，没有移植这些模块函数是无法使用的　　3.工程代码（只贴了一部分）：　　#define TEST_LENGTH_SAMPLES 320 　　#define BLOCK_SIZE 32 　　#define NUM_TAPS 29 　　/* -------------------- ----------------------------------------------- 　　* 输入信号和参考输出（用 MATLAB 计算）　　* 在 arm_fir_lpf_data.c 中外部定义。　　* ------------------------------------------------- ------------------ / 　　extern float32_t testInput_f32_1kHz_15kHz［TEST_LENGTH_SAMPLES］; 　　/ ------------------------------------------------ ----------- 　　* 声明测试输出缓冲区　　* ------------------------- ------------------------------------------ / 　　静态 float32_t testOutput［TEST_LENGTH_SAMPLES］ ; 　　/ ------------------------------------------------ ------------------- 　　* 声明大小为（numTaps + blockSize - 1）的状态缓冲区　　* ----------------- -------------------------------------------------- / 　　静态 float32_t firStateF32［BLOCK_SIZE + NUM_TAPS - 1］; 　　/ ------------------------------------------------ --------------- 　　 FIR 系数缓冲区使用 fir1（） MATLAB 函数生成。　　 fir1（28， 6/24）　　** --------------------------------- ----------------------------- / 　　const float32_t firCoeffs32［NUM_TAPS］ = { 　　-0.001355670276，-0.002235466149，-0.001963305054， 0.001088276156， 0.00700009428，　　0.01148389559，0.007103286684，-0.01001676917，-0.03250513598，-0.04089481756，　　-0.01371958014，0.05627118424，0.1515145153，0.2345061451，0.2674467266，　　0.2345061451，0.1515145153，0.05627118424，-0.01371958014，-0.04089481756，　　-0.03250513598，-0.01001676917，0.007103286684，0.01148389559，0.00700009428，　　0.001088276156，-0.001963305054，-0.002235466149，-0.001355670276 　　}; 　　/ ------------------------------------------------ ------------------ 　　* FIR LPF 示例的全局变量　　* ------------------------ ------------------------------------------- / 　　uint32_t blockSize = BLOCK_SIZE; 　　uint32_t numBlocks = TEST_LENGTH_SAMPLES/BLOCK_SIZE; 　　float32_t snr; 　　/ ------------------------------------------------ ---------------------- 　　* FIR LPF 示例　　* ----------------------- ---------------------------------------------------- / 　　int32_t main（void ）　　{ 　　uint32_t i; 　　浮动最小值，最大值；　　arm_fir_instance_f32 S; 　　// arm_status 状态；　　float32_t inputF32， outputF32; 　　/ 初始化输入和输出缓冲区指针 / 　　inputF32 = &testInput_f32_1kHz_15kHz［0］; 　　outputF32 = &testOutput［0］; 　　/ 调用FIR init函数初始化实例结构。/ 　　arm_fir_init_f32（&S， NUM_TAPS，（float32_t ）&firCoeffs32［0］， &firStateF32［0］， blockSize）; 　　Debug_USART_Config（）; 　　/* ------------------------------------------------ -------------- 　　为每个 blockSize 样本调用 FIR 处理函数　　 --------------- -------------------------------------------------- -- / 　　for（i=0; i 《 numBlocks; i++）　　{ 　　arm_fir_f32（&S， inputF32 + （i blockSize）， outputF32 + （i * blockSize）， blockSize）; 　　} 　　for（i=0;i《TEST_LENGTH_SAMPLES;i++）　　{ 　　printf（“%frn”，testOutput ）; 　　} 　　通过串口打印出来的数据在excel中表格描点绘制出来　　　　而且是可以清楚看到一个周期45个点，只剩下1K的正弦波（瑕疵是几个点必须舍去）　　三、同样方法得到的案例　　下面是用100K的采样频率输出的信号，1K的信号表面夹杂高次谐波　　　　通过增加的阶数，降低后的信号为　　　　遗憾，很可能造成信号的损失

2021-11-19 14:48:00 评论举报潘沪城