当前位置:首页 > 电源 > 数字电源
[导读]1 引言 在家庭影院、卡拉OK等音响系统中,实时显示音乐信号的频谱将为音响系统增添不少色彩。目前实际生产的音响系统产品,大多采用以下两种方法实现音频频谱显示:一是利用硬件滤波器和A/D转换器;二是利用DS


1 引言
   
在家庭影院、卡拉OK等音响系统中,实时显示音乐信号的频谱将为音响系统增添不少色彩。目前实际生产的音响系统产品,大多采用以下两种方法实现音频频谱显示:一是利用硬件滤波器和A/D转换器;二是利用DSP处理频谱显示。前者实现简单,但硬件成本高,后者软件和硬件实现都较复杂。这里针对单片机RAM资源少、运算速度慢的特点,提出一种切实可行的快速傅里叶变换算法实现频谱显示。

2 系统整体设计及原理
    该系统设计由单片机SST89V58RD2、音频数据采样电路、A/D转换电路、频谱显示电路等部分组成。图1为系统整体设计原理框图。

    该系统从功能上可划分成3部分:(1)音频数据采集电路实现模拟音频信号的采样保持和量化处理,包括音频采样电路和加转换电路;(2)频谱显示电路实现模拟音频信号频谱的分段显示,它将音频信号频谱划分成14段,每段按照14级量化,由VFD显示器件显示;(3)主控制器采用SST89V58RD2单片机。在完成系统其他控制任务的前提下,充分利用单片机剩余计算资源,采用优化FFT算法计算音频信号频谱,并将计算结果输出到频谱显示电路。


3 音频信号的采集和预处理
3.1 采样频率
    根据香农采样定理,一般采样频率至少应为所采样音频信号最高频率的2倍。由于人耳能够感受的频率为20 Hz~20 kHz,所以理论上采样频率最高取40 kHz。目前工业上广泛采用的采样频率大致有3种:44 kHz、16 bit的声音称作CD音质:22 kHz、16 bit的声音效果近似于立体声广播(FMStereo),称作广播音质;11 kHz、8 bit的声音称作电话音质。本文为提高频谱计算的精度,拟采用40 kHz的采样频率和8Bit的数据位长。
3.2 样本大小
    采样频率确定后,还需确定样本值,即完成一次FFT运算所需的采样点数。根据数字信号处理的基本原理,假设采样频率为Fs,采样点数为N,则FFT运算后,第n点所表示的频率为:Fn=[(n-1)×Fs]/N(1≤n≤N)。Fn若要精确到Hz,则需采样长度为(1/f)s的信号。提高频率分辨率,需增加采样点数,但这在一些实际应用中是不现实的,则采用有频率细分法,即采样比较短时间的信号,然后在后面补充一定数量的0,使其长度达到所需的点数,再作FFT,这在一定程度上能够提高频率分辨率。由于该系统是将音频信号频谱划分成14段显示,因此采用16点FFT运算,去掉第1点和第16点的结果即可。
3.3 音频数据采集电路
    A/D转换器选用TLC549,它具有8位转换结果,差分基准电压输入,3线数据串行输出接口,转换时间最大达17μs,每秒访问和转换次数达到40 000次,全部非校准误差为±0.5LSB,低功耗,最大为15 mW。因此,该器件完全满足系统40kHz采样频率的要求。
    图2为系统音频数据采集电路。音频信号来源于CD或DVD碟机输出的音乐信号,因而基本无噪音,考虑到音频信号是交流信号,电压范围在±0.1 V之间,而TLC549的输入必须是直流的,电压范围为0~5 V,所以还需整流和放大的预处理过程。A/D转换后得到8位数字值音频幅值的8位数字量。

 

4 音频频谱算法
    音频频谱值的计算采用快速傅里叶算法FFT(Fast Fourier Transform),为了提高显示器的刷新频率,系统每隔10 ms读取16次A/D转换值,得到16点实数序列,紧接着完成16点FFT运算得到16点复数序列。[!--empirenews.page--]
4.1 倒位序及其优化算法
    基2-FTT算法将原始数据倒位序存储,但运算后的结果则按正常顺序输出。原始采样数据放在数组float datalm[16]中,datalm[0]存放第1次读取的A/D转换值,datalm[1]存放第2次读取的A/D转换值,以此类推,可见第n(n=(b3b2b1b0)b)次读取的A/D转换值存放在dataIm[n]中。倒序操作后采样数据存储在float dataRe[16]中,原来第n次读取的A/D转换值存放在datalm[n](n=(b0b1b2b3)b)中。根据样本大小在系统代码段中建立倒序表数组,采用查表方式实现快速倒序操作,与移位操作等方法相比,可明显提高运算速度。
