From ControlTheoryPro.com

Contents
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
 LearnLaggedLinear
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
where
 is the roll, pitch, and yaw Euler angles,
 is the angular velocity of bodyfixed axis frame relative to earthfixed frame,
 is the linear velocity of the bodyfixed 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 pseudoinverse 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 inflight 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.
where
 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.
3.4 LearnLaggedLinear
3.5 Results
 General Least Squares: Fair results
 Coordinate Descent: Initial trials suggested that it was too numerically intensive for inflight use
 Recursive Estimation with Gaussian Noise: Performs well enough
 LearnLaggedLinear