15 lines
407 B
Mathematica
15 lines
407 B
Mathematica
|
function D = solveLinearSystem(Model,K,P)
|
|||
|
|
% Partition system of equations
|
|||
|
|
% [ Kff Kfs ] * [ Df ] = [ Pf ]
|
|||
|
|
% [ Ksf Kss ] [ Ds ] = [ Ps ]
|
|||
|
|
|
|||
|
|
Kff = K(1:Model.neqf,1:Model.neqf);
|
|||
|
|
Pf = P(1:Model.neqf); %<EFBFBD><EFBFBD><EFBFBD>ɽڵ<EFBFBD><EFBFBD><EFBFBD>Ӧ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ɶ<EFBFBD>
|
|||
|
|
Ds = zeros(Model.neqc,1);%<EFBFBD><EFBFBD>Լ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ɶ<EFBFBD>
|
|||
|
|
|
|||
|
|
% Solve for Df
|
|||
|
|
Df = Kff\Pf;
|
|||
|
|
|
|||
|
|
% Reconstruct the global unknown vector D
|
|||
|
|
D = [Df;Ds];
|
|||
|
|
end
|