2024.4.2 scipy.freqsによる周波数応答

アナログフィルタ（連続・一時遅れ系）のゲイン曲線をscipy.freqsを用いて求める。

$ G(s) = \frac{K}{Ts + 1}, ~~ K =1, T = 0.1

code:freqs01.py

from scipy import signal

import numpy as np

import matplotlib.pyplot as plt

T = 0.1

K = 3

w, h = signal.freqs(K, T, 1, worN=np.logspace(-2, 4, 1000))

#plt.semilogx(w, 20*np.log10(abs(h)))

https://scrapbox.io/files/660cf3bc4a7f9d0024fdaf71.png

code:freq02.py

from scipy import signal

import numpy as np

import matplotlib.pyplot as plt

omega_n = 10

zeta = 0.1

K = 1

w, h = signal.freqs(K, 1, 2*zeta*omega_n, omega_n**2, worN=np.logspace(-2, 4, 1000))

plt.semilogx(w, 20*np.log10(abs(h)))

plt.grid()

plt.show()