如何解决如何在python中将变量存储在矩阵中
一个快速的问题:如何在矩阵内存储变量或数学表达式?例如:
对于一维数组: T [1,2] =(T [0,1] + T [2,3])/ 2 T [2,3] =(T [1,2] + T [3,4])/ 2 ...
然后我解决这个方程组
P.S:我正在尝试创建一个2D矩阵来模拟和解决矩形板上的热扩散。
所以基本上我使用的是Central Eular方法,我需要构建网格所有点的方程组,显示它们,然后求解。我当然有边界条件和值。
这大约是我过去在Maple中执行的代码逻辑:
L = 15
H = 15
dx = 3
dy = 3
ndx = int(L/dx)
ndy = int(H/dy)
i_max = ndx-1
j_max = ndy-1
Tb = 10
Tl = 40
Tr = 60
Tt = 90
N = (ndx-1)*(ndy-1)
# T = np.ones(([ndx,ndy]))
T = np.empty(ndx,ndy,dtype=np)
# boundary conditions assignement:
T[:,0] = Tl
T[:,-1] = Tr
T[0,:] = Tt
T[-1,:] = Tb
T[0,0] = (Tt + Tl)/2
T[0,-1] = (Tt + Tr)/2
T[-1,0] = (Tb + Tl)/2
T[-1,-1] = (Tb + Tr)/2
# Equations:
for i in range(1,i_max):
for j in range(1,j_max):
T[i,j] = (T[i-1,j] + T[i+1,j] + T[i,j-1] + T[i,j+1])/(-4)
print(T[i,j])
# Solution:
for i in range(1,j_max):
s = np.solve(T[i,j],[i,j])
print(s)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。