118 | | <a name="l00067"></a>00067 <span class="comment">//</span> |
119 | | <a name="l00068"></a>00068 <span class="keywordflow">if</span> ( S.exists ( <span class="stringliteral">"profile"</span> ) ) { |
120 | | <a name="l00069"></a>00069 tmp-><a class="code" href="classpmsmDS.html#60d213a13fd9c3b4ac9e8202ca5d772e">set_profile</a> ( S[<span class="stringliteral">"tstep"</span>],getvec ( S[<span class="stringliteral">"profile"</span>] ) ); |
121 | | <a name="l00070"></a>00070 } |
122 | | <a name="l00071"></a>00071 <span class="keywordflow">else</span> { |
123 | | <a name="l00072"></a>00072 tmp-><a class="code" href="classpmsmDS.html#60d213a13fd9c3b4ac9e8202ca5d772e">set_profile</a> ( 1.0, vec ( <span class="stringliteral">"1"</span> ) ); |
124 | | <a name="l00073"></a>00073 } |
125 | | <a name="l00074"></a>00074 <span class="keywordflow">return</span> tmp; |
126 | | <a name="l00075"></a>00075 }; |
127 | | <a name="l00076"></a>00076 |
128 | | <a name="l00077"></a>00077 }; |
129 | | <a name="l00078"></a>00078 UIREGISTER ( <a class="code" href="classUIpmsmDS.html" title="UI for pmsmDS,.">UIpmsmDS</a> ); |
| 118 | <a name="l00067"></a>00067 |
| 119 | <a name="l00068"></a>00068 <span class="comment">// Default values of profiles for omega and Mz</span> |
| 120 | <a name="l00069"></a>00069 vec profW=vec(<span class="stringliteral">"1.0"</span>); |
| 121 | <a name="l00070"></a>00070 vec profM=vec(<span class="stringliteral">"0.0"</span>); |
| 122 | <a name="l00071"></a>00071 <span class="keywordtype">double</span> tstep=1.0; |
| 123 | <a name="l00072"></a>00072 |
| 124 | <a name="l00073"></a>00073 <span class="keywordflow">if</span> ( S.exists ( <span class="stringliteral">"tstep"</span> ) ) {tstep=S[<span class="stringliteral">"tstep"</span>];} |
| 125 | <a name="l00074"></a>00074 <span class="keywordflow">if</span> ( S.exists ( <span class="stringliteral">"profileW"</span> ) ) {profW=getvec ( S[<span class="stringliteral">"profileW"</span>] ) ;} |
| 126 | <a name="l00075"></a>00075 <span class="keywordflow">if</span> ( S.exists ( <span class="stringliteral">"profileM"</span> ) ) {profM=getvec ( S[<span class="stringliteral">"profileM"</span>] ) ;} |
| 127 | <a name="l00076"></a>00076 |
| 128 | <a name="l00077"></a>00077 tmp-><a class="code" href="classpmsmDS.html#bcea4f040063f65ae39a4ca21f6aff1f">set_profile</a> (tstep , profW, profM); |
| 129 | <a name="l00078"></a>00078 |
| 130 | <a name="l00079"></a>00079 <span class="keywordtype">string</span> opts=<span class="stringliteral">""</span>; |
| 131 | <a name="l00080"></a>00080 <span class="keywordflow">if</span> ( S.exists ( <span class="stringliteral">"options"</span> ) ) {opts=(<span class="keyword">const</span> <span class="keywordtype">char</span>*)S[<span class="stringliteral">"options"</span>];} |
| 132 | <a name="l00081"></a>00081 tmp-><a class="code" href="classpmsmDS.html#75dd22db70746bfd0c3e685a17bbf34a" title="parse options: &quot;modelu&quot; =&gt; opt_modu=true;">set_options</a>(opts); |
| 133 | <a name="l00082"></a>00082 |
| 134 | <a name="l00083"></a>00083 <span class="keywordflow">return</span> tmp; |
| 135 | <a name="l00084"></a>00084 }; |
| 136 | <a name="l00085"></a>00085 |
| 137 | <a name="l00086"></a>00086 }; |
| 138 | <a name="l00087"></a>00087 UIREGISTER ( <a class="code" href="classUIpmsmDS.html" title="UI for pmsmDS,.">UIpmsmDS</a> ); |
| 139 | <a name="l00088"></a>00088 |
| 140 | <a name="l00089"></a>00089 <span class="keyword">class </span>UIpmsmCRB: <span class="keyword">public</span> <a class="code" href="classbdm_1_1UIbuilder.html" title="Builds computational object from a UserInfo structure.">UIbuilder</a> { |
| 141 | <a name="l00090"></a>00090 <span class="keyword">public</span>: |
| 142 | <a name="l00091"></a>00091 UIpmsmCRB():<a class="code" href="classbdm_1_1UIbuilder.html" title="Builds computational object from a UserInfo structure.">UIbuilder</a>(<span class="stringliteral">"pmsmCRB"</span>){}; |
| 143 | <a name="l00092"></a>00092 <a class="code" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a>* <a class="code" href="classbdm_1_1UIbuilder.html#193cede6302669ff48dfd8ab3b7e31f8" title="Function building the computational object.">build</a> ( Setting &S )<span class="keyword"> const </span>{ |
| 144 | <a name="l00093"></a>00093 <a class="code" href="classbdm_1_1diffbifn.html" title="Class representing a differentiable function of two variables .">diffbifn</a>* IM; UIbuild(S[<span class="stringliteral">"IM"</span>],IM); |
| 145 | <a name="l00094"></a>00094 <a class="code" href="classbdm_1_1diffbifn.html" title="Class representing a differentiable function of two variables .">diffbifn</a>* OM; UIbuild(S[<span class="stringliteral">"OM"</span>],OM); |
| 146 | <a name="l00095"></a>00095 |
| 147 | <a name="l00096"></a>00096 <span class="comment">//parameters</span> |
| 148 | <a name="l00097"></a>00097 <a class="code" href="classpmsmCRB.html" title="This class behaves like BM but it is evaluating EKF.">pmsmCRB</a>* E; E=<span class="keyword">new</span> <a class="code" href="classpmsmCRB.html" title="This class behaves like BM but it is evaluating EKF.">pmsmCRB</a>; |
| 149 | <a name="l00098"></a>00098 E-><a class="code" href="classbdm_1_1EKFCh.html#50f9fbffad721f35e5ccb75d0f6b842a" title="Set nonlinear functions for mean values and covariance matrices.">set_parameters</a>(IM, OM, diag(getvec(S[<span class="stringliteral">"dQ"</span>])), diag(getvec(S[<span class="stringliteral">"dR"</span>]))); |
| 150 | <a name="l00099"></a>00099 |
| 151 | <a name="l00100"></a>00100 <span class="comment">//statistics</span> |
| 152 | <a name="l00101"></a>00101 <span class="keywordtype">int</span> dim=IM-><a class="code" href="classbdm_1_1fnc.html#083832294da9d1e40804158b979c4341" title="access function">dimension</a>(); |
| 153 | <a name="l00102"></a>00102 vec mu0; |
| 154 | <a name="l00103"></a>00103 mat P0; |
| 155 | <a name="l00104"></a>00104 <span class="keywordflow">if</span> (S.exists(<span class="stringliteral">"mu0"</span>)){mu0=getvec(S[<span class="stringliteral">"mu0"</span>]);}<span class="keywordflow">else</span>{mu0=zeros(dim);}; |
| 156 | <a name="l00105"></a>00105 <span class="keywordflow">if</span> (S.exists(<span class="stringliteral">"P0"</span>)){mu0=getmat(S[<span class="stringliteral">"P0"</span>],dim);}<span class="keywordflow">else</span>{P0=eye(dim);}; |
| 157 | <a name="l00106"></a>00106 E-><a class="code" href="classbdm_1_1KalmanCh.html#6e169272657ed101f3d128b49c59b890">set_statistics</a>(mu0,P0); |
| 158 | <a name="l00107"></a>00107 |
| 159 | <a name="l00108"></a>00108 <span class="comment">//connect</span> |
| 160 | <a name="l00109"></a>00109 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>* drv; UIbuild(S[<span class="stringliteral">"drv"</span>],drv); |
| 161 | <a name="l00110"></a>00110 E-><a class="code" href="classbdm_1_1BM.html#f135ae6dce7e9f30c9f88229c7930b96">set_drv</a>(*drv); |
| 162 | <a name="l00111"></a>00111 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>* rv; UIbuild(S[<span class="stringliteral">"rv"</span>],rv); |
| 163 | <a name="l00112"></a>00112 E-><a class="code" href="classbdm_1_1BM.html#b38d92f17620813ad872d86e01a26e5e">set_rv</a>(*rv); |
| 164 | <a name="l00113"></a>00113 <span class="keywordflow">return</span> E; |
| 165 | <a name="l00114"></a>00114 } |
| 166 | <a name="l00115"></a>00115 }; |
| 167 | <a name="l00116"></a>00116 UIREGISTER ( UIpmsmCRB ); |