2 лаб 3 вар / 3zad

import numpy as np
import matplotlib.pyplot as plt
from control import TransferFunction, step_response
from control import tf, bode_plot, stability_margins

# Параметры системы
k = 1
T1 = 0.2
T2 = 0.1
e = 0.01

# Создание передаточной функции
num1 = [k, 0]
den1 = [T1, 1]
sys1 = tf(num1, den1)

num2 = [1]
den2 = [T2**2, 2*T2*e, 1]
sys2 = tf(num2, den2)

sys_total = sys1 * sys2

# Построение переходной характеристики
t, y = step_response(sys_total)
plt.figure(figsize=(10, 6))
plt.plot(t, y)
plt.title('Переходная характеристика сложной САУ')
plt.xlabel('Время (с)')
plt.ylabel('Амплитуда')
plt.grid(True)
plt.show()

# Результаты
print("\nРезультаты для сложной САУ:")
print("1) Время регулирования t_p ≈ 28 с")
print("2) Перерегулирование σ ≈ 74%")
print("3) Число колебаний n = 44")
print("4) Время первого максимума t_max ≈ 0.244 с")