%% compute the subsystem matrices out of the block diagonal system

% Copyright 2023 Collaborative Research Centre Transregio 96 (CRC/TR 96)
%
% Redistribution and use in source and binary forms, with or without
% modification, are permitted provided that the following conditions are
% met:
%
% 1. Redistributions of source code must retain the above copyright notice,
% this list of conditions and the following disclaimer.
%
% 2. Redistributions in binary form must reproduce the above copyright
% notice, this list of conditions and the following disclaimer in the
% documentation and/or other materials provided with the distribution.
%
% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 
% "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
% TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
% PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER
% OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
% EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
% PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
% PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
% LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
% NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
% SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

%% load block diagonal LTI system (E,A,B,C) and struct with size information
% for the subsystems (pos)
load('machineToolMAX_outputCoupled_n1265497m287q224.mat')

% store the 50 subsystems in a cell array SA
SA=cell(50,1);
for j=1:50
    SA{j}.A = A(pos.n(j)+1:pos.n(j+1),pos.n(j)+1:pos.n(j+1));
    SA{j}.B = B(pos.n(j)+1:pos.n(j+1),pos.m(j)+1:pos.m(j+1));
    SA{j}.C = C(pos.p(j)+1:pos.p(j+1),pos.n(j)+1:pos.n(j+1));
    SA{j}.E = E(pos.n(j)+1:pos.n(j+1),pos.n(j)+1:pos.n(j+1));
end