拉普拉斯方程

拉普拉斯方程

-

2023年3月19日发(作者：世界地铁)

⼆维拉普拉斯⽅程的数值解法

importnumpyasnp

asplt

m=101

n=101

U=([m,n])

F=([m,n])

#定义边界条件

F[m-1,:]=100

e=0.001

count=0

foriinrange(10000):

count+=1

F[:,0]=0

F[:,n-1]=0

F[0,:]=0

F[m-1,:]=100

U[1:-1,1:-1]=(F[1:-1,:-2]+F[1:-1,2:]+F[:-2,1:-1]+F[2:,1:-1])/4

e=abs(U-F).max()#求助⽼师我这个值为何总是0

F=U

x=(0,m)

y=(0,n)

X,Y=id(x,y)

C=r(X,Y,U)

(C,inline=True,fontsize=10)

参考⽹址：