《古典法功率谱估计》由会员分享,可在线阅读,更多相关《古典法功率谱估计(19页珍藏版)》请在金锄头文库上搜索。
1、古典法功率谱估计一、 信号的产生(一)信号构成在本实验中,需要事先产生待估计的信号,为了使实验成果较为明显,我产生了由两个不同频率的正弦信号(频率差相对较大)和加性高斯白噪声构成的信号。(二)程序xn=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN信号figure(1)plot(n,xn);title(a) 两个正弦信号与白噪声叠加的时域波形)(三) 信号波形二、有关法功率谱估计(一)算法原理简介此措施以有关函数为媒介来计算功率谱,因此又叫间接法。它是1958年由Blackman 和Tukey提出。这
2、种措施的具体环节是: 第一步:从无限长随机序列x(n)中截取长度N的有限长序列 第二步:由N长序列求(2M-1)点的自有关函数序列。 第三步:由有关函数的傅式变换求功率谱。 以上过程中经历了两次截断,一次是将x(n)截成N长,称为加数据窗,一次是将想x(n)截成(2M-1)长,称为加延迟窗。因此所得的功率谱仅是近似值,也叫谱估计。一般取MN,由于只有当M较小时,序列傅式变换的点数才较小,功率谱的计算量才不至于大到难以实现,并且谱估计质量也较好。因此,在FFT问世之前,有关法是最常用的谱估计措施。当FFT问世后,状况有所变化。由于截断后的)(nxN可视作能量信号,由有关卷积定理可得这就将有关化为
3、线性卷积,而线性卷积又可以用迅速卷积来实现。我们可对上式两边取(2N-1)点DFT,则有于是将时域卷积变为频域乘积,用迅速有关求自有关函数估值的完整方案如下: 1. 对N长)(nxN的充(N-1)个零,成为(2N-1)长的。 2. 求(2N-1)点的FFT,得3、求 4、 求(2N-1)点的IFFT(二)运算简要框图X(n)迅速有关加窗截断(2M-1)点FFT输出矩形窗截断有关法谱估计运算简要框图图中迅速有关的输出时从-(N-1)到(N-1)的2N-1点,加窗后截取的是-(M-1)到(M-1)的 ,最后做(2M-1)点FFT,即可得到成果。(三)程序示例程序的重要思路就是按照运算框图一步一步进
4、行计算,下面附程序并进行简要解释:N=512,n=0:N-1; %N是FFT的变换区间xn=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN的信号Xk=fft(xn,1024); %进行2N-1点FFT,系统会自动补0Sk=abs(Xk).*(abs(Xk)./N; %取频谱幅度的平方,并除以N,以此作为对xn真实功率谱的估计Rn=ifft(Sk);Sk1=fft(Rn,512);figure(2)subplot(2,1,1);plot(n/N,Sk1);ylabel(Sk)title(b) 有关法估计
5、功率谱密度)Sk2=10*log(Sk1); %对估计出的Sk取对数,使画出的图更加突出特点subplot(2,1,2);plot(n/N,Sk2);ylabel(10log(PSD)(四)成果分析 下面是程序运营后的成果从上图中我们可以较为明显的看到信号中有两个频率分量,一种在0.2处,一种在0.4处,与产生的信号相一致。但是我们不难看出,估计出的功率谱谱线非常不平坦,有诸多起伏。三、 周期图法谱估计(一)算法原理简介周期图法又称直接法。它是从随机信号x(n)中截取N长的一段,把它视为能量有限x(n)真实功率谱的估计的抽样。其具体环节如下:第一步:由获得的N点数据构成有限长序列直接求傅里叶变
6、换,得频谱。第二步:取频谱幅度的平方,并除以N,以此作为对x(n)真实功率谱的估计。事实上,周期图法谱估计与自有关法谱估计的差别只是估计自有关函数的措施不同。(二)运算简要框图 矩形窗(长度N)截断N点FFTX(n) 图中用FFT来替代傅里叶变换(三)程序示例N=512,n=0:N-1;xn=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN的信号Xk1=fft(xn,512); %进行N点FFTSk3=abs(Xk1).*(abs(Xk1)./N;%取频谱幅度的平方,并除以N,以此作为对xn真实功率谱的
7、估计Sk4=10*log(Sk3); %对估计出的Sk取对数,使画出的图更加突出特点figure(3)subplot(2,1,1);plot(n/N,Sk3);title(c) 周期图法估计功率谱密度)ylabel(Sk)subplot(2,1,2);plot(n/N,Sk4);ylabel(10log(PSD)(四)成果分析下面是程序运营后的成果从上图中我们同样可以看到信号中有两个频率分量,一种在0.2处,一种在0.4处,与产生的信号相一致。但是我们不难看出,估计出的功率谱谱线与有关法功率谱估计同样非常不平坦,有诸多起伏。四、 Bartlett法功率谱估计(一)算法原理简介当我们用有关法或者
8、周期图法对信号的功率谱进行估计时,都不是对的一致估计,重要因素是方差大。于是就产生了周期图法的改善。改善的重要途径是平滑和平均。平滑是用一种合适的窗函数与计算的功率谱进行卷积,是谱线平滑。这种措施的出的谱估计是无偏的,方差也小,但辨别率下降。平均就是将截取的数据段再提成L个小段,分别计算功率谱后取功率谱的平均。由于L个平均的方差比随机变量的单独方差小L倍,因此当L趋于无穷时,L个平均的方差趋于零,可以达到一致估计的目的。(二)运算简要框图矩形窗截断对求得成果进行平均周期图法对每段进行计算X(n)提成L小段输出(三)程序示例%L=2时bartlett法N=256,n=0:255;x1n=2*co
9、s(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN的信号的前半段Xk1=fft(x1n,N); %进行N点FFTSk5=abs(Xk1).2./N;%取频谱幅度的平方,并除以N,以此作为对xn真实功率谱的估计n=256:511;x2n=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN的信号后半段Xk2=fft(x2n,N); %进行N点FFTSk6=abs(Xk2).2./N;%取频谱幅度的平方,并除以N,以此作为对xn真
10、实功率谱的估计Sk7=(Sk5+Sk6)/2; %相加求平均Sk8=10*log(Sk7);n=0:255;figure(4)subplot(2,1,1);plot(n/N,Sk7);title(d) Bartlett法估计功率谱密度 L=2)ylabel(Sk)subplot(2,1,2);plot(n/N,Sk8);ylabel(10log(PSD)%L=4时bartlett法N=128,n=0:127;x1n=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN的信号的1/4段Xk1=fft(x1n,
11、N);%进行N点FFTSk1=abs(Xk1).2./N;%取频谱幅度的平方,并除以N,以此作为对xn真实功率谱的估计n=128:255;x2n=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN的信号的1/4段Xk2=fft(x2n,N);%进行N点FFTSk2=abs(Xk2).*(abs(Xk2)./N;%取频谱幅度的平方,并除以N,以此作为对xn真实功率谱的估计N=128,n=256:383;x3n=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n
12、);%产生加有均值为0,方差为1的AWGN的信号的1/4段Xk3=fft(x3n,N);%进行N点FFTSk3=abs(Xk3).2./N;%取频谱幅度的平方,并除以N,以此作为对xn真实功率谱的估计n=384:511;x4n=2*cos(2*pi*0.2*n)+ cos(2*pi*0.4*n)+2*randn(size(n);%产生加有均值为0,方差为1的AWGN的信号的1/4段Xk4=fft(x4n,N);%进行N点FFTSk4=abs(Xk4).*(abs(Xk4)./N;%取频谱幅度的平方,并除以N,以此作为对xn真实功率谱的估计Sk5=(Sk1+Sk2+Sk3+Sk4)/4; %相加
13、求平均Sk6=10*log(Sk5);n=0:127;figure(5)subplot(2,1,1);plot(n/N,Sk5);title(e) Bartlett法估计功率谱密度 L=4)ylabel(Sk)subplot(2,1,2);plot(n/N,Sk6);ylabel(10log(PSD)(四)成果分析下面是程序运营后的成果上图分别是L=2和L=4时用bartlett法进行信号功率谱估计的波形。从上图中我们同样可以看到信号中有两个频率分量,一种在0.2处,一种在0.4处,与产生的信号相一致。但是我们不难看出,估计出的功率谱谱线与之前有关法功率谱估计和周期图法功率谱估计相比,波形相对
14、平坦了某些。L=2和L=4时用bartlett法进行信号功率谱估计的波形相比我们可以很明显的看出L大的那个波形更加平坦,这与之前在算法原理中简介的同样,L越大,平均后的方差就越小,越能达到一致估计的目的。五、 Welch法功率谱估计(一)算法原理简介目前比较常用的功率谱估计改善措施是Welch法,又叫加权交叠平均法。这种措施以加窗(加权)求取平滑,以分段重叠求得平均,因此集平均与平滑的长处于一体,同步也不可避免的带有两者的缺陷。其重要环节如下:第一步:将N长的数据段提成L个小段,每小段M点,相邻小段见交叠M/2点。第二步:对个小段加同样的品挂窗后求傅里叶变换第三步:求个小段功率谱的平均,得这里 (二)运算简要框图矩形窗截断相