数式の練習です。
伝達関数は、二次ローパス・フィルタの分子を$s^2$に変更すると、二次ハイパス・フィルタになります。
二次ハイパス・フィルタ
| 二次ハイパス・フィルタ 回路図 |
|---|
 |
| 伝達関数 |
| $T(s)=\LARGE \frac{s^2}{ω_0^2 + 2sζω_0 + s^2}$ |
| $( s はラプラス変数。 jω はsとする。 ω_0は角周波数(角振動数?)。 $ $ζは減衰係数。Q はフィルタのクオリティ・ファクタ。ζ = \frac{2}{Q}) $ |
| TeXの記述 |
$T(s)=\LARGE \frac{s^2}{ω_0^2 + 2sζω_0 + s^2}$ |
Qの取る値によって、フィルタの特性は大きく変わるのは、二次ローパス・フィルタと同じです。
matlabとLTspcieでも、$Q$（$ζ$）を2通りに変化させています。
matlab
time_constant = 0.16; % 時定数
natural_freq = 1 / time_constant; % 固有角振動数 ωn
damping_ratio = 1.0 % ζ（ゼータ）減衰係数（ダンピング比）。１＞；実数解、1；実数解、・・・1/√2；複素数解、1/√2＜；複素数解、ピークあり
G = tf([1, 0, 0 ], [1, 2 * damping_ratio * natural_freq, natural_freq^2])
h = bodeplot(G);
p = getoptions(h);
p.FreqUnits = 'Hz';
setoptions(h,p);
grid on
hold on
damping_ratio = 0.15 % ζ（ゼータ）減衰係数。１＞；実数解、1；実数解、・・・1/√2；複素数解、1/√2＜；複素数解、ピークあり
G2 = tf([1, 0, 0], [1, 2 * damping_ratio * natural_freq, natural_freq^2])
h = bodeplot(G2);
p = getoptions(h);
p.FreqUnits = 'Hz';
setoptions(h,p);
hold off
nyquistplot(G);
hold on
nyquistplot(G2);
grid on
実行例です。
LTspice
環境
matlab R2021b update1（ホームライセンス）、tr関数はControl System Toolboxが必要。ラプラス関数のSymbolic Math Toolboxはホームライセンスでは購入できない。