显式差分和隐式差分
s T i n 1 1 ( 2 2 s ) T i n 1 s T i n 1 1 s T i n 1 ( 2 2 s ) T i n s T i n 1
1
2
3
4
5
6
边界条件：
初始条件：
Ac(n1) Bc(n)
1 0
0
0
0
s
2 2s
s
0
0
0 s 2 2s s 0
A
0
0
s 2 2s s
0
0
0
0
c (n1)
T T
n 1
n 2
1 1
T T
n 3
n 4
1 1
0 0 0 0
0 s 2 2s s
0
0
0 1
T T
n 5
n 6
Ac(n1) Bc(n)
A = sparse(nx,nx);
1 0
0
0
0 0
s
2 2s
s
0
0
0
0 s 2 2s s 0 0
A
0
0
s 2 2s s
0
0 0 0 0
0 s 2 2s s
0
0
0 1
内部节点：
for i=2:nx-1 A(i,i-1) = -s; A(i,i ) = (2+2*s); A(i,i+1) = -s; end
j (i-1,j)
j-1 (i-1,j-1)
3
h2 0
(hi3,j)
h1
h44
(i,j)
(i+1,j+1)
1 (i+1,j) (i+1,j-1)
局部节点编号 总体节点编号
i-1
i i+1
U ( i 1 , j ) U ( i 1 , j ) U ( i , j 1 ) U ( i , j 1 ) 4 U ( i , j ) 0
1 (i+1,j) (i+1,j-1)
i-1
i i+1
U=0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 U=0
U ( i 1 , j ) U ( i 1 , j ) U ( i , j 1 ) U ( i , j 1 ) 4 U ( i , j ) 0
j+1 (i-1,j+1) (i,j+1) 2
Tj1Tj tf(Tj,tj)
Tj1Tj- tTj=T（ j 1- t）
Tj
T（0 1-
t
）j
现在改用向后差分形式进行近似，得到隐式差分格式：
显式差分格式
dT dt ttj
Tj Tj1O(t) t
Tj Tj1tf(Tj,tj)
f (T ,t) T
Tj Tj1(1t)1
隐式差分格式
Tj
T0(1
t
Cui, 2013
Cui, 2013
Cui, 2013
Cui, 2013
总结：
1、有限差分方法给出的数值解的精度取决于所用的差分形式（向 前、向后、中心）。
2、偏微分方程的显式有限差分格式通常是有条件稳定的，为了保 证得到精确的数值解，最关键的是需要根据稳定性条件选取正确的 空间和时间步长。
载荷项：
rhs = zeros(nx,1); rhs(2:nx-1) = Told(2:nx-1); rhs(1) = Tleft; rhs(nx) = Tright;
内部 边界
Crank-Nicolson 隐式差分格式的程序实现
s T i n 1 1 ( 2 2 s ) T i n 1 s T i n 1 1 s T i n 1 ( 2 2 s ) T i n s T i n 1
1 1
1 0
0
0
0 0
s
2 2s
s
0
0
0
0 s 2 2s s
0 0
B
0
0
s
2 2s
s
0
0 0
0
s 2 2s s
0 0
0
0
0 1
T T
n 1
n 2
c (n )
T T
n 3
n 4
T T
n 5
n 6
s T i n 1 1 ( 2 2 s ) T i n 1 s T i n 1 1 s T i n 1 ( 2 2 s ) T i n s T i n 1
B(1 ,1 ) = 1; B(nx,nx) = 1;
T T
n 1
n 2
c (n )
T T
n 3
n 4
T T
n 5
n 6
例子：牛顿冷却定律：温度高于周围环境的物体向周围媒质传递热量 逐渐冷却时所遵循的规律。当物体表面与周围存在温度差时，单位时 间从单位面积散失的热量与温度差成正比。
稳定性（stability）：如果偏微分方程的严格解析解有界，差分格式给出的 解也有界，称该差分格式是稳定的；如果差分格式给出的解是无界的，则 称该差分格式是不稳定的。
稳定性反映了差分格式在计算中控制误差传递的能力
收敛性（convergence）：如果当时间和空间步长趋于零时，FDE解趋于PDE
a(121,106)=-0.25; a(135,134)=-0.25; a(135,120)=-0.25; a(15,14)=-0.25; a(15,30)=-0.25; for i=2:14
a(i,i-1)=-0.25; a(i,i+1)=-0.25; a(i,i+15)=-0.25; end for i=122:134 a(i,i-1)=-0.25; a(i,i+1)=-0.25; a(i,i-15)=-0.25; end for i=1:7 for j=2:14; a(15*i+j,15*i+j-1)=-0.25; a(15*i+j,15*i+j+1)=-0.25; a(15*i+j,15*i+j+15)=-0.25; a(15*i+j,15*i+j-15)=-0.25; end end
b=a^(-1); c=zeros(135,1); for i=121:135
c(i,1)=25;end d=b*c; s=zeros(11,17); for i=2:16
s(11,i)=100; end for i=1:9 for j=1:15; s(i+1,j+1)=d(15*(i-1)+j,1); end end
一种隐式差分格式的程序实现
1
2
3
4
5
6
求解区域：
边界条件：
初始条件：
内部节点：
A = sparse(nx,nx); for i=2:nx-1 A(i,i-1) = -s; A(i,i ) = (1+2*s); A(i,i+1) = -s; end
边界节点：
A(1 ,1 ) = 1; A(nx,nx) = 1;
隐式差分格式（implicit difference scheme）
特点
G (Y(t),Y(t t))0
1. 时间步长和空间步长的选择不受限制；
2. 需要联立解方程组
例子： 1. 显式差分格式：
左端：n+1时刻的值； 右端：n时刻的值。
特点：结构简洁，直接求解，求解速度快。
但是，时间步长需满足：
显式差分格式才能得到稳定的数值解，否则，数值解将会不稳定而振荡。
显示差分格式示意图
2. 隐式差分格式：
时间一阶精度 空间二阶精度
隐式有限差分格式
Crank-Nicolson 隐式差分格式
T in 1 tT in 2 ( x )2(( T i n 1 1 2 T in 1 T i n 1 1 ) ( T i n 1 2 T in T i n 1 ))
a(i,i)=1; end; for i=1:7
a(15*i+1,15*i+2)=-0.25; a(15*i+1,15*i+16)=-0.25; a(15*i+1,15*i-14)=-0.25; end for i=1:7 a(15*i+15,15*i+14)=-0.25; a(15*i+15,15*i+30)=-0.25; a(15*i+15,15*i)=-0.25; End a(1,2)=-0.25; a(1,16)=-0.25; a(121,122)=-0.25;
解，称该差分格式是收敛的。
如果
h l0i,m t 0(Uim ,j,kuim ,j,k)0
则称该差分格式是收敛的。
收敛性描述的是当差分网格无限细化时，差分方程的解是否具有无限逼 近偏微分方程的解的能力
Lax等价定理（Lax equivalence theorem）：如果逼近一个给定问题的差 分格式是相容的，那么该差分格式的收敛性与稳定性互为充分必要条件。
subplot(1,2,1),mesh(s) axis([0,17,0,11,0,100]) subplot(1,2,2),contour(s,32)
11
10
100
9
80
8
60
7
40
6
20
5
4 0
10
3
15
5
2 10
5 00
1
显式算法与隐式算法的区别
显式算法与隐式算法的区别-CAL-FENGHAI.-(YICAI)-Company One1显式算法与隐式算法的区别所谓显式和隐式，是指求解方法的不同，即数学上的出发点不一样。
并不是说显式只能求动力学问题，隐式只能求静力学问题，只是求解策略不通。
显式求解是对时间进行差分，不存在迭代和收敛问题，最小时间步取决于最小单元的尺寸。
过多和过小的时间步往往导致求解时间非常漫长，但总能给出一个计算结果。
解题费用非常昂贵。
因此在建模划分网格时要非常注意。
隐式求解和时间无关，采用的是牛顿迭代法（线性问题就直接求解线性代数方程组），因此存在一个迭代收敛问题，不收敛就的不到结果。
两者求解问题所耗时间的长短理论上无法比较。
实际应用中一般感觉来说显式耗时多些。
由于两者解题的出发点，所以一般来说显式用于求解和时间相关的动力学问题。
隐式用来求解和时间无关的静力学问题。
但也不是绝对的。
比如，用隐式求解时，为了克服迭代不收敛，改用显式算，但是要多给点时间，这样虽然克服了不收敛的问题，但是求解的时间费用也是相当客观的。
另外，隐式也可以求解动力学问题。
这是ansys里面的两种求解方法。
大多数非线性动力学问题一般多是采用显式求解方法，特别是在求解大型结构的瞬时高度非线性问题时，显示求解方法有明显的优越性。
下面先简要对比一下隐式求解法和显示求解法。
动态问题涉及到时间域的数值积分方法问题。
在80年代中期以前，人们基本上采用纽曼法进行时间域的积分。
根据纽曼法，位移、速度和加速度有着如下关系：u(i+1)=u(i)+△t*v(i)[(1—2p)a(i)+2p*a(i+1)] (1)v(i+1)=V(i)+△t[(1-2q)a(i)+2qa(i+1)] (2)上面式子中 u(i+1),u(i)分别为当前时刻和前一时刻的位移，v(i+1)和V(i)为当前时刻和前一时刻的速度，a(i+1)和a(i)为当前时刻和前一时刻的加速度，p和q为两个待定参数，△t为当前时刻与前一时刻的时问差，符号 * 为乘号。
热传导方程的差分格式汇总
热传导方程的差分格式汇总1.显式差分格式：显式差分格式是最简单的一种方法，通过将导热方程时间和空间上的导数进行近似，引入差分算子，将方程转化为差分格式。
其中最常见的差分格式有：a. 前向差分法（Forward Difference Method）：利用当前节点和其相邻节点的温度值进行计算。
例如，在一维离散情况下，可以使用公式：u(i,j+1)=u(i,j)+α(u(i+1,j)-2u(i,j)+u(i-1,j))b. 后向差分法（Backward Difference Method）：利用当前节点和其相邻节点的温度值进行计算。
例如，在一维离散情况下，可以使用公式：u(i,j+1)=u(i,j)+α(u(i+1,j+1)-2u(i,j+1)+u(i-1,j+1))c. 中心差分法（Central Difference Method）：利用当前和其相邻节点的温度值进行计算。
例如，在一维离散情况下，可以使用公式：u(i,j+1)=u(i,j)+α(u(i+1,j)-2u(i,j)+u(i-1,j))+β(u(i+1,j)-u(i-1,j))其中α和β是时间和空间步长的比例因子。
2.隐式差分格式：显式差分格式具有较大的稳定性限制。
为了克服这个问题，可以使用隐式差分格式，其中使用下一个时间步长的温度值来求解当前时间步长。
常见的隐式差分格式有：a. C-N差分法（Crank-Nicolson Method）：利用前后两个时间步长的温度值进行计算。
例如，在一维离散情况下，可以使用公式：u(i,j+1)=u(i,j)+0.5α(u(i+1,j+1)-2u(i,j+1)+u(i-1,j+1))+0.5α(u(i+1,j)-2u(i,j)+u(i-1,j))b. 力学模拟法（Finite Element Method）：将空间离散化后，通过引入有限元方法，将热传导问题转化为线性方程组，再通过求解线性方程组得到温度分布。
第二章显式隐式格式
五、数值天气预报的概念和历 史回顾
numerical weather prediction uses numerical
methods to approximate a set of partially differential equations on discrete grid points in a finite area to predict the weather systems and processes in a finite area for a certain time in the future. In order to numerically integrate the partial differential equations, which govern the atmospheric motions and processes, with time, one needs to start the integration at certain time.
Fi n
收敛到 F ( xi , t n )
稳定性
差分近似的稳定性是指对于任意给定的初值， 当n无限增大时，任意时刻的数值是否有界的 问题。假如数值解是稳定有界的，则相应的数 值格式称为稳定的格式。 计算稳定性的分析方法： 冯纽曼方法（Von-Neumann方法，又称谐 波分析法）：通过测试差分格式近似解一个谐 波分量的稳定性，研究差分格式的稳定性。
收敛性convergence
设差分方程的解为 Fi n
微分方程的解为
n i n
F ( xi , t n )
Fi F ( xi , t n )
in 0
则称：差分方程的解
显式算法和隐式算法的介绍与比较
求解时间比较
使用显式方法，计算成本消耗与单元数量成正比，并且大致与最小 单元的尺寸成反比；分析中最消耗CPU的一项就是单元的处理。由于 CPU时间与积分的个数成正比，采用简化积分的单元便可极大地节省数 据存储量和运算次数，进而提高运算效率。但单点积分容易出现零能 模型且应力结果精确度与积分点相关。
（5）
U t t
Ut
Utt
[(1 2
)Ut Utt ]t2
（6）
将（5）（6）代入（1）中，得：
KˆUtt Rˆtt
其中：
Kˆ
1
t 2
M
t
C
K
（有效刚度矩阵）
Rˆtt
Rtt
1
M
[ t
2
Ut
1 t
Ut
(1 2
1)Ut ]
C[ t
U
t
(
1)Ut
( 2
1)tUt ]（有效载荷矢量）
可以看出，解当前 Utt，需要用到当前时刻 的 Rtt ，因此为隐式算法，需要用迭代实现。
特点：
当参数 0.5 ， 0.25(0.5 )2 时，Newmark法是无 条件稳定的，即 t 的大小不影响数值的稳定性。 此时时间步长的选择主要根据解的精度确定。
显式算法
总结
动态显式算法采用动力学方程的一些差分格式，不用直接求解切线 刚度，不需要进行平衡迭代，计算速度快，时间步长只要取的足够小， 一般不存在收敛性问题。因此需要的内存也比隐式算法要少。并且数值 计算过程可以很容易地进行并行计算，程序编制也相对简单。
LS-DYNA中中心差分法的步长控制
1D杆、梁单元：
显示算法和隐式算法单点积分算法和全积分算法
先说说显示算法和隐式算法：这是ansys里面的两种求解方法。
大多数非线性动力学问题一般多是采用显式求解方法，特别是在求解大型结构的瞬时高度非线性问题时，显示求解方法有明显的优越性。
下面先简要对比一下隐式求解法和显示求解法。
动态问题涉及到时间域的数值积分方法问题。
在80年代中期以前，人们基本上采用纽曼法进行时间域的积分。
根据纽曼法，位移、速度和加速度有着如下关系：u(i+1)=u(i)+△t*v(i)[(1—2p)a(i)+2p*a(i+1)] (1)v(i+1)=V(i)+△t[(1-2q)a(i)+2qa(i+1)] (2)上面式子中 u(i+1),u(i)分别为当前时刻和前一时刻的位移，v(i+1)和V(i)为当前时刻和前一时刻的速度，a(i+1)和a(i)为当前时刻和前一时刻的加速度，p和q为两个待定参数，△t为当前时刻与前一时刻的时问差，符号* 为乘号。
由式(1)和式(2)可知，在纽曼法中任一时刻的位移、速度、加速度都相互关联，这就使得运动方程的求解变成一系列相互关联的非线性方程的求解，这个求解过程必须通过迭代和求解联立方程组才能实现。
这就是通常所说的隐式求解法。
隐式求解法可能遇到两个问题。
一是迭代过程不一定收敛，二是联立方程组可能出现病态而无确定的解。
隐式求解法最大的优点是它具有无条件稳定性，即时间步长可以任意大。
如果采用中心差分法来进行动态问题的时域积分，则有如下位移、速度和加速度关系式：u(i+1)=2u(i)-u(i-1)+a(i)(△t)^2 (3)v(i+1)=[u(i+1)-u(i-1)]／2(△t) (4)式中u(i-1)，为i-1时刻的位移。
由式(3)可以看出，当前时刻的位移只与前一时刻的加速度和位移有关，这就意味着当前时刻的位移求解无需迭代过程。
另外，只要将运动过程中的质量矩阵和阻尼矩阵对角化，前一时刻的加速度求解无需解联立方程组，从而使问题大大简化，这就是所谓的显式求解法。
显式求解方法和隐式求解方法对比
采用有限元方法开展结构的动力学分析最终归结为求解离散后的常微分方程组tR KU U C U M =++ 。
在时域内求解该方程最常用的方法是直接积分法，而又根据求解过程中是否需要迭代求解线性方程组，将直接积分法分为隐式积分方法和显式积分方法两类。
隐式积分法认为ｔ＋Δｔ时刻系统的状态不仅与ｔ时刻状态有关，且与ｔ＋Δｔ时刻某些量有关。
因此隐式算法是根据t n 及t n-1...时刻体系的物理量值建立关于以t n+1时刻物理量为未知量的线性方程组，通过求解方程组确定t n+1时刻的物理量（常用的方法有线性加速度法、常平均加速度法、Newmark 方法、Wilson-θ法、Houbolt 方法等）。
而显式积分法认为ｔ＋Δｔ时刻系统的状态仅与ｔ时刻状态有关可，因此可由t n 及t n-1...时刻体系的物理量值直接外推t n+1时刻物理量值（如中心差分法），不需要求解线性方程组，实现了时间离散的解耦。
两种算法的比较 ：（1）隐式算法隐式算法基于虚功原理，要迭代计算。
隐式算法在每一增量步内都需要对静态平衡方程进行迭代求解，并且每次迭代都需要求解大型的线性方程组，这一过程需要占用相当数量的计算资源、磁盘空间和内存。
理论上在这个算法中的增量步可以很大，但是实际运算中上要受到接触以及摩擦等条件的限制。
随着单元数目的增加，计算时间几乎呈平方次增加。
由于需要矩阵求逆以及精确积分，对内存要求很高。
隐式算法的不利方面就是收敛问题不容易解决，且在开始起皱失稳时，在分叉点处刚度矩阵出现奇异。
（2）显式算法显示算法基于动力学方程，无需迭代，包括动态显式和静态显式算法。
动态显式算法采用动力学方程的中心差分格式，不用直接求解切线刚度，不需要进行平衡迭代，计算速度快，也不存在收敛控制问题。
该算法需要的内存也比隐式算法要少。
数值计算过程可以很容易地进行并行计算，程序编制也相对简单。
它也有一些不利方面。
显式算法要求质量矩阵为对角矩阵，而且只有在单元级计算尽可能少时速度优势才能发挥, 因而往往采用减缩积分方法，容易激发沙漏模式，影响应力和应变的计算精度。
显示积分和隐式积分法
显示积分和隐式积分法这是ansys里面的两种求解方法。
大多数非线性动力学问题一般多是采用显式求解方法，特别是在求解大型结构的瞬时高度非线性问题时，显示求解方法有明显的优越性。
下面先简要对比一下隐式求解法和显示求解法。
动态问题涉及到时间域的数值积分方法问题。
在80年代中期以前，人们基本上采用纽曼法进行时间域的积分。
根据纽曼法，位移、速度和加速度有着如下关系：u（i+1）=u（i）+△t*v（i）[（1—2p）a（i）+2p*a（i+1）]（1）v（i+1）=v（i）+△t[（1-2q）a（i）+2qa（i+1）]（2）一是迭代过程不一定收敛，二是联立方程组可能出现病态而无确定的解。
隐式求解法最大的优点是它具有无条件稳定性，即时间步长可以任意大。
如果采用中心差分法来进行动态问题的时域积分，则有如下位移、速度和加速度关系式：u（i+1）=2u（i）-u（i-1）+a（i）（△t）^2（3）v（i+1）=[u （i+1）-u（i-1）]／2（△t）（4）式中u（i-1），为i-1时刻的位移。
由式（3）可以看出，当前时刻的位移只与前一时刻的加速度和位移有关，这就意味着当前时刻的位移求解无需迭代过程。
另外，只要将运动过程中的质量矩阵和阻尼矩阵对角化，前一时刻的加速度求解无需解联立方程组，从而使问题大大简化，这就是所谓的显式求解法。
显式求解法的优点是它既没有收敛性问题，也不需要求解联立方程组，其缺点是时间步长受到数值积分稳定性的限制，不能超过系统的临界时间步长。
隐式求解法不考虑惯性效应[c]和[m]。
对于线性问题，无条件稳定，可以用大的时间步。
对于非线性问题，通过一系列线性逼近（newton-raphson）来求解；要求转置非线性刚度矩阵[k]，收敛时候需要小的时间步，对于高度非线性问题无法保证收敛。
因此，隐式求解一般用于线性分析和非线性结构静动力分析，包括结构固有频率和振型计算。
ansys使用的newmark时间积分法即为隐式求解法。
显式算法与隐式算法的区别
显式算法与隐式算法得区别1、显式算法最大优点就是有较好得稳定性。
动态显式算法采用动力学方程得一些差分格式(如广泛使用得中心差分法、线性加速度法、Newmark法与wilson法等),不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,时间步长只要取得足够小,一般不存在收敛性问题。
因此需要得内存也比隐式算法要少。
并且数值计算过程可以很容易地进行并行计算,程序编制也相对简单。
但显式算法要求质量矩阵为对角矩阵,而且只有在单元级计算尽可能少时速度优势才能发挥。
因而往往采用减缩积分方法,容易激发沙漏模式,影响应力与应变得计算精度。
静态显式法基于率形式得平衡方程组与Euler向前差分法,不需要迭代求解。
由于平衡方程式仅在率形式上得到满足,所以得出得结果会慢慢偏离正确值。
为了减少相关误差,必须每步使用很小得增量。
2、隐式算法隐式算法中,在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型得线性方程组,这以过程需要占用相当数量得计算资源、磁盘空间与内存。
该算法中得增量步可以比较大,至少可以比显式算法大得多,但就是实际运算中上要受到迭代次数及非线性程度得限制,需要取一个合理值。
3、求解时间t使用显式方法,计算成本消耗与单元数量成正比,并且大致与最小单元得尺寸成反比;应用隐式方法,经验表明对于许多问题得计算成本大致与自由度数目得平方成正比;因此如果网格就是相对均匀得,随着模型尺寸得增长,显式方法表明比隐式方法更加节省计算成本。
所谓显式与隐式,就是指求解方法得不同,即数学上得出发点不一样。
并不就是说显式只能求动力学问题,隐式只能求静力学问题,只就是求解策略不通。
显式求解就是对时间进行差分,不存在迭代与收敛问题,最小时间步取决于最小单元得尺寸。
过多与过小得时间步往往导致求解时间非常漫长,但总能给出一个计算结果。
解题费用非常昂贵。
因此在建模划分网格时要非常注意。
隐式求解与时间无关,采用得就是牛顿迭代法(线性问题就直接求解线性代数方程组),因此存在一个迭代收敛问题,不收敛就得不到结果。
期权定价的有限差分法
第9章期权定价的有限差分法【本章精粹】期权定价最终归结为一个二阶偏微分方程，而有限差分方法是计算偏微分方程的有效工具。
因此，本章介绍3种常见的有限差分方法，并且给出了相关程序，让读者了解有限差分方法计算的基本原理，熟悉用显式法和隐式法计算欧式看涨期权价格。
9.1 有限差分计算方法的基本原理偏微分方程在金融工程中占有重要位置，著名的Black-Scholes 方程就是以二阶偏微分方程形式给出的。
偏微分方程为求解复杂的金融衍生工具价格提供了有力手段，但是偏微分方程通常没有解析解，因此用数值计算方法求解衍生工具价格就成为金融工程的一项基本功。
求解金融衍生工具价格与求解通常偏微分方程的区别主要在于一般偏微分方程是给定初值求解终值，而衍生品定价问题是给定终值求初值，属于倒向随机偏微分方程求解。
有限差分方法的核心思想是对导数进行离散化，把偏微分方程转化为差分方程，然后利用迭代法求解。
根据对偏导数离散方法的不同，有限差分方法可分为显式差分法、隐式差分法和内含差分法，下面分别进行介绍。
假设,i j f 表示在i 时刻股票价格为第j 价位的期权价格，对f 一阶导数进行如下差分：,1,i j i j f f fS S+-∂=∂∆ (9-1) 1,,i j i j f f ft t+-∂=∂∆ (9-2) 式(9-1)、式(9-2)的差分方法称为显式差分法。
也可以对一阶导数做如下差分：,,1i j i j f f fS S--∂=∂∆ (9-3) ,1,i j i j f f ft t --∂=∂∆ (9-4) 式(9-3)、式(9-4)的差分方法称为隐式差分法。
也可以做如下差分：,1,12i j i j f f fS S+--∂=∂∆ (9-5) 1,1,2i j i j f f ft t +--∂=∂∆ (9-6) 式(9-5)、式(9-6)的差分方法称为内含差分法。
对二阶微分方程，用如下方法进行差分：2,1,,,12/i j i j i j i j f f f f f S S S S +---⎛⎫∂=-∆ ⎪∂∆∆⎝⎭ (9-7) 整理得：2,1,1,222i j i j i j f f f fS S +-+-∂=∂∆以上就是偏微分方程的常见离散方法。
一维热传导方程的差分法
一维热传导方程的差分法一维热传导方程是描述材料内部温度分布随时间变化的数学模型。
它在许多实际工程问题中起着重要的作用，比如热传导、材料加工、建筑设计等。
差分法是一种用于数值求解偏微分方程的常用方法，其原理是将偏微分方程中的导数项用差分近似代替，然后将求解区域划分为离散点，最终得到一个代数方程组。
本文将介绍一维热传导方程的差分法求解过程。
一维热传导方程可以写成如下形式：\[\frac{\partial u}{\partial t} = \alpha \frac{\partial^2 u}{\partial x^2}\]\(u(x, t)\)表示材料内部温度分布，\(x\)是空间坐标，\(t\)是时间，\(\alpha\)是热扩散系数。
为了使用差分法求解该方程，我们需要对空间和时间进行离散化。
假设求解区域为\(0 \leq x \leq L\)，时间区间为\(0 \leq t \leq T\)，将空间和时间分别划分成\(N_x\)和\(N_t\)个小区间，步长分别为\(\Delta x = \frac{L}{N_x}\)和\(\Delta t = \frac{T}{N_t}\)。
接下来，我们将使用显式差分格式对一维热传导方程进行离散化。
我们定义离散点\(u_i^n = u(i\Delta x, n\Delta t)\)，用\(u_i^n\)表示时间\(n\)、空间\(i\)处的温度。
那么热传导方程可以用差分格式表示为：\[\frac{u_i^{n+1} - u_i^n}{\Delta t} = \alpha \frac{u_{i+1}^n - 2u_i^n +u_{i-1}^n}{\Delta x^2}\]为了进行数值求解，我们需要给定初始条件和边界条件。
初始条件可以表示为：\[u_i^0 = f(i\Delta x)\]边界条件可以是温度固定或热传导定律，比如：\[u_0^n = g_1(t), u_{N_x}^n = g_2(t)\]或者\[\frac{\partial u}{\partial x}(0, t) = 0, \frac{\partial u}{\partial x}(L, t) = 0\]接下来，我们可以通过迭代计算离散点的温度值来求解一维热传导方程。