Hi,
In order to enforce that the two disturbance variances are identical in a time series state space model with a diffuse state, you can use the following approach:
- Set the state disturbance loading coefficient matrix B to [nan 0; 0 nan], as you mentioned.
- Create a new disturbance covariance matrix Q that has the desired variance for the cyclic term.
- Assign this new disturbance covariance matrix Q to both the cyclic term and the dummy cyclic term.
Here's an example of how you can implement this in MATLAB:
% Set up the state space model
A = ...; % State transition matrix
B = [nan 0; 0 nan]; % State disturbance loading coefficient matrix
C = ...; % Output matrix
D = ...; % Feedthrough matrix
% Create the disturbance covariance matrix
variance = ...; % Desired variance for the cyclic term
Q = diag([variance, variance]);
% Assign the disturbance covariance matrix to both terms
sys = ss(A, B, C, D);
sys.DisturbanceCovariance = Q;
By assigning the same disturbance covariance matrix Q to both the cyclic term and the dummy cyclic term, you ensure that their variances are identical.
This will probably resolve your issue, but if it doesn't, let me know and I'll be happy to look into it further.
