Неявно создайте рассеянную модель в пространстве состояний, содержащую компонент регрессии

В этом примере показано, как неявно создать рассеянную модель в пространстве состояний, которая содержит компонент регрессии в уравнении наблюдения. Модель состояния содержит ARMA (1,1) и случайный обход состояния.

Запишите функцию, которая задает как параметры в params сопоставьте с матрицами модели в пространстве состояний со значениями начального состояния, и к типу состояния. Задайте компонент регрессии путем выкачивания наблюдений в функции. Символически, модель:

$$\begin{array}{l}
\left[ {\begin{array}{*{20}{c}}
{{x_{1,t}}}\\
{{x_{2,t}}}\\
x_{3,t}
\end{array}} \right] = \left[ {\begin{array}{*{20}{c}}
{{\phi _1}}&{{\theta _1}}&0\\
0&0&0\\
0&0&1
\end{array}} \right]\left[ {\begin{array}{*{20}{c}}
{{x_{1,t - 1}}}\\
{{x_{2,t - 1}}}\\
x_{3,t - 1}
\end{array}} \right] + \left[ {\begin{array}{*{20}{c}}
{{\sigma _1}}& 0\\
1 & 0\\
0 & \sigma_3
\end{array}} \right]\left[ {\begin{array}{*{20}{c}}
{u_{1,t}}\\
{u_{3,t}}
\end{array}} \right]\\
{y_t} - \beta {z_t} = a_1{x_{1,t}} + a_2x_{3,t}  + {\sigma_2}{\varepsilon _t}.
\end{array}$$


% Copyright 2015 The MathWorks, Inc.

function [A,B,C,D,Mean0,Cov0,StateType,DeflateY] = diffuseRegressionParamMap(params,y,z)
% Diffuse state-space model with a regression component parameter mapping
% function example. This function maps the vector params to the state-space
% matrices (A, B, C, and D) and indicates the type of states (StateType).
% The state model contains an ARMA(1,1) model and a random walk.
    varu1 = exp(params(3)); % Positive variance constraint
    vare1 = exp(params(5));
    A = [params(1) params(2); 0 0];
    B = [sqrt(varu1) 0; 1 0]; 
    C = [varu1 0];
    D = sqrt(vare1);
    Mean0 = []; % Let software infer Mean0
    Cov0 = [];  % Let software infer Cov0
    StateType = [0 0 2];
    DeflateY = y - params(6)*z;
end

Сохраните этот код как файл с именем diffuseRegressionParamMap.m к папке на вашем пути MATLAB®.

Создайте рассеянную модель в пространстве состояний путем передачи diffuseRegressionParamMap как указатель на функцию к dssm.

Mdl = dssm(@(params)diffuseRegressionParamMap(params,y,z));

dssm неявно создает рассеянную модель в пространстве состояний. Обычно, вы не можете проверить неявно заданные модели в пространстве состояний.

Прежде, чем создать модель, гарантируйте что переменные y и z существуйте в своей рабочей области.

Смотрите также

Похожие темы

Для просмотра документации необходимо авторизоваться на сайте