Octaveでシンボリック行列を宣言するにはどうすればよいですか？
MatLabでは、シンボルを非常に簡単に宣言できます。
syms a,b
mat = [a,b]
ただし、Octaveでこれを複製しようとすると、エラーが発生します。これは私が使用しているコードです：
> symbols
> a = sym("a")
a =
a
> b = sym("b")
b =
b
> mat = [a,b]
error: octave_base_value::resize (): wrong type argument `ex'
error: octave_base_value::resize (): wrong type argument `<unknown type>'
octave-3.2.3.exe:4:C:\Octave\3.2.3_gcc-4.4.0\bin
オクターブでシンボリック行列をどのように宣言しますか？
this は役に立ちますか？
symbolic toolbox package が必要なようです。参照 here を参照してください。
シンボリックパッケージがない場合は、ダウンロードしてください。 Octaveコマンドライン、またはguiコマンドラインから。例えば.
octave> pkg install -forge symbolic
pythonとsympyがインストールされている場合は、オクターブフォージからパッケージがインストールされます。私はGoogleを使用してsympyをインストールする方法を見つけました。助けが必要な場合は私に連絡してください。
シンボリックパッケージがインストールされている場合は、「pkg load」を使用してパッケージ関数をインポートし、次にsyms関数を使用してシンボルを宣言します。
octave> pkg load symbolic
octave> syms a b
これは、シンボルaとbを定義しました。
octave> syms
Symbolic variables in current scope:
a
b
「syms」自体は、定義したすべてのシンボルを出力します。
octave> mat = [a,b]
mat = (sym) [a b] (1×2 matrix)
octave:34> mat * 2
ans = (sym) [2⋅a 2⋅b] (1×2 matrix)
このパッケージは、ロボットマニピュレータクラスの回転行列の計算に非常に役立ちました。お役に立てれば。
その他の例のスクリプトの一部を次に示します。
pkg load symbolic
syms psi phi theta psidot phidot thetadot
RzPsi = [[cos(psi), -sin(psi), 0]; [sin(psi), cos(psi), 0]; [0,0,1]]
RyTheta = [[cos(theta), 0, sin(theta)];[0,1,0];[-sin(theta), 0, cos(theta)]]
RzPhi = [[cos(phi), -sin(phi), 0]; [sin(phi), cos(phi), 0]; [0,0,1]]
RzPsi = (sym 3×3 matrix)
⎡cos(ψ) -sin(ψ) 0⎤
⎢ ⎥
⎢sin(ψ) cos(ψ) 0⎥
⎢ ⎥
⎣ 0 0 1⎦
RyTheta = (sym 3×3 matrix)
⎡cos(θ) 0 sin(θ)⎤
⎢ ⎥
⎢ 0 1 0 ⎥
⎢ ⎥
⎣-sin(θ) 0 cos(θ)⎦
RzPhi = (sym 3×3 matrix)
⎡cos(φ) -sin(φ) 0⎤
⎢ ⎥
⎢sin(φ) cos(φ) 0⎥
⎢ ⎥
⎣ 0 0 1⎦
octave> RzPhi * RyTheta
ans = (sym 3×3 matrix)
⎡cos(φ)⋅cos(θ) -sin(φ) sin(θ)⋅cos(φ)⎤
⎢ ⎥
⎢sin(φ)⋅cos(θ) cos(φ) sin(φ)⋅sin(θ)⎥
⎢ ⎥
⎣ -sin(θ) 0 cos(θ) ⎦
シンボリックツールボックスをインストールした後（Sudo apt-get install octave-symbolic）、次のことを行う必要があります。
symbols
x = sym('x')
ただし、シンボリック式を操作するOctaveの関数は、MATLABの関数よりもはるかに悪いことに注意してください。
Octave用のSymbolic Toolboxは多かれ少なかれ役に立たない。あなたの場合のように行列のサイズを変更することはできません。「-」演算子を使用することはできません。たとえば、単純なシンボリック演算を区別できます。
octave:1> symbols
octave:2> q1 = sym("q1");
octave:3> differentiate(Sin(q1)*Cos(q1),q1)
ans =
-sin(q1)^2+cos(q1)^2
しかし、これを行おうとすると：
octave:6> -Sin(q1)
error: unary operator `-' not implemented for `ex' operands
octave:6> -q1
error: unary operator `-' not implemented for `ex' operands
同じことがあなたの場合にも起こります、つまり、シンボリック値を含む行列のサイズを変更します
後世のもう一つの例。
このオクターブスクリプトを開発して実行するために http://octave-online.net/ を使用しました。
注意：結果を示すコメントとして出力を含めました。
disp("2-state markov chain symbolic analysis");
syms lambda mu
L = [lambda,0]
# L = (sym) [λ 0] (1×2 matrix)
U = [1;0]
#U =
# 1
# 0
C = [ [1,1]; [lambda,-mu]]
#C = (sym 2×2 matrix)
# ⎡1 1 ⎤
# ⎢ ⎥
# ⎣λ -μ⎦
C^-1
#ans = (sym 2×2 matrix)
# ⎡ λ -1 ⎤
# ⎢────── + 1 ──────⎥
# ⎢-λ - μ -λ - μ⎥
# ⎢ ⎥
# ⎢ -λ 1 ⎥
# ⎢ ────── ──────⎥
# ⎣ -λ - μ -λ - μ⎦
P = C^-1 * U
#P = (sym 2×1 matrix)
#
# ⎡ λ ⎤
# ⎢────── + 1⎥
# ⎢-λ - μ ⎥
# ⎢ ⎥
# ⎢ -λ ⎥
# ⎢ ────── ⎥
# ⎣ -λ - μ ⎦
lambda_sys = L * C^-1 * U
#lambda_sys = (sym)
#
# ⎛ λ ⎞
# λ⋅⎜────── + 1⎟
# ⎝-λ - μ ⎠
ハンドルの配列
Octave Struct Array を使用して、次のようなシンボリック行列を作成できます。
_b(1,1).vector = @sin;
b(1,2).vector = @cos;
b(2,1).vector = @sec;
b(2,2).vector = @csc;
b
b.vector
printf( "\n\nCalling each element:\n" )
b(1,1).vector
b(1,2).vector
b(2,1).vector
b(2,2).vector
printf( "\nCalculatin the sin of 1:\n" )
b(1,1).vector(1)
_これを実行すると、以下が返されます：
_b =
2x2 struct array containing the fields:
vector
ans = @sin
ans = @sec
ans = @cos
ans = @csc
Calling each element:
ans = @sin
ans = @cos
ans = @sec
ans = @csc
Calculatin the sin of 1:
ans = 0.841470984807897
_シンボルの配列
_@sin_、_@cos_などをsym("a")、sym("b")、sym("c")などに置き換えることができます。
_pkg load symbolic;
b(1,1).vector = sym("a");
b(1,2).vector = sym("b");
b(2,1).vector = sym("c");
b(2,2).vector = sym("d");
b
b.vector
printf( "\n\nCalling each element:\n" )
b(1,1).vector
b(1,2).vector
b(2,1).vector
b(2,2).vector
_これを実行すると、以下が返されます：
_b =
2x2 struct array containing the fields:
vector
ans = (sym) a
ans = (sym) c
ans = (sym) b
ans = (sym) d
Calling each element:
ans = (sym) a
ans = (sym) b
ans = (sym) c
ans = (sym) d
_参照：