1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
function [k,R0] = ac2rc(R)
// Convert autocorrelation sequence to reflection coefficients.
//
// Calling Sequence
// k = ac2rc(R)
// [k,R0] = ac2rc(R)
//
// Parameters
// R: The input autocorrelation sequence. If r is a matrix, each column of r is treated as a separate signal.
// k: Returns the reflection coefficients
// R0: the zero lag autocorrelation, R0, based on the autocorrelation sequence, R.
//
// Examples
// X = [7 6 5 8 3 6 8 7 5 2 4 7 4 3 2 5 4 9 5 3 5 7 3 9 4 1 2 0 5 4 8 6 4 6 5 3];
// [k,R0] = ac2rc(X)
// or t=[2 5 6; 8 6 5; 8 9 4]
// [k,R0] = ac2rc(t)
//
// Author
// Jitendra Singh
//
// call function "levin" before running this function
if or(type(R)==10) then
error ('Input arguments must be double.')
end
if isvector(R) then
R=R(:);
[x,y,z] = levin(R)
k=z;
R0=R;
else
n=size(R);
for i=1:n(2)
r=R(:,i);
[x,y, z] = levin(r)
kk(:,i)= z;
k=kk;
R0=R(1,:);
end
end
endfunction
|
