-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathUCdisturb.m
30 lines (28 loc) · 940 Bytes
/
UCdisturb.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
function sys = UCdisturb(sys)
% UCdisturb - Runs the Disturbance Smoother for UC models
%
% sys = UCdisturb(sys)
%
% Input:
% sys: structure of type UComp created with UCmodel
%
% Output:
% The same input structure with the appropiate fields filled in, in particular:
% yFit: Fitted values of output
% yFitV: Variance of fitted values of output
% a: State estimates
% P: Variance of state estimates (diagonal of covariance matrices)
% eta: State perturbations estimates
% eps: Observed perturbations estimates
%
% Authors: Diego J. Pedregal, Nerea Urbina
%
% Examples:
% load data/airpas
% m = UCmodel(log(y), 12)
% m = UCdisturb(m)
%
% See also UC, UCcomponents, UCsmooth, UCestim, UCfilter, UCmodel,
% UCsetup, UCvalidate
sys = filter_(sys,'disturb');
end