|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| In [[control theory]], a '''state observer''' is a system that provides an estimate of the [[state space (controls)|internal state]] of a given real system, from measurements of the [[Input/output|input]] and [[output]] of the real system. It is typically computer-implemented, and provides the basis of many practical applications.
| | It is very common to have a dental emergency -- a fractured tooth, an abscess, or severe pain when chewing. Over-the-counter pain medication is just masking the problem. Seeing an emergency dentist is critical to getting the source of the problem diagnosed and corrected as soon as possible.<br><br><br><br>Here are some common dental emergencies:<br>Toothache: The most common dental emergency. This generally means a badly decayed tooth. As the pain affects the tooth's nerve, treatment involves gently removing any debris lodged in the cavity being careful not to poke deep as this will cause severe pain if the nerve is touched. Next rinse vigorously with warm water. Then soak a small piece of cotton in oil of cloves and insert it in the cavity. This will give temporary relief until a dentist can be reached.<br><br>At times the pain may have a more obscure location such as decay under an old filling. As this can be only corrected by a dentist there are two things you can do to help the pain. Administer a pain pill (aspirin or some other analgesic) internally or dissolve a tablet in a half glass (4 oz) of warm water holding it in the mouth for several minutes before spitting it out. DO NOT PLACE A WHOLE TABLET OR ANY PART OF IT IN THE TOOTH OR AGAINST THE SOFT GUM TISSUE AS IT WILL RESULT IN A NASTY BURN.<br><br>Swollen Jaw: This may be caused by several conditions the most probable being an abscessed tooth. In any case the treatment should be to reduce pain and swelling. An ice pack held on the outside of the jaw, (ten minutes on and ten minutes off) will take care of both. If this does not control the pain, an analgesic tablet can be given every four hours.<br><br>Other Oral Injuries: Broken teeth, cut lips, bitten tongue or lips if severe means a trip to a dentist as soon as possible. In the mean time rinse the mouth with warm water and place cold compression the face opposite the injury. If there is a lot of bleeding, apply direct pressure to the bleeding area. If bleeding does not stop get patient to the emergency room of a hospital as stitches may be necessary.<br><br>Prolonged Bleeding Following Extraction: Place a gauze pad or better still a moistened tea bag over the socket and have the patient bite down gently on it for 30 to 45 minutes. The tannic acid in the tea seeps into the tissues and often helps stop the bleeding. If bleeding continues after two hours, call the dentist or take patient to the emergency room of the nearest hospital.<br><br>Broken Jaw: If you suspect the patient's jaw is broken, bring the upper and lower teeth together. Put a necktie, handkerchief or towel under the chin, tying it over the head to immobilize the jaw until you can get the patient to a dentist or the emergency room of a hospital.<br><br>Painful Erupting Tooth: In young children teething pain can come from a loose baby tooth or from an erupting permanent tooth. Some relief can be given by crushing a little ice and wrapping it in gauze or a clean piece of cloth and putting it directly on the tooth or gum tissue where it hurts. The numbing effect of the cold, along with an appropriate dose of aspirin, usually provides temporary relief.<br><br>In young adults, an erupting 3rd molar (Wisdom tooth), especially if it is impacted, can cause the jaw to swell and be quite painful. Often the gum around the tooth will show signs of infection. Temporary relief can be had by giving aspirin or some other painkiller and by dissolving an aspirin in half a glass of warm water and holding this solution in the mouth over the sore gum. AGAIN DO NOT PLACE A TABLET DIRECTLY OVER THE GUM OR CHEEK OR USE THE ASPIRIN SOLUTION ANY STRONGER THAN RECOMMENDED TO PREVENT BURNING THE TISSUE. The swelling of the jaw can be reduced by using an ice pack on the outside of the face at intervals of ten minutes on and ten minutes off.<br><br>If you liked this article and you would like to get much more information regarding [http://www.youtube.com/watch?v=90z1mmiwNS8 dentist DC] kindly visit our web page. |
| | |
| Knowing the system state is necessary to solve many [[control theory]] problems; for example, stabilizing a system using [[Full state feedback|state feedback]]. In most practical cases, the physical state of the system cannot be determined by direct observation. Instead, indirect effects of the internal state are observed by way of the system outputs. A simple example is that of vehicles in a tunnel: the rates and velocities at which vehicles enter and leave the tunnel can be observed directly, but the exact state inside the tunnel can only be estimated. If a system is [[Observability|observable]], it is possible to fully reconstruct the system state from its output measurements using the state observer.
| |
| | |
| == Typical observer model ==
| |
| | |
| The state of a linear, time-invariant physical discrete-time system is assumed to satisfy
| |
| | |
| : <math>x(k+1) = A x(k) + B u(k)</math>
| |
| | |
| : <math>y(k) = C x(k) + D u(k)</math>
| |
| | |
| where, at time <math>k</math>, <math>x(k)</math> is the plant's state; <math>u(k)</math> is its inputs; and <math>y(k)</math> is its outputs. These equations simply say that the plant's current outputs and its future state are both determined solely by its current states and the current inputs. (Although these equations are expressed in terms of [[discrete mathematics|discrete]] time steps, very similar equations hold for [[continuous function|continuous]] systems). If this system is [[Observability|observable]] then the output of the plant, <math>y(k)</math>, can be used to steer the state of the state observer.
| |
| | |
| The observer model of the physical system is then typically derived from the above equations. Additional terms may be included in order to ensure that, on receiving successive measured values of the plant's inputs and outputs, the model's state converges to that of the plant. In particular, the output of the observer may be subtracted from the output of the plant and then multiplied by a matrix <math>L</math>; this is then added to the equations for the state of the observer to produce a so-called ''[[David Luenberger|Luenberger]] observer'', defined by the equations below. Note that the variables of a state observer are commonly denoted by a "hat": <math>\hat{x}(k)</math> and <math>\hat{y}(k)</math> to distinguish them from the variables of the equations satisfied by the physical system.
| |
| <!-- insert plant with observer systems schematic -->
| |
| | |
| : <math>\hat{x}(k+1) = A \hat{x}(k) + L \left[y(k) - \hat{y}(k)\right] + B u(k)</math> | |
| | |
| : <math>\hat{y}(k) = C \hat{x}(k) + D u(k)</math> | |
| | |
| The observer is called asymptotically stable if the observer error <math>e(k) = \hat{x}(k) - x(k)</math> converges to zero when <math> k \rightarrow \infty </math>. For a Luenberger observer, the observer error satisfies <math> e(k+1) = (A - LC) e(k)</math>. The Luenberger observer for this discrete-time system is therefore asymptotically stable when the matrix <math> A - LC </math> has all the eigenvalues inside the unit circle. | |
| | |
| For control purposes the output of the observer system is fed back to the input of both the observer and the plant through the gains matrix <math>K</math>.
| |
| | |
| : <math>u(k)= -K \hat{x}(k)</math>
| |
| | |
| The observer equations then become:
| |
| | |
| : <math>\hat{x}(k+1) = A \hat{x}(k) + L \left(y(k) - \hat{y}(k)\right) - B K \hat{x}(k)</math>
| |
| | |
| : <math>\hat{y}(k) = C \hat{x}(k) - D K \hat{x}(k)</math>
| |
| | |
| or, more simply,
| |
| | |
| : <math>\hat{x}(k+1) = \left(A - B K \right) \hat{x}(k) + L \left(y(k) - \hat{y}(k)\right)</math>
| |
| | |
| : <math>\hat{y}(k) = \left(C - D K\right) \hat{x}(k)</math><!--
| |
| | |
| Substituting in the equation for <math>y(k)</math> from the plant system
| |
| | |
| : <math>\hat{x}(k+1) = \left(A - B K - L C) \right) \hat{x}(k) + L \hat{y}(k)</math>
| |
| | |
| : <math>\hat{y}(k) = \left(C - D K\right) \hat{x}(k)</math>
| |
| | |
| -->
| |
| | |
| Due to the [[separation principle]] we know that we can choose <math>K</math> and <math>L</math> independently without harm to the overall stability of the systems. As a rule of thumb, the poles of the observer <math>A-LC</math> are usually chosen to converge 10 times faster than the poles of the system <math>A-BK</math>.
| |
| | |
| === Continuous-time case ===
| |
| | |
| The previous example was for an observer implemented in a discrete-time LTI system. However, the process is similar for the continuous-time case; the observer gains <math>L</math> are chosen to make the continuous-time error dynamics converge to zero asymptotically (i.e., when <math>A-LC</math> is a [[Hurwitz matrix]]).
| |
| | |
| For a continuous-time linear system
| |
| | |
| : <math>\dot{x} = A x+ B u, </math>
| |
| | |
| : <math>y = C x, </math> ,
| |
| | |
| where <math>x \in \mathbb{R}^n, u \in \mathbb{R}^m ,y \in \mathbb{R}^r</math>, the observer looks similar to discrete-time case described above:
| |
| | |
| : <math>\dot{\hat{x}} = A \hat{x}+ B u + L \left(y - C \hat{x}\right) </math>.
| |
| | |
| The observer error <math>e=\hat{x}-x</math> satisfies the equation
| |
| | |
| : <math> \dot{e} = (A - LC) e</math>.
| |
| | |
| The eigenvalues of the matrix <math>A-LC</math> can be made arbitrarily by appropriate choice of the observer gain <math>L</math> when the pair <math>[A,C]</math> is observable, i.e. [[observability]] condition holds. In particular, it can be made Hurwitz, so the observer error <math>e(t) \rightarrow 0</math> when <math>t \rightarrow \infty</math>.
| |
| | |
| === Peaking and other observer methods ===
| |
| | |
| When the observer gain <math>L</math> is high, the linear Luenberger observer converges to the system states very quickly. However, high observer gain leads to a peaking phenomenon in which initial estimator error can be prohibitively large (i.e., impractical or unsafe to use).<ref name="Khalil02">{{Citation
| |
| | last = Khalil
| |
| | first = H.K.
| |
| | authorlink = Hassan K. Khalil
| |
| | year = 2002
| |
| | edition = 3rd
| |
| | url = http://www.egr.msu.edu/~khalil/NonlinearSystems/
| |
| | isbn = 0-13-067389-7
| |
| | title = Nonlinear Systems
| |
| | publisher = [[Prentice Hall]]
| |
| | location = Upper Saddle River, NJ}}</ref> As a consequence, nonlinear high gain observer methods are available that converge quickly without the peaking phenomenon. For example, [[sliding mode control]] can be used to design an observer that brings one estimated state's error to zero in finite time even in the presence of measurement error; the other states have error that behaves similarly to the error in a Luenberger observer after peaking has subsided. Sliding mode observers also have attractive noise resilience properties that are similar to a [[Kalman filter]].<ref name="UtkinGS99">{{citation|title=Sliding Mode Control in Electromechanical Systems|last1=Utkin|first1=Vadim|last2=Guldner|first2=Jürgen|last3=Shi|first3=Jingxin|year=1999|publisher=Taylor & Francis, Inc.|location=Philadelphia, PA|isbn=0-7484-0116-4}}</ref><ref name="Drakunov83">{{citation|title=An adaptive quasioptimal filter with discontinuous parameters|journal=Automation and Remote Control|last1=Drakunov|first1=S.V.|year=1983|volume=44|issue=9|pages=1167–1175}}</ref>
| |
| | |
| == State observers for nonlinear systems ==
| |
| | |
| Sliding mode observers can be designed for the non-linear systems as well. For simplicity, first consider the no-input non-linear system:
| |
| | |
| : <math>\dot{x} = f(x)</math>
| |
| | |
| where <math>x \in \mathbb{R}^n</math>. Also assume that there is a measurable output <math>y \in \mathbb{R}</math> given by
| |
| | |
| : <math>y = h(x).</math>
| |
| | |
| There are several non-approximate approaches for designing an observer. The two observers given below also apply to the case when the system has an input. That is,
| |
| | |
| : <math>\dot{x} = f(x) + B(x) u, </math>
| |
| : <math>y = h(x), </math> .
| |
| | |
| === Linearizable error dynamics ===
| |
| | |
| One suggested by Krener and Isidori<ref name="KrenerIsidori83">{{citation|doi=10.1016/0167−6911(83)90037−3|title=Linearization by output injection and nonlinear observers|journal=System and Control Letters|last1=Krener|first1=A.J.|year=1983|volume=3|pages=47–52|last2=Isidori|first2=Alberto}}</ref> and Krener and Respondek<ref name="KrenerRespondek85">{{citation|doi=10.1137/0323016|title=Nonlinear observers with linearizable error dynamics|journal=SIAM Journal on Control and Optimization|last1=Krener|first1=A.J.|last2=Respondek|first2=W.|year=1985|volume=23|pages=197–216|issue=2}}</ref> can be applied in a situation when there exists a linearizing transformation (i.e., a [[diffeomorphism]], like the one used in [[feedback linearization]]) <math>z=\Phi(x)</math> such that in new variables the system equations read
| |
| | |
| : <math>\dot{z} = A z+ \phi(y), </math>
| |
| : <math>y = Cz. </math>
| |
| | |
| The Luenberger observer is then designed as
| |
| | |
| : <math>\dot{\hat{z}} = A \hat{z}+ \phi(y) - L \left(C \hat{z}-y \right) </math>.
| |
| | |
| The observer error for the transformed variable <math>e=\hat{z}-z</math> satisfies the same equation as in classical linear case.
| |
| | |
| : <math> \dot{e} = (A - LC) e</math>.
| |
| | |
| As shown by Gauthier, Hammouri, and Othman<ref name="GauthierHammouriOthman92">{{citation|title=A simple observer for nonlinear systems applications to bioreactors |journal=IEEE Transactions on AC|last1=Gauthier|first1=J.P.|last2=Hammouri|first2=H.|last3=Othman|first3=S.|year=1992|doi=10.1109/9.256352|volume=37|issue=6|pages=875}}</ref>
| |
| and Hammouri and Kinnaert,<ref name="HammouriKinnaert96">{{citation|title=A New Procedure for Time-Varying Linearization up to Output Injection|journal=System and Control Letters|last1=Hammouri|first1=H.|last2=Kinnaert|year=1996|doi=10.1016/0167−6911(96)00022−9|first2=M.|volume=28|issue=3|pages=151}}</ref> if there exists transformation <math>z=\Phi(x)</math> such that the system can be transformed into the form | |
| | |
| : <math>\dot{z} = A(u(t)) z+ \phi(y,u(t) ), </math>
| |
| : <math>y = Cz, </math> | |
| | |
| then the observer is designed as
| |
| | |
| : <math>\dot{\hat{z}} = A(u(t)) \hat{z}+ \phi(y,u(t) ) - L(t) \left(C \hat{z}-y \right) </math>,
| |
| | |
| where <math>L(t)</math> is a time-varying observer gain.
| |
| | |
| === Sliding mode observer ===
| |
| | |
| As discussed for the linear case above, the peaking phenomenon present in Luenberger observers justifies the use of a [[Sliding mode control#Sliding mode observer|sliding mode observer]]. The sliding mode observer uses non-linear high-gain feedback to drive estimated states to a [[hypersurface]] where there is no difference between the estimated output and the measured output. The non-linear gain used in the observer is typically implemented with a scaled switching function, like the [[Sign function|signum]] (i.e., sgn) of the estimated – measured output error. Hence, due to this high-gain feedback, the vector field of the observer has a crease in it so that observer trajectories ''slide along'' a curve where the estimated output matches the measured output exactly. So, if the system is [[observable]] from its output, the observer states will all be driven to the actual system states. Additionally, by using the sign of the error to drive the sliding mode observer, the observer trajectories become insensitive to many forms of noise. Hence, some sliding mode observers have attractive properties similar to the [[Kalman filter]] but with simpler implementation.<ref name="UtkinGS99" /><ref name="Drakunov83" />
| |
| | |
| As suggested by Drakunov,<ref name="Drakunov92">{{citation|title=Sliding-Mode Observers Based on Equivalent Control Method|journal=Proceedings of the 31st IEEE Conference on Decision and Control (CDC), |last=Drakunov|first=S.V.|year=1992|issue=Tucson, Arizona, December 16–18|pages=2368–2370|url=http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=371368&isnumber=8509|isbn=0−7803−0872−7|doi=10.1109/CDC.1992.371368}}</ref> a [[Sliding mode control#Sliding mode observer|sliding mode observer]] can also be designed for a class of non-linear systems. Such an observer can be written in terms of original variable estimate <math>\hat{x}</math> and has the form
| |
| | |
| : <math> \dot{\hat{x}} = \left [ \frac{\partial H(\hat{x})}{\partial x}
| |
| | |
| \right]^{-1} M(\hat{x}) \, \operatorname{sgn}( V(t) - H(\hat{x}) )</math>
| |
| | |
| where:
| |
| | |
| * The <math>\operatorname{sgn}(\mathord{\cdot})</math> vector extends the scalar [[sign function|signum function]] to <math>n</math> dimensions. That is,
| |
| :: <math>\operatorname{sgn}(z) = \begin{bmatrix}
| |
| | |
| \operatorname{sgn}(z_1)\\
| |
| \operatorname{sgn}(z_2)\\
| |
| \vdots\\
| |
| \operatorname{sgn}(z_i)\\
| |
| \vdots\\
| |
| \operatorname{sgn}(z_n)
| |
| \end{bmatrix}</math>
| |
| | |
| : for the vector <math>z \in \mathbb{R}^n</math>.
| |
| * The vector <math>H(x)</math> has components that are the output function <math>h(x)</math> and its repeated Lie derivatives. In particular,
| |
| :: <math>H(x) \triangleq
| |
| | |
| \begin{bmatrix}
| |
| h_1(x)\\
| |
| h_2(x)\\
| |
| h_3(x)\\
| |
| \vdots\\
| |
| h_n(x)
| |
| \end{bmatrix}
| |
| \triangleq
| |
| \begin{bmatrix}
| |
| h(x)\\
| |
| L_{f}h(x)\\
| |
| L_{f}^2 h(x)\\
| |
| \vdots\\
| |
| L_{f}^{n-1}h(x)
| |
| \end{bmatrix}</math>
| |
| | |
| : where <math>L^i_f h</math> is the ''i''<sup>th</sup> [[Lie derivative]] of output function <math>h</math> along the vector field <math>f</math> (i.e., along <math>x</math> trajectories of the non-linear system). In the special case where the system has no input or has a [[feedback linearization|relative degree]] of ''n'', <math>H(x(t))</math> is a collection of the output <math>y(t)=h(x(t))</math> and its <math>n-1</math> derivatives. Because the inverse of the [[linearization|Jacobian linearization]] of <math>H(x)</math> must exist for this observer to be well defined, the transformation <math>H(x)</math> is guaranteed to be a local [[diffeomorphism]].
| |
| * The [[diagonal matrix]] <math>M(\hat{x})</math> of gains is such that
| |
| :: <math>M(\hat{x}) \triangleq
| |
| | |
| \operatorname{diag}( m_1(\hat{x}), m_2(\hat{x}), \ldots, m_n(\hat{x}) )
| |
| =
| |
| \begin{bmatrix}
| |
| m_1(\hat{x}) & & & & & \\
| |
| & m_2(\hat{x}) & & & & \\
| |
| & & \ddots & & & \\
| |
| & & & m_i(\hat{x}) & &\\
| |
| & & & & \ddots &\\
| |
| & & & & & m_n(\hat{x})
| |
| \end{bmatrix}</math>
| |
| | |
| : where, for each <math>i \in \{1,2,\dots,n\}</math>, element <math>m_i(\hat{x}) > 0</math> and suitably large to ensure reachability of the sliding mode.
| |
| * The observer vector <math>V(t)</math> is such that
| |
| :: <math>V(t) | |
| | |
| \triangleq
| |
| \begin{bmatrix}v_{1}(t)\\
| |
| v_2(t)\\
| |
| v_3(t)\\
| |
| \vdots\\
| |
| v_i(t)\\
| |
| \vdots\\
| |
| v_{n}(t)
| |
| \end{bmatrix}
| |
| \triangleq
| |
| \begin{bmatrix}
| |
| y(t)\\
| |
| \{ m_1(\hat{x}) \operatorname{sgn}( v_1(t) - h_1(\hat{x}(t)) ) \}_{\text{eq}}\\
| |
| \{ m_2(\hat{x}) \operatorname{sgn}( v_2(t) - h_2(\hat{x}(t)) ) \}_{\text{eq}}\\
| |
| \vdots\\
| |
| \{ m_{i-1}(\hat{x}) \operatorname{sgn}( v_{i-1}(t) - h_{i-1}(\hat{x}(t)) ) \}_{\text{eq}}\\
| |
| \vdots\\
| |
| \{ m_{n-1}(\hat{x}) \operatorname{sgn}( v_{n-1}(t) - h_{n-1}(\hat{x}(t)) ) \}_{\text{eq}}
| |
| \end{bmatrix}
| |
| </math>
| |
| | |
| : where <math>\operatorname{sgn}(\mathord{\cdot})</math> here is the normal [[signum]] function defined for scalars, and <math>\{ \ldots \}_{\text{eq}}</math> denotes an "equivalent value operator" of a discontinuous function in sliding mode.
| |
| | |
| The idea can be briefly explained as follows. According to the theory of sliding modes, in order to describe the system behavior, once sliding mode starts, the function <math>\operatorname{sgn}( v_{i}(t)\!-\! h_{i}(\hat{x}(t)) )</math> should be replaced by equivalent values (see ''equivalent control'' in the theory of [[sliding mode control|sliding mode]]s). In practice, it switches (chatters) with high frequency with slow component being equal to the equivalent value. Applying appropriate lowpass filter to get rid of the high frequency component on can obtain the value of the equivalent control, which contains more information about the state of the estimated system. The observer described above uses this method several times to obtain the state of the nonlinear system ideally in finite time.
| |
| | |
| The modified observation error can be written in the transformed states <math>e=H(x)-H(\hat{x})</math>. In particular,
| |
| | |
| : <math>\begin{align}
| |
| | |
| \dot{e}
| |
| &=
| |
| \frac{\operatorname{d}}{\operatorname{d}t} H(x)
| |
| -
| |
| \frac{\operatorname{d}}{\operatorname{d}t} H(\hat{x})\\
| |
| &=
| |
| \frac{\operatorname{d}}{\operatorname{d}t} H(x)
| |
| -
| |
| M(\hat{x}) \, \operatorname{sgn}( V(t) - H(\hat{x}(t)) ),
| |
| \end{align}</math>
| |
| and so
| |
| | |
| : <math>
| |
| | |
| \begin{align}
| |
| \begin{bmatrix}
| |
| \dot{e}_1\\
| |
| \dot{e}_2\\
| |
| \vdots\\
| |
| \dot{e}_i\\
| |
| \vdots\\
| |
| \dot{e}_{n-1}\\
| |
| \dot{e}_n
| |
| \end{bmatrix}
| |
| &=
| |
| \mathord{\overbrace{
| |
| \begin{bmatrix}
| |
| \dot{h}_1(x)\\
| |
| \dot{h}_2(x)\\
| |
| \vdots\\
| |
| \dot{h}_i(x)\\
| |
| \vdots\\
| |
| \dot{h}_{n-1}(x)\\
| |
| \dot{h}_n(x)
| |
| \end{bmatrix}
| |
| }^{\tfrac{\operatorname{d}}{\operatorname{d}t} H(x)}}
| |
| -
| |
| \mathord{\overbrace{
| |
| M(\hat{x}) \, \operatorname{sgn}( V(t) - H(\hat{x}(t)) )
| |
| }^{\tfrac{\operatorname{d}}{\operatorname{d}t} H(\hat{x})}}
| |
| =
| |
| \begin{bmatrix}
| |
| h_2(x)\\
| |
| h_3(x)\\
| |
| \vdots\\
| |
| h_{i+1}(x)\\
| |
| \vdots\\
| |
| h_n(x)\\
| |
| L_f^n h(x)
| |
| \end{bmatrix}
| |
| -
| |
| \begin{bmatrix}
| |
| m_1 \operatorname{sgn}( v_1(t) - h_1(\hat{x}(t)) )\\
| |
| m_2 \operatorname{sgn}( v_2(t) - h_2(\hat{x}(t)) )\\
| |
| \vdots\\
| |
| m_i \operatorname{sgn}( v_i(t) - h_i(\hat{x}(t)) )\\
| |
| \vdots\\
| |
| m_{n-1} \operatorname{sgn}( v_{n-1}(t) - h_{n-1}(\hat{x}(t)) )\\
| |
| m_n \operatorname{sgn}( v_n(t) - h_n(\hat{x}(t)) )
| |
| \end{bmatrix}\\
| |
| &=
| |
| \begin{bmatrix}
| |
| h_2(x) - m_1(\hat{x}) \operatorname{sgn}( \mathord{\overbrace{ \mathord{\overbrace{v_1(t)}^{v_1(t) = y(t) = h_1(x)}} - h_1(\hat{x}(t)) }^{e_1}} )\\
| |
| h_3(x) - m_2(\hat{x}) \operatorname{sgn}( v_2(t) - h_2(\hat{x}(t)) )\\
| |
| \vdots\\
| |
| h_{i+1}(x) - m_i(\hat{x}) \operatorname{sgn}( v_i(t) - h_i(\hat{x}(t)) )\\
| |
| \vdots\\
| |
| h_n(x) - m_{n-1}(\hat{x}) \operatorname{sgn}( v_{n-1}(t) - h_{n-1}(\hat{x}(t)) )\\
| |
| L_f^n h(x) - m_n(\hat{x}) \operatorname{sgn}( v_n(t) - h_n(\hat{x}(t)) )
| |
| \end{bmatrix}.
| |
| \end{align}
| |
| </math> | |
| So:
| |
| | |
| # As long as <math>m_1(\hat{x}) \geq |h_2(x(t))|</math>, the first row of the error dynamics, <math>\dot{e}_1 = h_2(\hat{x}) - m_1(\hat{x}) \operatorname{sgn}( e_1 )</math>, will meet sufficient conditions to enter the <math>e_1 = 0</math> sliding mode in finite time.
| |
| # Along the <math>e_1 = 0</math> surface, the corresponding <math>v_2(t) = \{m_1(\hat{x}) \operatorname{sgn}( e_1 )\}_{\text{eq}}</math> equivalent control will be equal to <math>h_2(x)</math>, and so <math>v_2(t) - h_2(\hat{x}) = h_2(x) - h_2(\hat{x}) = e_2</math>. Hence, so long as <math>m_2(\hat{x}) \geq |h_3(x(t))|</math>, the second row of the error dynamics, <math>\dot{e}_2 = h_3(\hat{x}) - m_2(\hat{x}) \operatorname{sgn}( e_2 )</math>, will enter the <math>e_2 = 0</math> sliding mode in finite time.
| |
| # Along the <math>e_i = 0</math> surface, the corresponding <math>v_{i+1}(t) = \{\ldots\}_{\text{eq}}</math> equivalent control will be equal to <math>h_{i+1}(x)</math>. Hence, so long as <math>m_{i+1}(\hat{x}) \geq |h_{i+2}(x(t))|</math>, the <math>(i+1)</math><sup>th</sup> row of the error dynamics, <math>\dot{e}_{i+1} = h_{i+2}(\hat{x}) - m_{i+1}(\hat{x}) \operatorname{sgn}( e_{i+1} )</math>, will enter the <math>e_{i+1} = 0</math> sliding mode in finite time.
| |
| | |
| So, for sufficiently large <math>m_i</math> gains, all observer estimated states reach the actual states in finite time. In fact, increasing <math>m_i</math> allows for convergence in any desired finite time so long as each <math>|h_i(x(0))|</math> function can be bounded with certainty. Hence, the requirement that the map <math>H:\mathbb{R}^n \rightarrow \mathbb{R}^n </math> is a [[diffeomorphism]] (i.e., that its [[linearization|Jacobian linearization]] is invertible) asserts that convergence of the estimated output implies convergence of the estimated state. That is, the requirement is an observability condition.
| |
| | |
| In the case of the sliding mode observer for the system with the input, additional conditions are needed for the observation error to be independent of the input. For example, that
| |
| | |
| : <math> \frac{\partial H(x)}{\partial x} B(x)</math>
| |
| | |
| does not depend on time. The observer is then
| |
| | |
| : <math>
| |
| | |
| \dot{\hat{x}} = \left[ \frac{\partial H(\hat{x})}{\partial x}
| |
| \right]^{-1} M(\hat{x}) \operatorname{sgn}(V(t) - H(\hat{x}))+B(\hat{x})u.
| |
| </math>
| |
| | |
| == Bounding observers ==
| |
| | |
| The Bounding <ref>http://www.nt.ntnu.no/users/skoge/prost/proceedings/ecc03/pdfs/437.pdf</ref> or Interval observers <ref>http://www.nt.ntnu.no/users/skoge/prost/proceedings/cdc-2008/data/papers/1446.pdf</ref> constitute a class of observers that provide two estimations
| |
| of the state simultaneously: one of the estimation provides an upper bound on the real value of the state,
| |
| whereas the second one provides a lower bound. The real value of the state is then known to be always within these two
| |
| estimations.
| |
| | |
| These bounds are very important in practical applications,<ref>http://www.iaeng.org/publication/WCE2010/WCE2010_pp656-661.pdf</ref><ref>http://www.sciencedirect.com/science/article/pii/S0959152499000748</ref> as they make possible to know at each time the precision of the estimation.
| |
| | |
| Mathematically, two Luenberger observers can be used, if <math> L </math> is properly selected, using, for example, [[positive systems]] properties:<ref>Ait Rami, M., Tadeo, F., Helmke, U. (2011),
| |
| "Positive observers for linear positive systems, and their implications",
| |
| International Journal of Control 84
| |
| </ref> one for the upper bound <math> \hat{x}_U(k) </math> (that ensures that <math> e(k) = \hat{x}_U(k) - x(k) </math> converges to zero from above when <math> k \rightarrow \infty </math>, in the absence of noise and [[uncertainty]]), and a lower bound <math> \hat{x}_L(k) </math> (that ensures that <math> e(k) = \hat{x}_L(k) - x(k) </math> converges to zero from below). That is, always <math> \hat{x}_U(k) \ge x(k) \ge \hat{x}_L(k) </math>
| |
| | |
| == See also ==
| |
| | |
| * [[Kalman filter]]
| |
| * [[Extended Kalman filter]]
| |
| * [[Positive systems]]
| |
| | |
| == References ==
| |
| | |
| ; In-line references
| |
| | |
| {{reflist}}
| |
| | |
| ; General references
| |
| | |
| {{refbegin}}
| |
| | |
| * {{Citation
| |
| | |
| | last = Sontag
| |
| | first = Eduardo
| |
| | authorlink = Eduardo D. Sontag
| |
| | year = 1998
| |
| | title = Mathematical Control Theory: Deterministic Finite Dimensional Systems. Second Edition
| |
| | publisher = Springer
| |
| | isbn = 0-387-98489-5
| |
| }}
| |
| | |
| {{refend}}
| |
| | |
| [[Category:Control theory]]
| |
| [[Category:Signal estimation]]
| |