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の位まで表示したい)ため周波数(横軸)とゲイン(縦軸)の分解能を上げる方法があれば教えてください.

 採用された回答

covao
covao 2023 年 11 月 20 日

0 投票

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 件のコメント

和希
和希 2023 年 11 月 23 日
ありがとうございます。うまくボード線図を表示できました。

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

ヘルプ センター および File ExchangeControl System Toolbox 入門 についてさらに検索

製品

リリース

R2023a

タグ

質問済み:

2023 年 7 月 7 日

コメント済み:

2023 年 11 月 23 日

Community Treasure Hunt

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

Start Hunting!