1

線形システム AX = B を解きたいと思います。ここで、A は定数要素を持つ nxn 行列で、行列 B は nx1 型です。ただし、行列 B の各要素は nx1 型のベクトルです (これは、各要素 bij が時変であるために発生します)。

%%% Example


t = 0:0.002:0.5;    %% Time

A = [1 0 -1 0 0 0; ...
      0 -1 0 0 1 0; ...
     r12y, r32y-r12x r32x 0 0; ...
      0 0 -1 0 1 0; ...
      0 -1 0 0 1 0; ...
      0 r23y-7-r43y r23x r43x];


  %% Where rij is constant
 
% Construction 6x1 matrix C

C = [m2.*A2x ; ...
    m2.*FG2-a2y; ...
    ICM2.*Alpha2; ...
    m3.*A3X ; ...
    m3.*a3y-FG3; ...
    Icm3.*Alph8a3];


%% Where A2x, a2y, A3X, a3y, alpha2, Alpha3 are elements of the matrix C that are time-varying.

セグンテ形式を解決しようとしました：

C = rand (6,1,251);
A = rand (6,6);%


X = zeros (6, size (C, 3));
for i = 1: size (C, 3)
     X (:, i) = A \ C (:,:, i);
end

しかし、それが最善の方法かどうかはわかりません。

4

1 に答える 1

3

次の操作を行うだけです。

C = reshape(rand(6,1,251), 6, 251); % Or just create rand(6, 251);
A = rand(6,6);
X = A \ C;

これにより、同じ結果が得られ、for ループよりも高速になります。

于 2012-12-31T17:54:35.257 に答える