LPCAA2RF LPC: Convert vocal tract areas to reflection coefficients RF=(AA)
0001 function rf=lpcaa2rf(aa) 0002 %LPCAA2RF LPC: Convert vocal tract areas to reflection coefficients RF=(AA) 0003 0004 % Copyright (C) Mike Brookes 1998 0005 % 0006 % Last modified Tue May 12 15:58:43 1998 0007 % 0008 % VOICEBOX home page: http://www.ee.ic.ac.uk/hp/staff/dmb/voicebox/voicebox.html 0009 % 0010 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0011 % This program is free software; you can redistribute it and/or modify 0012 % it under the terms of the GNU General Public License as published by 0013 % the Free Software Foundation; either version 2 of the License, or 0014 % (at your option) any later version. 0015 % 0016 % This program is distributed in the hope that it will be useful, 0017 % but WITHOUT ANY WARRANTY; without even the implied warranty of 0018 % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 0019 % GNU General Public License for more details. 0020 % 0021 % You can obtain a copy of the GNU General Public License from 0022 % ftp://prep.ai.mit.edu/pub/gnu/COPYING-2.0 or by writing to 0023 % Free Software Foundation, Inc.,675 Mass Ave, Cambridge, MA 02139, USA. 0024 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0025 0026 [nf,p2]=size(aa); 0027 rf = (aa(:,2:p2)-aa(:,1:p2-1))./(aa(:,2:p2)+aa(:,1:p2-1));