信号处理与数据分析第10章作业题目与答案.pdfVIP

1.(书稿 10.20)一段数据包含N 点采样，其抽样频率为f s 1000Hz 。用平均法改进周期图法 估计时将数据分成了互不交叠的K 段，每段数据长度M N / K 。假定在频谱中有两个相距 为0.04π(rad) 的谱峰，为要分辨它们，M 应该取多大？ 解答：谱峰间隔为0.04π(rad) ，即为20Hz，也就是要求频率分辨率小于20Hz，根据傅里叶 变换可知，其频率分辨率为F0 fs / M ，所以M 50 。 2.(编程题) 已知AR(4)模型为： x(n) 2.7607x(n 1) 3.8106x(n 2) 2.6535x(n 3) 0.9238x(n 4) w(n) 其中，w(n) 是零均值、单位方差的白噪声，用此模型分别产生256 点和1024 点数据，然后 分别利用周期图法，Welch 方法以及AR 参数模型方法计算功率谱。将不同数据点下的功率 谱进行画图比较。 解答：MATLAB 参考程序 （见FTP 中“Power_spectrum.m”）如下： %%功率谱估计 %% 以数据点数等于1024 为例，类似地更改N 值大小即可 clear all; close all; clc; N=1024;%数据点数 m=4;%AR 模型阶数 wn=randn(1,N);%零均值、单位方差白噪声 a1=-2.7607;a2=3.8106;a3=-2.6535;a4=0.9238;%系数 xn=zeros(1,N); xn(1)=wn(1); xn(2)=wn(2)-a1*xn(1); xn(3)=wn(3)-a1*xn(2)-a2*xn(1); xn(4)=wn(4)-a1*xn(3)-a2*xn(2)-a3*xn(1); for n=5:N xn(n)=wn(n)-a1*xn(n-1)-a2*xn(n-2)-a3*xn(n-3)-a4*xn(n-4); end [Pxx1,W1]=periodogram(xn);%周期图法 figure; subplot(3,1,1); plot(W1/max(W1),10*log10(Pxx1));grid on; xlabel(归一化频率); ylabel(功率谱密度/dB); title(周期图法功率谱估计（N=1024)); axis([0 1 -30 50]); set(gca,XTick,0:0.2:1); set(gca,YTick,-30:10:50); [Pxx2,W2]=pwelch(xn);%Welch 法 subplot(3,1,2); plot(W2/max(W2),10*log10(Pxx2));grid on; xlabel(频率); ylabel(功率谱密度/dB); title(Welch 法功率谱估计（N=1024)); axis([0 1 -30 50]); set(gca,XTick,0:0.2:1); set(gca,YTick,-30:10:50); [Pxx3,W3]=pyulear(xn,4);%AR 模型法 subplot(3,1,3); plot(W3/max(W3),10*log10(Pxx3));grid on; xlabel(归一化频率); ylabel(功率谱密度/dB); title(AR 模型法功率谱估计（N=1024)); axis([0 1 -30 50]); set(gca,XTick,0:0.2:1); set(gca,YTick,-30:10:50); 数据点数为1024 下的功率谱比较如下图所示： 周期图法功率谱估计（ N=1024) 50 40 B 30 d / 度 20 密 10 谱 0 率 功 -1