python扩散模拟_扩散模拟 - 问答

我试图模拟扩散过程，并有以下代码来模拟扩散方程：dx = 0.1

dt = 0.1

t = np.arange(0, 10, dt)

x = np.arange(0, 10, dx)

D = 1/20

k = 1

# We have an empty array

Cxt = np.tile(np.nan, (len(t), len(x)))

# Definition of concentration profile at t = 0.

Cxt[0] = np.sin(k*2*np.pi*x/10)+1

for j in range(len(t) - 1):

# Second derivative to x: C_xx

C_xx = (np.roll(Cxt[j], -1) + np.roll(Cxt[j], 1) - 2*Cxt[j]) / dx**2

# Concentrationprofile in the next time step

Cxt[j+1] = Cxt[j] + dt * D * C_xx

# Plot the concentration profiles in qt

%matplotlib qt

plt.waitforbuttonpress()

for i in range(len(t)):

ti = t[i]

Ci = Cxt[i]

plt.cla()

plt.plot(x, Ci, label='t={}'.format(ti))

plt.xlabel('x')

plt.ylabel('C(x)')

plt.axis([0, 10, 0, 2])

plt.title('t={0:.2f}'.format(ti))

plt.show()

plt.pause(0.01)

%matplotlib inline

我想看看正弦最大值消失得有多快。为了做到这一点，我想把振幅(最大值和平均值之间的距离)绘制成时间的函数，但是我怎么做呢？在

我怎么知道振幅在什么时候比开始的时候小呢？在

python扩散模拟_扩散模拟 - 问答 - Python中文网相关推荐

python编程基础_月隐学python第2课
python编程基础_月隐学python第2课学习目标掌握变量的输入和输出掌握数据类型的基本概念掌握算数运算 1.变量的输入和输出 1.1 变量输入使用input输入 input用于输入数据 ...
查看Python的版本_查看当前安装Python的版本
一.查看Python的版本_查看当前安装Python的版本具体方法: 首先按[win+r]组合键打开运行: 然后输入cmd,点击[确定]: 最后执行[python --version]命令即可. 特 ...
Matlab 李萨如图模拟_动态模拟李萨如图形形成过程保存为.GIF动图
九层妖塔起于垒土 Matlab 李萨如图模拟_动态模拟李萨如图形形成过程保存为.GIF动图一.数学分析二.程序设计及使用说明三.动图四.代码一.数学分析二.程序设计及使用说明三. ...
随机模拟_随机模拟可帮助您掌握统计概念
随机模拟模拟有助于提炼概念 (Simulation helps distilling concepts) 掌握与统计相关的概念可能很困难 (Grasping statistics-related c ...
python 对象锁_也许你对 Python GIL 锁的理解是错的。
摄影:产品经理甜白与草莓更配~ 我刚到现在这个公司时,听到当时一个高级工程师(现已离职)大声地跟他旁边的同事说:  Python 有 GIL 锁,所以它的多线程实际上是单线程,所以写多线程代码不用考 ...
适合小学生python的书_小学生开始学Python,最接近AI的编程语言:安利一波Python书单...
原标题:小学生开始学Python,最接近AI的编程语言:安利一波Python书单 AlphaGo 都在使用的 Python 语言,是最接近 AI 的编程语言. 教育部考试中心近日发布了"关于 ...
python基本原理概论_怎样开始自学Python？
本人才疏学浅,学识大多浅尝辄止,故文章若有错误,不论是文字笔误还是理解有错,烦请您留言以告知,本人必定感激不尽! **Python分类下的系列文章,不断更新中,如果你迫不及待地想要看看写得如何可以先试 ...
数据结构python课后答案_数据结构与算法:Python语言描述 1~5章课后习题
数据结构与算法:Python语言描述 1~5章课后习题发布时间:2018-07-19 20:42, 浏览次数:1885 , 标签: Python MarkDown语法写的,不知道为啥上传到CSDN不 ...
python opencv手册_教你用Python实现5毛钱特效（给你的视频来点料）
一.前言请务必看到最后.Python牛已经不是一天两天的事了,但是我开始也没想到,Python能这么牛.前段时间接触了一个批量抠图的模型库,而后在一些视频中找到灵感,觉得应该可以通过抠图的方式,给视 ...
python并行运算库_最佳并行绘图Python库简介：“ HiPlot”
python并行运算库 HiPlot is Facebook's Python library to support visualization of high-dimensional data ta ...

python扩散模拟_扩散模拟 - 问答 - Python中文网

python扩散模拟_扩散模拟 - 问答 - Python中文网相关推荐

最新文章

热门文章