分类信息网站手机版,长沙市住房和建设局官方网站,设计小程序多少钱,建筑工程网站建设方案1、方程及其离散
1.1一维热传导方程 1.2离散化
设定步长#xff0c;依据上述方程得到递推关系#xff1a; 2、python求解实现
import numpy as np
import matplotlib.pyplot as plth 0.1#空间步长
N 30#空间步数
dt 0.0001#时间步长
M 10000#时间的步数
A dt/(h**2)…1、方程及其离散
1.1一维热传导方程 1.2离散化
设定步长依据上述方程得到递推关系 2、python求解实现
import numpy as np
import matplotlib.pyplot as plth 0.1#空间步长
N 30#空间步数
dt 0.0001#时间步长
M 10000#时间的步数
A dt/(h**2) #lambda*tau/h^2
U np.zeros([N1,M1])#建立二维空数组
Space np.arange(0,(N1)*h,h)#建立空间等差数列从0到3公差是h#边界条件
for k in np.arange(0,M1):U[0,k] 0.0U[N,k] 0.0#初始条件
for i in np.arange(0,N):U[i,0]4*i*h*(3-i*h)#递推关系
for k in np.arange(0,M):for i in np.arange(1,N):U[i,k1]A*U[i1,k](1-2*A)*U[i,k]A*U[i-1,k]
#不同时刻的温度随空间坐标的变化
plt.plot(Space,U[:,0], g-, labelt0,linewidth1.0)
plt.plot(Space,U[:,3000], b-, labelt3/10,linewidth1.0)
plt.plot(Space,U[:,6000], k-, labelt6/10,linewidth1.0)
plt.plot(Space,U[:,9000], r-, labelt9/10,linewidth1.0)
plt.plot(Space,U[:,10000], y-, labelt1,linewidth1.0)
plt.ylabel(u(x,t), fontsize20)
plt.xlabel(x, fontsize20)
plt.xlim(0,3)
plt.ylim(-2,10)
plt.legend(locupper right)
plt.show()
