bt algorithm.m

近代信号处理功率谱估计的BT算法

%序列的间接法功率谱估计
%to test BT.m;

clear all;
Fs=1000;
%产生含有噪声的序列
var=sqrt(1/exp(1.0));
n=0:1/Fs:1;
N=length(n);
e=var*randn(1,N);
w0=100*pi;
w1=50*pi;
xn=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e;

%绘制信号波形
subplot(311)
plot(n,abs(xn))
xlabel('n')
title('x(n)=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e(n)')

m=-500:500;
[r,lag]=xcorr(xn,500,'biased');
subplot(312)
hndl=stem(m,r);
set(hndl,'Marker','.')
set(hndl,'Markersize',2);
ylabel('自相关函数R(m)')

%利用间接法计算功率谱
k=0:1000;
w=(pi/500)*k;
M=k/500;
X=r*(exp(-j*pi/500).^(m'*k));
magX=abs(X);
subplot(313)
plot(M,10*log10(magX));
axis([0 2 -50 50])
xlabel('功率谱的BT法估计')