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

38 vues (au cours des 30 derniers jours)
和希
和希 le 7 Juil 2023
Commenté : 和希 le 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の位まで表示したい)ため周波数(横軸)とゲイン(縦軸)の分解能を上げる方法があれば教えてください.

Réponse acceptée

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

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Control System Toolbox 入門 dans Help Center et File Exchange

Tags

Produits


Version

R2023a

Community Treasure Hunt

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

Start Hunting!