2020数学建模a炉温曲线python代码

import numpy as np
import matplotlib.pyplot as plt

# 定义模型的参数
a = 0.02 # 动力系统的系数
b = 0.01 # 非线性扰动项的系数
c = 0.05 # 热泵系统的系数
d = 0.03 # 辐射损失项的系数
T_env = 25 # 外界温度

# 定义模型的初始条件
T0 = 50 # 初始温度
t_start = 0 # 起始时间
t_end = 100 # 结束时间
dt = 0.1 # 时间步长

# 定义函数来计算温度的变化率
def temperature_derivative(T):
return a * (T_env - T) - b * (T - T0) + c * (T - T_env) - d * (T - T_env) ** 2

# 使用欧拉方法求解微分方程
t = np.arange(t_start, t_end, dt)
T = np.zeros_like(t)
T[0] = T0

for i in range(1, len(t)):
T[i] = T[i-1] + dt * temperature_derivative(T[i-1])

# 绘制温度曲线
plt.plot(t, T)
plt.xlabel('时间')
plt.ylabel('温度')
plt.title('2020数学建模A题-炉温曲线')
plt.grid(True)
plt.show()