parzen.m

来自「matlab谱分析工具箱」· M 代码 · 共 28 行

function w=parzen(N)
% function w=parzen(N)
% N should be odd
%
% Programmed by: Dimitris Manolakis, 1996
%
%-----------------------------------------------------------
% Copyright 2000, by Dimitris G. Manolakis, Vinay K. Ingle,
% and Stephen M. Kogon.  For use with the book
% "Statistical and Adaptive Signal Processing"
% McGraw-Hill Higher Education.
%-----------------------------------------------------------



w=zeros(N,1);
N1=floor(N/2);
w1=zeros(N1+1,1);
N2=floor(N1/2);
n1=(0:N2)';
n2=(N2+1:N1)';

w1(n1+1)=1-6*(n1/N1).^2+6*(abs(n1)/N1).^3;
w1(n2+1)=2*(1-abs(n2)/N1).^3;

w(1:N1+1)=flipud(w1);
w(N1+2:N)=w1(2:N1+1);