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

Power Spectrum Estimation上机实验 设输入音频信号，取f=1KHz，fs=20KHz，N=128，W(n)为三角窗序列。用计算机求出功率谱估值及不分段的，打印出曲线图、列出程序。 测量流程图为： 取样 取样 fs=20KHz 三角窗加权 规范运算 N/2点ODFT R2+I2 除NU或除N 补奇数谱线 内存U K2 Xa(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 功率谱估值程序 clear all; %清内存 clc; %清屏 N=128; %采样点数为128 f=1000; %原始音频信号频率为1000HZ fs=20000; %采样频率为20000HZ t=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]); for m=1:N/2 u(m)=(Xt(m)-j*Xt(m+N/2))*exp(-j*(m-1)*pi/N); %生成Xt(n)的规范运算u(n) end W=fft(u,N/2); %通过规范运算u(n)的N/2 FFT运算来得到ODFT的偶数谱线 k=0:N/2-1; %取一个和后面离散20行W1点数一样的自变量 subplot(3,1,2); W1=abs(W).^2/N; %按题中流程图——先求模的平方，再除N stem(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; %按题中流程图—先求模的平方，再除N subplot(3,1,3); L=0:N-1; %取一个和后面34行离散W1点数一样的自变量 stem(L,W3); title('功率谱估值'); xlabel('自变量'); ylabel('幅度'); 1.2 不分段功率谱估值程序 clear all ; %清内存 clc ; %清屏 f=1000 ; %原始音频信号频率为1000HZ fs=20000 ; %采样频率为20000HZ N=128 ; %采样点数为128 t=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的窗