daniell_zfj.m

来自「使用matlab编写的8个现代信号处理程序」· M 代码 · 共 28 行

%*************************************%
%            Daniell Methed           % 
%*************************************%
%产生信号序列
n=linspace(-pi,pi,128);
x=2*sin(n*40)+cos(n*100);
figure,plot(x);
title('Input Signal');
%信号补零并计算Ppw
y1=fft(x,1024);
y2=conj(y1);
Ppw=y1.*y2/1024;
figure,plot(Ppw);
grid on
%对Ppw补零并计算PDw
Pp=zeros(1,64);
Ppw=[Pp,Ppw,Pp];
for k=1+64:1024+64
    PD1(k-64)=0;
    for j=k-64:k+64
        PD1(k-64)=PD1(k-64)+Ppw(j);
    end
    PDw(k-64)=PD1(k-64)/129;
end
%绘图
figure,plot(PDw);
grid on
title('Daniell');