基于matlab的功率譜估計(jì).doc

《基于matlab的功率譜估計(jì).doc》由會(huì)員上傳分享，免費(fèi)在線(xiàn)閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。

1、題目基于matlab的功率譜估計(jì)學(xué)院通信工程學(xué)院專(zhuān)業(yè)通信與信息系統(tǒng)第一章功率譜估計(jì)分析及比較1.實(shí)驗(yàn)?zāi)康?1)掌握Welch算法的概念、應(yīng)用及特點(diǎn)；(2)了解譜估計(jì)在信號(hào)分析中的作用；(3)能夠利用Welch法對(duì)信號(hào)作譜估計(jì)，對(duì)信號(hào)的特點(diǎn)加以分析。2.實(shí)驗(yàn)內(nèi)容(1)讀入實(shí)驗(yàn)數(shù)據(jù)。(2)編寫(xiě)一利用Welch法作估計(jì)的算法程序。(3)將計(jì)算結(jié)果表示成圖形的形式，給出信號(hào)譜的分布情況圖。3.譜估計(jì)方法簡(jiǎn)介(1)周期圖法周期圖法是直接建立在功率譜的定義式上的，也稱(chēng)之為直接法。原理計(jì)算如下：a)取N點(diǎn)數(shù)據(jù)的DTFT（DFT）；b)求模之平方并除以N；(2)自相關(guān)法自相關(guān)法的原理是由維納-

2、辛欽公式，經(jīng)自相關(guān)函數(shù)間接獲得的。原理計(jì)算如下：a)x(n)，N點(diǎn)T2N-1點(diǎn)，得Rxx(m)；b)按2N-1點(diǎn)對(duì)Rxx(m)；作DFT，Rxx(m)能推出SBT(k)；(3)加窗平滑法（BT法）加窗平滑法的原理是先做自相關(guān)估計(jì)，在選擇合適的窗函數(shù)相乘，也即截?cái)?，然后作DFT。(4)平均周期圖法（Bartlett法）平均周期圖法的原理是個(gè)獨(dú)立同分布的隨機(jī)變量的均值之方差，等于單個(gè)變量方差的1/k。具體的方法步驟是長(zhǎng)數(shù)據(jù)N分成k段，每段M=N/k，針對(duì)每段分別用周期圖法求譜，然后k段平均后求的的就是譜估計(jì)。ww(5)韋爾奇（Welch）譜估計(jì)法韋爾奇譜估計(jì)法將加窗平滑法和平均周期

3、法二者相結(jié)合，其原理在于：a)把N個(gè)數(shù)據(jù)分成k段，每段可以互相獨(dú)立(如平均周期圖法)a)再把每段數(shù)據(jù)乘上窗函數(shù)w(n)(如加窗平滑法)后作DFT。第二部分仿真結(jié)果圖未加窗時(shí)周期圖：取平均后的圖形為：加窗周期圖法：加矩形窗：取平均后的圖形為：加漢明窗后的圖形為：取平均后的圖形為：下圖為采樣頻率為100，窗函數(shù)分別取矩形窗、漢寧窗和漢明窗時(shí)的圖形：未加窗時(shí)自相關(guān)法：取平均后的圖形為：B-T法：加漢明窗時(shí)：取平均后的圖形為：welch法：L=512,加漢明窗，50%重疊取平均后的圖形為：L=512,加矩形窗，50%重疊取平均后的圖形為：L=256,加漢明窗，50%重疊取平均后的圖形為

4、：L=256,加矩形窗，50%重疊取平均后的圖形為：L=128,加漢明窗，50%重疊取平均后的圖形為：Bartlett法：L=512，加漢明窗取平均后的圖形為：L=256，加漢明窗取平均后的圖形為：L=128,加漢明窗取平均后的圖形為：第三章結(jié)果分析（1）周期圖法周期圖法所求得的功率譜振蕩劇烈，信號(hào)方差較大，不利于對(duì)功率信號(hào)的分析。缺少了統(tǒng)計(jì)平均時(shí)，記錄的信號(hào)序列長(zhǎng)度一定的條件下，要保證足夠高的譜分辨率，譜估計(jì)的方差就會(huì)很大，譜的正確性會(huì)很差。當(dāng)數(shù)據(jù)長(zhǎng)度N太大時(shí)，譜曲線(xiàn)呈現(xiàn)較大的起伏；當(dāng)數(shù)據(jù)長(zhǎng)度N太小時(shí)，譜的分辨率又不好。(2)自相關(guān)法先根據(jù)實(shí)驗(yàn)所給數(shù)據(jù)求解出自相關(guān)函數(shù)，然后對(duì)

5、自相關(guān)函數(shù)進(jìn)行傅里葉變換，從而得到功率譜估計(jì)。從圖中可以看出，自相關(guān)法得到的結(jié)果與周期圖法相似，同樣是方差較大，信號(hào)振蕩大。(3)加窗平滑法（BT）根據(jù)加窗平滑法的求解步驟進(jìn)行編程，取窗函數(shù)為矩形窗。但是可以看出在L取較大時(shí)性能不是很好。(4)平均周期圖法（Bartlett）將數(shù)據(jù)平均分為K段，Bartlett法很好地改善了直接法的方差特性，但是它是以犧牲偏差和分辨率為代價(jià)的。(5)Welch法Welch法各段允許交疊，從而增大了段數(shù)L，這樣可以更好地改善方差特性。但是，數(shù)據(jù)的交疊又減小了每一段的不相關(guān)性，使方差的減小不會(huì)達(dá)到理論計(jì)算的程度。另外，選擇合適的窗函數(shù)可以減小頻譜的

6、泄漏，改善分辨率。另外，由矩形窗處理的譜估計(jì)的主瓣寬度最窄，分辨率最好，但是其旁瓣比其他窗函數(shù)的旁瓣要高，因此其正弦譜線(xiàn)附近的旁瓣泄漏比較嚴(yán)重，而且其起伏性較大，所以其方差特性最差。由海明窗處理譜估計(jì)的主瓣寬度最寬，因此其分辨率相對(duì)較差，但其旁瓣較小，大大改善了由矩形窗處理的譜估計(jì)旁瓣較大所產(chǎn)生的譜失真。究其原因，選擇不同的窗函數(shù)其主瓣寬度不一樣，造成譜估計(jì)的分辨率也不相同；另外，選擇不同的窗函數(shù)旁瓣的衰減速度也不相同，因而譜估計(jì)旁瓣的泄漏程度也不一樣。附錄：未加窗時(shí)周期圖：clear;Fs=1000;n=0:1/Fs:1;nfft=1024;x=cos(100*pi*n)+c

7、os(200*pi*n)+0.25*cos(400*pi*n)+randn(size(n));y=fft(x);Rx=abs(y).^2/length(n);a=length(Rx);sum=zeros(1,a);fori=0:49x=cos(100*pi*n)+cos(200*pi*n)+0.25*cos(400*pi*n)+randn(size(n));y=fft(x);Rx=abs(y).^2/length(n);sum=sum+Rx;endxpsd=10*log10(sum);plo