P(w)二—X(ejw)peFN―其中 X(ejw)=SxNN(n)e-jwn现代信号处理作业实验题目:设信号 x(n)二 cos(0
35 兀 n+0)+cos(0
47 兀 n+0)+0
25cos(0
8 兀 n+0)+v(n),123其中 0,0,0 是[—兀,兀]内的独立随机变量,v(n)是单位高斯白噪声
利用周期图法对序列进行功率谱估计
数据窗采用汉明窗
利用 BT 法对序列进行功率谱估计,自相关函数的最大相关长度为 M=64,128,256,512 采用 BARTLETT 窗
利用 Welch 法对序列进行功率谱估计,50%重叠,采用汉明窗,L=256,128,64
利用 Burg 法对序列进行 AR 模型功率谱估计,阶数分别为 10,13
要求每个实验都取 1024 个点,fft 作为谱估计,取 50 个样本序列的算术平均,画出平均的功率谱图
实验原理:1)
周期图法:又称间接法,它把随机信号的 N 个观察值 xN(n)直接进行傅里叶变换,得到 XJejw),然后取其幅值的平方,再除以 N,作为对 x(n)真实功率谱的估计
BT 法:对于 N 个观察值 x(0),x(1),
,x(N-1),令 xN(n)=a(n)x(n)
计算 r(m)为xAiP1MU Xi(W)2M其中 Xi(w)=EXi(n)e-jwn
由此得到平均周期图为MMn=0r(m)=—NEX(n)x(n+m),|m|