最近学现代控制,需要求解 G T P G − P = − Q {G^T}PG - P = - Q GTPG−P=−Q这种形式的矩阵方程,如下图。维数一高,手算简直要命。 下面给出利用MATLAB进行求解的方法,需要用到solve函数。
输入系数矩阵: >> 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 求解 >> 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] 到工作区查看求解结果