function x = vqvectors(s,fs)
% function x = vqvectors(s,fs)
%   Return the LSF-vectors of each frame of a signal 's' with sampling
%   frequency 'fs'.

m = 10;                         % model order
winlen_ms = 20;                 % window length in ms
winlen = fs*winlen_ms/1000;     % window length in samples
frames = floor(length(s)/winlen);   % number of frames
win = hamming(winlen);          % windowing function
s = s(:);                       % make sure the signal is a column vector

% iterate for each frame
for k=1:frames
   h = win.*s(1+(k-1)*winlen+(0:winlen-1));
   a = lpc(h,m);
   x(k,:) = poly2lsf(a);
end