如何在python中使用传递函数反转矩阵

Question

我知道在matlab中我可以执行以下操作：

s = tf('s')
G11 = (s + 1)/(s + 2)
G12 = 1/(2*s + 1)
G21 = 1/(3*s + 1)
G22 = 1/(4*s + 1)
A = [G11 G12; G21, G22]
Ai = inv(A)
bode(A)

它会很好用的。在python中，我试图做一些类似的事情：

import control as co
import numpy as np
s = co.tf('s')
G11 = (s + 1)/(s + 2)
G12 = 1/(2*s + 1)
G21 = 1/(3*s + 1)
G22 = 1/(4*s + 1)
A = np.array([[G11, G12], [G21, G22]])
Ai = np.linalg.inv(A)
co.bode(A)

但这不起作用——numpy不知道如何反转这个矩阵。

在python中有一个好方法可以做到这一点吗？我知道我可以使用以s为符号的scipy，但我认为这对我使用控制工具箱中的其他工具没有帮助。

编辑：

numpy返回以下错误：

---------------------------------------------------------------------------
UFuncTypeError                            Traceback (most recent call last)
<ipython-input-1-ec46afd90eb6> in <module>
10 
11 A = np.array([[G11, G12], [G21, G22]])
---> 12 Ai = np.linalg.inv(A)
13 co.bode(A)
<__array_function__ internals> in inv(*args, **kwargs)
/usr/local/lib/python3.7/dist-packages/numpy/linalg/linalg.py in inv(a)
543     signature = 'D->D' if isComplexType(t) else 'd->d'
544     extobj = get_linalg_error_extobj(_raise_linalgerror_singular)
--> 545     ainv = _umath_linalg.inv(a, signature=signature, extobj=extobj)
546     return wrap(ainv.astype(result_t, copy=False))
547 
UFuncTypeError: Cannot cast ufunc 'inv' input from dtype('O') to dtype('float64') with casting rule 'same_kind'

Answer 1

Numpy(提示：它就在名称中(只是一个数字库；它不做符号运算。Sympy(也在名字中(做符号数学，所以使用它：

import sympy
s = sympy.Symbol('s', imaginary=True)
g11 = (s + 1)/(s + 2)
g12 = 1/(2*s + 1)
g21 = 1/(3*s + 1)
g22 = 1/(4*s + 1)
A = sympy.Matrix((
(g11, g12),
(g21, g22),
))
sympy.pprint(A.inv())

带输出

⎡    3       2                       3       2              ⎤
⎢ 6⋅s  + 17⋅s  + 11⋅s + 2      - 12⋅s  - 31⋅s  - 15⋅s - 2   ⎥
⎢ ───────────────────────      ──────────────────────────   ⎥
⎢     3      2                      3      2                ⎥
⎢  6⋅s  + 7⋅s  - 3⋅s - 1         6⋅s  + 7⋅s  - 3⋅s - 1      ⎥
⎢                                                           ⎥
⎢     3       2                 4       3       2           ⎥
⎢- 8⋅s  - 22⋅s  - 13⋅s - 2  24⋅s  + 50⋅s  + 35⋅s  + 10⋅s + 1⎥
⎢─────────────────────────  ────────────────────────────────⎥
⎢     3      2                      3      2                ⎥
⎣  6⋅s  + 7⋅s  - 3⋅s - 1         6⋅s  + 7⋅s  - 3⋅s - 1      ⎦

Answer 2

看起来control.tf返回了control.TransferFunction类的对象。这与返回符号函数对象的MATLAB版本不同。

通过查看文档，我没有看到将control.TransferFunction对象转换为符号函数对象的内置方法，但我确实看到了num和den方法，您可以使用这些值构造符号函数。然后你可以应用Reinderien的答案。