dsp课程设计-基于STM的高速频谱分析仪系统.docx

目录摘要2一.设计目的和要求3二.设计原理3三.设计内容3（1）快速傅里叶变换(FFT)3（2）.窗函数的分析与计算4（3）. A/D采样的分析与计算5(4)快速傅立叶变换以及相关原理7（5）音频的频率范围及表现力度：8(6).部分代码10结论与致谢12参考文献13附录14摘要本系统是以STM32F103为核心，主要采用FIFO来做高速缓存。高速信号先通过AD采样，然后先将采样后的数据给FIFO先缓存处理，然后再通过STM32F407进行加Blackman预处理，再做1024个点FFT进行频谱分析，最后将数据显示在LCD上，以便进行人机交互！该系统可实现任意波形信号的频谱显示，以及可以自动寻找各谐波分量的幅值，频率以及相位并进行8位有效数据显示。关键词：STM32F103 数字示波器 高速缓存设计目的和要求 基于STM的高速频谱分析仪系统（1）完成STM32的寄存器开发和库函数开发。 （2）掌握嵌入式串口通信及上位机操作；实现LCD屏界面设计与开发；以及AD信号采集和实现高速频谱分析仪系统项目的开发。（3）实现项目功能并调试程序二.设计原理以STM32F103为核心，主要采用FIFO来做高速缓存。高速信号先通过AD采样，然后先将采样后的数据给FIFO先缓存处理，然后再通过STM32F407进行加Blackman预处理，再做1024个点FFT进行频谱分析，最后将数据显示在LCD上，以便进行人机交互！该系统可实现任意波形信号的频谱显示，以及可以自动寻找各谐波分量的幅值，频率以及相位并进行8位有效数据显示。三.设计内容（1）快速傅里叶变换(FFT)DFT(离散傅里叶变换)是数字信号分析与处理中的一种重要变换，它可以使数字信号处理可以在频域内采用数值运算的方法进行，大大增加了数字信号处理的灵活性，但直接计算DFT的计算量与变换区间长度N的平方成正比，当N较大时，计算量太大，所以在快速傅里叶变换(FFT)出现以前，直接用DFT算法进行谱分析和信号的实时处理是不切实际的。DFT的定义：设是一个长度为N的有限长序列，定义的N点离散傅里叶变换为 能提高DFT速度的唯一可利用的是因子。称为旋转因子，可表示为，具有以下两个重要性质。①对称性②周期性利用的周期性和对称性可把DFT的计算次数大大减小。 而本方案采用的就是ST公司官方固件库里DSP库中已经写好的关于FFT的库函数，调用的是基4复数序列1024个点的FFT函数库。（2）.窗函数的分析与计算在频谱分析过程中，加窗是提高信号分析精度的一个重要措施，对加窗函数的基本要求：时域为改善截断处的不连续状态（由于吉布斯现象造成的振荡）；频域为窗谱的主瓣窄而高，以提高分辨率，旁瓣幅值应小，正负交替接近相等，以减小泄露和假频。因此在选择窗函数的时候，应考虑被分析信号的性质和处理要求，如果仅要求精确求出主瓣频率，而不考虑幅值精度，则可选择用主瓣宽度比较窄而便于分辨的矩形窗，例如测量物体的自振频率等；如果分析窄带信号，且有较强的干扰噪声，则应选用旁瓣幅度较小的窗函数，如汉宁窗（Hanning）等。 为了使加窗函数后的功率谱和幅值谱不受窗函数的影响，必须根据一定的原则推导出恢复系数。加窗后的恢复系数一般遵守两个原则之一：幅值相等或能量相等的原则。 因此加窗模块要实现的功能：提高采样信号分析的幅值精度。频谱分析中恢复系数的使用原则 在频谱分析中，根据不同用途采用不同恢复系数，在进行倍频程和三分之一倍频程分析时，为了使频带内总能量不变，一定要采用能量相等的恢复系数；而进行谱分析时，更关心的是各峰值频率对应的幅值，此时只能采用幅值相等的恢复系数。常用的各种窗函数的恢复系数窗函数数学表达式幅值相等恢复系数功率相等恢复系数矩形窗11汉宁窗21.633海明窗1.8521.586三角窗21.732高斯窗2.3961.840指数窗1.5821.521布莱克曼窗2.3811.812平顶窗1.1101.069（3）. A/D采样的分析与计算本系统采用的AD转换芯片为TLC5510A，TLC5510A是采用高速CMOS技术，8位的，最大转换速率为20MSPS的AD转换芯片。支持+5V电源供电，内部包含采样保持电路，输出带有高阻态模式，以及带有内部参考电阻。输出数据在时钟的下降沿有效，数据流水线结构导致了2.5个时钟的延时。虽然AD转换数据输出有2.5个时钟延时，但是只要量化的结果是正确的和数据是顺序输出，这个数据输出延时我们可以不用关心。外部还需要接一个4V的参考电压，这样才能量化到0V~4V的输入信号。本设计采用STM32F1作为核心处理器，该处理器内核架构 ARM Cortex-M3，具有高性能、低成本、低功耗等特点。 主控板包括电源模块、红外通信模