如何在Octave中声明符号矩阵？

Question

在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

如何在八度音程中声明符号矩阵？

Answer 1

如果您还没有符号包,请下载它.从Octave命令行或gui命令行.例如

octave> pkg install -forge symbolic

如果您安装了python和sympy,那么将从octave forge为您安装软件包.我用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(?)    ?

Answer 2

这会有帮助吗？

看起来您可能需要符号工具箱包,请参考此处.

Answer 3

安装符号工具箱(您可以通过发布在某些环境中执行此操作sudo apt-get install octave-symbolic)后,您必须执行以下操作:

symbols
x = sym('x')

但请注意,Octave的操作符号表达式的功能比MATLAB更糟糕.