《MATLAB仿真 實(shí)現(xiàn)經(jīng)典譜估計(jì)(采用周期圖法)》由會(huì)員分享，可在線(xiàn)閱讀，更多相關(guān)《MATLAB仿真 實(shí)現(xiàn)經(jīng)典譜估計(jì)(采用周期圖法)（4頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、數(shù)字信號(hào)處理 課程實(shí)驗(yàn)報(bào)告
實(shí)驗(yàn)指導(dǎo)教師：黃啟宏
實(shí)驗(yàn)名稱(chēng)
MATLAB仿真實(shí)現(xiàn)經(jīng)典譜估計(jì)（采用周期圖法）
專(zhuān)業(yè)、班級(jí)
電子與通信工程
姓 名
張帥
實(shí)驗(yàn)地點(diǎn)
仿古樓301
實(shí)驗(yàn)日期
2013.11.17
1、 實(shí)驗(yàn)內(nèi)容
采用周期圖法（直接法）實(shí)現(xiàn)經(jīng)典譜估計(jì)。
二、實(shí)驗(yàn)?zāi)康?
（1）掌握周期圖法（直接法）估計(jì)出功率譜的步驟和方法；
（2）在實(shí)驗(yàn)的過(guò)程中找到影響經(jīng)典譜估計(jì)的因素；
（3）了解周期圖法（直接法）估計(jì)功率譜的缺陷。
三、實(shí)驗(yàn)原理
把隨機(jī)信號(hào)的N點(diǎn)觀察數(shù)據(jù)視為一能量有限信號(hào)，直接取得傅里葉變換，得，然后再取其幅值的平方，并除
2、以N，作為對(duì)真實(shí)的功率譜的估計(jì)。即為：
四、涉及實(shí)驗(yàn)的相關(guān)情況介紹（包含使用軟件或?qū)嶒?yàn)設(shè)備等情況）
一臺(tái)安裝MATLAB軟件的電腦
5、 實(shí)驗(yàn)記錄（程序、相關(guān)的圖形、相關(guān)數(shù)據(jù)記錄及分析）
%采用直接法（周期圖法）估計(jì)功率譜；
clear
Fs = 1000;%采樣頻率
n = 0:1 /Fs: .3;%產(chǎn)生含有噪聲的序列
xn = cos(200*pi*n)+0.1*randn(size(n));
subplot(311);%輸出隨機(jī)信號(hào)xn；
plot(n,xn);xlabel('時(shí)間');ylabel('幅度');title('輸入信號(hào)x（n）');
ax
3、is([0 0.3 -2 2]);
grid on;
window = boxcar( length( xn) ) ;%矩形窗
nfft = 512;
[Pxx f]= periodogram( xn,window,nfft,Fs) ;%直接法
subplot(312)
plot( f,10* log10( Pxx) ) ;
title('直接法經(jīng)典譜估計(jì),512點(diǎn)');
xlabel('頻率（Hz）');
ylabel('功率譜密度');
grid on;
window = boxcar( length( xn) );
nfft = 1024;
[Pxx,f]= periodogram( xn,window,nfft,Fs) ;
subplot(313) ;
plot( f,10*log10( Pxx) ) ;
title('直接法經(jīng)典譜估計(jì),1024點(diǎn)');
xlabel('頻率（Hz）');
ylabel('功率譜密度');
grid on;
6、 實(shí)驗(yàn)總結(jié)
從上圖我們可以得到這樣的結(jié)論：在增加數(shù)據(jù)長(zhǎng)度N時(shí)，就會(huì)使互不相關(guān)的點(diǎn)數(shù)增加，提高譜曲線(xiàn)的分辨力，但是加劇譜曲線(xiàn)的起伏。經(jīng)典功率譜估計(jì)不是一致估計(jì)，這是周期圖法（直接法）的一個(gè)嚴(yán)重的缺點(diǎn)。
MATLAB仿真 實(shí)現(xiàn)經(jīng)典譜估計(jì)(采用周期圖法)
