00001 #include <math.h>
00002 #include "d_lpc10.h"
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015 VOID irc2pc( FLOAT rc[], FLOAT pc[], FLOAT gprime, FLOAT *g2pass )
00016 {
00017 INDEX i, j;
00018 FLOAT temp[ORDER], t;
00019
00020 t = 1.;
00021
00022 for (i = 0; i < ORDER; i++)
00023 t *= (1 - rc[i] * rc[i]);
00024
00025 *g2pass = (FLOAT)(gprime * sqrt(t));
00026 pc[0] = rc[0];
00027
00028 for (i = 1; i < ORDER; i++) {
00029 for (j = 0; j < i; j++)
00030 temp[j] = pc[j] - rc[i] * pc[i-j-1];
00031
00032 for (j = 0; j < i; j++)
00033 pc[j] = temp[j];
00034
00035 pc[i] = rc[i];
00036 }
00037 }
00038
00039