dz.py - master - PollyM/Group_Homework_Python

dz.py
79 строк2 KB
Михаил Ершов
Загрузка dz.py
07 ноя 2025, 22:53
07 ноя 2025, 22:532e755cd
# Задаём переменные со значениями по умолчанию

default_m1 = 1
default_m2 = 2
default_c1 = 200
default_c2 = 300
default_c3 = 500


def err01(arg_m1):
    if abs(float(arg_m1)) <= 10 ** (-6):
        arg_m1 = default_m1
    return float(arg_m1)


def err02(arg_m2):
    if abs(float(arg_m2)) <= 10 ** (-6):
        arg_m2 = default_m2
    return float(arg_m2)


def err03(arg_c1):
    if abs(float(arg_c1)) <= 10 ** (-6):
        arg_c1 = default_c1
    return float(arg_c1)


def err04(arg_c2):
    if abs(float(arg_c2)) <= 10 ** (-6):
        arg_c2 = default_c2
    return float(arg_c2)


def err05(arg_c3):
    if abs(float(arg_c3)) <= 10 ** (-6):
        arg_c3 = default_c3
    return float(arg_c3)


m1 = err01(input())
m2 = err02(input())
c1 = err03(input())
c2 = err04(input())
c3 = err05(input())

print(f"Масса тела 1 составляет {m1} кг")
print(f"Масса тела 2 составляет {m2} кг")
print(f"Жёсткость пружины 1 составляет {c1} Н/м")
print(f"Жёсткость пружины 2 составляет {c2} Н/м")
print(f"Жёсткость пружины 3 составляет {c3} Н/м")


# коэффициенты a11, a12, a22, c11, c12, c22
a11 = 1.5 * m1
a12 = 0
a22 = m2
c11 = c1
c12 = c1
c22 = c1 + c2 + c3

print(f"Коэффициент a11 = {a11}")
print(f"Коэффициент a12 = {a12}")
print(f"Коэффициент a22 = {a22}")
print(f"Коэффициент c11 = {c11}")
print(f"Коэффициент c12 = {c12}")
print(f"Коэффициент c22 = {c22}")

# Численное интегрирование с помощью решателя Рунге-Кутта 45

t = 0.000  # начальный момент времени
dt = 0.1  # шаг интегрирования
t_end = 5.000  # время интегрирования

while t < t_end:
    print(f"t = {t}; ")
    t += dt

# y(1) - перемещение тела 1, y(2) - скорость тела 1
# y(3) - перемещение тела 2, y(4) - скорость тела 2