數(shù)字信號(hào)處理研討_DFT近似計(jì)算信號(hào)頻譜

1、beijing jiaotong university數(shù)字信號(hào)處理研討dft近似計(jì)算信號(hào)頻譜學(xué) 院：電子信息工程學(xué)院小組成員：指導(dǎo)教師:2013.06.24dft近似計(jì)算信號(hào)頻譜專(zhuān)題研討【目的】(1) 掌握利用dft近似計(jì)算不同類(lèi)型信號(hào)頻譜的原理和方法。(2) 理解誤差產(chǎn)生的原因及減小誤差的方法。(3) 培養(yǎng)學(xué)生自主學(xué)習(xí)能力，以及發(fā)現(xiàn)問(wèn)題、分析問(wèn)題和解決問(wèn)題的能力?！狙杏戭}目】基本題1. 己知一離散序列為xk = sin(0.2 兀眾)，眾=0，1，31 用l=32點(diǎn)dft計(jì)算該序列的頻譜，求出頻譜中譜峰的頻率；(2) 對(duì)序列進(jìn)行補(bǔ)零，然后分別用l=64、128、256、512點(diǎn)dft u算該

2、序列的頻譜，求出頻譜中譜峰 的頻率；(3) 討論所獲得的結(jié)果，給出你的結(jié)論。該結(jié)論對(duì)序列的頻譜計(jì)算有何指導(dǎo)意義？【題目分析】本題討論補(bǔ)零對(duì)離散序列頻譜計(jì)算的影響?！緶仨嗵崾尽吭谟?jì)算離散非周期序列頻譜吋常用x2/7t作為橫坐標(biāo)，稱(chēng)/2/k為歸一化頻率(normalized frequency)。 在畫(huà)頻譜時(shí)需給出橫坐標(biāo)。每幅圖下都需給出簡(jiǎn)要的文字說(shuō)明。由于離散非周期序列頻譜是周期的，所以在計(jì)算時(shí)不必用ffishift函數(shù)對(duì)fft計(jì)算的結(jié)果進(jìn)行重 新排列?！拘蛄蓄l譜計(jì)算的基本方法】連續(xù)信號(hào)，通過(guò)其抽樣的離散信號(hào)，和離散信號(hào)的dft變換存在如下關(guān)系：n-x(ejq) n =dft氏) = x aot

3、cq=vzwa-=o通過(guò)如上關(guān)系，我們就可以通過(guò)dft來(lái)求信號(hào)的頻譜?！痉抡娼Y(jié)果】 dd k 、s os 矚0fm =0.9375fm =0.9219fm =0.9141fm =0.9102fm =0.9082【結(jié)果分析】增加dft的點(diǎn)數(shù)可以使頻譜更容易觀察，即減輕了柵欄效應(yīng)帶來(lái)的影響。頻譜的橫坐標(biāo)為歸一 化頻率，所以原信號(hào)的峰值第一次應(yīng)該出現(xiàn)在0.2處，隨著dft點(diǎn)數(shù)的增大，頻譜表示也越來(lái)越精 確。【自主學(xué)習(xí)內(nèi)容】1. 歸一化頻率相關(guān)知識(shí)。2. 通過(guò)matlab計(jì)算dft和matlab的繪圖操作。【閱讀文獻(xiàn)】1. 數(shù)字信號(hào)處理2. 補(bǔ)零對(duì)有限長(zhǎng)序列頻譜及dft的影響?！景l(fā)現(xiàn)問(wèn)題】(專(zhuān)題研討或

4、相關(guān)知識(shí)點(diǎn)學(xué)習(xí)屮發(fā)現(xiàn)的問(wèn)題)：dft點(diǎn)數(shù)的增多是否能提高頻譜分辨率？【問(wèn)題探究】雖然dft點(diǎn)數(shù)增加使圖像更加細(xì)致，但是因?yàn)椴徽揹ft點(diǎn)數(shù)是多少，抽樣點(diǎn)數(shù)都是相同的，所 以每個(gè)頻譜所包含的信息相同，頻譜分辨率只與抽樣點(diǎn)數(shù)相關(guān)，與點(diǎn)數(shù)無(wú)關(guān)，頻譜分辨率相同。 所以不能通過(guò)增大dft點(diǎn)數(shù)而減少信息損失。dft點(diǎn)數(shù)的增多不能提過(guò)頻率分辨率?！痉抡娉绦颉?k = 0:31;n = 16;x = sin (0.2*pi*k);for i = 1:5n = 2*n;x = f ft (x, n); k = 0:n-1; subplot (5, 1, i); plot (k/n,abs(x); xlabel (

5、f/hz; ylabel (丨magnitude *); fm = find( x = max(x)/n; fmend2已知一離散序列為x k=acosfi)k+bcos (z)+a/2)z:)o用長(zhǎng)度ams4的哈明窗對(duì)信號(hào)截短后近似計(jì)算其頻譜。試川不同的和z?的值(如和打近似相等，4和z?差距較大)，確定用哈明窗能分辯的最小的譜峰間隔a/2w =十c的值。n【題目分析】本題討論用哈明窗h算序列頻譜時(shí)的頻率分辨率。【仿真結(jié)果】c=4 a/b=16c=4 ab=420100201002010c=4 >vb=120100201002010020100152025c=3 ab=12c=4 ab

