资源描述：
《功率谱密度估计方法的MATLAB实现.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、.word格式.功率谱密度估计方法的MATLAB实现在应用数学和物理学中，谱密度、功率谱密度和能量谱密度是一个用于信号的通用概念，它表示每赫兹的功率、每赫兹的能量这样的物理量纲。在物理学中，信号通常是波的形式，例如电磁波、随机振动或者声波。当波的频谱密度乘以一个适当的系数后将得到每单位频率波携带的功率，这被称为信号的功率谱密度（powerspectraldensity,PSD）或者谱功率分布（spectralpowerdistribution,SPD）。功率谱密度的单位通常用每赫兹的瓦特数（W/Hz）表示，
2、或者使用波长而不是频率，即每纳米的瓦特数（W/nm）来表示。信号的功率谱密度当且仅当信号是广义的平稳过程的时候才存在。如果信号不是平稳过程，那么自相关函数一定是两个变量的函数，这样就不存在功率谱密度，但是可以使用类似的技术估计时变谱密度。信号功率谱的概念和应用是电子工程的基础，尤其是在电子通信系统中，例如无线电和微波通信、雷达以及相关系统。因此学习如何进行功率谱密度估计十分重要，借助于Matlab工具可以实现各种谱估计方法的模拟仿真并输出结果。下面对周期图法、修正周期图法、最大熵法、Levinson递推法和
3、Burg法的功率谱密度估计方法进行程序设计及仿真并给出仿真结果。以下程序运行平台：MatlabR2015a（8.5.0.197613）一、周期图法谱估计程序1、源程序Fs=100000;%采样频率100kHzN=1024;%数据长度N=1024n=0:N-1;t=n/Fs;xn=sin(2000*2*pi*t);%正弦波,f=2000HzY=awgn(xn,10);%加入信噪比为10db的高斯白噪声subplot(2,1,1);plot(n,Y)title('信号').专业.专注..word格式.xlabe
4、l('时间');ylabel('幅度');gridon;window=boxcar(length(xn));%矩形窗nfft=N/4;%采样点数[Pxxf]=periodogram(Y,window,nfft,Fs);%直接法subplot(2,1,2);plot(f,10*log10(Pxx));gridon;title(['周期图法谱估计,',int2str(N),'点']);xlabel('频率（Hz）');ylabel('功率谱密度');1、仿真结果一、修正周期图法（加窗）谱估计程序1、源程序Fs=
5、100000;%采样频率100kHzN=512;%数据长度M=32;%汉明窗宽度.专业.专注..word格式.n=0:N-1;t=n/Fs;xn=sin(2000*2*pi*t);%正弦波,f=2000HzY=awgn(xn,10);%加入信噪比为10db的高斯白噪声subplot(2,1,1);subplot(2,1,1);plot(n,Y)title('信号')xlabel('时间');ylabel('幅度');gridon;window=hamming(M);%汉明窗[Pxxf]=pwelch(Y,w
6、indow,10,256,Fs);subplot(2,1,2);plot(f,10*log10(Pxx));gridon;title(['修正周期图法谱估计N=',int2str(N),'M=',int2str(M)]);xlabel('频率（Hz）');ylabel('功率谱密度');2、仿真结果.专业.专注..word格式.一、最大熵法谱估计程序1、源程序fs=1;%设采样频率N=128;%数据长度改变数据长度会导致分辨率的变化；f1=0.2*fs;%第一个sin信号的频率,f1/fs=0.2f2=0.
7、3*fs;%第二个sin信号的频率，f2/fs=0.2或者0.3P=10;%滤波器阶数n=1:N;s=sin(2*pi*f1*n/fs)+sin(2*pi*f2*n/fs);%s为原始信号x=awgn(s,10);%x为观测信号，即对原始信号加入白噪声，信噪比10dBfigure(1);%画出原始信号和观测信号subplot(2,1,1);plot(s,'b'),xlabel('时间'),ylabel('幅度'),title('原始信号s');grid;subplot(2,1,2);plot(x,'r'),
8、xlabel('时间'),ylabel('幅度'),title('观测信号x');.专业.专注..word格式.[Pxx1,f]=pmem(x,P,N,fs);%最大熵谱估计figure(2);plot(f,10*log10(Pxx1));xlabel('频率(Hz)');ylabel('功率谱(dB)');title(['最大熵法谱估计模型阶数P=',int2str(P),'数据长度N=',int2str(N)
