资源描述：
《基于Matlab系统的信号FFT频谱分析与显示.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第16卷第4期科技通报Vol.16No.42000年7月BULLETINOFSCIENCEANDTECHNOLOGYJuly2000基于Matlab系统的信号FFT频谱分析与显示121曾尚璀,沈华,俞振利(1.浙江大学信息与电子工程学系,浙江杭州310028;2.杭州职业技术学院,浙江杭州310015)摘要:给出一种用Matlab系统实现信号频谱分析与显示的方法.Matlab是具有很强的科学计算和图形显示界面的软件系统.该法可对语音信号(或其它类型的似平稳信号)进行基于FFT的短时频谱分析,频谱图的伪彩色映射及显示.频谱图的类型有宽带和窄带两种,伪彩色显示的映
2、射可按需求设置.文中介绍了基于FFT频谱分析与频谱图显示的原理,Matlab的相关功能函数,并给出了一个程序实例及其实验结果.实验结果表明,该法所需的编程简单,对设备硬件要求不高,易于实现,能快速地在普通通用PC机上完成信号的频谱分析和频谱图显示.关键词:信号处理;FFT频谱分析;伪彩色显示;Matlab中图分类号:TN912.3文献标识码:A文章编号:1001-7119(2000)04-0241-060引言语音信号(或其它类型的似平稳信号)的频谱分析是语音识别、合成与编码的最基本任务之一;而灵活和适应人类视觉效果的动态频谱显示则是语音信号处理研究和技术开发的
3、重要手段.对于许多研究者来说,如果不是把频谱分析与显示技术和手段作为研究的终极目标,而仅仅作为一种研究手段,都希望手头有一个现成的频谱分析与显示工具.目前,可用于此目的[1]的优秀工具确实不少,如美国Entropic公司的基于Unix系统的ESPS系统中的“Xwaves”,美国Kay公司的基于PC机软硬一体化的“CSLModel4300B”(计算机语音工作站)或其软件[2]版“Model3700Multi-Speech”.但是,此类系统价格昂贵,如“ESPS”价格在10万元以上,“CSL4300B”价格在7万元以上,即使“Model3700”也需1万元以上.M
4、atlab(英文MatrixLaboratory的缩写)是70年代末期MathWorksInc.公司的首席科学家CleverMoler博士为进行矩阵理论、线性代数和数值计算研究而设计的一个软件包.后来被学术界的研究者和工程界工程师们广泛采用.现在,Matlab已经发展到5.1版本,功能变得非常强大和灵活.特别由于它具有非常规范的数据组织形式、广泛的数值计算能力、专业工具箱(Toolbox)、现成的功能子程序和优秀的图形界面,使它成为科学研究和技术开发不可多[3]得的一个软件工具.收稿日期:1999-06-22基金项目:国家自然科学基金(69972064)、浙江
5、省自然科学基金(698076)和浙江省教委科研基金资助作者简介:曾尚璀,女,1958年生,广东揭阳人,实验师;沈华,男,1945年生,浙江杭州人,讲师;俞振利,男,1956年生,浙江余杭人,副教授,博士.242科技通报16卷本文作者根据语音频谱分析和频谱图显示的原理与算法,应用Matlab实现了一种可在PC机上实现的信号频谱分析与频谱图伪彩色映射显示的方法.此方法所需的编程简单,对设备要求不高,容易实现.1频谱分析与显示原理1.1离散信号的短时FFT频谱分析[4,5]现代信号频谱分析是基于离散时域的短时Fourier分析.设离散时域采样信号为x(n),n=0,
6、1,…,N-1是时域采样点序号,N是信号长度.在数字信号处理中,用加窗方法将信号分段(帧).此时x(n)表示成xm(n),n=0,1,…,N-1,m是帧的序号,n是帧同步的时间序号,N则是一帧内的采样点数(帧长).信号xm(n)的加窗离散时域Fourier变换(DTFT)为N-1jk-jknX(m,e)=∑wm(n)xm(n)e(1)n=0为了便于离散计算,采用wm(n)·xm(n)的离散Fourier变换(DFT):N-1X(m,k)=w-j2πnk/2N∑m(n)xm(n)e,k=0,…,N-1(2)n=0
7、X(m,k)
8、便是xm(n)的短时幅度谱估计.把
9、m当作时间变量,k当作频率变量,则
10、X(m,k)
11、就是信号x(n)的动态频谱.由20·log10(
12、X(m,k)
13、)就是以dB表示的动态频谱.[6]为了高效计算
14、X(m,k)
15、,可以运用快速Fourier变换算法,即FFT算法.1.2伪彩色编码的动态频谱图显示在二维平面坐标上,将m作为横坐标,k作为纵坐标,把
16、X(m,k)
17、的值作为灰度作图,则得到x(n)的动态频谱图.对
18、X(m,k)
19、的值实施伪彩色映射,可得到伪彩色的动态频谱图,以获得更好的分辨率和视觉效果,并增强频谱图的可读性.把
20、X(m,k)
21、的最小值Xmin映射为归一化的0电平,最大值Xmax映射为归
22、一化的1电平,而将
23、X(m,k)
24、线性
