diff options
Diffstat (limited to 'macros/arburg.sci~')
| -rw-r--r-- | macros/arburg.sci~ | 55 |
1 files changed, 0 insertions, 55 deletions
diff --git a/macros/arburg.sci~ b/macros/arburg.sci~ deleted file mode 100644 index 1ca6e01..0000000 --- a/macros/arburg.sci~ +++ /dev/null @@ -1,55 +0,0 @@ -function varargout = arburg( x, poles, criterion ) -//This function calculates coefficients of an autoregressive (AR) model of complex data. -//Calling Sequence -//a = arburg(x, poles) -//a = arburg(x, poles, criterion) -//[a, v] = arburg(...) -//[a, v, k] = arburg(...) -//Parameters -//x: vector of real or complex numbers, of length > 2 -//poles: positive integer value < length(x) - 2 -//criterion: string value, takes in "AKICc", "KIC", "AICc", "AIC" and "FPE", default it not using a model-selection criterion -//a, v, k: Output variables -//Description -//This is an Octave function. -//This function calculates coefficients of an autoregressive (AR) model of complex data x using the whitening lattice-filter method of Burg. -//The first argument is the data sampled. The second argument is the number of poles in the model (or limit in case a criterion is supplied). -//The third parameter takes in the criterion to limit the number of poles. The acceptable values are "AIC", "AKICc", "KIC", "AICc" which are based on information theory. -//Output variable a is a list of P+1 autoregression coefficients. -//Output variable v is the mean square of residual noise from the whitening operation of the Burg lattice filter. -//Output variable k corresponds to the reflection coefficients defining the lattice-filter embodiment of the model. -//Examples -//arburg([1,2,3,4,5],2) -//ans = -// 1.00000 -1.86391 0.95710 - -funcprot(0); -rhs = argn(2) -lhs = argn(1) -if(lhs>3) -error("Wrong number of output arguments.") -elseif(rhs<2) -error("Wrong number of input arguments.") -end - - select(lhs) - case 1 then - if(rhs==2) - a = callOctave("arburg",x,poles) - elseif(rhs==3) - a = callOctave("arburg",x,poles,criterion) - end - case 2 then - if(rhs==2) - [a,v] = callOctave("arburg",x,poles) - elseif(rhs==3) - [a,v] = callOctave("arburg",x,poles,criterion) - end - case 3 then - if(rhs==2) - [a,v,k] = callOctave("arburg",x,poles) - elseif(rhs==3) - [a,v,k] = callOctave("arburg",x,poles,criterion) - end - end -endfunction |