public void testFactorFomelExample() { float[] r = {24.0f,242.0f,867.0f,1334.0f,867.0f,242.0f,24.0f}; int[] lag1 = {0,1,2,3}; CausalFilter cf = new CausalFilter(lag1); cf.factorWilsonBurg(10,0.0f,r); float[] a = cf.getA(); assertEquals(24.0f,a[0],10*FLT_EPSILON); assertEquals(26.0f,a[1],10*FLT_EPSILON); assertEquals( 9.0f,a[2],10*FLT_EPSILON); assertEquals( 1.0f,a[3],10*FLT_EPSILON); }
}; CausalFilter cf = new CausalFilter(lag1,lag2); cf.factorWilsonBurg(100,FLT_EPSILON,r); float[][] s = new float[3][3]; float[][] t = new float[3][3];
}; CausalFilter cf = new CausalFilter(lag1,lag2,lag3); cf.factorWilsonBurg(100,FLT_EPSILON,r); float[][][] s = new float[3][3][3]; float[][][] t = new float[3][3][3];
{ -p12*p12, -2.0f*m12*p12, -m12*m12} }; cf.factorWilsonBurg(maxiter,epsilon,r); dump(r); zero(s);
{ 0.0f, -m23*m23, 0.0f} }}; cf.factorWilsonBurg(maxiter,epsilon,r); float[] a = cf.getA(); for (int j=0; j<m; ++j) {