%*************************************%
%              BT Method              % 
%*************************************%
%產生信號序列
n=linspace(-pi,pi,1024);
xn=2*sin(n*40)+cos(n*100);
figure,plot(xn);
title('Input Signal');
%計算序列的自相關函數
cxn=xcorr(xn,'unbiased');
%截短相關函數
m=round(length(n)/10);
cxn1=cxn(length(n)-m:length(n)+m);
%計算序列的PSD
CXk=fft(cxn1);
Pbt=abs(CXk);
%繪圖
figure,plot(Pbt);
grid on
title('BT Method');