function [Y,A] = whiten(X,P)
% y = whiten(X,P)
%  Fit a P'th order LPC to the whole of X, and inverse-filter by
%  it.  Y is the whitened version.
%  A returns the actual whitening FIR coefficients.
% 2010-11-27 Dan Ellis dpwe@ee.columbia.edu

if nargin < 2;  P = 40; end

A = lpc(X,P);
Y = filter(A,1,X);