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

38 visualizzazioni (ultimi 30 giorni)
和希
和希 il 7 Lug 2023
Commentato: 和希 il 23 Nov 2023
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の位まで表示したい)ため周波数(横軸)とゲイン(縦軸)の分解能を上げる方法があれば教えてください.

Risposta accettata

covao
covao il 20 Nov 2023
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 Commento
和希
和希 il 23 Nov 2023
ありがとうございます。うまくボード線図を表示できました。

Accedi per commentare.

Più risposte (0)

Categorie

Scopri di più su Control System Toolbox 入門 in Help Center e File Exchange

Tag

Prodotti


Release

R2023a

Community Treasure Hunt

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

Start Hunting!