資源描述：
《88段功率譜估值的流程圖matlab程序仿真結(jié)果》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、PowerSpectrumEstimation上機實驗設(shè)輸入音頻信號X“（r）=cos（2/）,取f=KHz,Q20KHZ,N=128,W何為三角窗序列。用計算機求出功率譜估值尺,伙）及不分段的3；伙），打印出曲線圖、列出程序。測量流程圖為：注：K1斷開，K2合上得硝伙）；K1合上，K2斷開得耳伙）用計算機求丘伙）及昭伙）時的編程流程圖如下（供參考）。1?用MATLAB軟件求功率譜估值丘（幻及不分段的3：：伙）,程序及注釋分別如下:stem(2*k,Wl);titleC偶數(shù)譜線J;xlabelf自變量')；ylabel(*幅度J;axis(LO,140,0,40]);k=k+l;
2、WW(2*k)=W(k);WW(N+1-2*k)=conj(W);W3=abs(WW).A2/N;1.1功率譜估值程序clearall;%清內(nèi)存clc;%清屏N=128;%采樣點數(shù)為128f=1000;%原始咅頻信號頻率為I000HZfs=20000;%采樣頻率為20000HZt=0:l/fs:N/fs;%采樣時間的范圍Xt=cos(2*pi*f*t);subplot(3,l,l);stem(t,Xt);xlabel('時間t(s)');ylabel(幗度')；titleC原始音頻信號axis([0,N/fs,-2,2]);form=l:N/2%生成原始音頻采樣信號u(m)=(Xt
3、(m)-j*Xt(m+N/2))*exp(-j*(m-l)*pi/N);endW=fft(u,N/2);k=0:N/2-l;subplot(3,l,2);Wl=abs(W).A2/N;%生成Xt(n)的規(guī)范運算u(n)%通過規(guī)范運算u(n)的N/2FFT運算來得到ODFT的偶數(shù)譜線%取一個和后面離散20行W1點數(shù)一樣的自變量%按題中流程圖一一先求模的平方，再除N%MATLAB的下標不能從0開始%WW(n)是一個總序列，該步是將偶W(n)譜加入該總序列%求偶譜線的共轆為奇譜，根據(jù)安老師課件公式(52)得此，補奇譜線，補全總序列WW(n)%按題中流程圖一先求模的平方，再除Nsubplo
4、t(3,l,3)；L=0:N-l;%取一個和后面34行離散W1點數(shù)一樣的自變量stem(L,W3);title。功率譜估值xlabelC自變量')；ylabel('幅度);clearall;1.2不分段功率譜估值程序%清內(nèi)存clc;f=1000;fs=20000;%清屏%原始咅頻信號頻率為1000HZ%采樣頻率為20000HZN=128;%釆樣點數(shù)為128t=O:l/fs:(N-l)/fs;Xt=cos(2*pi*f*t);subplot(4,l,l);stem(t,Xt);%采樣時間的范圍%生成原始音頻采樣信號%生成原始音頻采樣信號的圖title(原始音頻信號');axis([0
5、,(N-l)/fs,-2,2]);xlabel(時間t(s)');ylabel('幅度xt');wnd=triang(N);U=sum(wnd.*wnd);%生成點數(shù)為N的窗函數(shù)%按題中所給取U為三角窗函數(shù)的模平方序列Xt_wnd=Xt.*(wnd');subplot(4,l,2);stem(t,Xt_wnd);%生成原始音頻采樣信號通過三角窗函數(shù)后的Xt(n)%牛成原始音頻采樣信號通過三角窗函數(shù)后的Xt(n)-wnd的圖title(加三角窗后的音頻信號J;axis([0,(N-l)/fs,-2,2]);xlabelC時間(⑸)；ylabel('幅度xt_wnd);form=l:N
6、/2u(m)=(Xt_wnd(m)-j*Xt_wnd(m+N/2))*exp(-j*(m-1)*pi/N);endk=0:N/2-l;W=fft(u,N/2);W1二abs(W).人2/N/U;subplot(4,l,3)；stem(2*k,Wl);titleC偶數(shù)譜J;xlabelC自變量')；%生成Xt(n)-wnd的規(guī)范運算u(n)%取一個和78行離散W1點數(shù)一樣的自變量%快速FFT得到偶數(shù)譜線%按題中流程圖一一先求模的平方，再除NU%生成偶譜線W1的圖ylabelC幅度);k=k+l;WW(2*k)=W(k);WW(N+1-2*k)=conj(W);%MATLAB的下標不能
7、從0開始%WW(n)是一個總序列，該步是將偶W(n)譜加入該總序列%求偶譜線的共轆為奇譜，根據(jù)安老師課件公式(52)得此，補奇譜線,W3=abs(WW).A2/N/U;補全總序列WW(n)%按題中流程圖——先求模的平方，再除NUsubplot(4,l,4);L=O:N-1;%取一個和前面87行離散W3點數(shù)一樣的自變量stem(L,W3);title。功率譜估值J;xlabeK'g變量')；ylabel('幅度);2?程序?qū)崿F(xiàn)的曲線圖下圖1(a),(b)所示分別為功率譜
