00001 /* 00002 * SPUC - Signal processing using C++ - A DSP library 00003 * 00004 * This program is free software; you can redistribute it and/or modify 00005 * it under the terms of the GNU General Public License as published by 00006 * the Free Software Foundation; either version 2, or (at your option) 00007 * any later version. 00008 * 00009 * This program is distributed in the hope that it will be useful, 00010 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00012 * GNU General Public License for more details. 00013 * 00014 * You should have received a copy of the GNU General Public License 00015 * along with this program; if not, write to the Free Software 00016 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 00017 */ 00018 #ifndef AUTOC 00019 #define AUTOC 00020 namespace SPUC { 00021 // \brief Compute the autocorrelation of the Vector 00024 //: <font color="red"><i>Under construction!</i></font> 00025 template <class T> Vector<T> auto_corr(Vector<T> x) { 00026 int i,j; 00027 T d; 00028 long N = x.size(); 00029 Vector<T> autoc(N); 00030 00031 for (j=0;j<N;j++) { 00032 for (i = j, d = 0; i < N; i++) d += x[i] * x[i-j]; 00033 autoc[j] = d; 00034 } 00035 return(autoc); 00036 } 00037 } 00038 #endif