% The function bl_filt.m creates a unit-area lowpass filter
% with a narrow transition band (truncated sinc) in order
% to generate bandlimited, white disturbances to Ra or Qvo.
%
% Function arguments:
%	fco - desired cutoff frequency
%	Sgran - desired granularity
%
% Function outputs:
%	h - vector representing discrete-time, autonomic filter
%

function h = bl_filt(fco,Sgran)

% Forming truncated sinc filter.
hra = 2*fco*sinc(2*fco*(-60:Sgran:60));

% Normalizing the filter to unity area.
h = hra/sum(hra);