ЛР_0026-2 (Анализ рекурсивных цифровых фильтров 1-го и 2-го порядка) / cos_26_2_ach_dlya_punkta_3_FIX

import matplotlib.pyplot as plt
import numpy as np
###

### АЧХ В ФИЛЬТРАХ, С 5А ПО 10А.

a1 = 0
a2 = 0.9
b0 = 1 # У ВСЕХ ОДИНАКОВО. НЕ МЕНЯЕМ

f = np.linspace(0, 8000, 8000) # НЕ ТРОШ, ЗАРУИНИШЬ
T = 1/8000 # НЕ ТРОШ, ЗАРУИНИШЬ

im = a1 * np.sin(2 * 3.14 * f * T) + a2 * np.sin(4 * 3.14 * f * T) # НЕ ТРОШ, ЗАРУИНИШЬ
re = 1 - a1 * np.cos(2 * 3.14 * f * T) - a2 * np.cos(4 * 3.14 * f * T) # НЕ ТРОШ, ЗАРУИНИШЬ
ach = b0/(np.power(re, 2) + np.power(im, 2)) # НЕ ТРОШ, ЗАРУИНИШЬ

###

if __name__ == '__main__':

plt.plot(f, ach, color="blue")
plt.axhline(0, color='black', linestyle='-')
plt.xticks([0, 1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000])
plt.grid(True)

plt.xlabel('f, Частота [Гц]', fontsize=9, fontweight='bold')
plt.ylabel('Значение', fontsize=9, fontweight='bold')
plt.title('АЧХ')

plt.show()