Hi Valeriy,
As per my understanding, you calculated the process noise covariance matrix ('Q') for a discrete linear system identified using the MATLAB System Identification Toolbox and want to make sure the results are accurate.
The formula you used attempts to account for two primary sources of uncertainty in the model:
- Inherent process noise represented by 'NoiseVariance', which is directly scaled by the identity matrix ('I') to match the dimensions of the state space. This noise component is assumed to be white noise affecting all states equally and independently.
- Input-driven uncertainty represented by the term (B*(inputVariance)*B'). considers how variations in the input (with a known variance 'inputVariance') propagate through the system dynamics (modeled by (B)) to affect the states.
Your calculation of 'Q' results in a matrix that is almost diagonal with (0.0026) on the diagonal. The off-diagonal elements are very close to zero, which is expected given the magnitude of elements in 'B'. This suggests that the process noise and input-driven uncertainties are assumed to affect each state independently, with minimal cross-correlation. This is a common simplification that makes the Kalman filter easier to implement and understand. The small magnitude of 'B' means the input-driven uncertainty's contribution to 'Q' is negligible compared to the inherent process noise 'NoiseVariance'. This is why the resulting 'Q' matrix values are very close to NoiseVariance and the matrix appears almost as an identity matrix scaled by 'NoiseVariance '.
Your approach and the resulting calculation are conceptually correct within the assumptions made (e.g., white noise, independence of state noises, linear propagation of input variance through 'B'). The almost identity appearance of 'Q' is a direct consequence of the small magnitude of 'B' in your system, which minimizes the impact of input variance on the process noise covariance.I would suggest you validate the assumptions leading to the formulation of 'Q' against the specific characteristics and requirements of your application. Adjustments may be necessary if, for example, there is significant cross-correlation between the noises affecting different states or if the input variance significantly affects the system dynamics in a non-linear manner.
Hope that helps!
