MATLABでシンボリック式を評価するときにゼロ除算

Question

clear all 

syms s w 

G = 1/((s)*(s+1)*(s+2)); %transfer function 
G_w = subs(G,s,j*w); 
W= [-100:0.01:100]; %[min_range:step size:max_range] 
nyq = eval(subs(G_w,w,W)); 

x = real(nyq) 
y = imag(nyq) 

plot(x,y) 

このコードを実行できないようで、20行未満しかない行100 ++でエラーが表示され続けます。

Error using symengine (line 59) 
Division by zero. 
Error in sym/subs>mupadsubs (line 139) 
G = mupadmex('symobj::fullsubs',F.s,X2,Y2); 
Error in sym/subs (line 124) 
    G = mupadsubs(F,X,Y); 
Error in nyquist2 (line 8) 
nyq = eval(subs(G_w,w,W)); %replace W with w in equation G_w 

これは表示されているエラーですが、どの専門家がこれを助けることができますか？

Answer 1

エラーは、配列Wを使用してG_wを計算していて、この配列に値が含まれているためエラーが発生し、エラーが発生したためです。 symengineを用い

エラー（線59）
除算ゼロによる。

あなたはこれを回避するために何ができるか、epsでWで0を交換することです。サイドノートとして

% Replace zeros with epsilon 
W(W == 0) = eps; 

nyq = eval(subs(G_w,w,W)); 

x = real(nyq) 
y = imag(nyq) 

plot(x,y) 

は、エラーが あなたコードの行100+の問題に文句されるのではなく、スタックトレースは、エラーが実際に関数内 に由来することを述べてあなたは、エラーは、あなたがそれ

Error using symengine (line 59)    <--- WHERE THE ERROR HAPPENED 
Division by zero.       <--- THIS IS THE ERROR MESSAGE 

Error in sym/subs>mupadsubs (line 139)  <--- THIS FUNCTION CALLED symengine 
G = mupadmex('symobj::fullsubs',F.s,X2,Y2); <--- THIS IS THE LINE THAT CALLED symengine 

Error in sym/subs (line 124)    <--- THIS FUNCTION CALLED mupadsubs 
G = mupadsubs(F,X,Y);      <--- THIS IS THE LINE THAT CALLED mupadsubs 

Error in nyquist2 (line 8)     <--- THIS FUNCTION (YOURS) CALLED subs 
nyq = eval(subs(G_w,w,W))     <--- THIS IS THE LINE THAT CALLED subs