powerp.m

这里几种功率谱估计方法

% -------------------------------------------------------------------------
% 功能：通过仿真，研究平稳信号的功率谱估计的方法和特点
% 作者：
% 日期：2007年4月19日
% 状态：调试未完成
% -------------------------------------------------------------------------
% function powerp
clear all,close all
fs=1000;
T=1;
t=0:1/fs:T;
sn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(size(t));
% sn=exp(i*2*pi*50*t)+2*exp(i*2*pi*120*t)+randn(size(t));
figure
plot(t,sn);
% ------------经典周期法---------------------------------------------------%
nFFT1=256;
S1=fft(sn,nFFT1);
k1=0:round(nFFT1/2-1);
f1=k1*fs/nFFT1;
Pxx1=abs(S1(k1+1)).^2/nFFT1;
% Pxx=abs(S(k+1)).^2/length(t);
Pxx1=10*log10(Pxx1);
figure(2)
plot(f1,Pxx1);
% % ------------平均周期法---------------------------------------------------%
% nFFT2=256;
% Pxx2=zeros(1,256);
% % i=1;
% for i=0:2
%     S2=zeros(1,256);
%     S2=fft(sn((1+nFFT2*i):(nFFT2*(i+1))),nFFT2);
%     S2=abs(S2).^2/(nFFT2*3);
%     Pxx2=Pxx2+S2;
% end
% k2=0:round(nFFT2/2-1);
% f2=k2*fs/nFFT2;
% Pxx2=10*log10(Pxx2(k2+1));
% figure(3)
% plot(f2,Pxx2);
% % ------------平均修正周期法-----------------------------------------------%
% nFFT3=256;
% noverlap=128;
% window=hanning(256);
% [Pxx3,f3]=psd(sn,nFFT3,fs,window,noverlap);
% Pxx3=10*log10(Pxx3);
% figure(4)
% plot(f3,Pxx3);
% % ------------多正交窗口法-------------------------------------------------%
% nFFT4=256;
% nw=4;
% [Pxx4,f4]=pmtm(sn,nw,nFFT4,fs);
% Pxx4=10*log10(Pxx4);
% figure(5)
% plot(f4,Pxx4);
% % ------------最大墒法-----------------------------------------------------%
% nFFT5=256;
% order=100;
% [Pxx5,f5]=pmem(sn,order,nFFT5,fs);
% Pxx5=10*log10(Pxx5);
% figure(6)
% plot(f5,Pxx5);
% % ------------多信号分类特征向量法------------------------------------------%
% nFFT6=256;
% p=6;
% [Pxx6,f6]=pmusic(sn,p,nFFT6,fs);
% Pxx6=10*log10(Pxx6);
% figure(7)
% plot(f6,Pxx6);
% % ------------自相关函数法-------------------------------------------------%