【MATLAB】求解矩阵方程

最近学现代控制，需要求解$G^{T}PG-P=-Q$这种形式的矩阵方程，如下图。维数一高，手算简直要命。

下面给出利用MATLAB进行求解的方法，需要用到solve函数。

1. 输入系数矩阵：

>> p=[1 4 0;-3 -2 -3;2 0 0]

p =

     1     4     0
    -3    -2    -3
     2     0     0

>> pt=transpose(p)

pt =

     1    -3     2
     4    -2     0
     0    -3     0

2. 求解

>> solve(pt*x*p-x==-eye(3))

ans = 

  包含以下字段的 struct:

    x1_1: [1×1 sym]
    x1_2: [1×1 sym]
    x1_3: [1×1 sym]
    x2_1: [1×1 sym]
    x2_2: [1×1 sym]
    x2_3: [1×1 sym]
    x3_1: [1×1 sym]
    x3_2: [1×1 sym]
    x3_3: [1×1 sym]

3. 到工作区查看求解结果