1
<?xml version="1.0" encoding="ISO-8859-1" standalone="no"?>
2
<!DOCTYPE MAN SYSTEM "../../manrev.dtd">
4
<LANGUAGE>eng</LANGUAGE>
6
<TYPE>Scilab Function</TYPE>
7
<DATE>April 1993</DATE>
8
<SHORT_DESCRIPTION name="srfaur"> square-root algorithm</SHORT_DESCRIPTION>
10
<CALLING_SEQUENCE_ITEM>[p,s,t,l,rt,tt]=srfaur(h,f,g,r0,n,p,s,t,l) </CALLING_SEQUENCE_ITEM>
15
<PARAM_NAME>h, f, g</PARAM_NAME>
17
<SP>: convenient matrices of the state-space model.</SP>
21
<PARAM_NAME>r0</PARAM_NAME>
27
<PARAM_NAME>n</PARAM_NAME>
29
<SP>: number of iterations.</SP>
33
<PARAM_NAME>p</PARAM_NAME>
35
<SP>: estimate of the solution after n iterations.</SP>
39
<PARAM_NAME>s, t, l</PARAM_NAME>
41
<SP>: intermediate matrices for successive iterations;</SP>
45
<PARAM_NAME>rt, tt</PARAM_NAME>
47
<SP>: gain matrices of the filter model after <VERB>n</VERB> iterations.</SP>
51
<PARAM_NAME>p, s, t, l</PARAM_NAME>
53
<SP>: may be given as input if more than one recursion is desired (evaluation of intermediate values of <VERB>p</VERB>).</SP>
60
square-root algorithm for the algebraic Riccati equation.</P>
65
x=%pi/10:%pi/10:102.4*%pi;
66
rand('seed',0);rand('normal');
67
y=[1;1]*sin(x)+[sin(2*x);sin(1.9*x)]+rand(2,1024);
68
//COMPUTE CORRELATIONS
69
c=[];for j=1:2,for k=1:2,c=[c;corr(y(k,:),y(j,:),64)];end;end
71
//FINDING H,F,G with 6 states
76
[P,s,t,l,Rt,Tt]=srfaur(H,F,G,r0,200);
77
//Make covariance matrix exactly symetric
89
<LINK>lindquist</LINK>