www.pudn.com > g729_audio_encode.rar > ld8k.h
/*
ITU-T G.729 Speech Coder with Annex B ANSI-C Source Code
Version 1.3 Last modified: August 1997
Copyright (c) 1996,
AT&T, France Telecom, NTT, Universite de Sherbrooke, Lucent Technologies,
Rockwell International
All rights reserved.
*/
/*--------------------------------------------------------------*
* LD8K.H *
* ~~~~~~ *
* Function prototypes and constants use in G.729 *
* *
*--------------------------------------------------------------*/
/*--------------------------------------------------------------------------*
* Codec constant parameters (coder, decoder, and postfilter) *
*--------------------------------------------------------------------------*/
#define L_TOTAL 240 /* Total size of speech buffer. */
#define L_WINDOW 240 /* Window size in LP analysis. */
#define L_NEXT 40 /* Lookahead in LP analysis. */
#define L_FRAME 80 /* Frame size. */
#define L_SUBFR 40 /* Subframe size. */
#define M 10 /* Order of LP filter. */
#define MP1 (M+1) /* Order of LP filter + 1 */
#define MM1 (M-1) /* Order of LP filter - 1 */
#define PIT_MIN 20 /* Minimum pitch lag. */
#define PIT_MAX 143 /* Maximum pitch lag. */
#define L_INTERPOL (10+1) /* Length of filter for interpolation. */
#define PRM_SIZE 11 /* Size of vector of analysis parameters. */
#define SERIAL_SIZE (80+2) /* Bits/frame + bfi+ number of speech bits */
#define SHARPMAX 13017 /* Maximum value of pitch sharpening 0.8 Q14 */
#define SHARPMIN 3277 /* Minimum value of pitch sharpening 0.2 Q14 */
#define L_SUBFRP1 (L_SUBFR + 1)
#define GPCLIP 15564 /* Maximum pitch gain if taming is needed Q14*/
#define GPCLIP2 481 /* Maximum pitch gain if taming is needed Q9 */
#define GP0999 16383 /* Maximum pitch gain if taming is needed */
#define L_THRESH_ERR 983040000L /* Error threshold taming 16384. * 60000. */
/*--------------------------------------------------------------------------*
* Mathematic functions. *
*--------------------------------------------------------------------------*/
Word32 Inv_sqrt( /* (o) Q30 : output value (range: 0<=val<1) */
Word32 L_x /* (i) Q0 : input value (range: 0<=val<=7fffffff) */
);
void Log2(
Word32 L_x, /* (i) Q0 : input value */
Word16 *exponent, /* (o) Q0 : integer part of Log2. (range: 0<=val<=30)*/
Word16 *fraction /* (o) Q15 : fractional part of Log2. (range: 0<=val<1) */
);
Word32 Pow2( /* (o) Q0 : result (range: 0<=val<=0x7fffffff) */
Word16 exponent, /* (i) Q0 : integer part. (range: 0<=val<=30) */
Word16 fraction /* (i) Q15 : fractional part. (range: 0.0<=val<1.0) */
);
/*--------------------------------------------------------------------------*
* Pre and post-process. *
*--------------------------------------------------------------------------*/
void Init_Pre_Process(void);
void Init_Post_Process(void);
void Pre_Process(
Word16 signal[], /* (i/o) : input/output signal */
Word16 lg /* (i) : length of signal */
);
void Post_Process(
Word16 signal[], /* (i/o) : input/output signal */
Word16 lg /* (i) : length of signal */
);
/*--------------------------------------------------------------------------*
* Main coder and decoder functions *
*--------------------------------------------------------------------------*/
void Init_Coder_ld8k(void);
void Coder_ld8k(
Word16 ana[], /* (o) : analysis parameters */
Word16 synth[], /* (o) : local synthesis */
Word16 frame, /* (i) : frame counter */
Word16 vad_enable /* (i) : VAD enable flag */
);
void Init_Decod_ld8k(void);
void Decod_ld8k(
Word16 parm[], /* (i) : vector of synthesis parameters
parm[0] = bad frame indicator (bfi) */
Word16 voicing, /* (i) : voicing decision from previous frame */
Word16 synth[], /* (o) : synthesized speech */
Word16 A_t[], /* (o) : decoded LP filter for 2 subframes */
Word16 *T0_first,/* (o) : decoded pitch lag in first subframe */
Word16 *Vad /* (o) : frame type */
);
/*--------------------------------------------------------------------------*
* taming functions *
*--------------------------------------------------------------------------*/
void Init_exc_err(void);
Word16 test_err(Word16 t0, Word16 t0_frac);
void update_exc_err(Word16 gain_pit, Word16 t0);
/*--------------------------------------------------------------------------*
* LPC analysis and filtering *
*--------------------------------------------------------------------------*/
void Autocorr(
Word16 x[], /* (i) : input signal */
Word16 m, /* (i) : LPC order */
Word16 r_h[], /* (o) : autocorrelations (msb) */
Word16 r_l[], /* (o) : autocorrelations (lsb) */
Word16 *exp_R0
);
void Lag_window(
Word16 m, /* (i) : LPC order */
Word16 r_h[], /* (i/o) : autocorrelations (msb) */
Word16 r_l[] /* (i/o) : autocorrelations (lsb) */
);
void Levinson(
Word16 Rh[], /* (i) : Rh[m+1] autocorrelation coefficients (msb) */
Word16 Rl[], /* (i) : Rl[m+1] autocorrelation coefficients (lsb) */
Word16 A[], /* (o) Q12 : A[m] LPC coefficients (m = 10) */
Word16 rc[], /* (o) Q15 : rc[M] Reflection coefficients. */
Word16 *Err /* (o) : Residual energy */
);
void Az_lsp(
Word16 a[], /* (i) Q12 : predictor coefficients */
Word16 lsp[], /* (o) Q15 : line spectral pairs */
Word16 old_lsp[] /* (i) : old lsp[] (in case not found 10 roots) */
);
void Lsp_Az(
Word16 lsp[], /* (i) Q15 : line spectral frequencies */
Word16 a[] /* (o) Q12 : predictor coefficients (order = 10) */
);
void Lsf_lsp(
Word16 lsf[], /* (i) Q15 : lsf[m] normalized (range: 0.0<=val<=0.5) */
Word16 lsp[], /* (o) Q15 : lsp[m] (range: -1<=val<1) */
Word16 m /* (i) : LPC order */
);
void Lsp_lsf(
Word16 lsp[], /* (i) Q15 : lsp[m] (range: -1<=val<1) */
Word16 lsf[], /* (o) Q15 : lsf[m] normalized (range: 0.0<=val<=0.5) */
Word16 m /* (i) : LPC order */
);
void Int_lpc(
Word16 lsp_old[], /* (i) : LSP vector of past frame */
Word16 lsp_new[], /* (i) : LSP vector of present frame */
Word16 lsf_int[], /* (o) : interpolated lsf coefficients */
Word16 lsf_new[], /* (o) : new lsf coefficients */
Word16 Az[] /* (o) : interpolated Az() for the 2 subframes */
);
void Int_qlpc(
Word16 lsp_old[], /* (i) : LSP vector of past frame */
Word16 lsp_new[], /* (i) : LSP vector of present frame */
Word16 Az[] /* (o) : interpolated Az() for the 2 subframes */
);
/*--------------------------------------------------------------------------*
* PWF constant parameters *
*--------------------------------------------------------------------------*/
#define A1 4567 /* 2.23 in Q11 */
#define L_B1 3271557L /* 0.78 in Q22 */
#define A2 11776 /* 5.75 in Q11 */
#define L_B2 16357786L /* 3.90 in Q22 */
#define A3 27443 /* 13.40 in Q11 */
#define L_B3 46808433L /* 11.16 in Q22 */
#define SEG1 1299 /* 0.6341 in Q11 */
#define SEG2 1815 /* 0.8864 in Q11 */
#define SEG3 1944 /* 0.9490 in Q11 */
#define THRESH_L1 -3562 /* -1.74 in Q11 */
#define THRESH_L2 -3116 /* -1.52 in Q11 */
#define THRESH_H1 1336 /* 0.65 in Q11 */
#define THRESH_H2 890 /* 0.43 in Q11 */
#define GAMMA1_0 32113 /* 0.98 in Q15 */
#define GAMMA1_1 30802 /* 0.94 in Q15 */
#define GAMMA2_0_L 13107 /* 0.40 in Q15 */
#define GAMMA2_0_H 22938 /* 0.70 in Q15 */
#define GAMMA2_1 19661 /* 0.60 in Q15 */
#define ALPHA 19302 /* 6*pi in Q10 */
#define BETA 1024 /* 1 in Q10 */
void perc_var (
Word16 *gamma1, /* bandwidth expansion parameter */
Word16 *gamma2, /* bandwidth expansion parameter */
Word16 *lsfint, /* Interpolated LSP vector : 1st subframe */
Word16 *lsfnew, /* New LSP vector : 2nd subframe */
Word16 *r_c /* Reflection coefficients */
);
void Weight_Az(
Word16 a[], /* (i) Q12 : a[m+1] LPC coefficients */
Word16 gamma, /* (i) Q15 : Spectral expansion factor. */
Word16 m, /* (i) : LPC order. */
Word16 ap[] /* (o) Q12 : Spectral expanded LPC coefficients */
);
void Residu(
Word16 a[], /* (i) Q12 : prediction coefficients */
Word16 x[], /* (i) : speech (values x[-m..-1] are needed (m=10) */
Word16 y[], /* (o) : residual signal */
Word16 lg /* (i) : frame size */
);
void Syn_filt(
Word16 a[], /* (i) Q12 : a[m+1] prediction coefficients (m=10) */
Word16 x[], /* (i) : input signal */
Word16 y[], /* (o) : output signal */
Word16 lg, /* (i) : size of filtering */
Word16 mem[], /* (i/o) : memory associated with this filtering. */
Word16 update /* (i) : 0=no update, 1=update of memory. */
);
void Convolve(
Word16 x[], /* (i) : input vector */
Word16 h[], /* (i) Q12 : impulse response */
Word16 y[], /* (o) : output vector */
Word16 L /* (i) : vector size */
);
/*--------------------------------------------------------------------------*
* LTP constant parameters *
*--------------------------------------------------------------------------*/
#define THRESHPIT 27853 /* Threshold to favor small pitch 0.85 Q15 */
#define UP_SAMP 3 /* upsampling factor for fractional pitch */
#define L_INTER4 4 /* length/2 for interpolation filter */
#define FIR_SIZE_ANA (UP_SAMP*L_INTER4+1)
#define L_INTER10 10 /* length/2 for interpolation filter */
#define FIR_SIZE_SYN (UP_SAMP*L_INTER10+1)
/*--------------------------------------------------------------------------*
* Pitch functions. *
*--------------------------------------------------------------------------*/
Word16 Pitch_ol( /* (o) : open loop pitch lag */
Word16 signal[], /* (i) : signal used to compute the open loop pitch */
/* signal[-pit_max] to signal[-1] should be known */
Word16 pit_min, /* (i) : minimum pitch lag */
Word16 pit_max, /* (i) : maximum pitch lag */
Word16 L_frame /* (i) : length of frame to compute pitch */
);
Word16 Pitch_fr3( /* (o) : pitch period. */
Word16 exc[], /* (i) : excitation buffer */
Word16 xn[], /* (i) : target vector */
Word16 h[], /* (i) Q12 : impulse response of filters. */
Word16 L_subfr, /* (i) : length of subframe */
Word16 t0_min, /* (i) : minimum value in the searched range. */
Word16 t0_max, /* (i) : maximum value in the searched range. */
Word16 i_subfr, /* (i) : indicator for first subframe. */
Word16 *pit_frac /* (o) : chosen fraction. */
);
Word16 G_pitch( /* (o) Q14 : Gain of pitch lag saturated to 1.2 */
Word16 xn[], /* (i) : Pitch target. */
Word16 y1[], /* (i) : Filtered adaptive codebook. */
Word16 g_coeff[], /* (i) : Correlations need for gain quantization. */
Word16 L_subfr /* (i) : Length of subframe. */
);
Word16 Enc_lag3( /* (o) : Return index of encoding */
Word16 T0, /* (i) : Pitch delay */
Word16 T0_frac, /* (i) : Fractional pitch delay */
Word16 *T0_min, /* (i/o) : Minimum search delay */
Word16 *T0_max, /* (i/o) : Maximum search delay */
Word16 pit_min, /* (i) : Minimum pitch delay */
Word16 pit_max, /* (i) : Maximum pitch delay */
Word16 pit_flag /* (i) : Flag for 1st subframe */
);
void Dec_lag3( /* (o) : return integer pitch lag */
Word16 index, /* (i) : received pitch index */
Word16 pit_min, /* (i) : minimum pitch lag */
Word16 pit_max, /* (i) : maximum pitch lag */
Word16 i_subfr, /* (i) : subframe flag */
Word16 *T0, /* (o) : integer part of pitch lag */
Word16 *T0_frac /* (o) : fractional part of pitch lag */
);
Word16 Interpol_3( /* (o) : interpolated value */
Word16 *x, /* (i) : input vector */
Word16 frac /* (i) : fraction */
);
void Pred_lt_3(
Word16 exc[], /* (i/o) : excitation buffer */
Word16 T0, /* (i) : integer pitch lag */
Word16 frac, /* (i) : fraction of lag */
Word16 L_subfr /* (i) : subframe size */
);
Word16 Parity_Pitch( /* (o) : parity bit (XOR of 6 MSB bits) */
Word16 pitch_index /* (i) : index for which parity to compute */
);
Word16 Check_Parity_Pitch( /* (o) : 0 = no error, 1= error */
Word16 pitch_index, /* (i) : index of parameter */
Word16 parity /* (i) : parity bit */
);
/*--------------------------------------------------------------------------*
* fixed codebook excitation. *
*--------------------------------------------------------------------------*/
/*--------------------------------------------------------------------------*
* FCB constant parameters *
*--------------------------------------------------------------------------*/
#define DIM_RR 616 /* size of correlation matrix */
#define NB_POS 8 /* */
#define STEP 5 /* spacing for individual pulse */
#define MSIZE 64
#define THRESHFCB 13107 /* 0.4 in Q15 */
#define MAX_TIME 75 /* maximum number of iterations */
/*--------------------------------------------------------------------------*
* FCB functions. *
*--------------------------------------------------------------------------*/
Word16 ACELP_Codebook( /* (o) : index of pulses positions */
Word16 x[], /* (i) : Target vector */
Word16 h[], /* (i) Q12 : Impulse response of filters */
Word16 T0, /* (i) : Pitch lag */
Word16 pitch_sharp,/* (i) Q14: Last quantized pitch gain */
Word16 i_subfr, /* (i) : Indicator of 1st subframe, */
Word16 code[], /* (o) Q13 : Innovative codebook */
Word16 y[], /* (o) Q12 : Filtered innovative codebook */
Word16 *sign /* (o) : Signs of 4 pulses */
);
void Decod_ACELP(
Word16 sign, /* (i) : signs of 4 pulses. */
Word16 index, /* (i) : Positions of the 4 pulses. */
Word16 cod[] /* (o) Q13 : algebraic (fixed) codebook excitation */
);
/*--------------------------------------------------------------------------*
* LSP quantizer constant parameters *
*--------------------------------------------------------------------------*/
#define NC 5 /* NC = M/2 */
#define MA_NP 4 /* MA prediction order for LSP */
#define MODE 2 /* number of modes for MA prediction */
#define NC0_B 7 /* number of first stage bits */
#define NC1_B 5 /* number of second stage bits */
#define NC0 (1< -2 */
/* , -2, 2 */
Word16 exp_coeff[],/* (i) : Q-Format g_coeff[] */
Word16 L_subfr, /* (i) : Subframe length. */
Word16 *gain_pit, /* (o) Q14 : Pitch gain. */
Word16 *gain_cod, /* (o) Q1 : Code gain. */
Word16 tameflag /* (i) : flag set to 1 if taming is needed */
);
void Dec_gain(
Word16 index, /* (i) : Index of quantization. */
Word16 code[], /* (i) Q13 : Innovative vector. */
Word16 L_subfr, /* (i) : Subframe length. */
Word16 bfi, /* (i) : Bad frame indicator */
Word16 *gain_pit, /* (o) Q14 : Pitch gain. */
Word16 *gain_cod /* (o) Q1 : Code gain. */
);
void Gain_predict(
Word16 past_qua_en[],/* (i) Q10 :Past quantized energies */
Word16 code[], /* (i) Q13 : Innovative vector. */
Word16 L_subfr, /* (i) : Subframe length. */
Word16 *gcode0, /* (o) Qxx : Predicted codebook gain */
Word16 *exp_gcode0 /* (o) : Q-Format(gcode0) */
);
void Gain_update(
Word16 past_qua_en[],/* (i) Q10 :Past quantized energies */
Word32 L_gbk12 /* (i) Q13 : gbk1[indice1][1]+gbk2[indice2][1] */
);
void Gain_update_erasure(
Word16 past_qua_en[]/* (i) Q10 :Past quantized energies */
);
void Corr_xy2(
Word16 xn[], /* (i) Q0 :Target vector. */
Word16 y1[], /* (i) Q0 :Adaptive codebook. */
Word16 y2[], /* (i) Q12 :Filtered innovative vector. */
Word16 g_coeff[], /* (o) Q[exp]:Correlations between xn,y1,y2 */
Word16 exp_g_coeff[] /* (o) :Q-format of g_coeff[] */
);
/*--------------------------------------------------------------------------*
* Postfilter constants and functions *
*--------------------------------------------------------------------------*/
/* short term pst parameters : */
#define GAMMA1_PST 22938 /* denominator weighting factor (Q15) */
#define GAMMA2_PST 18022 /* numerator weighting factor (Q15) */
#define LONG_H_ST 20 /* impulse response length */
#define GAMMA3_PLUS 6554 /* tilt weighting factor when k1>0 (Q15) */
#define GAMMA3_MINUS 29491 /* tilt weighting factor when k1<0 (Q15) */
/* long term pst parameters : */
#define F_UP_PST 8 /* resolution for fractional delay */
#define LH2_S 4 /* length of short interp. subfilters */
#define L2_LH2_L 4 /* log2(LH2_L) */
#define LH2_L (1 << L2_LH2_L)
#define MIN_GPLT 21845 /* LT gain minimum (Q15) */
#define LH_UP_S (LH2_S/2)
#define LH_UP_SM1 (LH_UP_S-1)
#define LH_UP_L (LH2_L/2)
#define LH2_L_P1 (LH2_L + 1)
/* gain adjustment parameters */
#define AGC_FAC 32358 /* gain adjustment factor 0.9875 (Q15) */
#define AGC_FAC1 (Word16)(32768L - AGC_FAC)
/* Array sizes */
#define MEM_RES2 (PIT_MAX + 1 + LH_UP_L)
#define SIZ_RES2 (MEM_RES2 + L_SUBFR)
#define SIZ_Y_UP ((F_UP_PST-1) * L_SUBFRP1)
#define SIZ_TAB_HUP_L ((F_UP_PST-1) * LH2_L)
#define SIZ_TAB_HUP_S ((F_UP_PST-1) * LH2_S)
void Init_Post_Filter( void);
void Post(
Word16 t0, /* (i) : 1st subframe delay given by coder */
Word16 *signal_ptr, /* (i) : input signal (pointer to current subframe */
Word16 *coeff, /* (i) : LPC coefficients for current subframe */
Word16 *sig_out, /* (o) : postfiltered output */
Word16 *vo, /* (o) : voicing decision 0 = uv, > 0 delay */
Word16 Vad /* (i) : frame type */
);
/*--------------------------------------------------------------------------*
* Bitstream constants and functions *
*--------------------------------------------------------------------------*/
#define BIT_0 (short)0x007f /* definition of zero-bit in bit-stream */
#define BIT_1 (short)0x0081 /* definition of one-bit in bit-stream */
#define SYNC_WORD (short)0x6b21 /* definition of frame erasure flag */
#define SIZE_WORD (short)80 /* number of speech bits */
void prm2bits_ld8k(
Word16 prm[], /* (i) : coder parameters */
Word16 bits[] /* (o) : bit stream */
);
void bits2prm_ld8k(
Word16 bits[], /* (i) : bit stream */
Word16 prm[] /* (o) : coder parameters */
);
/*--------------------------------------------------------------------------*
* Prototypes for auxiliary functions. *
*--------------------------------------------------------------------------*/
void Copy(
Word16 x[], /* (i) : input vector */
Word16 y[], /* (o) : output vector */
Word16 L /* (i) : vector length */
);
void Set_zero(
Word16 x[], /* (o) : vector to clear */
Word16 L /* (i) : length of vector */
);
Word16 Random(Word16 *seed);