伝達行列の行列式を計算し、ナイキスト プロットを作成してナイキスト解析を行うつもりですが、行列式コマンドが伝達行列を認識しないという問題があります。コードを以下に示します
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」コマンドは避けたいと思います。その場合、ナイキスト プロットを実行できないからです。同じのナイキストプロットを計算する別の方法はありますか?