私は、例えば次のように記述される運動方程式を解いています。
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);
それが効率的な方法かどうかはわかりません。
私は、例えば次のように記述される運動方程式を解いています。
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);
それが効率的な方法かどうかはわかりません。
あなたの方程式はまだ標準形ではありませんが、それに還元できます。
コンピューティング
>> 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