0

私は、例えば次のように記述される運動方程式を解いています。

syms u1 u2 f3 f4;
K=magic(4);
f=[2 3 f3 f4];
u=[u1 u2 9 7];

ここで、f と K の両方に変数が含まれます。

これを行うための最も簡単な、または最も効率的な方法はどれですか? 私は何かが欲しい

solution=solve(f==K*u);

それが効率的な方法かどうかはわかりません。

4

1 に答える 1

2

あなたの方程式はまだ標準形ではありませんが、それに還元できます。

コンピューティング

>> syms u1 u2 f3 f4;
>> K = magic(4);
>> f = [2 3 f3 f4].';
>> u = [u1 u2 9 7].';

>> K*u - f

ans =
      16*u1  +  2*u2            + 116
       5*u1  + 11*u2            + 143
       9*u1  +  7*u2 - f3       + 138
       4*u1  + 14*u2      - f4  + 142

わかりやすくするために書式を変更しました。これから、定義すると

b = [-116 -143 -138 -142].';
A  = [16  2   0   0 
       5 11   0   0 
       9  7  -1   0
       4 14   0  -1];

次に、標準形式の方程式があります

A*X = b

どこ

X = [u1 u2 f3 f4].'

となると、簡単なことですが、

>> X = A\b
ans =
   -5.963855421686747e+000    % u1
   -1.028915662650603e+001    % u2
   +1.230120481927711e+001    % f3
   -2.590361445783132e+001    % f4
于 2013-07-01T15:37:08.557 に答える