root/library/doc/html/classbdm_1_1KalmanCh.html @ 477

Revision 472, 36.6 kB (checked in by mido, 15 years ago)

obnovena dokumentace, vcetne pridani mnoha novych doc-souboru do svn

Line 
1<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
2<html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
3<title>mixpp: bdm::KalmanCh Class Reference</title>
4<link href="tabs.css" rel="stylesheet" type="text/css">
5<link href="doxygen.css" rel="stylesheet" type="text/css">
6</head><body>
7<!-- Generated by Doxygen 1.5.9 -->
8<script type="text/javascript">
9<!--
10function changeDisplayState (e){
11  var num=this.id.replace(/[^[0-9]/g,'');
12  var button=this.firstChild;
13  var sectionDiv=document.getElementById('dynsection'+num);
14  if (sectionDiv.style.display=='none'||sectionDiv.style.display==''){
15    sectionDiv.style.display='block';
16    button.src='open.gif';
17  }else{
18    sectionDiv.style.display='none';
19    button.src='closed.gif';
20  }
21}
22function initDynSections(){
23  var divs=document.getElementsByTagName('div');
24  var sectionCounter=1;
25  for(var i=0;i<divs.length-1;i++){
26    if(divs[i].className=='dynheader'&&divs[i+1].className=='dynsection'){
27      var header=divs[i];
28      var section=divs[i+1];
29      var button=header.firstChild;
30      if (button!='IMG'){
31        divs[i].insertBefore(document.createTextNode(' '),divs[i].firstChild);
32        button=document.createElement('img');
33        divs[i].insertBefore(button,divs[i].firstChild);
34      }
35      header.style.cursor='pointer';
36      header.onclick=changeDisplayState;
37      header.id='dynheader'+sectionCounter;
38      button.src='closed.gif';
39      section.id='dynsection'+sectionCounter;
40      section.style.display='none';
41      section.style.marginLeft='14px';
42      sectionCounter++;
43    }
44  }
45}
46window.onload = initDynSections;
47-->
48</script>
49<div class="navigation" id="top">
50  <div class="tabs">
51    <ul>
52      <li><a href="main.html"><span>Main&nbsp;Page</span></a></li>
53      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
54      <li class="current"><a href="annotated.html"><span>Classes</span></a></li>
55      <li><a href="files.html"><span>Files</span></a></li>
56    </ul>
57  </div>
58  <div class="tabs">
59    <ul>
60      <li><a href="annotated.html"><span>Class&nbsp;List</span></a></li>
61      <li><a href="classes.html"><span>Class&nbsp;Index</span></a></li>
62      <li><a href="hierarchy.html"><span>Class&nbsp;Hierarchy</span></a></li>
63      <li><a href="functions.html"><span>Class&nbsp;Members</span></a></li>
64    </ul>
65  </div>
66  <div class="navpath"><b>bdm</b>::<a class="el" href="classbdm_1_1KalmanCh.html">KalmanCh</a>
67  </div>
68</div>
69<div class="contents">
70<h1>bdm::KalmanCh Class Reference</h1><!-- doxytag: class="bdm::KalmanCh" --><!-- doxytag: inherits="Kalman&lt; chmat &gt;" --><a class="el" href="classbdm_1_1Kalman.html" title="Kalman filter with covariance matrices in square root form.">Kalman</a> filter in square <a class="el" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a> form. 
71<a href="#_details">More...</a>
72<p>
73<code>#include &lt;<a class="el" href="kalman_8h_source.html">kalman.h</a>&gt;</code>
74<p>
75
76<p>
77<a href="classbdm_1_1KalmanCh-members.html">List of all members.</a><table border="0" cellpadding="0" cellspacing="0">
78<tr><td></td></tr>
79<tr><td colspan="2"><br><h2>Public Member Functions</h2></td></tr>
80<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="24ce65bdaa538d4d5153d709a929b996"></a><!-- doxytag: member="bdm::KalmanCh::_copy_" ref="24ce65bdaa538d4d5153d709a929b996" args="() const " -->
81<a class="el" href="classbdm_1_1BM.html">BM</a> *&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1KalmanCh.html#24ce65bdaa538d4d5153d709a929b996">_copy_</a> () const </td></tr>
82
83<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">copy constructor <br></td></tr>
84<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="20a4d4c664e8ac8a3f1bb7b0d11c6d87"></a><!-- doxytag: member="bdm::KalmanCh::set_parameters" ref="20a4d4c664e8ac8a3f1bb7b0d11c6d87" args="(const mat &amp;A0, const mat &amp;B0, const mat &amp;C0, const mat &amp;D0, const chmat &amp;Q0, const chmat &amp;R0)" -->
85void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1KalmanCh.html#20a4d4c664e8ac8a3f1bb7b0d11c6d87">set_parameters</a> (const mat &amp;A0, const mat &amp;B0, const mat &amp;C0, const mat &amp;D0, const <a class="el" href="classchmat.html">chmat</a> &amp;Q0, const <a class="el" href="classchmat.html">chmat</a> &amp;R0)</td></tr>
86
87<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Set parameters with check of relevance. <br></td></tr>
88<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="6e169272657ed101f3d128b49c59b890"></a><!-- doxytag: member="bdm::KalmanCh::set_statistics" ref="6e169272657ed101f3d128b49c59b890" args="(const vec &amp;mu0, const chmat &amp;P0)" -->
89void&nbsp;</td><td class="memItemRight" valign="bottom"><b>set_statistics</b> (const vec &amp;mu0, const <a class="el" href="classchmat.html">chmat</a> &amp;P0)</td></tr>
90
91<tr><td class="memItemLeft" nowrap align="right" valign="top">void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1KalmanCh.html#b41fe5540548100b08e1684c3be767b6">bayes</a> (const vec &amp;dt)</td></tr>
92
93<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Here dt = [yt;ut] of appropriate dimensions.  <a href="#b41fe5540548100b08e1684c3be767b6"></a><br></td></tr>
94<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="9264fc6b173ecb803d2684b883f32c68"></a><!-- doxytag: member="bdm::KalmanCh::set_est" ref="9264fc6b173ecb803d2684b883f32c68" args="(const vec &amp;mu0, const chmat &amp;P0)" -->
95void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#9264fc6b173ecb803d2684b883f32c68">set_est</a> (const vec &amp;mu0, const <a class="el" href="classchmat.html">chmat</a> &amp;P0)</td></tr>
96
97<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Set estimate values, used e.g. in initialization. <br></td></tr>
98<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="f75e487ff6c129d7012d702030f8c890"></a><!-- doxytag: member="bdm::KalmanCh::posterior" ref="f75e487ff6c129d7012d702030f8c890" args="() const" -->
99const <a class="el" href="classbdm_1_1epdf.html">epdf</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#f75e487ff6c129d7012d702030f8c890">posterior</a> () const</td></tr>
100
101<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">access function <br></td></tr>
102<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c34989b1e53c7d4ecdaea63a95ddbd77"></a><!-- doxytag: member="bdm::KalmanCh::_e" ref="c34989b1e53c7d4ecdaea63a95ddbd77" args="() const" -->
103const <a class="el" href="classbdm_1_1enorm.html">enorm</a>&lt; <a class="el" href="classchmat.html">chmat</a> &gt; *&nbsp;</td><td class="memItemRight" valign="bottom"><b>_e</b> () const</td></tr>
104
105<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c788ec6e6c6f5f5861ae8a56d8ede277"></a><!-- doxytag: member="bdm::KalmanCh::__K" ref="c788ec6e6c6f5f5861ae8a56d8ede277" args="()" -->
106mat &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#c788ec6e6c6f5f5861ae8a56d8ede277">__K</a> ()</td></tr>
107
108<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">access function <br></td></tr>
109<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="a250d1dbe7bba861dba2a324520cfa48"></a><!-- doxytag: member="bdm::KalmanCh::_dP" ref="a250d1dbe7bba861dba2a324520cfa48" args="()" -->
110vec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#a250d1dbe7bba861dba2a324520cfa48">_dP</a> ()</td></tr>
111
112<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">access function <br></td></tr>
113<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="af0e6941ba5919cae1bc530e2181b109"></a><!-- doxytag: member="bdm::KalmanCh::to_string" ref="af0e6941ba5919cae1bc530e2181b109" args="()" -->
114virtual string&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1root.html#af0e6941ba5919cae1bc530e2181b109">to_string</a> ()</td></tr>
115
116<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">This method returns a basic info about the current instance. <br></td></tr>
117<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="0551e3121091c5199bf4413b50522176"></a><!-- doxytag: member="bdm::KalmanCh::from_setting" ref="0551e3121091c5199bf4413b50522176" args="(const Setting &amp;set)" -->
118virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1root.html#0551e3121091c5199bf4413b50522176">from_setting</a> (const Setting &amp;set)</td></tr>
119
120<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">This method arrange instance properties according the data stored in the Setting structure. <br></td></tr>
121<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="67d954d255ede776eade7334d4895790"></a><!-- doxytag: member="bdm::KalmanCh::to_setting" ref="67d954d255ede776eade7334d4895790" args="(Setting &amp;set) const " -->
122virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1root.html#67d954d255ede776eade7334d4895790">to_setting</a> (Setting &amp;set) const </td></tr>
123
124<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">This method save all the instance properties into the Setting structure. <br></td></tr>
125<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="1c314bd6d6dacb8ba78ea5eb88fd9516"></a><!-- doxytag: member="bdm::KalmanCh::validate" ref="1c314bd6d6dacb8ba78ea5eb88fd9516" args="()" -->
126virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1root.html#1c314bd6d6dacb8ba78ea5eb88fd9516">validate</a> ()</td></tr>
127
128<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">This method TODO. <br></td></tr>
129<tr><td colspan="2"><div class="groupHeader">Mathematical operations</div></td></tr>
130<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="1dee3fddaf021e62d925289660a707dc"></a><!-- doxytag: member="bdm::KalmanCh::bayesB" ref="1dee3fddaf021e62d925289660a707dc" args="(const mat &amp;Dt)" -->
131virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#1dee3fddaf021e62d925289660a707dc">bayesB</a> (const mat &amp;Dt)</td></tr>
132
133<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Batch Bayes rule (columns of Dt are observations). <br></td></tr>
134<tr><td class="memItemLeft" nowrap align="right" valign="top">virtual double&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#50257e0c1e5b5c73153ea6e716ad8ae0">logpred</a> (const vec &amp;dt) const </td></tr>
135
136<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="0e8ebe61fb14990abe1254bd3dda5fae"></a><!-- doxytag: member="bdm::KalmanCh::logpred_m" ref="0e8ebe61fb14990abe1254bd3dda5fae" args="(const mat &amp;dt) const " -->
137vec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#0e8ebe61fb14990abe1254bd3dda5fae">logpred_m</a> (const mat &amp;dt) const </td></tr>
138
139<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix version of logpred. <br></td></tr>
140<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="688d7a2aced1e06aa1c468d73a9e5eba"></a><!-- doxytag: member="bdm::KalmanCh::epredictor" ref="688d7a2aced1e06aa1c468d73a9e5eba" args="() const " -->
141virtual <a class="el" href="classbdm_1_1epdf.html">epdf</a> *&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#688d7a2aced1e06aa1c468d73a9e5eba">epredictor</a> () const </td></tr>
142
143<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Constructs a predictive density <img class="formulaInl" alt="$ f(d_{t+1} |d_{t}, \ldots d_{0}) $" src="form_14.png">. <br></td></tr>
144<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="598b25e3f3d96a5bc00a5faeb5b3c912"></a><!-- doxytag: member="bdm::KalmanCh::predictor" ref="598b25e3f3d96a5bc00a5faeb5b3c912" args="() const " -->
145virtual <a class="el" href="classbdm_1_1mpdf.html">mpdf</a> *&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#598b25e3f3d96a5bc00a5faeb5b3c912">predictor</a> () const </td></tr>
146
147<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Constructs a conditional density 1-step ahead predictor. <br></td></tr>
148<tr><td colspan="2"><div class="groupHeader">Access to attributes</div></td></tr>
149<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="ff2d8755ba0b3def927d31305c03b09c"></a><!-- doxytag: member="bdm::KalmanCh::_drv" ref="ff2d8755ba0b3def927d31305c03b09c" args="() const " -->
150const <a class="el" href="classbdm_1_1RV.html">RV</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><b>_drv</b> () const </td></tr>
151
152<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="f135ae6dce7e9f30c9f88229c7930b96"></a><!-- doxytag: member="bdm::KalmanCh::set_drv" ref="f135ae6dce7e9f30c9f88229c7930b96" args="(const RV &amp;rv)" -->
153void&nbsp;</td><td class="memItemRight" valign="bottom"><b>set_drv</b> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &amp;rv)</td></tr>
154
155<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="b38d92f17620813ad872d86e01a26e5e"></a><!-- doxytag: member="bdm::KalmanCh::set_rv" ref="b38d92f17620813ad872d86e01a26e5e" args="(const RV &amp;rv)" -->
156void&nbsp;</td><td class="memItemRight" valign="bottom"><b>set_rv</b> (const <a class="el" href="classbdm_1_1RV.html">RV</a> &amp;rv)</td></tr>
157
158<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5be65d37dedfe33a3671e7065f523a70"></a><!-- doxytag: member="bdm::KalmanCh::_ll" ref="5be65d37dedfe33a3671e7065f523a70" args="() const " -->
159double&nbsp;</td><td class="memItemRight" valign="bottom"><b>_ll</b> () const </td></tr>
160
161<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="236b3abbcc93594fc97cd86d82c1a83f"></a><!-- doxytag: member="bdm::KalmanCh::set_evalll" ref="236b3abbcc93594fc97cd86d82c1a83f" args="(bool evl0)" -->
162void&nbsp;</td><td class="memItemRight" valign="bottom"><b>set_evalll</b> (bool evl0)</td></tr>
163
164<tr><td colspan="2"><br><h2>Protected Attributes</h2></td></tr>
165<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="48611c8582706cfa62e832be0972e75d"></a><!-- doxytag: member="bdm::KalmanCh::preA" ref="48611c8582706cfa62e832be0972e75d" args="" -->
166mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1KalmanCh.html#48611c8582706cfa62e832be0972e75d">preA</a></td></tr>
167
168<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">pre array (triangular matrix) <br></td></tr>
169<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="bcbd68f51d4b57246e7784ca5900171f"></a><!-- doxytag: member="bdm::KalmanCh::postA" ref="bcbd68f51d4b57246e7784ca5900171f" args="" -->
170mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1KalmanCh.html#bcbd68f51d4b57246e7784ca5900171f">postA</a></td></tr>
171
172<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">post array (triangular matrix) <br></td></tr>
173<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="3fe475a1e920b20b63bb342c0e1571f7"></a><!-- doxytag: member="bdm::KalmanCh::rvy" ref="3fe475a1e920b20b63bb342c0e1571f7" args="" -->
174<a class="el" href="classbdm_1_1RV.html">RV</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#3fe475a1e920b20b63bb342c0e1571f7">rvy</a></td></tr>
175
176<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Indetifier of output rv. <br></td></tr>
177<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="149e27424fd1a7cc1c998ea088618a94"></a><!-- doxytag: member="bdm::KalmanCh::rvu" ref="149e27424fd1a7cc1c998ea088618a94" args="" -->
178<a class="el" href="classbdm_1_1RV.html">RV</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#149e27424fd1a7cc1c998ea088618a94">rvu</a></td></tr>
179
180<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Indetifier of exogeneous rv. <br></td></tr>
181<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="ba7699cdb3b1382a54d3e28b9b7517fa"></a><!-- doxytag: member="bdm::KalmanCh::dimx" ref="ba7699cdb3b1382a54d3e28b9b7517fa" args="" -->
182int&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#ba7699cdb3b1382a54d3e28b9b7517fa">dimx</a></td></tr>
183
184<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of rv.count() <br></td></tr>
185<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="d2c36ba01760bf207b985bf321b7817f"></a><!-- doxytag: member="bdm::KalmanCh::dimy" ref="d2c36ba01760bf207b985bf321b7817f" args="" -->
186int&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#d2c36ba01760bf207b985bf321b7817f">dimy</a></td></tr>
187
188<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of rvy.count() <br></td></tr>
189<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c5136ef617f6ac0e426bea222755d92b"></a><!-- doxytag: member="bdm::KalmanCh::dimu" ref="c5136ef617f6ac0e426bea222755d92b" args="" -->
190int&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#c5136ef617f6ac0e426bea222755d92b">dimu</a></td></tr>
191
192<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of rvu.count() <br></td></tr>
193<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="0a2072e2090c10fac74ad30a023a4ace"></a><!-- doxytag: member="bdm::KalmanCh::A" ref="0a2072e2090c10fac74ad30a023a4ace" args="" -->
194mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#0a2072e2090c10fac74ad30a023a4ace">A</a></td></tr>
195
196<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix A. <br></td></tr>
197<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5977b2c81857948a35105f0e7840203c"></a><!-- doxytag: member="bdm::KalmanCh::B" ref="5977b2c81857948a35105f0e7840203c" args="" -->
198mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#5977b2c81857948a35105f0e7840203c">B</a></td></tr>
199
200<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix B. <br></td></tr>
201<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="818eba63a23972786a4579ad30294177"></a><!-- doxytag: member="bdm::KalmanCh::C" ref="818eba63a23972786a4579ad30294177" args="" -->
202mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#818eba63a23972786a4579ad30294177">C</a></td></tr>
203
204<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix C. <br></td></tr>
205<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="7b56ac423d0654b5755e4f852a870456"></a><!-- doxytag: member="bdm::KalmanCh::D" ref="7b56ac423d0654b5755e4f852a870456" args="" -->
206mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#7b56ac423d0654b5755e4f852a870456">D</a></td></tr>
207
208<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix D. <br></td></tr>
209<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="70f8bf19e81b532c60fd3a7a152425ee"></a><!-- doxytag: member="bdm::KalmanCh::Q" ref="70f8bf19e81b532c60fd3a7a152425ee" args="" -->
210<a class="el" href="classchmat.html">chmat</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#70f8bf19e81b532c60fd3a7a152425ee">Q</a></td></tr>
211
212<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix Q in square-root form. <br></td></tr>
213<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="475b088287cdfbba4dc60a3d027728b7"></a><!-- doxytag: member="bdm::KalmanCh::R" ref="475b088287cdfbba4dc60a3d027728b7" args="" -->
214<a class="el" href="classchmat.html">chmat</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#475b088287cdfbba4dc60a3d027728b7">R</a></td></tr>
215
216<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix R in square-root form. <br></td></tr>
217<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="383f329ff18bbe219254c8b3b916f40d"></a><!-- doxytag: member="bdm::KalmanCh::est" ref="383f329ff18bbe219254c8b3b916f40d" args="" -->
218<a class="el" href="classbdm_1_1enorm.html">enorm</a>&lt; <a class="el" href="classchmat.html">chmat</a> &gt;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#383f329ff18bbe219254c8b3b916f40d">est</a></td></tr>
219
220<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">posterior density on $x_t$ <br></td></tr>
221<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="ba555c394c429f6831c9bbabfa2c944c"></a><!-- doxytag: member="bdm::KalmanCh::fy" ref="ba555c394c429f6831c9bbabfa2c944c" args="" -->
222<a class="el" href="classbdm_1_1enorm.html">enorm</a>&lt; <a class="el" href="classchmat.html">chmat</a> &gt;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#ba555c394c429f6831c9bbabfa2c944c">fy</a></td></tr>
223
224<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">preditive density on $y_t$ <br></td></tr>
225<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="bd69dfb802465f22dd84d73a180d5c92"></a><!-- doxytag: member="bdm::KalmanCh::_K" ref="bd69dfb802465f22dd84d73a180d5c92" args="" -->
226mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#bd69dfb802465f22dd84d73a180d5c92">_K</a></td></tr>
227
228<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">placeholder for Kalman gain <br></td></tr>
229<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c249d45258c8578b13858ad3e7b729b1"></a><!-- doxytag: member="bdm::KalmanCh::_yp" ref="c249d45258c8578b13858ad3e7b729b1" args="" -->
230vec &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#c249d45258c8578b13858ad3e7b729b1">_yp</a></td></tr>
231
232<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of fy.mu <br></td></tr>
233<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="2dd268f2d7fbe6382cb8825a1114192a"></a><!-- doxytag: member="bdm::KalmanCh::_Ry" ref="2dd268f2d7fbe6382cb8825a1114192a" args="" -->
234<a class="el" href="classchmat.html">chmat</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#2dd268f2d7fbe6382cb8825a1114192a">_Ry</a></td></tr>
235
236<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of fy.R <br></td></tr>
237<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="fa172078091e45561343fa513dd573b0"></a><!-- doxytag: member="bdm::KalmanCh::_mu" ref="fa172078091e45561343fa513dd573b0" args="" -->
238vec &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#fa172078091e45561343fa513dd573b0">_mu</a></td></tr>
239
240<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of est.mu <br></td></tr>
241<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="00c27b0bf324f0018497921ca23c71ed"></a><!-- doxytag: member="bdm::KalmanCh::_P" ref="00c27b0bf324f0018497921ca23c71ed" args="" -->
242<a class="el" href="classchmat.html">chmat</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#00c27b0bf324f0018497921ca23c71ed">_P</a></td></tr>
243
244<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of est.R <br></td></tr>
245<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c400357e37d27a4834b2b1d9211009ed"></a><!-- doxytag: member="bdm::KalmanCh::drv" ref="c400357e37d27a4834b2b1d9211009ed" args="" -->
246<a class="el" href="classbdm_1_1RV.html">RV</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed">drv</a></td></tr>
247
248<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Random variable of the data (optional). <br></td></tr>
249<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="4064b6559d962633e4372b12f4cd204a"></a><!-- doxytag: member="bdm::KalmanCh::ll" ref="4064b6559d962633e4372b12f4cd204a" args="" -->
250double&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a">ll</a></td></tr>
251
252<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Logarithm of marginalized data likelihood. <br></td></tr>
253<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="faff0ad12556fe7dc0e2807d4fd938ee"></a><!-- doxytag: member="bdm::KalmanCh::evalll" ref="faff0ad12556fe7dc0e2807d4fd938ee" args="" -->
254bool&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee">evalll</a></td></tr>
255
256<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">If true, the filter will compute likelihood of the data record and store it in <code>ll</code> . Set to false if you want to save computational time. <br></td></tr>
257<tr><td colspan="2"><br><h2>Extension to conditional BM</h2></td></tr>
258<tr><td colspan="2">This extension is useful e.g. in Marginalized Particle Filter (<a class="el" href="classbdm_1_1MPF.html">bdm::MPF</a>). Alternatively, it can be used for automated connection to <a class="el" href="classbdm_1_1DS.html" title="Abstract class for discrete-time sources of data.">DS</a> when the condition is observed <br><br></td></tr>
259<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="b5e8e016be16e4029c789063cb3f3b33"></a><!-- doxytag: member="bdm::KalmanCh::_rvc" ref="b5e8e016be16e4029c789063cb3f3b33" args="() const " -->
260const <a class="el" href="classbdm_1_1RV.html">RV</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#b5e8e016be16e4029c789063cb3f3b33">_rvc</a> () const </td></tr>
261
262<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">access function <br></td></tr>
263<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="6799f4b16a6a59ed58b1d0d6e17116f4"></a><!-- doxytag: member="bdm::KalmanCh::condition" ref="6799f4b16a6a59ed58b1d0d6e17116f4" args="(const vec &amp;val)" -->
264virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#6799f4b16a6a59ed58b1d0d6e17116f4">condition</a> (const vec &amp;val)</td></tr>
265
266<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Substitute <code>val</code> for <code>rvc</code>. <br></td></tr>
267<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="de95c590d650b2db423d53992f0906cb"></a><!-- doxytag: member="bdm::KalmanCh::rvc" ref="de95c590d650b2db423d53992f0906cb" args="" -->
268<a class="el" href="classbdm_1_1RV.html">RV</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#de95c590d650b2db423d53992f0906cb">rvc</a></td></tr>
269
270<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Name of extension variable. <br></td></tr>
271<tr><td colspan="2"><br><h2>Logging of results</h2></td></tr>
272<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="adb539fa2585c8a2a91c54e5aebf330e"></a><!-- doxytag: member="bdm::KalmanCh::set_options" ref="adb539fa2585c8a2a91c54e5aebf330e" args="(const string &amp;opt)" -->
273virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#adb539fa2585c8a2a91c54e5aebf330e">set_options</a> (const string &amp;opt)</td></tr>
274
275<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Set boolean options from a string, recognized are: "logbounds,logll". <br></td></tr>
276<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="2298ca6af9b13a78a4c9e18ab43f1827"></a><!-- doxytag: member="bdm::KalmanCh::log_add" ref="2298ca6af9b13a78a4c9e18ab43f1827" args="(logger &amp;L, const string &amp;name=&quot;&quot;)" -->
277virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#2298ca6af9b13a78a4c9e18ab43f1827">log_add</a> (<a class="el" href="classbdm_1_1logger.html">logger</a> &amp;L, const string &amp;name=&quot;&quot;)</td></tr>
278
279<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Add all logged variables to a <a class="el" href="classbdm_1_1logger.html" title="Class for storing results (and semi-results) of an experiment.">logger</a>. <br></td></tr>
280<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="b517e1679eaa94e803ea4cd0b8efbcd7"></a><!-- doxytag: member="bdm::KalmanCh::logit" ref="b517e1679eaa94e803ea4cd0b8efbcd7" args="(logger &amp;L)" -->
281virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><b>logit</b> (<a class="el" href="classbdm_1_1logger.html">logger</a> &amp;L)</td></tr>
282
283<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="109c1a626a69031658e3a44e9e500cca"></a><!-- doxytag: member="bdm::KalmanCh::LIDs" ref="109c1a626a69031658e3a44e9e500cca" args="" -->
284ivec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca">LIDs</a></td></tr>
285
286<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll]. <br></td></tr>
287<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="398b7e7a162cf37d0022d6126495f207"></a><!-- doxytag: member="bdm::KalmanCh::LFlags" ref="398b7e7a162cf37d0022d6126495f207" args="" -->
288ivec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207">LFlags</a></td></tr>
289
290<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Flags for logging - same size as LIDs, each entry correspond to the same in LIDs. <br></td></tr>
291</table>
292<hr><a name="_details"></a><h2>Detailed Description</h2>
293<a class="el" href="classbdm_1_1Kalman.html" title="Kalman filter with covariance matrices in square root form.">Kalman</a> filter in square <a class="el" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a> form.
294<p>
295Trivial example: <div class="fragment"><pre class="fragment"><span class="preprocessor">#include "<a class="code" href="kalman_8h.html" title="Bayesian Filtering for linear Gaussian models (Kalman Filter) and extensions.">estim/kalman.h</a>"</span>
296<span class="keyword">using namespace </span>bdm;
297       
298<span class="comment">// estimation of AR(0) model</span>
299<span class="keywordtype">int</span> main() {
300        <span class="comment">//dimensions</span>
301        <span class="keywordtype">int</span> dx=3, dy=3, du=1;
302        <span class="comment">// matrices</span>
303        mat A = eye(dx);
304        mat B = zeros(dx,du);
305        mat C = eye(dx);
306        mat D = zeros(dy,du);
307        mat Q = eye(dx);
308        mat R = 0.1*eye(dy);
309        <span class="comment">//prior</span>
310        mat P0 = 100*eye(dx);
311        vec mu0 = zeros(dx);
312        <span class="comment">// Estimator</span>
313        KalmanCh KF;
314        KF.set_parameters(A,B,C,D,<span class="comment">/*covariances*/</span> Q,R);
315        KF.set_statistics(mu0,P0);
316        <span class="comment">// Estimation loop</span>
317        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0;i&lt;100;i++){
318                KF.bayes(randn(dx+du));
319        }
320        <span class="comment">//print results</span>
321        cout &lt;&lt; <span class="stringliteral">"Posterior estimate of x is: "</span>  &lt;&lt; endl;
322        cout &lt;&lt; <span class="stringliteral">"mean: "</span>&lt;&lt; KF.posterior().mean()&lt;&lt; endl;
323        cout &lt;&lt; <span class="stringliteral">"variance: "</span>&lt;&lt; KF.posterior().variance()&lt;&lt; endl;
324}
325</pre></div> <hr><h2>Member Function Documentation</h2>
326<a class="anchor" name="b41fe5540548100b08e1684c3be767b6"></a><!-- doxytag: member="bdm::KalmanCh::bayes" ref="b41fe5540548100b08e1684c3be767b6" args="(const vec &amp;dt)" -->
327<div class="memitem">
328<div class="memproto">
329      <table class="memname">
330        <tr>
331          <td class="memname">void bdm::KalmanCh::bayes           </td>
332          <td>(</td>
333          <td class="paramtype">const vec &amp;&nbsp;</td>
334          <td class="paramname"> <em>dt</em>          </td>
335          <td>&nbsp;)&nbsp;</td>
336          <td><code> [virtual]</code></td>
337        </tr>
338      </table>
339</div>
340<div class="memdoc">
341
342<p>
343Here dt = [yt;ut] of appropriate dimensions.
344<p>
345The following equality hold::<p class="formulaDsp">
346<img class="formulaDsp" alt="\[ \left[\begin{array}{cc} R^{0.5}\\ P_{t|t-1}^{0.5}C' &amp; P_{t|t-1}^{0.5}CA'\\ &amp; Q^{0.5}\end{array}\right]<\mathrm{orth.oper.}>=\left[\begin{array}{cc} R_{y}^{0.5} &amp; KA'\\ &amp; P_{t+1|t}^{0.5}\\ \\\end{array}\right]\]" src="form_31.png">
347<p>
348<p>
349Thus this object evaluates only predictors! Not filtering densities.
350<p>Reimplemented from <a class="el" href="classbdm_1_1Kalman.html#4a39330c14eff8d13179e868a1d1aa8c">bdm::Kalman&lt; chmat &gt;</a>.</p>
351
352<p>Reimplemented in <a class="el" href="classbdm_1_1EKFCh.html#4c8609c37290b158f88a31dae4047225">bdm::EKFCh</a>.</p>
353
354<p>References <a class="el" href="chmat_8h_source.html#l00073">chmat::_Ch()</a>, <a class="el" href="kalman_8h_source.html#l00099">bdm::Kalman&lt; chmat &gt;::_K</a>, <a class="el" href="kalman_8h_source.html#l00105">bdm::Kalman&lt; chmat &gt;::_mu</a>, <a class="el" href="kalman_8h_source.html#l00107">bdm::Kalman&lt; chmat &gt;::_P</a>, <a class="el" href="kalman_8h_source.html#l00103">bdm::Kalman&lt; chmat &gt;::_Ry</a>, <a class="el" href="kalman_8h_source.html#l00101">bdm::Kalman&lt; chmat &gt;::_yp</a>, <a class="el" href="kalman_8h_source.html#l00081">bdm::Kalman&lt; chmat &gt;::A</a>, <a class="el" href="kalman_8h_source.html#l00083">bdm::Kalman&lt; chmat &gt;::B</a>, <a class="el" href="kalman_8h_source.html#l00085">bdm::Kalman&lt; chmat &gt;::C</a>, <a class="el" href="kalman_8h_source.html#l00087">bdm::Kalman&lt; chmat &gt;::D</a>, <a class="el" href="kalman_8h_source.html#l00079">bdm::Kalman&lt; chmat &gt;::dimu</a>, <a class="el" href="kalman_8h_source.html#l00075">bdm::Kalman&lt; chmat &gt;::dimx</a>, <a class="el" href="kalman_8h_source.html#l00077">bdm::Kalman&lt; chmat &gt;::dimy</a>, <a class="el" href="bdmbase_8h_source.html#l00741">bdm::BM::evalll</a>, <a class="el" href="exp__family_8h_source.html#l00051">bdm::eEF::evallog()</a>, <a class="el" href="kalman_8h_source.html#l00096">bdm::Kalman&lt; chmat &gt;::fy</a>, <a class="el" href="bdmbase_8h_source.html#l00739">bdm::BM::ll</a>, <a class="el" href="kalman_8h_source.html#l00148">postA</a>, and <a class="el" href="kalman_8h_source.html#l00146">preA</a>.</p>
355
356</div>
357</div><p>
358<a class="anchor" name="50257e0c1e5b5c73153ea6e716ad8ae0"></a><!-- doxytag: member="bdm::KalmanCh::logpred" ref="50257e0c1e5b5c73153ea6e716ad8ae0" args="(const vec &amp;dt) const " -->
359<div class="memitem">
360<div class="memproto">
361      <table class="memname">
362        <tr>
363          <td class="memname">virtual double bdm::BM::logpred           </td>
364          <td>(</td>
365          <td class="paramtype">const vec &amp;&nbsp;</td>
366          <td class="paramname"> <em>dt</em>          </td>
367          <td>&nbsp;)&nbsp;</td>
368          <td> const<code> [inline, virtual, inherited]</code></td>
369        </tr>
370      </table>
371</div>
372<div class="memdoc">
373
374<p>
375Evaluates predictive log-likelihood of the given data record I.e. marginal likelihood of the data with the posterior integrated out.
376<p>Reimplemented in <a class="el" href="classbdm_1_1ARX.html#080a7e531e3aa06694112863b15bc6a4">bdm::ARX</a>, <a class="el" href="classbdm_1_1MixEF.html#da724da464a75e07521941e430929efa">bdm::MixEF</a>, and <a class="el" href="classbdm_1_1multiBM.html#e157b607c1e3fa91d42aeea44458e2bf">bdm::multiBM</a>.</p>
377
378<p>Referenced by <a class="el" href="bdmbase_8h_source.html#l00770">bdm::BM::logpred_m()</a>.</p>
379
380</div>
381</div><p>
382<hr>The documentation for this class was generated from the following files:<ul>
383<li><a class="el" href="kalman_8h_source.html">kalman.h</a><li>kalman.cpp</ul>
384</div>
385<hr size="1"><address style="text-align: right;"><small>Generated on Wed Aug 5 00:06:58 2009 for mixpp by&nbsp;
386<a href="http://www.doxygen.org/index.html">
387<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.9 </small></address>
388</body>
389</html>
Note: See TracBrowser for help on using the browser.