Difference between revisions of "ControlTheory"

From PaparazziUAV
Jump to navigation Jump to search
 
(6 intermediate revisions by one other user not shown)
Line 1: Line 1:
== Dynamic Systems ==
== Dynamic Systems ==
* Notion of State
  The state of a dynamic system represents the minimal quantity of information required in order to be able to predict its future,
knowing its equation of dynamic and its inputs.
You can think of the state of a system as the frontier between past and future : two systems of identical dynamic and identical state will have the same trajectory in the future when applied the same inputs, whatever their past can be.
* State Space Representation


== State Estimation ==
== State Estimation ==
Line 12: Line 20:
* Direct Differentiation
* Direct Differentiation


A first idea would be to take the difference between two consecutive altitude measurements and divide that quantity by the time between the two measurements. This is indeed a bad idea. Let's see one of the reasons why.
A first idea would be to take the difference between two consecutive altitude measurements and divide that quantity by the duration between the two measurements. This is indeed a bad idea. Let's see one of the reasons why.


Our barometer being a real world instrument, it will give us noisy measurements ( think electrical noise, etc...). Let's use a little [[media:se_ex1.sce|simulation]] to see what happens when you differentiate noisy measurements:
Our barometer being a real world instrument, it will output noisy measurements ( think electrical noise, etc...). Let's use a little [[media:se_ex1.sce|simulation]] to see what happens when you differentiate noisy measurements:


[[Image:se_ex1.png|1024px]]
[[Image:se_ex1.png|1024px]]


The noise gets amplified by the division by the sampling interval ( and we did not even introduce inaccuracy in that sampling interval ).
The noise gets amplified by the division by the sampling interval ( and we did not even introduce inaccuracy in that sampling interval ).
In the above simulation, starting with a very good sensor having 1cm noise standard deviation and sampling at 100Hz, we obtain an estimate for the vertical speed having 1.4m of standard deviation. Pretty unusable...
In the above simulation, starting with a very good sensor having 1cm noise standard deviation and sampling at 100Hz, we obtain an estimate for the vertical speed having 1.4m/s of standard deviation. Pretty unusable...


* Let's filter that noise out
* Let's filter that noise out


A second idea would be to attempt and filter that noise. Let's run for example our estimate through a first order low pass filter. Here is the result of the [[media:se_ex2.sce|simulation]
A second idea would be to attempt and filter that noise. Let's run for example our estimate through a first order low pass filter. Here is the result of the [[media:se_ex2.sce|simulation]]


[[Image:se_ex2.png|1024px]]
[[Image:se_ex2.png|1024px]]


Yeah !  We reduced noise on our vertical speed estimate to 0.15m/s
Yay !  We reduced noise on our vertical speed estimate to 6cm/s. Did we achieve something interesting when doing so ? sadly no. As shown on the next plot where we have run the same algorithm on a moving vehicle, we have introduced lag in our estimate.
 
[[Image:se_ex3.png|1024px]]
 
And lag is bad for control as we'll show later. What would have happened if we had filtered the measurement before differentiating ? same punishment, as shown in the next plot produced by the following [[media:se_ex4.sce|simulation]]...
 
[[Image:se_ex4.png|1024px]]
 
so, is there a better way ?
 
* Let's introduce Kinematics
 
<math>X=\begin{pmatrix}z\\\dot{z}\end{pmatrix}</math>
 
<math>\dot{X}(t)=\begin{pmatrix}0&1\\0&0\end{pmatrix}X(t) + \begin{pmatrix}0\\1\end{pmatrix}a(t)</math>
 
We don't really have an idea of the acceleration of our system, so let's asume it's 0 for now.
When using this hypothesis and switching to discrete time (more on this later), we obtain the following equation:
 
<math>X_{t+dt}=\begin{pmatrix}1&dt\\0&1\end{pmatrix}X_{t}</math>
 
 
 
 






--[[User:Poine|Poine]] 14:43, 13 November 2009 (CET)
--[[User:Poine|Poine]] 14:43, 13 November 2009 (CET)
[[Category:Developer_Documentation]]

Latest revision as of 18:27, 9 March 2011

Dynamic Systems

  • Notion of State
 The state of a dynamic system represents the minimal quantity of information required in order to be able to predict its future,
knowing its equation of dynamic and its inputs.

You can think of the state of a system as the frontier between past and future : two systems of identical dynamic and identical state will have the same trajectory in the future when applied the same inputs, whatever their past can be.

  • State Space Representation

State Estimation

In order to be able to control a vehicle, ie to steer it along a trajectory, a first prerequisite is to know its state. Typically, sensors are used to obtain this information. Sadly, sensors won't give you a straightforward knowledge of the state of the vehicle. The information is likely to be partial ( we have an altitude sensor, but nothing to measure vertical speed) and will also likely contain "errors" (bias, noise, latency, etc... ). You can always get more sensors and/or better sensors, but there's no such thing as a perfect sensor. The process of reconstructing the state of a dynamic system using a number of incomplete and erroneous measurements is called state estimation or state observation.

Let's start with an example :

We are interested in controlling the altitude of a vehicle. We measure this altitude using for example a barometer. As we will see later, the vertical speed of the vehicle is a mandatory information in order to be able to achieve a satisfactory control. So, how do we obtain the vertical speed information ?

  • Direct Differentiation

A first idea would be to take the difference between two consecutive altitude measurements and divide that quantity by the duration between the two measurements. This is indeed a bad idea. Let's see one of the reasons why.

Our barometer being a real world instrument, it will output noisy measurements ( think electrical noise, etc...). Let's use a little simulation to see what happens when you differentiate noisy measurements:

Se ex1.png

The noise gets amplified by the division by the sampling interval ( and we did not even introduce inaccuracy in that sampling interval ). In the above simulation, starting with a very good sensor having 1cm noise standard deviation and sampling at 100Hz, we obtain an estimate for the vertical speed having 1.4m/s of standard deviation. Pretty unusable...

  • Let's filter that noise out

A second idea would be to attempt and filter that noise. Let's run for example our estimate through a first order low pass filter. Here is the result of the simulation

Se ex2.png

Yay ! We reduced noise on our vertical speed estimate to 6cm/s. Did we achieve something interesting when doing so ? sadly no. As shown on the next plot where we have run the same algorithm on a moving vehicle, we have introduced lag in our estimate.

Se ex3.png

And lag is bad for control as we'll show later. What would have happened if we had filtered the measurement before differentiating ? same punishment, as shown in the next plot produced by the following simulation...

Se ex4.png

so, is there a better way ?

  • Let's introduce Kinematics

We don't really have an idea of the acceleration of our system, so let's asume it's 0 for now. When using this hypothesis and switching to discrete time (more on this later), we obtain the following equation:




--Poine 14:43, 13 November 2009 (CET)