Matlab实现一维和二维扩散方程

收录时间:2014-04-15
资源分类:Matlab 工具:MATLAB 7.10 (R2010a)

该程序采用有限差分方法隐式和显式)仿真了一维和二维扩散方程

扩散方程是一类偏微分方程, 用来描述扩散现象中的物质密度的变化. 通常也用来和扩散类似的现象。

有限差分法是微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替, 这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似, 积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组 , 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。
在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题,即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下:
1、区域离散化,即把所给偏微分方程的求解区域细分成由有限个格点组成的网格;
2、近似替代,即采用有限差分公式替代每一个格点的导数;
3、逼近求解。换而言之,这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程

Diffusion in 1D and 2D

The diffusion equation is simulated using finite differencing methods (both implicit and explicit) in both 1D and 2D domains. In both cases central difference is used for spatial derivatives and an upwind in time.

文件下载列表
Diffusion.zip (3.31KB)  
附件内容(只显示3中的2个)
Diffusion_1D.m  Diffusion_2D.m  
标签: 扩散方程 有限差分法 
更多

目前尚无评论

用户反馈   关于我们
Copyright (©) ZHIHUISHI.COM 2013 All Rights Reserved.
京ICP备14040384号