6、=1220100152025c=22 a/b=1220100152025201002010020100201002010v152025c=3 a/b=16152520101520c=22 ab=16252010201015252025152025152025c=2 a«=12c=2 a/b=1gc=2 ab=41520c=19 a/b=122010020100251520c=19 ab=16152025c:18 ./vb=122010025201002010020100152025c:1861520c:18 abm2010251520c=15 ab=161520c=15，a/b=12

7、2010020100«:20252520100z/15202s2010252015file editinsert took desktopwindow help、【結(jié)果分析】第一排是c=4,第二排c=3,遠(yuǎn)大于教材定義的標(biāo)準(zhǔn)c=2。無(wú)論a/b的比值是多少，區(qū)分譜峰毫 無(wú)壓力。第3排巾，c=2是教材定義的標(biāo)準(zhǔn)，a與b接近的時(shí)候，區(qū)分較為容易，當(dāng)a/b為1.6的時(shí)候, 幾乎是無(wú)法區(qū)分，而a/b=4的時(shí)候，完企無(wú)法區(qū)分。所以說(shuō)在沒(méi)計(jì)使用哈明窗的時(shí)候，如果不同頻 譜的信號(hào)所占比例相差較大，哈明窗的c要比2大。從第一列可以看岀，如果a/b象近的時(shí)候，即 使c稍小于2,也是可以區(qū)分的。每一列中，

8、a/b的比值是相同的，隨著c的減小，譜峰的區(qū)別能力很容易看出是降低的。【自主學(xué)習(xí)內(nèi)容】不同窗函數(shù)的屏幕分辨率。【閱讀文獻(xiàn)】數(shù)字信號(hào)處理【發(fā)現(xiàn)問(wèn)題】（專(zhuān)題研討或相關(guān)知識(shí)點(diǎn)學(xué)習(xí)中發(fā)現(xiàn)的問(wèn)題）：哈明窗所起的作川？哈明窗與矩形窗誰(shuí)的分辨率高？在仿真的時(shí)候c該選用怎樣的值？【問(wèn)題探宄】哈明窗主瓣較寬，分辨率比矩形窗差，但是哈明窗可以使能呈集中在主瓣處。所以哈明窗可以 起到減小泄漏的作用。在仿真的時(shí)候c該選用怎樣的值？在開(kāi)始的時(shí)候選用的c為1248 16,對(duì)比效果不明顯，發(fā)現(xiàn)c的取值過(guò)于大，反s嘗試確定在2附近多取值，其他值也不可取得過(guò)于大，因?yàn)閏稍大實(shí)驗(yàn)現(xiàn)象基本相同，所以沒(méi)有意義?！痉抡娉绦颉縩 = 6