4.2 蝶形运算及其优化算法
    根据基2-FFT算法,N点FFT运算可以分成log2N级,每一级都有N/2个蝶形运算,如图3所示。 

 

    蝶形运算公式的推导过程如下:


    将式(1)化简成实部和虚部的形式,得到:


    可见每个蝶形运算的输出都是由其输入值与某一正弦函数和余弦函数的乘积累加得到的。由式(3)~式(6)编制正弦和余弦表,每次做蝶形运算时可查表加快运算速度。
    基2-FFT算法的基本思想是用3层循环完成全部N点FFT运算:(1)最里层循环处理单独的一个蝶形运算,采用查表方法实现乘法运算;(2)中间层循环完成每一级的N/2个蝶形运算;(3)最外层循环完成log2N级蝶形运算。
    由此可看出:在每一级中,最里层循环完成N/2L个蝶形运算;中间层循环控制最里层循环进行2L-1次运算。因此,中间层循环完成时,共进行2L-1xN/2L=N/2个蝶形运算。实际上最里层和中间层循环完成了第L级计算,最外层则最终完成log2N级蝶形运算。
    需要加以说明的数据是:(1)在第L级中,每个蝶形的两个输入端相距b=2L-1一个点;(2)同一乘数对应着相邻间隔为2L个点的N/2L个蝶形;(3)第L级的2L-1个蝶形因子WPN中的P,可表示为P=jx25-L,其中j=0,1,2,…(2L-1-1)。
    完成16点FFT运算的RAM需求量是128字节,而单片机SST89V58RD2的RAM共1 K字节:显示器每10 ms刷新一次,而单片机SST89V58RD2的时钟频率是40 MHz,完成一次16点FFT运算实际所需时间不到6 ms,因此该系统完全满足FFT运算的时间复杂度和空间复杂度要求。


5 频谱值在VFD上的显示
    系统要求将音频信号频谱划分成14段,每段按14级量化,再使用VFD显示器显示,因此对于FFT运算结果还要作一定转换才能输出到显示器。第n点的FFT运算结果是复数,实部是dataRe[n],虚部是datalm[i]。该点的模值除以2/N就是对应该频率下信号的幅度(对于第1个点则是除以N);该点的相位即是对应该频率下信号的相位。最后的结果保存在dataRe[i]中,因为音频信号频谱被划分成14段,所以dataRe[0]和dataRe[15]的值应该舍去。同时,dataRe[i]可能不是整数,而VFD显示器要求每个频段按照14级量化,因此还需将dataRe[i]的值量化成0~14整数,最后输出到VFD电路上显示。


6 结束语
    讨论了单片机实现音响系统频谱显示的快速傅里叶变换算法,针对SST89V58RD2单片机进行算法优化,并详细论述系统的实现方法,结果证明该方法具有可行性。

本站声明: 本文章由作者或相关机构授权发布,目的在于传递更多信息,并不代表本站赞同其观点,本站亦不保证或承诺内容真实性等。需要转载请联系该专栏作者,如若文章内容侵犯您的权益,请及时联系本站删除。
换一批
延伸阅读

9月2日消息,不造车的华为或将催生出更大的独角兽公司,随着阿维塔和赛力斯的入局,华为引望愈发显得引人瞩目。

关键字: 阿维塔 塞力斯 华为

加利福尼亚州圣克拉拉县2024年8月30日 /美通社/ -- 数字化转型技术解决方案公司Trianz今天宣布,该公司与Amazon Web Services (AWS)签订了...

关键字: AWS AN BSP 数字化

伦敦2024年8月29日 /美通社/ -- 英国汽车技术公司SODA.Auto推出其旗舰产品SODA V,这是全球首款涵盖汽车工程师从创意到认证的所有需求的工具,可用于创建软件定义汽车。 SODA V工具的开发耗时1.5...

关键字: 汽车 人工智能 智能驱动 BSP

北京2024年8月28日 /美通社/ -- 越来越多用户希望企业业务能7×24不间断运行,同时企业却面临越来越多业务中断的风险,如企业系统复杂性的增加,频繁的功能更新和发布等。如何确保业务连续性,提升韧性,成...

关键字: 亚马逊 解密 控制平面 BSP

8月30日消息,据媒体报道,腾讯和网易近期正在缩减他们对日本游戏市场的投资。

关键字: 腾讯 编码器 CPU

8月28日消息,今天上午,2024中国国际大数据产业博览会开幕式在贵阳举行,华为董事、质量流程IT总裁陶景文发表了演讲。

关键字: 华为 12nm EDA 半导体

8月28日消息,在2024中国国际大数据产业博览会上,华为常务董事、华为云CEO张平安发表演讲称,数字世界的话语权最终是由生态的繁荣决定的。

关键字: 华为 12nm 手机 卫星通信

要点: 有效应对环境变化,经营业绩稳中有升 落实提质增效举措,毛利润率延续升势 战略布局成效显著,战新业务引领增长 以科技创新为引领,提升企业核心竞争力 坚持高质量发展策略,塑强核心竞争优势...

关键字: 通信 BSP 电信运营商 数字经济

北京2024年8月27日 /美通社/ -- 8月21日,由中央广播电视总台与中国电影电视技术学会联合牵头组建的NVI技术创新联盟在BIRTV2024超高清全产业链发展研讨会上宣布正式成立。 活动现场 NVI技术创新联...

关键字: VI 传输协议 音频 BSP

北京2024年8月27日 /美通社/ -- 在8月23日举办的2024年长三角生态绿色一体化发展示范区联合招商会上,软通动力信息技术(集团)股份有限公司(以下简称"软通动力")与长三角投资(上海)有限...

关键字: BSP 信息技术
关闭
关闭