热传导问题的数值求解
📁 python, 问题, 错误, 示例, 人工智能
需要求解区间上的热传导问题
求解时使用了显式格式
N = 10 # 沿x的最大步数K = 10 # 沿t的最大步数l = 1 # 右边界x的值h = l / N # 沿x的网格步长T = 1 # 右边界时间t的最大值t = T / K # 时间网格步长# 定义网格x_i = np.arange(0, N, h) # 沿x的节点值t_j = np.arange(0, K, t) # 沿t的节点值r_j = len(t_j) # 沿t的节点数r_i = len(x_i) # 沿x的节点数w_h_t = np.zero