import matplotlib.pyplot as plt
import numpy as np
class diffusion:
def __init__(self,constant=0.3,l=500,update_times=50):
self.l=l
self.c=constant
self.ss=[[0]*self.l]
self.n=update_times
def diff(self):
self.ss[-1][int(self.l/2)]=1
temp=[0]*self.l
counter=1
while(1):
for i in range(len(self.ss[-1])-2):
temp[i+1]=self.ss[-1][i+1]+self.c*(self.ss[-1][i+2]+self.ss[-1][i]-2*self.ss[-1][i+1])
self.ss.append(temp)
counter+=1
if counter>self.n:
break
def show(self):
x=np.arange(0,self.l,1)
plt.plot(x,self.ss[-1],'.',label='time=%.2f'%(len(self.ss)-1))
plt.legend(frameon=True)
plt.grid(True)
plt.xlabel('x')
plt.ylabel('density')
plt.title('diffusion in one dimension')
a=diffusion(update_times=100)
a.diff()
a.show()
a=diffusion(update_times=500)
a.diff()
a.show()
a=diffusion(update_times=1000)
a.diff()
a.show()
a=diffusion(update_times=10000)
a.diff()
a.show()
a=diffusion(update_times=50000)
a.diff()
a.show()
代码5
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 审查日期：2017-08-11 审查人员：鸣人 代码文件：SystemSettingAboutActivity（关...
- Alamofire.request request 函数签名 request 函数实现 SessionManage...
