最新版功率谱密度估计方法的MATLAB实现解析.docxVIP

精品文档 精品文档 PAGE PAGE #欢迎。下载 精品文档 精品文档 PAGE PAGE #欢迎下载 功率谱密度估计方法的 MATLA实现 在应用数学和物理学中， 谱密度、功率谱密度和能量谱密度是一个用于信号的通用概念， 它表示每赫兹的功率、每赫兹的能量这样的物理量纲。在物理学中，信号通常是波的形式， 例如电磁波、 随机振动或者声波。 当波的频谱密度乘以一个适当的系数后将得到每单位频率 波携带的功率，这被称为信号的功率谱密度( power spectral density, PSD )或者谱功率 分布( spectral power distribution, SPD )。功率谱密度的单位通常用每赫兹的瓦特数 (W/Hz)表示，或者使用波长而不是频率，即每纳米的瓦特数( W/nm来表示。信号的功率 谱密度当且仅当信号是广义的平稳过程的时候才存在。 如果信号不是平稳过程， 那么自相关 函数一定是两个变量的函数， 这样就不存在功率谱密度， 但是可以使用类似的技术估计时变 谱密度。 信号功率谱的概念和应用是电子工程的基础， 尤其是在电子通信系统中， 例如无线 电和微波通信、雷达以及相关系统。因此学习如何进行功率谱密度估计十分重要，借助于 Matlab 工具可以实现各种谱估计方法的模拟仿真并输出结果。下面对周期图法、修正周期 图法、最大熵法、 Levinson 递推法和 Burg 法的功率谱密度估计方法进行程序设计及仿真并 给出仿真结果。 以下程序运行平台： Matlab R2015a ( 8.5.0.197613 ) 周期图法谱估计程序 采样频率 100kHz数据长度 采样频率 100kHz 数据长度 N=1024 正弦波 ,f=2000Hz 加入信噪比为 10db 的高斯白噪声 Fs=100000; % N=1024; % n=0:N-1; t=n/Fs; xn=sin(2000*2*pi*t); % Y=awgn(xn,10); % subplot(2,1,1); plot(n,Y) title( 信号 ) xlabel( 时间);ylabel( 幅度）; grid on; win dow=boxcar(le ngth(x n)); % 矩形窗 n fft=N/4; % 采样点数 [Pxx f]=periodogram(Y,w in dow ,n fft,Fs); % 直接法 subplot(2,1,2); plot(f,10*log10(Pxx)); grid on; title([ 周期图法谱估计,,int2str(N), 点]); xlabel(频率(Hz) );ylabel( 功率谱密度); 2、仿真结果 12000-20■80400 600 SOO时间周期罔法谱估计,1024点-4Q -2 200 1000 1200 0 -20 ■80 400 600 SOO 时间 周期罔法谱估计,1024点 -4Q -60 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 频率 Hz) x104 修正周期图法（加窗）谱估计程序1 修正周期图法（加窗）谱估计程序 1、源程序 Fs=100000; % Fs=100000; % 采样频率100kHz N=512; N=512; % 数据长度 M=32; % 汉明窗宽度 n=0:N-1; t=n/Fs; xn=sin(2000*2*pi*t); % 正弦波 ,f=2000Hz Y=awgn(xn,10); % 加入信噪比为 10db 的高斯白噪声 subplot(2,1,1);subplot(2,1,1); plot(n,Y) title( 信号 ) xlabel( 时间 );ylabel( 幅度 ); grid on; window=hamming(M); % 汉明窗 [Pxx f]=pwelch(Y,window,10,256,Fs); subplot(2,1,2); plot(f,10*log10(Pxx)); grid on; title([ 修正周期图法谱估计 N=,int2str(N), M=,int2str(M)]); xlabel( 频率( Hz) );ylabel( 功率谱密度 ); 2、 仿真结果 精品文档 精品文档 PAGE PAGE #欢迎下载 精品文档 精品文档 PAGE PAGE #欢迎下载 时间 修正周期图法谱估计N-512 M=32 ^40-50 ^40 -50 詣0 0 0.5 1 1.5 2 25 3 3.5 4 4.S 5 频＜HZ) x104 三、最大熵法谱估计程序 1、源程序 fs=1; % 设采样频率 N=128; % 数据长度 改变数据长度会导致分辨率的变化；