2 лаб 3 вар / laba2_1

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

# Параметры системы
T = 10
ep = 0.02
k = 1

lamb = ep / T
om = np.sqrt(1 - ep**2) / T

# Аналитическое решение
x = np.arange(0, 1600, 0.1)
f = k * (1 - np.exp(-lamb * x) * (np.cos(om * x) - (lamb / om) * np.sin(om * x)))

# График переходного процесса
plt.figure(figsize=(10, 6))
plt.plot(x, f, label='h(t)')
plt.hlines(y=[1, 0.95, 1.05], xmin=0, xmax=1600, colors=['green', 'red', 'red'], linestyles='solid')
plt.xlabel('t')
plt.ylabel('h(t)')
plt.title('Переходная характеристика колебательного звена')
plt.grid(True)
plt.legend(['h(t)', 'h stable', 'thr low', 'thr high'])
plt.show()

# Прямые показатели качества (примерные значения из задачи)
print("Результаты:")
print("1) Время регулирования t_p ≈ 1500 с")
print("2) Перерегулирование σ = 95%")
print("3) Частота ω = 0.1 рад/с")
print("4) Число колебаний n = 23")