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="faurre"> filter computation by simple Faurre algorithm</SHORT_DESCRIPTION>
10
<CALLING_SEQUENCE_ITEM>[P,R,T]=faurre(n,H,F,G,R0) </CALLING_SEQUENCE_ITEM>
15
<PARAM_NAME>n</PARAM_NAME>
17
<SP>: number of iterations.</SP>
21
<PARAM_NAME>H, F, G</PARAM_NAME>
23
<SP>: estimated triple from the covariance sequence of <VERB>y</VERB>.</SP>
27
<PARAM_NAME>R0</PARAM_NAME>
33
<PARAM_NAME>P</PARAM_NAME>
35
<SP>: solution of the Riccati equation after n iterations.</SP>
39
<PARAM_NAME>R, T</PARAM_NAME>
41
<SP>: gain matrix of the filter.</SP>
48
This function computes iteratively the minimal solution of the algebraic
49
Riccati equation and gives the matrices <VERB>R</VERB> and <VERB>T</VERB> of the
51
The algorithm tries to compute the solution P as the growing limit of a
52
sequence of matrices Pn such that</P>
56
Pn+1=F*Pn*F'+(G-F*Pn*h')*(R0-H*Pn*H') *(G'-H*Pn*F')
62
Note that this method may not converge,especially when F has poles
63
near the unit circle. Use preferably the srfaur function.</P>
70
<LINK>lindquist</LINK>
76
<AUTHOR>G. Le V. </AUTHOR>