view general/numerical/array/gauss_win.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

function w=gauss_win(n,sigma)
% gauss_win - Gaussian window of given length and std deviation
%
% gauss_win :: N:natural, nonneg~'std dev as multiple of N/2' -> [[N]].

s=(1:n)'-(n+1)/2;
w=expquadf(2*s,sigma*n);