- 1 Introduction to System Identification of Cessna 182 Model UAV
- 2 Vehicle Dynamical Model - the basis for System Identification
- 3 Machine Learning for System Identification
- 4 References
1 Introduction to System Identification of Cessna 182 Model UAV
Typically system identification is employed in order to determine model parameters for a system whose parameters are unknown or whose system model is complex. Complex models can often be approximated accurately using test data.
This particular article is on the use of system identification with test data to create a system model for a quarter scale model of a Cessna 182 UAV. The article will briefly discuss system identification via the following machine learning algorithms
- General Least Squares
- Coordinate Descent
- Recursive Estimation with Gaussian Noise
2 Vehicle Dynamical Model - the basis for System Identification
The determination of system identification parameters requires a model form. In other words, the control designer or modeler needs to have a basic idea of the form the model will take. In the case of a SISO system this a matter of deciding how many poles and zeros the system model will have.
The student authors of the project this is based on an assumption of straight and steady flight. This allows the lateral and longitudinal dynamics to be decoupled with each having 4 states. After linearizing based on the small angle approximation ad zero cross product term , the lateral equations of motion become
and the longitudinal dynamics are
- is the roll, pitch, and yaw Euler angles,
- is the angular velocity of body-fixed axis frame relative to earth-fixed frame,
- is the linear velocity of the body-fixed axis frame, and
- is the aileron, rudder, elevator and throttle inputs.
Remember that these equations can be written in the form of
3 Machine Learning for System Identification
3.1 General Least Squares
For system identification, least squares regression can be used to find the system parameters. This is accomplished by minimizing the following objective function
The data can be rearranged into the conventional form
wher the state data and unknown parameters are reordered as . The pseudo-inverse is used to find .
3.2 Coordinate Descent
The performance of the dynamic model is based on its simulated trajectory relative to the actual trajectory. Minimizing the mean squared error between simulated states and actual states at each time step. However, this appears to be farr too computationally expensive.
3.3 Recursive Estimation with Gaussian Noise
Recursion methods are ideal for in-flight applications since they are formulated to use only present data thus minimizing calculation time at each step. For this formulation all past data is contained in the current state and covariance estimates.
- is the optimal state and
- is the covariance estimate given at k measurements to .
Application of this method to system identification of the model requires an estimate of the prior covariances and parameters. The initial estimates on the parameters were chosen to be 0. The covariances are large and decoupled giving more weight to current data than prior data. This leads to the observation that these initial values are most likely not optimal.
- General Least Squares: Fair results
- Coordinate Descent: Initial trials suggested that it was too numerically intensive for in-flight use
- Recursive Estimation with Gaussian Noise: Performs well enough