PowerSpectrumEstimation上機實驗設(shè)輸入音頻信號，取f=1KHz，fs=20KHz，N=128，W(n)為三角窗序列。用計算機求出功率譜估值及不分段的，打印出曲線圖、列出程序。測量流程圖為：取樣取樣fs=20KHz三角窗加權(quán)規(guī)范運算N/2點ODFTR2+I2除NU或除N補奇數(shù)譜線內(nèi)存UK2Xa(t)K1或D(2k)注：K1斷開，K2合上得；K1合上，K2斷開得 用計算機求及時的編程流程圖如下（供參考）。開始開始輸入N及Ts=1/fs對音視頻信號取N=128個樣點為X(k)三角窗加權(quán)求：規(guī)范運算取實部Y(N)，虛部Z(N)N/2點ODFT，得偶數(shù)譜線，其實部為A(k)，虛部B(k)求模的平方R2+I2除N由ODFT譜線的共軛奇偶對稱補奇數(shù)譜線結(jié)束規(guī)范運算取實部Y1(N)，虛部Z1(N)N/2點ODFT，得偶數(shù)譜線，其實部為C(k)，虛部D(k)求模的平方R2+I2除NU由ODFT譜線的共軛奇偶對稱補奇數(shù)譜線1.用MATLAB軟件求功率譜估值及不分段的，程序及注釋分別如下：1.1功率譜估值程序clearall;%清內(nèi)存clc;%清屏N=128;%采樣點數(shù)為128f=1000;%原始音頻信號頻率為1000HZfs=20000;%采樣頻率為20000HZt=0:1/fs:N/fs;%采樣時間的范圍Xt=cos(2*pi*f*t);%生成原始音頻采樣信號subplot(3,1,1);stem(t,Xt);xlabel('時間t(s)');ylabel('幅度');title('原始音頻信號');axis([0,N/fs,-2,2]);form=1:N/2u(m)=(Xt(m)-j*Xt(m+N/2))*exp(-j*(m-1)*pi/N);%生成Xt(n)的規(guī)范運算u(n)endW=fft(u,N/2);%通過規(guī)范運算u(n)的N/2FFT運算來得到ODFT的偶數(shù)譜線k=0:N/2-1;%取一個和后面離散20行W1點數(shù)一樣的自變量subplot(3,1,2);W1=abs(W).^2/N;%按題中流程圖——先求模的平方，再除Nstem(2*k,W1);title('偶數(shù)譜線');xlabel('自變量');ylabel('幅度');axis([0,140,0,40]);k=k+1;%MATLAB的下標(biāo)不能從0開始WW(2*k)=W(k);%WW(n)是一個總序列，該步是將偶W(n)譜加入該總序列WW(N+1-2*k)=conj(W);%求偶譜線的共軛為奇譜，根據(jù)安老師課件公式（52）得此，補奇譜線，補全總序列WW(n)W3=abs(WW).^2/N;%按題中流程圖—先求模的平方，再除Nsubplot(3,1,3);L=0:N-1;%取一個和后面34行離散W1點數(shù)一樣的自變量stem(L,W3);title('功率譜估值');xlabel('自變量');ylabel('幅度');1.2不分段功率譜估值程序clearall;%清內(nèi)存clc;%清屏f=1000;%原始音頻信號頻率為1000HZfs=20000;%采樣頻率為20000HZN=128;%采樣點數(shù)為128t=0:1/fs:(N-1)/fs;%采樣時間的范圍Xt=cos(2*pi*f*t);%生成原始音頻采樣信號subplot(4,1,1);stem(t,Xt);%生成原始音頻采樣信號的圖title('原始音頻信號');axis([0,(N-1)/fs,-2,2]);xlabel('時間t(s)');ylabel('幅度xt');wnd=triang(N);%生成點數(shù)為N的窗函數(shù)U=sum(wnd.*wnd);%按題中所給取U為三角窗函數(shù)的模平方序列Xt_wnd=Xt.*(wnd');%生成原始音頻采樣信號通過三角窗函數(shù)后的Xt(n)subplot(4,1,2);stem(t,Xt_wnd);%生成原始音頻采樣信號通過三角窗函數(shù)后的Xt(n)-wnd的圖title('加三角窗后的音頻信號');axis([0,(N-1)/fs,-2,2]);xlabel('時間t(s)');ylabel('幅度xt_wnd');form=1:N/2u(m)=(Xt_wnd(m)-j*Xt_wnd(m+N/2))*exp(-j*(m-1)*pi/N);%生成Xt(n)-wnd的規(guī)范運算u(n)endk=0:N/2-1;%取一個和78行離散W1點數(shù)一樣的自變量W=fft(u,N/2);%快速FFT得到偶數(shù)譜線W1=abs(W).^2/N/U;%按題中流程圖——先求模的平方，再除NUsubplot(4,1,3);stem(2*k,W1);%生成偶譜線W1的圖title('偶數(shù)譜');xlabel('自變量');ylabel('幅度');k=k+1;%MATLAB的下標(biāo)不能從0開始WW(2*k)=W(k);%WW(n)是一個總序列，該步是將偶W(n)譜加入該總序列WW(N+1-2*k)=conj(W);%求偶譜線的共軛為奇譜，根據(jù)安老師課件公式（52）得此，補奇譜線，補全總序列WW(n)W3=abs(WW).^2/N/U;%按題中流程圖——先求模的平方，再除NUsubplot(4,1,4);L=0:N-1;%取一個和前面87行離散W3點數(shù)一樣的自變量