levinson durbin algo (again)

Started by John McDermick in comp.dsp10 years ago 4 replies

Hi, I keep on getting a divide-by-zero exception when I execute the following code: Can somebody please tell me where my fixed-point math...

Hi, I keep on getting a divide-by-zero exception when I execute the following code: Can somebody please tell me where my fixed-point math is "broken" ??? Thank you! void levinsonDurbin(SINT32 *a, SINT32 *k, SINT32 *r, SINT32 n) { SINT32 a_temp[10+1], alpha, epsilon; int i, j; SINT32 tmpval; k[0] = 0; a[0] = 32768; a_temp[0] = 32768; alpha =...


Levinson durban recursion

Started by mobi in comp.dsp14 years ago 4 replies

Hi, I am trying to get forward/backward predictor coefficients from the ladder relection coefficients using the following recursion....

Hi, I am trying to get forward/backward predictor coefficients from the ladder relection coefficients using the following recursion. (From haykins book) a_m(n) = [a_{m-1}(n) 0] + k_{f, m} [0 c_{m-1}(n-1)] --------- (1) c_m(n) = [0 c_{m-1}(n-1) + k_{b, m} [a_{m-1}(n) 0] --------- (2) where a_m and c_m correspond to forward/backward predictor coefficients in this manner. a_m = [1 -w_...