Connecting a model to a Kalman filter or extended Kalman filter

2 views (last 30 days)
Hi, how do I connect a model, e.g a polynomial model identified from system identification toolbox to a Kalman filter in Simulink? What are the statetransitionfcn and the measurementfcn? Any example? Thank you!

Accepted Answer

Rajiv Singh
Rajiv Singh on 7 Nov 2022
Suppose your starting point is a linear model, SYS, identified using an offline identification routine such as TFEST, SSEST, ARMAX, etc. Follow these steps:
  1. Extract state-space data from the model using the IDSSDATA command, [A,B,C,D,K] = idssdata(SYS)
  2. In the Kalman filter block, check the "Use G and H matrices" checkbox.
  3. In the Kalman filter block, set A, B, C, D matrices. Also set Q = SYS.NoiseVariance, R = N = 0. Also set G = K and H = eye(ny), where ny = size(C,1).
If you are using the identified linear model is an Extended Kalman Filter, Unscented Kalman Filter or a Particle Filter block:
Create a MATLAB function file that returns the state update x(t+1) as a function of the current state x(t) and input u(t) as follows:
%-------------------------------
function dx = stateFcn(x,u,A,B)
% state transition fcn
dx = A*x + B*u;
end
%-------------------------------
Similarly, write a function that returns the output as a function of the current state and input values:
%-------------------------------
function y = measurementFcn(x,u,C,D)
% measurement function
y = C*x + D*u;
end
%-------------------------------
In the block, set the state transition function to stateFcn, and the measurement function to measurementFcn. Make sure these functions are on MATLAB path.
  1 Comment
Learning
Learning on 10 Nov 2022
Thank you for your response! Question: from a block layout standpoint in Simulink, how would it look like? Do I connect the input block as measurements1 in the Extended Kalman filter? Can I just pass the input through a block containing just the ss model and use the output of the ss model as measurement1 in the Extended Kalman filter? In this case what would the measurement function and state transition functions be?
Thank you!

Sign in to comment.

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by