dtawin.m

来自「A MATLAB tool for analysis of Transient 」· M 代码 · 共 32 行

function w=dtawin(fs,start,rise,len,fall)

%DTAWIN makes cosine taper window for ILO88 data
%   w = dtawin(fs, start, rise, len, fall)
%
%   in: fs      sampling frequency (default 25000 Hz)
%       start   start time [ms] (default 2.5 ms)
%       rise    rise time of the window [ms] (default 2.5 ms)
%       len     length of the window (default 512)
%       fall    fall time in the end [ms] (default 2.0 ms)    
%
%  out: w       window

% Pekka Kumpulainen 17.8.1993 (22.9.1993 P.K)

if nargin < 1; fs = 25000; end
if nargin < 2; start = 2.5; end
if nargin < 3; rise = 2.5; end
if nargin < 4; len = 512; end
if nargin < 5; fall = 2; end

start = floor(fs/1000*start);
rise = floor(fs/1000*rise);
if rise > len; error('Rise time cannot be longer than window'); end
fall = floor(fs/1000*fall);


w = [zeros(start-1,1);...
  .5 - 1/2*cos(pi*(0:rise)'/rise);...
  ones(len-start-rise-fall-1,1);...
  .5 + 1/2*cos(pi*(0:fall)'/fall)];