88段功率谱估值的流程图 MATLAB程序 仿真结果

PowerSpectrumEstimation上机实验设输入音频信号，取f=1KHz，fs=20KHz，N=128，W(n)为三角窗序列。用计算机求出功率谱估值及不分段的，打印出曲线图、列出程序。测量流程图为：取样取样fs=20KHz三角窗加权规范运算N/2点ODFTR2+I2除NU或除N补奇数谱线内存UK2Xa(t)K1或D(2k)注：K1断开，K2合上得；K1合上，K2断开得 用计算机求及时的编程流程图如下（供参考）。开始开始输入N及Ts=1/fs对音视频信号取N=128个样点为X(k)三角窗加权求：规范运算取实部Y(N)，虚部Z(N)N/2点ODFT，得偶数谱线，其实部为A(k)，虚部B(k)求模的平方R2+I2除N由ODFT谱线的共轭奇偶对称补奇数谱线结束规范运算取实部Y1(N)，虚部Z1(N)N/2点ODFT，得偶数谱线，其实部为C(k)，虚部D(k)求模的平方R2+I2除NU由ODFT谱线的共轭奇偶对称补奇数谱线1.用MATLAB软件求功率谱估值及不分段的，程序及注释分别如下：1.1功率谱估值程序clearall;%清内存clc;%清屏N=128;%采样点数为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)的规范运算u(n)endW=fft(u,N/2);%通过规范运算u(n)的N/2FFT运算来得到ODFT的偶数谱线k=0:N/2-1;%取一个和后面离散20行W1点数一样的自变量subplot(3,1,2);W1=abs(W).^2/N;%按题中流程图——先求模的平方，再除Nstem(2*k,W1);title('偶数谱线');xlabel('自变量');ylabel('幅度');axis([0,140,0,40]);k=k+1;%MATLAB的下标不能从0开始WW(2*k)=W(k);%WW(n)是一个总序列，该步是将偶W(n)谱加入该总序列WW(N+1-2*k)=conj(W);%求偶谱线的共轭为奇谱，根据安老师课件公式（52）得此，补奇谱线，补全总序列WW(n)W3=abs(WW).^2/N;%按题中流程图—先求模的平方，再除Nsubplot(3,1,3);L=0:N-1;%取一个和后面34行离散W1点数一样的自变量stem(L,W3);title('功率谱估值');xlabel('自变量');ylabel('幅度');1.2不分段功率谱估值程序clearall;%清内存clc;%清屏f=1000;%原始音频信号频率为1000HZfs=20000;%采样频率为20000HZN=128;%采样点数为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);%生成点数为N的窗函数U=sum(wnd.*wnd);%按题中所给取U为三角窗函数的模平方序列Xt_wnd=Xt.*(wnd');%生成原始音频采样信号通过三角窗函数后的Xt(n)subplot(4,1,2);stem(t,Xt_wnd);%生成原始音频采样信号通过三角窗函数后的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的规范运算u(n)endk=0:N/2-1;%取一个和78行离散W1点数一样的自变量W=fft(u,N/2);%快速FFT得到偶数谱线W1=abs(W).^2/N/U;%按题中流程图——先求模的平方，再除NUsubplot(4,1,3);stem(2*k,W1);%生成偶谱线W1的图title('偶数谱');xlabel('自变量');ylabel('幅度');k=k+1;%MATLAB的下标不能从0开始WW(2*k)=W(k);%WW(n)是一个总序列，该步是将偶W(n)谱加入该总序列WW(N+1-2*k)=conj(W);%求偶谱线的共轭为奇谱，根据安老师课件公式（52）得此，补奇谱线，补全总序列WW(n)W3=abs(WW).^2/N/U;%按题中流程图——先求模的平方，再除NUsubplot(4,1,4);L=0:N-1;%取一个和前面87行离散W3点数一样的自变量