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

Revision 538, 36.4 kB (checked in by smidl, 15 years ago)

Documentation regenerated

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.8 -->
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="classbdm_1_1chmat.html">chmat</a> &amp;Q0, const <a class="el" href="classbdm_1_1chmat.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="classbdm_1_1chmat.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="classbdm_1_1chmat.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="c788ec6e6c6f5f5861ae8a56d8ede277"></a><!-- doxytag: member="bdm::KalmanCh::__K" ref="c788ec6e6c6f5f5861ae8a56d8ede277" args="()" -->
103mat &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#c788ec6e6c6f5f5861ae8a56d8ede277">__K</a> ()</td></tr>
104
105<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">access function <br></td></tr>
106<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="a250d1dbe7bba861dba2a324520cfa48"></a><!-- doxytag: member="bdm::KalmanCh::_dP" ref="a250d1dbe7bba861dba2a324520cfa48" args="()" -->
107vec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#a250d1dbe7bba861dba2a324520cfa48">_dP</a> ()</td></tr>
108
109<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">access function <br></td></tr>
110<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="af0e6941ba5919cae1bc530e2181b109"></a><!-- doxytag: member="bdm::KalmanCh::to_string" ref="af0e6941ba5919cae1bc530e2181b109" args="()" -->
111virtual string&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1root.html#af0e6941ba5919cae1bc530e2181b109">to_string</a> ()</td></tr>
112
113<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">This method returns a basic info about the current instance. <br></td></tr>
114<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)" -->
115virtual 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>
116
117<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>
118<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 " -->
119virtual 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>
120
121<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">This method save all the instance properties into the Setting structure. <br></td></tr>
122<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="1c314bd6d6dacb8ba78ea5eb88fd9516"></a><!-- doxytag: member="bdm::KalmanCh::validate" ref="1c314bd6d6dacb8ba78ea5eb88fd9516" args="()" -->
123virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1root.html#1c314bd6d6dacb8ba78ea5eb88fd9516">validate</a> ()</td></tr>
124
125<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">This method TODO. <br></td></tr>
126<tr><td colspan="2"><div class="groupHeader">Mathematical operations</div></td></tr>
127<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)" -->
128virtual 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>
129
130<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Batch Bayes rule (columns of Dt are observations). <br></td></tr>
131<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>
132
133<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 " -->
134vec&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>
135
136<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix version of logpred. <br></td></tr>
137<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="688d7a2aced1e06aa1c468d73a9e5eba"></a><!-- doxytag: member="bdm::KalmanCh::epredictor" ref="688d7a2aced1e06aa1c468d73a9e5eba" args="() const " -->
138virtual <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>
139
140<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>
141<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="598b25e3f3d96a5bc00a5faeb5b3c912"></a><!-- doxytag: member="bdm::KalmanCh::predictor" ref="598b25e3f3d96a5bc00a5faeb5b3c912" args="() const " -->
142virtual <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>
143
144<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Constructs conditional density of 1-step ahead predictor <img class="formulaInl" alt="$ f(d_{t+1} |d_{t+h-1}, \ldots d_{t}) $" src="form_131.png">. <br></td></tr>
145<tr><td colspan="2"><div class="groupHeader">Access to attributes</div></td></tr>
146<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="ff2d8755ba0b3def927d31305c03b09c"></a><!-- doxytag: member="bdm::KalmanCh::_drv" ref="ff2d8755ba0b3def927d31305c03b09c" args="() const " -->
147const <a class="el" href="classbdm_1_1RV.html">RV</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><b>_drv</b> () const </td></tr>
148
149<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)" -->
150void&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>
151
152<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)" -->
153void&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>
154
155<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5be65d37dedfe33a3671e7065f523a70"></a><!-- doxytag: member="bdm::KalmanCh::_ll" ref="5be65d37dedfe33a3671e7065f523a70" args="() const " -->
156double&nbsp;</td><td class="memItemRight" valign="bottom"><b>_ll</b> () const </td></tr>
157
158<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)" -->
159void&nbsp;</td><td class="memItemRight" valign="bottom"><b>set_evalll</b> (bool evl0)</td></tr>
160
161<tr><td colspan="2"><br><h2>Protected Attributes</h2></td></tr>
162<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="48611c8582706cfa62e832be0972e75d"></a><!-- doxytag: member="bdm::KalmanCh::preA" ref="48611c8582706cfa62e832be0972e75d" args="" -->
163mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1KalmanCh.html#48611c8582706cfa62e832be0972e75d">preA</a></td></tr>
164
165<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">pre array (triangular matrix) <br></td></tr>
166<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="bcbd68f51d4b57246e7784ca5900171f"></a><!-- doxytag: member="bdm::KalmanCh::postA" ref="bcbd68f51d4b57246e7784ca5900171f" args="" -->
167mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1KalmanCh.html#bcbd68f51d4b57246e7784ca5900171f">postA</a></td></tr>
168
169<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">post array (triangular matrix) <br></td></tr>
170<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="3fe475a1e920b20b63bb342c0e1571f7"></a><!-- doxytag: member="bdm::KalmanCh::rvy" ref="3fe475a1e920b20b63bb342c0e1571f7" args="" -->
171<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>
172
173<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Indetifier of output rv. <br></td></tr>
174<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="149e27424fd1a7cc1c998ea088618a94"></a><!-- doxytag: member="bdm::KalmanCh::rvu" ref="149e27424fd1a7cc1c998ea088618a94" args="" -->
175<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>
176
177<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Indetifier of exogeneous rv. <br></td></tr>
178<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="ba7699cdb3b1382a54d3e28b9b7517fa"></a><!-- doxytag: member="bdm::KalmanCh::dimx" ref="ba7699cdb3b1382a54d3e28b9b7517fa" args="" -->
179int&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#ba7699cdb3b1382a54d3e28b9b7517fa">dimx</a></td></tr>
180
181<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of rv.count() <br></td></tr>
182<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="d2c36ba01760bf207b985bf321b7817f"></a><!-- doxytag: member="bdm::KalmanCh::dimy" ref="d2c36ba01760bf207b985bf321b7817f" args="" -->
183int&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#d2c36ba01760bf207b985bf321b7817f">dimy</a></td></tr>
184
185<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of rvy.count() <br></td></tr>
186<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c5136ef617f6ac0e426bea222755d92b"></a><!-- doxytag: member="bdm::KalmanCh::dimu" ref="c5136ef617f6ac0e426bea222755d92b" args="" -->
187int&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#c5136ef617f6ac0e426bea222755d92b">dimu</a></td></tr>
188
189<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of rvu.count() <br></td></tr>
190<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="0a2072e2090c10fac74ad30a023a4ace"></a><!-- doxytag: member="bdm::KalmanCh::A" ref="0a2072e2090c10fac74ad30a023a4ace" args="" -->
191mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#0a2072e2090c10fac74ad30a023a4ace">A</a></td></tr>
192
193<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix A. <br></td></tr>
194<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="5977b2c81857948a35105f0e7840203c"></a><!-- doxytag: member="bdm::KalmanCh::B" ref="5977b2c81857948a35105f0e7840203c" args="" -->
195mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#5977b2c81857948a35105f0e7840203c">B</a></td></tr>
196
197<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix B. <br></td></tr>
198<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="818eba63a23972786a4579ad30294177"></a><!-- doxytag: member="bdm::KalmanCh::C" ref="818eba63a23972786a4579ad30294177" args="" -->
199mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#818eba63a23972786a4579ad30294177">C</a></td></tr>
200
201<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix C. <br></td></tr>
202<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="7b56ac423d0654b5755e4f852a870456"></a><!-- doxytag: member="bdm::KalmanCh::D" ref="7b56ac423d0654b5755e4f852a870456" args="" -->
203mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#7b56ac423d0654b5755e4f852a870456">D</a></td></tr>
204
205<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix D. <br></td></tr>
206<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="70f8bf19e81b532c60fd3a7a152425ee"></a><!-- doxytag: member="bdm::KalmanCh::Q" ref="70f8bf19e81b532c60fd3a7a152425ee" args="" -->
207<a class="el" href="classbdm_1_1chmat.html">chmat</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#70f8bf19e81b532c60fd3a7a152425ee">Q</a></td></tr>
208
209<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix Q in square-root form. <br></td></tr>
210<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="475b088287cdfbba4dc60a3d027728b7"></a><!-- doxytag: member="bdm::KalmanCh::R" ref="475b088287cdfbba4dc60a3d027728b7" args="" -->
211<a class="el" href="classbdm_1_1chmat.html">chmat</a>&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#475b088287cdfbba4dc60a3d027728b7">R</a></td></tr>
212
213<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Matrix R in square-root form. <br></td></tr>
214<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="383f329ff18bbe219254c8b3b916f40d"></a><!-- doxytag: member="bdm::KalmanCh::est" ref="383f329ff18bbe219254c8b3b916f40d" args="" -->
215<a class="el" href="classbdm_1_1enorm.html">enorm</a>&lt; <a class="el" href="classbdm_1_1chmat.html">chmat</a> &gt;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#383f329ff18bbe219254c8b3b916f40d">est</a></td></tr>
216
217<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">posterior density on $x_t$ <br></td></tr>
218<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="ba555c394c429f6831c9bbabfa2c944c"></a><!-- doxytag: member="bdm::KalmanCh::fy" ref="ba555c394c429f6831c9bbabfa2c944c" args="" -->
219<a class="el" href="classbdm_1_1enorm.html">enorm</a>&lt; <a class="el" href="classbdm_1_1chmat.html">chmat</a> &gt;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#ba555c394c429f6831c9bbabfa2c944c">fy</a></td></tr>
220
221<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">preditive density on $y_t$ <br></td></tr>
222<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="bd69dfb802465f22dd84d73a180d5c92"></a><!-- doxytag: member="bdm::KalmanCh::_K" ref="bd69dfb802465f22dd84d73a180d5c92" args="" -->
223mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#bd69dfb802465f22dd84d73a180d5c92">_K</a></td></tr>
224
225<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">placeholder for Kalman gain <br></td></tr>
226<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c249d45258c8578b13858ad3e7b729b1"></a><!-- doxytag: member="bdm::KalmanCh::_yp" ref="c249d45258c8578b13858ad3e7b729b1" args="" -->
227vec &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#c249d45258c8578b13858ad3e7b729b1">_yp</a></td></tr>
228
229<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of fy.mu <br></td></tr>
230<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="2dd268f2d7fbe6382cb8825a1114192a"></a><!-- doxytag: member="bdm::KalmanCh::_Ry" ref="2dd268f2d7fbe6382cb8825a1114192a" args="" -->
231<a class="el" href="classbdm_1_1chmat.html">chmat</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#2dd268f2d7fbe6382cb8825a1114192a">_Ry</a></td></tr>
232
233<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of fy.R <br></td></tr>
234<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="fa172078091e45561343fa513dd573b0"></a><!-- doxytag: member="bdm::KalmanCh::_mu" ref="fa172078091e45561343fa513dd573b0" args="" -->
235vec &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#fa172078091e45561343fa513dd573b0">_mu</a></td></tr>
236
237<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of est.mu <br></td></tr>
238<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="00c27b0bf324f0018497921ca23c71ed"></a><!-- doxytag: member="bdm::KalmanCh::_P" ref="00c27b0bf324f0018497921ca23c71ed" args="" -->
239<a class="el" href="classbdm_1_1chmat.html">chmat</a> &amp;&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1Kalman.html#00c27b0bf324f0018497921ca23c71ed">_P</a></td></tr>
240
241<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">cache of est.R <br></td></tr>
242<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="c400357e37d27a4834b2b1d9211009ed"></a><!-- doxytag: member="bdm::KalmanCh::drv" ref="c400357e37d27a4834b2b1d9211009ed" args="" -->
243<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>
244
245<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Random variable of the data (optional). <br></td></tr>
246<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="4064b6559d962633e4372b12f4cd204a"></a><!-- doxytag: member="bdm::KalmanCh::ll" ref="4064b6559d962633e4372b12f4cd204a" args="" -->
247double&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a">ll</a></td></tr>
248
249<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Logarithm of marginalized data likelihood. <br></td></tr>
250<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="faff0ad12556fe7dc0e2807d4fd938ee"></a><!-- doxytag: member="bdm::KalmanCh::evalll" ref="faff0ad12556fe7dc0e2807d4fd938ee" args="" -->
251bool&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee">evalll</a></td></tr>
252
253<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>
254<tr><td colspan="2"><br><h2>Extension to conditional BM</h2></td></tr>
255<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>
256<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="b5e8e016be16e4029c789063cb3f3b33"></a><!-- doxytag: member="bdm::KalmanCh::_rvc" ref="b5e8e016be16e4029c789063cb3f3b33" args="() const " -->
257const <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>
258
259<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">access function <br></td></tr>
260<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)" -->
261virtual 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>
262
263<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Substitute <code>val</code> for <code>rvc</code>. <br></td></tr>
264<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="de95c590d650b2db423d53992f0906cb"></a><!-- doxytag: member="bdm::KalmanCh::rvc" ref="de95c590d650b2db423d53992f0906cb" args="" -->
265<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>
266
267<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Name of extension variable. <br></td></tr>
268<tr><td colspan="2"><br><h2>Logging of results</h2></td></tr>
269<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)" -->
270virtual 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>
271
272<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Set boolean options from a string, recognized are: "logbounds,logll". <br></td></tr>
273<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;)" -->
274virtual 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>
275
276<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>
277<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)" -->
278virtual 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>
279
280<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="109c1a626a69031658e3a44e9e500cca"></a><!-- doxytag: member="bdm::KalmanCh::LIDs" ref="109c1a626a69031658e3a44e9e500cca" args="" -->
281ivec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca">LIDs</a></td></tr>
282
283<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>
284<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="398b7e7a162cf37d0022d6126495f207"></a><!-- doxytag: member="bdm::KalmanCh::LFlags" ref="398b7e7a162cf37d0022d6126495f207" args="" -->
285ivec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207">LFlags</a></td></tr>
286
287<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>
288</table>
289<hr><a name="_details"></a><h2>Detailed Description</h2>
290<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.
291<p>
292Trivial 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>
293<span class="keyword">using namespace </span>bdm;
294
295<span class="comment">// estimation of AR(0) model</span>
296<span class="keywordtype">int</span> main() {
297        <span class="comment">//dimensions</span>
298        <span class="keywordtype">int</span> dx = 3, dy = 3, du = 1;
299        <span class="comment">// matrices</span>
300        mat A = eye ( dx );
301        mat B = zeros ( dx, du );
302        mat C = eye ( dx );
303        mat D = zeros ( dy, du );
304        mat Q = eye ( dx );
305        mat R = 0.1 * eye ( dy );
306        <span class="comment">//prior</span>
307        mat P0 = 100 * eye ( dx );
308        vec mu0 = zeros ( dx );
309        <span class="comment">// Estimator</span>
310        KalmanCh KF;
311        KF.set_parameters ( A, B, C, D,<span class="comment">/*covariances*/</span> Q, R );
312        KF.set_statistics ( mu0, P0 );
313        <span class="comment">// Estimation loop</span>
314        <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i = 0; i &lt; 100; i++ ) {
315                KF.bayes ( randn ( dx + du ) );
316        }
317        <span class="comment">//print results</span>
318        cout &lt;&lt; <span class="stringliteral">"Posterior estimate of x is: "</span>  &lt;&lt; endl;
319        cout &lt;&lt; <span class="stringliteral">"mean: "</span> &lt;&lt; KF.posterior().mean() &lt;&lt; endl;
320        cout &lt;&lt; <span class="stringliteral">"variance: "</span> &lt;&lt; KF.posterior().variance() &lt;&lt; endl;
321}
322</pre></div> <hr><h2>Member Function Documentation</h2>
323<a class="anchor" name="b41fe5540548100b08e1684c3be767b6"></a><!-- doxytag: member="bdm::KalmanCh::bayes" ref="b41fe5540548100b08e1684c3be767b6" args="(const vec &amp;dt)" -->
324<div class="memitem">
325<div class="memproto">
326      <table class="memname">
327        <tr>
328          <td class="memname">void bdm::KalmanCh::bayes           </td>
329          <td>(</td>
330          <td class="paramtype">const vec &amp;&nbsp;</td>
331          <td class="paramname"> <em>dt</em>          </td>
332          <td>&nbsp;)&nbsp;</td>
333          <td><code> [virtual]</code></td>
334        </tr>
335      </table>
336</div>
337<div class="memdoc">
338
339<p>
340Here dt = [yt;ut] of appropriate dimensions.
341<p>
342The following equality hold::<p class="formulaDsp">
343<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">
344<p>
345<p>
346Thus this object evaluates only predictors! Not filtering densities.
347<p>Reimplemented from <a class="el" href="classbdm_1_1Kalman.html#4a39330c14eff8d13179e868a1d1aa8c">bdm::Kalman&lt; chmat &gt;</a>.</p>
348
349<p>Reimplemented in <a class="el" href="classbdm_1_1EKFCh.html#4c8609c37290b158f88a31dae4047225">bdm::EKFCh</a>.</p>
350
351<p>References <a class="el" href="chmat_8h-source.html#l00084">bdm::chmat::_Ch()</a>, <a class="el" href="kalman_8h-source.html#l00096">bdm::Kalman&lt; chmat &gt;::_K</a>, <a class="el" href="kalman_8h-source.html#l00102">bdm::Kalman&lt; chmat &gt;::_mu</a>, <a class="el" href="kalman_8h-source.html#l00104">bdm::Kalman&lt; chmat &gt;::_P</a>, <a class="el" href="kalman_8h-source.html#l00100">bdm::Kalman&lt; chmat &gt;::_Ry</a>, <a class="el" href="kalman_8h-source.html#l00098">bdm::Kalman&lt; chmat &gt;::_yp</a>, <a class="el" href="kalman_8h-source.html#l00078">bdm::Kalman&lt; chmat &gt;::A</a>, <a class="el" href="kalman_8h-source.html#l00080">bdm::Kalman&lt; chmat &gt;::B</a>, <a class="el" href="kalman_8h-source.html#l00082">bdm::Kalman&lt; chmat &gt;::C</a>, <a class="el" href="kalman_8h-source.html#l00084">bdm::Kalman&lt; chmat &gt;::D</a>, <a class="el" href="kalman_8h-source.html#l00076">bdm::Kalman&lt; chmat &gt;::dimu</a>, <a class="el" href="kalman_8h-source.html#l00072">bdm::Kalman&lt; chmat &gt;::dimx</a>, <a class="el" href="kalman_8h-source.html#l00074">bdm::Kalman&lt; chmat &gt;::dimy</a>, <a class="el" href="bdmbase_8h-source.html#l00852">bdm::BM::evalll</a>, <a class="el" href="exp__family_8h-source.html#l00049">bdm::eEF::evallog()</a>, <a class="el" href="kalman_8h-source.html#l00093">bdm::Kalman&lt; chmat &gt;::fy</a>, <a class="el" href="bdmbase_8h-source.html#l00850">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>
352
353</div>
354</div><p>
355<a class="anchor" name="50257e0c1e5b5c73153ea6e716ad8ae0"></a><!-- doxytag: member="bdm::KalmanCh::logpred" ref="50257e0c1e5b5c73153ea6e716ad8ae0" args="(const vec &amp;dt) const " -->
356<div class="memitem">
357<div class="memproto">
358      <table class="memname">
359        <tr>
360          <td class="memname">virtual double bdm::BM::logpred           </td>
361          <td>(</td>
362          <td class="paramtype">const vec &amp;&nbsp;</td>
363          <td class="paramname"> <em>dt</em>          </td>
364          <td>&nbsp;)&nbsp;</td>
365          <td> const<code> [inline, virtual, inherited]</code></td>
366        </tr>
367      </table>
368</div>
369<div class="memdoc">
370
371<p>
372Evaluates predictive log-likelihood of the given data record I.e. marginal likelihood of the data with the posterior integrated out.
373<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>
374
375<p>Referenced by <a class="el" href="bdmbase_8h-source.html#l00886">bdm::BM::logpred_m()</a>.</p>
376
377</div>
378</div><p>
379<hr>The documentation for this class was generated from the following files:<ul>
380<li><a class="el" href="kalman_8h-source.html">kalman.h</a><li>kalman.cpp</ul>
381</div>
382<hr size="1"><address style="text-align: right;"><small>Generated on Sun Aug 16 17:58:19 2009 for mixpp by&nbsp;
383<a href="http://www.doxygen.org/index.html">
384<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.8 </small></address>
385</body>
386</html>
Note: See TracBrowser for help on using the browser.