%% Extract Zero-Pole-Gain Models from Identified Model
% For this example, extract the measured and noise components of an identified polynomial model into two separate zero-pole-gain models.
% Load the Box-Jenkins polynomial model ltiSys in identifiedModel.mat.

load('identifiedModel.mat','ltiSys');

% ltiSys is an identified discrete-time model of the form: y(t)=\frac{B}{F}u(t)+\frac{C}{D}e(t), where \frac{B}{F} represents the measured component and \frac{C}{D} the noise component.
%
% Extract the measured and noise components as zero-pole-gain models.

sysMeas = zpk(ltiSys,'measured') 
sysNoise = zpk(ltiSys,'noise')

% The measured component can serve as a plant model, while the noise component can be used as a disturbance model for control system design.
% Copyright 2020 The MathWorks, Inc.