bartlett_zfj.m

来自「文件包含现代信号处理的常用算法的matlab源代码」· M 代码 · 共 22 行

%产生信号序列
n=linspace(-pi,pi,1024);
xn=2*sin(n*40)+cos(n*100);
figure,plot(xn);
title('Input Signal');
%计算序列的DFT及PSD
M=64;
L=length(n)/M;
Pbw=zeros(1,M);
for i=1:L
    for j=1:M
        x(j)=xn(M*(i-1)+j);
    end
    y1=fft(x);
    y2=conj(y1);
    Plw=y1.*y2/M;
    Pbw=Plw+Pbw;
end
%绘图
figure,plot(Pbw/L);
grid on
title('Bartlett Method');