kaiser.m

来自「matlabDigitalSigalProcess内有文件若干」· M 代码 · 共 21 行

function w = kaiser(n_est,beta)
%KAISER Kaiser window.
%   W = KAISER(N,beta) returns the BETA-valued N-point Kaiser window.

%   Author(s): L. Shure, 3-4-87
%   Copyright (c) 1988-98 by The MathWorks, Inc.
%   $Revision: 1.11 $  $Date: 1997/11/26 20:13:36 $

[nn,w,trivalwin] = check_order(n_est);
if trivalwin, return, end;

nw = round(nn);
bes = abs(besseli(0,beta));
odd = rem(nw,2);
xind = (nw-1)^2;
n = fix((nw+1)/2);
xi = (0:n-1) + .5*(1-odd);
xi = 4*xi.^2;
w = besseli(0,beta*sqrt(1-xi/xind))/bes;
w = abs([w(n:-1:odd+1) w])';