1

Im trying to get MATLAB to create a state space model with initial values, which i later have to fine tune using the system identification toolbox.

The problem is that i am doing this with "white-box" models, which means that the A-matrix in my model isn't just a matrix of numeric values, but consists of several other parameters.

Since each of the parameters that makes up the matrix has a physical meaning, i'd like to estimate that seperately instead of just the value of the A-matrix.

Is this possible in any way? Have looked at both genss, idss and several other tools, but haven't been able to make any of them work as of yet.

A version of the system i'll be working with is as follows (C and D matrices are not relevant, since this is just testing really):

A = [-H_infil/Ca+Hsa_i^2/(Ca*(Hms_i+Hsa_i))-Hsa_i/Ca-H_win/Ca+Hsa_e^2/(Ca*(Hsa_e+Hms_e))-Hsa_e/Ca-H_vent/Ca Hsa_i*Hms_i/((Hms_i+Hsa_i)*Ca) Hsa_e*Hms_e/((Hsa_e+Hms_e)*Ca);
     Hms_i*Hsa_i/((Hms_i+Hsa_i)*Ci) Hms_i^2/((Hms_i+Hsa_i)*Ci)-Hms_i/Ci 0;
     Hms_e*Hsa_e/((Hsa_e+Hms_e)*Ce) 0 Hms_e^2/((Hsa_e+Hms_e)*Ce)-Hms_e/Ce-Hem/Ce];

Bc = [Hsa_i*S*wi/(Ca*(Hms_i+Hsa_i))+S*wa/Ca Hsa_i*fi/(Ca*(Hms_i+Hsa_i))+fa/Ca+Hsa_e*fe/(Ca*(Hsa_e+Hms_e)) H_infil/Ca+H_win/Ca H_vent/Ca;
    Hms_i*S*wi/((Hms_i+Hsa_i)*Ci) Hms_i*fi/((Hms_i+Hsa_i)*Ci) 0 0;
    0 Hms_e*fe/((Hsa_e+Hms_e)*Ce) Hem/Ce 0];

Cc = [0 0 0];

Dc = [0 0 0 0]

I realize that i may run into a problem with several solutions, but then i figured that the system identification would just choose one.

Any ideas?

Thanks in advance!

oh, and initial values goes as follows:

%% Q_Heating ratio factors are defined
    % Factors
   ki = 0.6;
   ke = 0;
   ka = 1 - ki - ke;

%% Internal heat gain factors are defined
   fe = 0;             % Internal heat gain ratio to external constructions
   fi = 0;             % Internal heat gain ratio to internal constructions
   fa = 1 - fe - fi;   % Internal heat gain ratio to room air

%% Solar heat gain factors are defined
   S = 0.558; % Shading factor
   wi = 0.6; % Solar heat gain ratio to internal constructions
   wa = 1 - wi; % Solar heat gain ratio to room air

%% Capacities
   % Air node
   Ca = 60708.707;
   % Internal node
   Ci = 6.160*10^6*1.3;
   % External node
   Ce = 5.489*10^5*1.3;

%% Heat transfer coefficients
    % Air to surface heat transfer coefficients
    Hsa_i = 581.286;
    Hsa_e = 52.094;
    % Surface to mass heat transfer coefficients
    Hms_i = 3235.007;
    Hms_e = 284.432;
    % Mass to external air/ground heat transfer coefficient
    Hem = 0.986*1.3;
    % Room air to external air heat transfer coefficients
    H_infil = 0;
    H_win = 1.627;
    % Ventilation heat transfer coefficient
    H_vent = 0;

0 Answers0