You could do it brute-force by using a pseudo-inverse:
%From state space
C_tilde=C*pinv(A);
state_space = ssm(A,B,C_tilde,D-C_tilde*B) ;
compare_irfs.m (483 Bytes)
ABCD_test.m (4.6 KB)
parameterfile.mat (3.2 KB)
rentalhousing_steadystate.m (5.7 KB)
rentalhousing.mod (17.7 KB)