自数字信号处理实验8音频频谱分析仪设计与实现.docx
《自数字信号处理实验8音频频谱分析仪设计与实现.docx》由会员分享，可在线阅读，更多相关《自数字信号处理实验8音频频谱分析仪设计与实现.docx（36页珍藏版）》请在冰豆网上搜索。
自数字信号处理实验8音频频谱分析仪设计与实现
信号处理实验
实验八:
音频频谱分析仪设计与实现
一、实验名称:
音频频谱分析仪设计与实现
二、实验原理：
MATLAB是一个数据信息和处理功能十分强大的工程实用软件,其数据采集工具箱为实现数据的输入和输出提供了十分方便的函数和命令。
本实验可以用ＭATＬAB进行音频信号频谱分析仪的设计与实现。
1、信号频率、幅值和相位估计
（１）频率（周期）检测
对周期信号来说,可以用时域波形分析来确定信号的周期,也就是计算相邻的两个信号波峰的时间差、或过零点的时间差。
这里采用过零点（ti）的时间差T（周期）。
频率即为ｆ ＝1/T,由于能够求得多个Ｔ值（tｉ有多个），故采用它们的平均值作为周期的估计值。
（２）幅值检测
在一个周期内,求出信号最大值yｍａx与最小值ｙmin的差的一半，即A=（ymａx-ｙmin）／2,同样，也会求出多个A值，但第1个A值对应的ｙmax和yｍiｎ不是在一个周期内搜索得到的，故以除第１个以外的A值的平均作为幅值的估计值。
（3）相位检测
采用过零法，即通过判断与同频零相位信号过零点时刻，计算其时间差,然后换成相应的相位差。
φ＝2π（1-tｉ/Ｔ）,{x｝表示x的小数部分,同样，以φ的平均值作为相位的估计值。
频率、幅值和相位估计的流程如图所示。
其中tin表示第n个过零点，yi为第i个采样点的值，Fs为采样频率。
2、数字信号统计量估计
（1）峰值P的估计
在样本数据ｘ中找出最大值与最小值,其差值为双峰值,双峰值的一半即为峰值。
P＝０.5［max（yi）-mｉｎ（ｙi）]
（２）均值估计
式中，N为样本容量,下同。
（3）均方值估计
（4）方差估计
2、频谱分析原理
时域分析只能反映信号的幅值随时间的变化情况，除单频率分量的简单波形外,很难明确提示信号的频率组成和各频率分量大小，而频谱分析能很好的解决此问题。
（1）DＦT与ＦFT
对于给定的时域信号y,可以通过Fｏｕrｉer变换得到频域信息Ｙ。
Ｙ可按下式计算
式中，N为样本容量,Δt=1／Fs为采样间隔。
采样信号的频谱是一个连续的频谱，不可能计算出所有的点的值,故采用离散Ｆourｉｅr变换（DFＴ）,即
式中,Δf=Fｓ／N。
但上式的计算效率很低,因为有大量的指数（等价于三角函数）运算,故实际中多采用快速Ｆoｕrier变换（ＦFＴ）。
其原理即是将重复的三角函数算计的中间结果保存起来,以减少重复三角函数计算带来的时间浪费。
由于三角函数计算的重复量相当大,故ＦFT能极大地提高运算效率。
（2）频率、周期的估计
对于Ｙ（kΔｆ）,如果当kΔf= 时,Y（kΔｆ）取最大值，则为频率的估计值,由于采样间隔的误差,也存在误差，其误差最大为Δf ／２。
周期Ｔ=1/f。
从原理上可以看出,如果在标准信号中混有噪声，用上述方法仍能够精确地估计出原标准信号的频率和周期，这个将在下一章做出验证
3、频谱图
为了直观地表示信号的频率特性,工程上常常将Ｆouｒier变换的结果用图形的方式表示，即频谱图。
以频率ｆ为横坐标，|Y（f）|为纵坐标,可以得到幅值谱；
以频率f为横坐标，aｒgＹ（f）为纵坐标，可以得到相位谱；
以频率ｆ为横坐标,ReY（ｆ）为纵坐标，可以得到实频谱;
以频率f为横坐标,ImＹ（f）为纵坐标,可以得到虚频谱。
根据采样定理,只有频率不超过Fｓ／2的信号才能被正确采集，即Ｆｏuｒier变换的结果中频率大于Ｆs/2的部分是不正确的部分,故不在频谱图中显示。
即横坐标f∈[0,Fs/2]
4、模块划分
模块化就是把程序划分成独立命名且可独立访问的模块,每个模块完成一个子功能,把这些模块集成起来构成一个整体,可以完成指定的功能满足用户需求。
根据人类解决一般问题的经验，如果一个问题由两个问题组合而成，那么它的复杂程度大于分别考虑每个问题时的复杂程度之和,也就是说把复杂的问题分解成许多容易解决的小问题,原来的问题也就容易解决了。
这就是模块化的根据。
在模块划分时应遵循如下规则:
改进软件结构提高模块独立性；模块规模应该适中；深度、宽度、扇出和扇入都应适当；模块的作用域应该在控制域之内;力争降低模块接口的复杂程度;设计单入口单出口的模块;模块功能应该可以预测。
三、实验要求:
参考以上原理，查阅相关资料,构建交互界面,设计一个音频频谱分析仪，实现一下功能:
（1）音频信号信号输入,从声卡输入、从WAＶ文件输入、从标准信号发生器输入;
（2）信号波形分析,包括幅值、频率、周期、相位的估计,以及统计量峰值、均值、均方值和方差的计算；
（3）信号频谱分析,频率、周期的估计，图形显示幅值谱、相位谱、实频谱、虚频谱和功率谱的曲线。
四、实验界面及程序
程序：
fuｎctionvarargout =yinｐin（ｖarａrgin）
％YＩNPIN Ｍ－ｙinpｉn.fig
% YINPIN,byitsｅｌf,cｒeaｔesaｎeｗYＩNPINoｒrａises the eｘisｔing
％singｌｅｔon＊.
%
% H= ＹINＰINreturｎｓtheｈaｎdletoaneｗYINPINorthｅhａnｄｌe ｔo
% the ｅxiｓtingｓingleｔｏｎ*.
%
% YINPＩN（＇CALLBACK'，hOｂｊｅｃt,eventData，hａndlｅｓ,．．.）ｃalｌstheloｃal
% fuｎctionnamｅd CALLＢACKinYINPIN．Ｍwithｔhegｉｖenｉnputａrｇumeｎts.
％
％YINPIN（'Ｐｒoｐｅｒｔy'，'Valuｅ＇,..．）crｅatesaｎewYＩNPIN ｏｒｒaiseｓ thｅ
％ exisｔｉng ｓingleton＊.Startinｇｆrｏm thelｅft，pｒｏpｅrｔｙｖａlｕepairs ａｒe
% applｉed tｏthｅGUＩｂｅｆore yinpin_OpenｉｎgFunctｉongｅts callｅｄ．An
% ｕｎrecognｉzeｄpｒopｅrty name orinｖalｉｄ valｕemakesproperｔyａppｌication
％ stop.Ａlｌinpuｔsareｐasｓｅdtoｙinpiｎ_OｐenｉngFcn ｖiａvarａｒｇin.
％
% *SeeGＵIOptionsｏｎGUIDE'sTools menu. Choose"ＧUIallowsonlyonｅ
%insｔaｎｃｅtｏｒun（ｓiｎgleton）＂.
％
%Ｓee alｓo:
GUＩDE,ＧUIDATＡ,GUIHＡNＤＬＥS
%Edit ｔheａbovｅ tｅxt ｔｏｍoｄifytheｒesｐonsetoｈelpyinｐｉn
%LaｓtMoｄiｆieｄbyGUIDEｖ２.5 07-Ｎｏv-20１21７:
08:
53
% Beginiｎｉｔialiｚatｉonｃode－DＯNOTＥDＩT
gｕi_Singleton=1;
gui_Statｅ=struct（'gui_Name'， m， ．．.
'gui_Sinｇleｔon',guｉ_Sｉngleton,．..
'gui＿OpeninｇFcn'， @ｙinpｉn＿ＯpｅningFcn， ...
'ｇui＿OuｔpuｔFcn'，@yinｐin_ＯutputＦcｎ, ..．
'ｇuｉ_LayｏutＦｃｎ', ［] ，...
'ｇui_Callback＇,[]）;
ｉf nargin&&iscｈar（ｖararｇin{1｝）
gui_Statｅ.ｇuｉ_Callbａck=ｓｔr２ｆｕnｃ（vａrargin｛1}）;
end
if nａrｇｏuｔ
[varａrgoｕt{1：
naｒgｏｕt｝］=gui＿ｍａiｎｆｃｎ（ｇuｉ_Sｔａte,varargｉn｛:
}）;
else
gui＿ｍaｉnｆｃn（gui＿Sｔaｔe,vaｒａrgin{：
｝）；
end
%Endinｉtiaｌizaｔｉon codｅ －DOＮOTEDＩT
%－-－Ｅｘｅcutｅｓ jusｔbeｆore ｙinpｉnｉs made ｖisｉbｌe.
functionyinｐin_OpenｉngFcｎ（hObject, eventｄａｔa,hanｄles, ｖararｇｉn）
%Thisｆunｃtiｏｎ has nooutpｕｔaｒｇs, seeOutｐutＦcｎ．
% hObjecthａｎdle toｆｉgure
％ eventdatareserved －ｔobｅ definｅdiｎafutureversioｎｏｆＭATLAＢ
％handｌeｓstruｃtuｒeｗｉthhandles anｄ userdａta（seeGUIＤATＡ）
%ｖarargincomｍａndline aｒｇumentｓtoyinpｉｎ（seeVAＲＡRGＩN）
%Cｈoosedefault coｍｍandlinｅｏutput ｆoｒyinｐin
ｈandles.oｕtｐut =hObｊect；
%Ｕpdatehandleｓstｒucｔｕre
ｇuidａtａ（hObjｅcｔ,hａndles）;
%ＵIWＡIＴmａkes ｙiｎｐｉnwａitfｏruｓｅｒ ｒesponｓe（ｓeeUＩREＳUME）
% ｕiwaｉt（handles.figuｒe１）;
%-－-Ouｔｐutｓｆrｏmthisfunctioｎaｒeｒｅtｕｒnedto thecomｍaｎdline.
fｕnctionｖaｒaｒgouｔ＝ｙinpｉn_OutpｕｔＦcn（hOｂjeｃt， ｅvenｔdａta, handles）
% ｖaｒarｇouｔ ｃell arｒａy forretuｒniｎｇoutｐutａｒｇs（sｅｅVＡRＡRＧＯＵT）;
%ｈObject handｌeｔo fｉguｒe
％ｅventｄataｒesｅｒｖed－tｏ bedeｆinedｉna ｆuｔure version ｏｆＭＡTLＡＢ
% ｈａndleｓ ｓｔruｃtｕrewithhaｎdlesandｕseｒdaｔａ（sｅｅＧUIＤＡTA）
%Gｅtdｅｆaultcｏmｍａndlineoutｐutfromhanｄles structure
varargout{１}= hanｄles.outｐuｔ;
% －-－ Ｅxｅcutesdurｉｎgoｂjecｔcreation, ａftersettingallpropertｉes.
funcｔionedｉt1_CrｅateFcn（ｈObｊecｔ,eｖeｎtdata， handles）
%hObｊecthandletoｅdiｔ１（see GCBO）
%eｖｅntdataresｅｒｖｅｄ- tｏbe deｆiｎeｄiｎafutureveｒsionｏf MＡTLAB
％handｌeｓeｍpｔy-handlｅs notcrｅated untilａfｔeraｌlCreateFｃnｓcalled
%Hinｔ:
editｃｏntroｌsｕsｕａｌｌyhave awhiteｂaｃkｇrouｎdｏnＷindｏws.
% SeeIＳPCａndCＯMPUTER.
ifisｐc&＆isｅqual（get（hＯbject,'BackgrouｎdColｏr'）, get（0,'defaultＵiｃonｔｒolBacｋgｒoｕｎｄＣｏlor'））
set（ｈObject,＇BackgroundColor'，＇whitｅ＇）;
end
％---Ｅxecuｔｅｓduｒiｎｇ objｅｃtcｒeation,after settingalｌ proｐertｉes.
functioncaiyaｎgpｉnlv＿CreateFｃｎ（hObject, eventｄata,hａndles）
%hOｂｊect haｎdlｅtocａｉyaｎｇpiｎlv（seeGＣＢＯ）
%evｅntdatａreｓeｒｖeｄ- ｔobeｄeｆiｎｅｄinaｆuｔuｒe versｉon of MＡＴLＡB
%hanｄlｅｓemｐｔｙ -ｈandlｅｓnoｔcｒeated untilａftｅrallCｒeａteFcｎscaｌled
%Hiｎt:
ｅdiｔcontrｏls usuallyhavｅ awhitｅbaｃkｇｒoｕnｄoｎＷindows.
% Sｅｅ ＩSPC and CＯＭPUTER.
ｉfiｓpc＆&iseｑuａｌ（get（ｈObｊect,'BaｃｋgｒouｎdColor'）,get（0,'defａｕltＵiｃonｔrolBaｃｋgrouｎdColoｒ'））
set（hObject，'BacｋgrouｎｄColｏr'，'wｈitｅ'）;
eｎd
%--－Executesｄurｉngobｊectcreaｔion，aｆｔersettinｇ all prｏpertｉes.
ｆunｃｔｉonｃaiyaｎｇdiａnshｕ_ＣreateFｃn（hＯbject,eventｄａｔa,haｎdlｅs）
％ｈObjecｔ haｎｄｌe to caiｙanｇdｉaｎｓhｕ（sｅeGＣBO）
%eveｎtdatａｒeserｖed-tobｅｄeｆiｎeｄ ｉnａｆuturｅｖeｒｓion ofMＡＴLAB
%handles eｍpty-hanｄles notcreaｔeduntilaftｅrallCreａteFcnsｃalled
% Ｈint:
eｄitcｏｎtｒols usualｌyhavｅaｗhiteｂackgｒoundon Ｗinｄｏws.
% SeeＩSPCandCOMPＵTＥR.
ifispｃ&&ｉseqｕａｌ（geｔ（hＯbjeｃt,＇BａckgｒｏuｎdＣoloｒ'）,ｇet（０,'ｄｅfａultUｉcontrolＢaｃｋgrｏuｎdCoｌor'））
seｔ（ｈOｂｊect,'BａｃkgrouｎdColoｒ','wｈｉｔｅ＇）；
enｄ
％ -－- Execuｔｅsonbuｔtｏnpresｓinluyin．
fｕnｃtｉonluyｉn＿Ｃａｌlbacｋ（hObjeｃt，eventdata，hanｄlｅs）
%hObjecｔｈaｎdleto ｌｕyin（seeＧCＢO）
%eventdata rｅｓｅrｖed -to be dｅfineｄin aｆutｕre ｖｅrｓiｏn ofMAＴＬAB
%hanｄlｅsstructurｅ witｈhanｄleｓandusｅrdata（seeＧUIDATA）
%Hint:
geｔ（hObjeｃt,＇Vａlue'） reｔｕrns togglesｔatｅofluyiｎ
set（handｌｅs.luyin,'value',1）;
set（ｈａndles.ｗavfiｌe,＇value',0）；
seｔ（haｎｄｌes.xinhaｏfasｈｅng,'vａlue＇,0）；
h=fｉndｏbｊ（＇Tａg',＇recordtime'）；ﻩ
ｓet（ｈ,'enａbｌｅ＇,'oｎ'）;
h=fｉndｏbｊ（'Ｔag'，'ｋａｉshiｌｕｙｉn'）;
ｓｅt（ｈ,＇enａble＇,'on'）;
h=findｏbj（'Tａg','wavnａme'）；
set（h,'eｎａbｌe'，'off'）；
h=ｆｉｎdoｂj（'Tag'，＇liｕlａｎ'）;
sｅt（h，＇eｎａble＇,'off＇）;
h=ｆindobｊ（'Tag＇，＇inwave'）;
set（ｈ,'enａbｌｅ',＇off'）;
h=ｆiｎdobｊ（'Tａg＇，'ｉnfudu'）;
ｓｅt（h,'enａble＇,'off'）；
h=fiｎdobj（'Tag'，＇iｎfｒｅ'）;
ｓｅt（h,＇enable＇,'off＇）；
h＝fｉｎdｏｂｊ（'Tａg',＇ｉnｐhase＇）；
ｓeｔ（h，'ｅｎable'，'off'）;
ｈ＝ｆiｎdｏbj（'Tag'，'ｍiｘ'）;
set（h,'enaｂｌe',＇off＇）；
h＝fiｎｄobｊ（'Tag',＇sｈengｃｈeｎgｂｏxing'）；
set（ｈ,＇eｎａｂle'，'off'）;
％ －--Exeｃuteｓoｎ buttonｐreｓsinwavfile.
functiｏｎ wav（hＯbjｅcｔ, eveｎtdatａ，handles）
% hObject hａndleｔo wａvfｉｌe（seeＧＣBO）
％ eveｎtdata rｅservｅｄ-ｔoｂedｅｆineｄiｎaｆｕtureverｓｉｏｎoｆＭＡTLAB
% hanｄlessｔrucｔureｗiｔh ｈaｎｄlesanduser datａ （ｓeｅＧUIDATＡ）
% Hint:
get（hＯbjecｔ,'Valｕe'）returnｓtogｇlestatｅofwａvｆilｅ
set（handles.luyｉｎ,＇vaｌue',０）；
seｔ（hａndles.wavfiｌe，'valｕｅ',1）;
ｓet（ｈａndleｓ.xｉnhaofａshｅnｇ,'value＇，0）;
h=fｉｎdobｊ（'Ｔag＇，'reｃoｒdｔimｅ＇）;ﻩ
ｓet（ｈ,'enaｂｌe',＇off'）；
h=fiｎdobj（'Ｔag＇,'kaｉshiluyｉn'）;
set（h,'eｎaｂle＇,'ｏff＇）；
h=ｆｉnｄobj（'Ｔag','wａvｎame'）;
ｓet（h,'enａble','ｏn'）;
h=finｄｏｂj（＇Tag','liulａn＇）；
set（ｈ,＇eｎａblｅ','on'）;
h=findobｊ（＇Tａg＇,'inwaｖe'）;
sｅｔ（ｈ,＇enable'，'off'）；
h=findobｊ（'Tag',＇ｉnfｕｄu'）;
sｅt（ｈ,'enａble',＇ｏｆf＇）;
h＝fｉndobj（＇Taｇ','ｉｎfre'）；
set（h，＇ｅnable','off'）;
h=ｆindobj（'Taｇ＇,＇inpｈaｓe＇）;
seｔ（h，'eｎable'，＇oｆｆ'）;
ｈ=findoｂｊ（'Tag'，'mix'）;
set（ｈ,＇ｅnａble','off'）;
h=fiｎｄoｂj（'Tａg','sｈengｃhengboｘing＇）;
sｅｔ（ｈ,'enaｂle',＇ofｆ'）；
%--－Executesonbuttonｐresｓinxｉnhａofａshｅnｇ.
funcｔioｎxiｎhaofasheng＿Cａllback（hObｊeｃｔ,ｅventdatａ,handleｓ）
％ hObject haｎdletｏxｉnhaofasheng（seｅGCＢO）
％ ｅveｎtdatａreserveｄ-to bedefiｎｅｄ inafuturｅｖｅrsioｎｏfMATLAB
%handles ｓtrｕcturｅwiｔh hａnｄlesanｄuserdata（seeGUIDATA）
％ Hｉnｔ:
get（hＯｂjｅct,'Valuｅ'）rｅturｎstoｇglesｔatｅof xｉnｈaoｆasheng
ｓet（handles.luyiｎ,'ｖaｌue'，0）;
set（handｌeｓ.wavfiｌe,＇vaｌue',0）;
set（hanｄles.xinhａofasheｎg,＇vaｌue＇,1）;
ｈ=findobj（＇Tag'，'reｃordｔimｅ'）；
seｔ（h，'enabｌｅ＇，＇ofｆ'）；
ｈ=findoｂj（＇Ｔag'，'kaｉshiｌuyin'）;
set（h，'ｅｎａble','ofｆ＇）；
h=fｉndobj（＇Tag','wａvnａｍe'）；
ｓet（h,'eｎable＇，'off'）；
h＝fiｎｄobj（'Ｔaｇ',＇liｕlan'）;
ｓet（h,＇eｎable','off＇）;
ｈ＝findoｂｊ（'Tag',＇inwave'）;
ｓet（h,＇eｎable','on'）;
ｈ＝findobj（＇Taｇ',＇infuｄｕ＇）;
ｓet（h，'ｅnable＇，'ｏn'）;
h=findobｊ（'Ｔaｇ',＇ｉｎｆre＇）;
ｓｅt（h,'eｎabｌe',＇on'）；
h=fiｎdobj（'Tag','inｐhａsｅ'）;
set（h,'enable＇,'ｏn'）;
h=ｆindｏbｊ（'Ｔaｇ','miｘ'）;
set（ｈ，'enaｂlｅ＇,＇on＇）;
h=ｆiｎｄoｂj（'Taｇ','ｓhｅngchｅｎgboｘing'）;
seｔ（h，＇enabｌe＇，'ｏn'）;
%－－-Execｕｔeｓ duｒｉｎgｏbjectcｒeaｔｉｏn， ａftｅr settingａｌlpｒoｐｅrties.
fuｎcｔionrｅｃｏｒdｔｉｍｅ＿CreatｅＦcn（hOｂjｅｃt, ｅveｎｔdaｔa,hanｄleｓ）
％ｈObject handle tｏ ｒecｏrdｔime （sｅeGCBＯ）
%ｅvenｔdａtａ rｅsｅrveｄ-tｏｂedｅfinediｎ afuturevｅrsionｏfＭATＬAB
%handｌesempｔy -hａndlesnoｔ created untiｌaftｅｒａl
升级会员 