1

我打算计算传递矩阵的行列式，然后通过制作奈奎斯特图进行奈奎斯特分析，但问题是行列式命令无法识别传递矩阵。代码如下所示

clc
clear all;
close all;

g11 = tf(12.8,[16.7 1],'InputDelay',1)
g12 = tf(-18.9,[21 1],'InputDelay',3)
g21 = tf(6.6,[10.9 1],'InputDelay',7)
g22 = tf(-19.4,[14.4 1],'InputDelay',3)

G=[g11 g12 ; g21 g22]

[re,im,w] = nyquist(G)

F=2.55;


s=tf('s');
%syms s;

ggc11 = g11*(0.96*(1+3.25*F*s)/(3.25*F^2*s))
ggc12 = g12*(0.534*(1+3.31*F*s)/(3.31*F^2*s))
ggc21 = g21*(0.96*(1+3.25*F*s)/(3.25*F^2*s))
ggc22 = g22*(0.534*(1+3.31*F*s)/(3.31*F^2*s))

GGc=[ggc11 ggc12 ; ggc21 ggc22];

L=eye(2)+ GGc;

W= -1 + det(L)

nyquist(W)

出现的错误如下

'ss' 类型的输入参数的未定义函数 'det'。

BLT_code 中的错误（第 30 行）W= -1 + det(L)

我想避免使用“syms”命令，因为那时我将无法进行奈奎斯特图。是否有任何替代方法来计算相同的奈奎斯特图？

4

1 回答 1

0

我被困在同一条船上，试图计算传递函数矩阵的行列式以检查 MIMO 奈奎斯特稳定性标准，请参阅MIMO Stability ETH Zurich Lecture slides (pg 10)。不幸的是，似乎没有一个简单的 MATLAB 命令。我认为它可以手动评估。

如果您有以下形式的 TF 矩阵 G(s)：

G = [g_11 g_12; g_21 g_22];

您可以通过根据其原始定义对其进行评估来获得行列式

det_G = g_11*g_22 - g_12*g_21;

这将产生一个 1x1 TF 变量。当然，这种方法对于 2x2 系统以上的任何东西来说都太复杂了。

于 2019-11-15T11:16:57.507 回答