diff options
Diffstat (limited to 'macros/rc2poly.sci')
| -rw-r--r-- | macros/rc2poly.sci | 18 |
1 files changed, 11 insertions, 7 deletions
diff --git a/macros/rc2poly.sci b/macros/rc2poly.sci index 19eed2a..a3665f2 100644 --- a/macros/rc2poly.sci +++ b/macros/rc2poly.sci @@ -1,28 +1,32 @@ function [a, efinal] = rc2poly(kr, R0) - + //rc2poly function convert reflection coefficients to prediction polynomial. // Calling Sequence // a = rc2poly(kr) // [a, efinal] = rc2poly(kr,R0) // Parameters -// kr: Refelection coefficient. +// kr: Refelection coefficient. // R0: the zero lag autocorrelation, R0. // a: Return the prediction polynomial. // efinal: Return the final prediction error. // Examples //X = [7 6 5 8 3 6] -// [a, efinal] = rc2poly(X) +// [a, efinal] = rc2poly(X) //error as only one input parameter is specified.(R0 is not mentioned) +// +//k = [0.3090 0.9800 0.0031 0.0082 -0.0082]; +//a = rc2poly(k) +//EXPECTED OUTPUT:a= 1. 0.6148162 0.9898814 0.0000243 0.0031580 - 0.0082 // // See also // // Author // Jitendra Singh -// +// if or(type(kr)==10) then error ('Input arguments must be double.') -end +end if (size(kr,1) > 1) & (size(kr,2) > 1) error ('The reflection coefficients must be stored in a vector.') @@ -47,10 +51,10 @@ a=kr(1); for i=2:length(kr) a = [a + a(i-1:-1:1)*kr(i) kr(i)] - ee=ee.*(1-kr(i)'.*kr(i)); + ee=ee.*(1-kr(i)'.*kr(i)); end a=[1 a]; -efinal=ee; +efinal=ee; endfunction |