function eval_vcf_lut()
N = 1000;
fs = 48000;
fmin = 18;
fmax = 18000;
sweepbase = 10;

filterParam = struct('sweepbase', 10, 'omega_min', 2*fmin/fs, 'omega_max', 2*fmax/fs, 'num_octaves', log(fmax/fmin)/log(sweepbase));
filterParam

cv = (0:N-1)/N;
plot (c
function cv = omega2cv(filterParam, omega)
	cv = log(omega/filterParam.omega_min)/(filterParam.num_octaves * log(filterParam.sweepbase));
	cv = min(1, max(0, cv));

function omega = cv2omega(filterParam, cv)
	omega = filterParam.omega_min*filterParam.sweepbase ^ (filterParam.num_octaves*min(1, max(0, cv))); 
