bode関数を用いて​出力したボード線図の​周波数とゲインの分解​能を高める方法

bode関数を用いて伝達関数G(s)=1/(s^2 + T1s+T2)のボード線図を出力しようと考えています.そこで以下のようなコードを書きました.
%横軸Hz表示
opts=bodeoptions();
opts.FreqUnits='Hz';
%パラメータ
T1=1;
T2=300;
num = 1;%分子
den = [1 T1 T2];%分母
sysc = tf(num, den);%伝達関数
bode(sysc,opts);
上記を実行すると下記のボード線図が得られます.横軸の周波数はHzで表しています.ボード線図をクリックすることで周波数やゲインを見ることができますがその値をより細かく見たい(例えば周波数を1/1000の位まで表示したい)ため周波数(横軸)とゲイン(縦軸)の分解能を上げる方法があれば教えてください.

 Respuesta aceptada

covao
covao el 20 de Nov. de 2023

0 votos

bode関数の周波数[rad/s]はベクトルで設定することができます。
logspace関数を使うと対数的に等間隔なベクトルを生成できます。
下記の例は、周波数10^-1から10^2 [Hz]、1000点をプロットしています。
%横軸Hz表示
opts=bodeoptions();
opts.FreqUnits='Hz';
%パラメータ
T1=1;
T2=300;
num = 1;%分子
den = [1 T1 T2];%分母
sysc = tf(num, den);%伝達関数
n=1000;
w = logspace(-1,2,n)*2*pi();
bode(sysc,w,opts,"+-");

1 comentario

和希
和希 el 23 de Nov. de 2023
ありがとうございます。うまくボード線図を表示できました。

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Control System Toolbox 入門 en Centro de ayuda y File Exchange.

Productos

Versión

R2023a

Etiquetas

Preguntada:

el 7 de Jul. de 2023

Comentada:

el 23 de Nov. de 2023

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!