9、4;k = 0:n-1;c_set = 4 3 2.2 2 1.9 1.8 1.5;a_set = 1 1.2 1.6 4;b = 1;wo = 0.5*pi;for c_index = 1:length(c_set) c = c_set(c_index); dw = c*2*pi/n;for a_index = 1:length(a_set)a = a_set(a_index); x = a*cos(wo*k)+b*cos(wo+dw)*k);xh = fft (x.*hamming(n) 1);subplot(length(c_set),length(a_set),length(a_set

10、)*(c_index-l)+a_index); plot(k,abs(xh); hold on;title ( 1c=1 num2str(c);a/b=1 num2str(a);axis ( 12 25 0 25);endend3 己知一離散序列為 x*=cos(/)+0.75cos(/2|z:)，0<*<63 其屮從二0.471，,=+71/64(1) 對(duì)x幻做64點(diǎn)fft,畫(huà)山此時(shí)信號(hào)的頻譜。(2) 如果(1)中顯示的譜不能分辨兩個(gè)譜峰，是否可對(duì)(1)中的64點(diǎn)信號(hào)補(bǔ)零而分辨出兩個(gè)譜 峰。通過(guò)編程進(jìn)行證實(shí)，并解釋其原因。(3) 給出一種能分辨出信號(hào)中兩個(gè)譜峰的計(jì)算方案，并進(jìn)行仿

11、真實(shí)驗(yàn)?！绢}目分析】分析影響譜峰分辨率的主要因數(shù)，進(jìn)一步認(rèn)識(shí)補(bǔ)苓在在頻譜計(jì)算屮的作用。【仿真結(jié)果】 figure 1rue toix v»ewioo«5 uewcxopwindow ne<p> jj ji k c © 3 x -q on=4 l=45.380.390.40410 420 430 44l=256401 0 380 39040 410 420 430 44n=128 l=266恥 256 l=256n=256 l=32768【結(jié)果分析】頻譜的橫坐標(biāo)是歸一化頻率，圖示只觀察要區(qū)別的部分。每一個(gè)列的抽樣點(diǎn)相同，每一行dft 點(diǎn)相同。但看第一列，

12、雖然增加dft點(diǎn)數(shù)，甚至是增加到32768,也不能夠區(qū)分兩個(gè)譜峰，所以要想增 加頻率的分辨率，可以增加抽樣點(diǎn)數(shù)，第二列和第三列抽樣點(diǎn)數(shù)分別是128和256,都可以區(qū)分兩 個(gè)譜峰?！咀灾鲗W(xué)習(xí)內(nèi)容】如何增加頻譜分辨率?！鹃喿x文獻(xiàn)】數(shù)字信號(hào)處理。【發(fā)現(xiàn)問(wèn)題】(專(zhuān)題研討或相關(guān)知識(shí)點(diǎn)學(xué)習(xí)中發(fā)現(xiàn)的問(wèn)題)：如何計(jì)算所需抽樣點(diǎn)數(shù)？可以通過(guò)信號(hào)頻率的差值計(jì)算。【問(wèn)題探究】1、2、3題討論的是離散信號(hào)頻譜的計(jì)算問(wèn)題。與連續(xù)信號(hào)頻譜計(jì)算問(wèn)題相比較，其計(jì)算誤差 有何不同？連續(xù)信號(hào)的頻率是非周期的，離散信號(hào)的頻譜是連續(xù)信號(hào)頻譜的周期話，可能會(huì)有混疊誤差?！痉抡娉绦颉縲0 = 0.4*pi;dw = pi/64;wl =

13、 wo+dw;n_set = 64 128 256;l_set = 64 128 256 512 32768;for n_index = 1:length(n_set);n = n_set(n_index); k = 0:n-l;x = cos (w0*k)+0.75*cos(wl*k); for l_index =1:length(l_set)l = l_set(l_index); if l<ncontinueendx = f ft (x, l); m = (o:l-1)*2/l;subplot(length(l一set),length(n_set),(lindex1*length(n

14、set)+nindex); plot(m,abs(x); axis (0.38 0.44 0 40);title(1n=1 num2str(n) l= num2str(l);endend4試用dft近似計(jì)算島斯信號(hào)以z) = cxp(-6/z2)的頻譜抽樣值。島斯信號(hào)頻譜的理論值為 g( j co)=通過(guò)與理論伉比較，討論信號(hào)的時(shí)域截取長(zhǎng)度和抽樣頻率對(duì)計(jì)算誤差的影響?！绢}目分析】連續(xù)非周期信號(hào)頻譜計(jì)算的基本方法。計(jì)算中出現(xiàn)誤差的主要原因及減小誤差的方法?！痉抡娼Y(jié)果】filt ea vitwtook desktop window h«4p d k %b© a os cutli

15、m»=0 125cutlim«so 0625cimtime=0 25 fsam=16【結(jié)果分析】如圖所示，橫軸坐標(biāo)為角頻率，左斜線上的閣截取時(shí)間相同，垂線上的抽樣頻率相同，很容易能否看出，增加截取時(shí)間和抽樣頻率都可以減小誤差【自主學(xué)習(xí)內(nèi)容】連續(xù)非周期信號(hào)的頻譜計(jì)算方法?！痉抡娉绦颉縟 = pi;fsam_set =1248 16;n_set = 1248;l = 64 ;for fsam_index = 1:length(fsam_set); fsam = fsam_set(fsam_index);tsam = 1/fsam;for nindex = 1:length(ns

16、et);n = n_set(n_index); cuttime = n*tsam; k = -n/2:n/2-l; t = k*tsam; x = exp(-d*t.a2);x = tsam*fftshift(fft(xz l); w = 1inspace(-pi*fsam,pi*fsam,l);subplot(length(n_set),length(fsam_set),fsam_index+length(fsam_set)*(n_index -1);plot(w,abs(x); hold on;g = exp(-w.*w/(4*d); plot(w,g,'r)；title(cutt

17、ime=1 num2str(cuttime) fsam=1 num2str(fsam);endend擴(kuò)展題5本題研究連續(xù)周期信號(hào)頻譜的近似計(jì)算w題。周期為7；的連續(xù)時(shí)間周期信號(hào);v（0可用fourier級(jí)數(shù)表示為zi=oo其屮稱(chēng)力連續(xù)時(shí)間周期信號(hào)a;的頻譜函數(shù)。叫=2兀/7； = 2兀乂稱(chēng)力信號(hào)的基頻（基波），n叫稱(chēng)力 信號(hào)的諧波。如果信號(hào)x(z)函數(shù)表達(dá)式已知，則可由積分得出信號(hào)的頻譜。如果信號(hào)x(0函數(shù)表達(dá)式未知，或 者40闌數(shù)表達(dá)式非常復(fù)雜，則很難巾積分得信號(hào)的頻譜。本題的fi的就是研宂如何利用dft近似 計(jì)算連續(xù)時(shí)間周期信號(hào)的頻譜。(1)若在信號(hào)的一個(gè)周期r0閃抽樣tv個(gè)點(diǎn)，即tq=

18、nt, r為抽樣周期(間隔)，可獲得序列對(duì)%/: = x(f)| z=at; z:二 0,1，tv -1試分析序列對(duì)幻的dft與連續(xù)時(shí)間周期信號(hào)x(0的頻譜汾似的關(guān)系；(2) 由(1)的結(jié)論，給出由dft近似計(jì)算周期信號(hào)頻譜的方案；(3) 周期信號(hào)x的周期7>1，x在區(qū)間0, 1的表達(dá)式為x(0=20r( 1 -z)4cos( 1 2tu)(a) 試畫(huà)出信號(hào)x(/)在區(qū)間0, 1的波形；(b) 若要用6次以a的諧波近似表示xw，試給出計(jì)算方案，并計(jì)算出近似表示的誤差。討論出現(xiàn)誤 差的原因及減小誤差的方法。【理論推導(dǎo)】dft計(jì)算所得結(jié)果xm與連續(xù)周期信號(hào)頻譜x(rtob)的關(guān)系x0= ix

19、cziq,)/=-00l=kt"=-oo2/r由于 = nt為t = -，于是得co.2/rk=nnk今n=m+rn;m=0，l,2,n 1，r為整數(shù)，上式可化為(m+rn)k'v-1<»y4眾=x+ /w)690)-c a，'nt=q /i=oo化簡(jiǎn)可得n-«>刺zx(m + raqty0).e vzw=o /!=«>對(duì)比idft 2/r1 amaw=exm，e，amkm=0xmoo可得x (m + zw)6y0); m = 0,1，2,/v -1"=-oo【計(jì)算方案】根據(jù)理論推導(dǎo)結(jié)果設(shè)計(jì)近似計(jì)算方案。分析產(chǎn)

20、生誤差的主要原因。dft是連續(xù)信號(hào)頻譜的周期化，如果頻譜信號(hào)高頻信號(hào)較多，就會(huì)產(chǎn)生混疊，所以有誤差。如 果主要是低頻信號(hào)，周期話后互補(bǔ)產(chǎn)生影響，就沒(méi)有誤差?！緮U(kuò)展分析】如果周期信號(hào)x（0是帶限信號(hào)，即信號(hào)的最高頻率分量為m吻（是正整數(shù)），試確定在一個(gè)周期閃 的最少抽樣點(diǎn)m使得在頻譜的計(jì)算過(guò)程當(dāng)屮不存在混疊誤差。與抽樣定理給出的結(jié)論比較，發(fā)表 你的看法。很容易看出，要想不產(chǎn)生混疊，必須7v22m+1。該結(jié)論與抽樣定理相似。因?yàn)橹芷谛盘?hào)是一種特殊的連續(xù)信號(hào)，所以抽樣定理對(duì)周期信號(hào)也是適用的?！痉抡娼Y(jié)果】n=16k=6diviation:mean:0.07187 max:0.1339diviatio

21、n:n=12k=6diviation:mean:0.1564 max:0.3297diviation:【結(jié)果分析】在n=16，k=6時(shí)，誤差最小，最大誤差力0.1229，平均誤差為0.07187。無(wú)論是減小n和減小k都很大的增大了誤差。為了分別測(cè)試n和k對(duì)誤差的影響，分別設(shè)k=6,求不同n吋候的誤差，n=16,不同k對(duì)誤差 的影響，分別得到了下圖（在原程序上稍加改動(dòng)，所以在仿真程序中不再另寫(xiě)了）在k=6的時(shí)候，誤差有明顯的減小?？梢?jiàn)原函數(shù)最大的高頻分量為6*0,根據(jù)前面推出的理論，n取6*2+1=13點(diǎn)的吋侯誤差減小。 實(shí)驗(yàn)結(jié)果符合結(jié)論 【自主學(xué)習(xí)內(nèi)容】周期信號(hào)的頻譜與dft的關(guān)系。【發(fā)現(xiàn)問(wèn)題】原函數(shù)最高頻分量大約是多少？【問(wèn)題探究】根據(jù)前面推理，最高頻分量為6*w0。【仿真程序】n = 1024;t = 1inspace (0,1,n);n_set = 12 16;h_set = 4 6;to = 1;count = 0;for n_index = 1:length(n_set)n = n_set(n_index); dt = to