第六讲_2.3高阶有限差分(6)
1 T j 1 T j (k1 2k2 2k3 k4 ) 6 k1 dtf (t j , T j ) k dt ,Tj 1 ) 2 2 k2 dt k3 dtf (t j , T j ) 2 2 k4 dtf (t j dt , T j k3 ) k2 dtf (t j
y f ( x, y ) x
c=1, 得到Euler公式
h yn 1 yn f ( xn h), yn hf ( xn , yn ) 2
1 T j 1 T j (k1 2k2 2k3 k4 ) 6 k1 dtf (t j , T j ) k dt ,Tj 1 ) 2 2 k dt k3 dtf (t j , T j 2 ) 2 2 k4 dtf (t j dt , T j k3 ) k2 dtf (t j
回顾：
显式差分方法
隐式差分方法
n1 n1 n1 n n n sTi sTi 1 (2 2s)T i 1 sT i 1 (2 2s)T i sT i 1
Crank-Nicolson隐式差分
显式
最简单隐式
Crank-Nicolson 隐式
牛顿冷却问题不同差分格式的matlab程序
1 1 , w2 2dx 2dx
插值权重：
w1 0.5, w2 0.5
w1
x-2dx x-dx 1 a 2 b
x
x+dx 3 c 4 d
x+2dx
四点三阶
x-2dx x-dx 1 a 2 b
x
x+dx 3 c 4 d
x+2dx
四点插值
f ( x)
1 2 2 1 f ( x 2dx) f ( x dx) f ( x dx) f ( x 2dx) 6 3 3 6
%Malab-1D clear; clc; figure('color','w'); t0=1; % initial temperature tau=0.7; % time constant dt=0.01; % time interval t_total=10; nt=round(t_total/dt); % total time steps T_e(1)=t0; T_i(1)=t0; T_m(1)=t0; for i=1:nt; xt(i)=(i-1)*dt; T_e(i+1)=T_e(i)*(1-dt/tau); %explicit T_i(i+1)=T_i(i)/(1+dt/tau); % implicit T_m(i+1)=T_m(i)*(1-dt/2/tau)/(1+dt/2/tau); %mix T_a(i)=t0*exp(-xt(i)/tau); %analytical results end
小结
• 1、高阶精度的有限差分算子可以利用Taylor级数展 开推导出来。 • 2、对于双侧算子，权重系数随着算子长度的增加 迅速减小。 • 3、对于单侧算子，权重系数随着算子长度的增加 迅速增大。 • 4、以上给出的高阶算子提高了空间微分近似的精 度，我们还需要给出提高时间域外推的精度的方法。
算子长度
微分
微分 函数f(x)在某一位置x的近似值以及导数（一阶、二阶、高阶等）如何计算？ 办法：充分利用更多相邻点的信息。 问题：如何计算相邻点的权重（系数）？
i-3
i-2
i-1
i
i+1
i+2
i+3
Taylor算子
构造如下形式的表达式：
在不同相邻点进行Taylor级数展开
f (i ) ( x) w(ji ) f ( x j )
显式结构
隐式结构
几类显式Runge-Kutta方法
对于L=2，则
其局部截断误差是
令
1
四阶R-K方法是最常见的高精度方法：
1 T j 1 T j (k1 2k2 2k3 k4 ) 6 k1 dtf (t j , T j ) k dt ,Tj 1 ) 2 2 k2 dt k3 dtf (t j , T j ) 2 2 k4 dtf (t j dt , T j k3 ) k2 dtf (t j
为什么需要高阶算子
以上的差分形式又称为微分算子(derivative operator)。在实际应用中，基于一阶 和二阶差分形式的差分格式的精度不能满足要求，在以上差分形式中，我们只 用到相邻网格节点上的值。 是否可以利用更多的信息提高微分算子的精度？ 用待定系数法构造精度较高的差分格式。
Taylor算子
牛顿冷却问题不同方法结果的对比
四阶R-K方法精度最高！
牛顿冷却问题不同方法结果的对比
四阶R-K方法精度最高！
j 1,l
f(x-dx)
f(x)
f(x+dx)
一阶情况：
f ( x dx) f x f ( x)dx
f ( x dx) f x f ( x)dx
f(x-dx)
f(x)
f(x+dx)
1 两点一阶情况：
2
af af af dx
bf bf bf dx
董良国等，2000，地球物理学报
实例： 一阶弹性波方程交错网格高阶差分解法
高阶Taylor外推
高阶外推
for n=0,1,2,3,…,N-1 xn+1=xn+dx; k1=dx*f(xn,yn); k2=dx*f(xn+1,yn+k1); yn+1=yn+1/2*(k1+k2) end
四点一阶导数
f ( x)
1 ( f ( x 2dx) 8 f ( x dx) 8 f ( x dx) f ( x 2dx)) 12dx
四点三阶导数
f ( x)
1 ( f ( x 2dx) 2 f ( x dx) 2 f ( x dx) f ( x 2dx)) 3 2(dx)
x-2dx x-dx 1 a 2 b
x 3 c
x+dx 4 d
x+2dx 5 e
五点二阶导数
f ( x)
1 ( f ( x 2dx) 16 f ( x dx) 30 f ( x) 16 f ( x dx) f ( x 2dx)) 2 12(dx)
单侧算子
dT f (T , t ) dt
四阶R-K有限差分解牛顿冷却问题
f (T , t )
T
T_rk(1)=t0; for i=1:nt; xt(i)=(i-1)*dt; k1=-dt*T_rk(i)/tau; k2=-dt*(T_rk(i)+k1/2)/tau; k3=-dt*(T_rk(i)+k2/2)/tau; k4=-dt*(T_rk(i)+k3)/tau; T_rk(i+1)=T_rk(i)+1/6*(k1+2*k2+2*k3+k4); end plot(xt,T_rk,'k.-');hold on % R_K_4 difference scheme
两式相加：
af bf (a b) f (a b) f dx
注意：上式对任意a，b都成立！
插值：
令
a b 1 a b 0
a 0.5 微分： a b 0 令 a b 1/ dx b 0.5
插值权重：
1 1 a ,b ຫໍສະໝຸດ 2dx 2dx预估校正
for n=0,1,2,3,...,n-1 tn 1 tn dt ; k1 k2 dt
预估
dt
Tn (Tn k 1);
校正
1 Tn 1 Tn (k 1 k 2); 2 End
高阶外推
y f ( x, y ) x
显式Euler方法是最简单的单步法，它是一阶的，它可以看作Talylor展 开后取前两项。因此，得到高阶方法的一个直接想法是用Talylor展开， 如果能计算的高阶导数，则可写出p阶方法的计算方法
不同差分格式计算结果对比
混合差分格式精度最高！
不同差分格式计算结果对比
混合差分格式精度最高！
不同差分格式计算结果对比
混合差分格式精度最高！
不同差分格式计算结果对比
混合差分格式精度最高！
不同差分格式计算结果对比
混合差分格式精度最高！
高阶差分算子
主讲人：胡才博 中国科学院大学地球科学学院 中国科学院计算地球动力学重点实验室
xt(nt+1)=nt*dt; plot(xt,T_e,'b.-', xt,T_i,'g.-', xt,T_m,'m.-', xt,T_a,'r.-',);hold on % set(gca,'DataAspectRatio',[(max(xt)min(xt))/(max(T_e)-min(T_e))/3 1 1]); xlabel('Time (s)','Fontname','times new roman','FontSize',14); ylabel('Temperature','Fontname','times new roman','FontSize',14); title('dt=0.01 tau=0.7');
合集下载
高阶差分方程
2013-7-25 经济管理学院财务与投资系 刘亚娟 7
• 第一种情况（不同的实根）：当a12＞4a2时， b1和b2为不同的实根。在这种情况下，b1t 和b2t线性无关，余函数可以简单地写成b1t 和b2t的线性组合，即yc=A1b1t+A2b2t。
2013-7-25
经济管理学院财务与投资系 刘亚娟
2013-7-25 经济管理学院财务与投资系 刘亚娟 2
• 因此，yt的二阶差分可以转换为包含两期时滞的 项的和。因为像Δ2yt和Δyt这样的表达式写起来很 麻烦，所以我们将二阶差分方程重新定义为包含 变量的两期时滞的方程。类似地，三阶差分方程 为包含三期时滞的方程；等等。 • 我们首先集中讨论二阶差分方程的解法，然后再 在后面的章节中将其推广至高阶差分方程。为控 制讨论的范围，在本章，我们仅讨论常系数线性 差分方程。但对常数项和可变项两种形式，均作 考察。
• 所以，yc A1 (1) A2 (5) A1 A2 (5) • 因此，方程的通解为
t t t
y A A (5) 3t
c 1 2
2013-7-25 经济管理学院财务与投资系 刘亚娟 13
t
• (3)该方程的特别积分为： y
p
8 2 2 4t 2t
• 该方程的特征方程为：b2-2b+1=0，所以特 征根为： b 2 4 4 1 2 1 b
1 1 2 1 2
1
1
• 联立方程求解A1＝5和A2=3，最后把它代入 通解中可得特解：
y 5 (2) 3(8) 2
t
2013-7-25 经济管理学院财务与投资系 刘亚娟 12
t
t
• （2）该方程的特别积分为：
• 第一种情况（不同的实根）：当a12＞4a2时， b1和b2为不同的实根。在这种情况下，b1t 和b2t线性无关，余函数可以简单地写成b1t 和b2t的线性组合，即yc=A1b1t+A2b2t。
2013-7-25
经济管理学院财务与投资系 刘亚娟
2013-7-25 经济管理学院财务与投资系 刘亚娟 2
• 因此，yt的二阶差分可以转换为包含两期时滞的 项的和。因为像Δ2yt和Δyt这样的表达式写起来很 麻烦，所以我们将二阶差分方程重新定义为包含 变量的两期时滞的方程。类似地，三阶差分方程 为包含三期时滞的方程；等等。 • 我们首先集中讨论二阶差分方程的解法，然后再 在后面的章节中将其推广至高阶差分方程。为控 制讨论的范围，在本章，我们仅讨论常系数线性 差分方程。但对常数项和可变项两种形式，均作 考察。
• 所以，yc A1 (1) A2 (5) A1 A2 (5) • 因此，方程的通解为
t t t
y A A (5) 3t
c 1 2
2013-7-25 经济管理学院财务与投资系 刘亚娟 13
t
• (3)该方程的特别积分为： y
p
8 2 2 4t 2t
• 该方程的特征方程为：b2-2b+1=0，所以特 征根为： b 2 4 4 1 2 1 b
1 1 2 1 2
1
1
• 联立方程求解A1＝5和A2=3，最后把它代入 通解中可得特解：
y 5 (2) 3(8) 2
t
2013-7-25 经济管理学院财务与投资系 刘亚娟 12
t
t
• （2）该方程的特别积分为：
数值方法课件_有限差分法
土
工
数
值
计
算
方
法
法分差限有 讲二第
2.5 基本原理
1)问题的提出 1)问题的提出 2)稳定性和收敛性 2)稳定性和收敛性 3)差分法的求解步骤 3)差分法的求解步骤
土
工
数
值
计
算
方
法
法分差限有 讲二第
2.5 基本原理
3)差分法的求解步骤 3)差分法的求解步骤
①对求解区域进 ①对求解区域进 行网格划分； 行网格划分； ②选择逼近微分 ②选择逼近微分 方程定解问题 方程定解问题 的差分格式； 的差分格式；
¶ 2V c ¶V 2 c -1 ¶V = ( ) + V ¶T ¶ 2 z V ¶z v
1.7
土
工
数
值
计
算
方
法
2.6 例子
法分差限有 讲二第
软粘土地基非线性一维固结分析 软粘土地基非线性一维固结分析
定解条件变为： 定解条件变为： V = 1 ; （1） T v = 0 : （2） Z = 0 : V = b ; （3） Z = 1 :
算
方
法
法分差限有 讲二第
2.1 优点与局限性
不 适 应
规则边界的问题 规则边界的问题
应 适
简便、易编程 简便、易编程 不规则边界的问题 不规则边界的问题
土
工
数
值
计
算
方
法
法分差限有 讲二第
2.2 基本思路
1
3
2
将求解区域 将求解区域 划分成网格 划分成网格
差分方程解 差分方程解 作为微分方 作为微分方 程近似解。 程近似解。
有限差分法PPT课件
有限差分法在求解导热微分方程中的应用
1
有限差分方法是一种微分方法,广泛用于计算机求解偏微分方程 。
为求解由偏微分方程定解问题所构造的数学模型，有限差分法 是将定解区域（场区）离散化为网格离散节点的集合。并以各离 散点上函数的差商来近似该点的偏导数，使待求的偏微分方程定 解问题转化为一组相应的差分方程。根据差分方程组解出各离散 点处的待求函数值——离散解。
Q c hc (T Ta )
Qr (T4Ta4)
代 入
C pz T t kz 2 T 2 h c T 2T 4 2 h c T a 2T a 4
上 式Leabharlann 边界条件： x=0m ,x=1m, y=1m ； q=0 w/m2
y=1m
； T=300 K
12
（2）利用matlab中的pdetool工具箱，首先绘出空间区域，并以0.1m为 步长对其进行网格划分。 （3）输入已知的参数并设定边界条件
2
建立控制方程及定解条件
确定节点（区域离散化）
建立节点物理量的代数方程
设立迭代初值
求解代数方程组 否
收敛？ 是
解的分析
改进初场
3
1. 建立控制方程及定解条件
根据实际问题建立偏微分方程，同时给出边界条件。
2. 区域离散化
理论上可以通过任意的网格划分把求解区域划分成许多求解区域，以网格 线的交点作为需要确定的物理量的空间位置。实际应用中根据边界的形状采用 最简单、最有规律，和边界拟合程度最佳的方法来分割。
建立节点物理量的离散方程节点类型内节点边界节点泰勒级数展开法热平衡法泰勒级数展开法热平衡法热平衡法多运用于非均分网格划分下离散方程的建立其物理概念清晰推导过程简洁我们以二维稳态无内热源矩形均分下的温度场为例先用泰勒级数展开法对内节点由ab两个式子即可推出一阶导数和二阶导数的差分一般取中心差分更为精确一阶导数的中心差分
1
有限差分方法是一种微分方法,广泛用于计算机求解偏微分方程 。
为求解由偏微分方程定解问题所构造的数学模型，有限差分法 是将定解区域（场区）离散化为网格离散节点的集合。并以各离 散点上函数的差商来近似该点的偏导数，使待求的偏微分方程定 解问题转化为一组相应的差分方程。根据差分方程组解出各离散 点处的待求函数值——离散解。
Q c hc (T Ta )
Qr (T4Ta4)
代 入
C pz T t kz 2 T 2 h c T 2T 4 2 h c T a 2T a 4
上 式Leabharlann 边界条件： x=0m ,x=1m, y=1m ； q=0 w/m2
y=1m
； T=300 K
12
（2）利用matlab中的pdetool工具箱，首先绘出空间区域，并以0.1m为 步长对其进行网格划分。 （3）输入已知的参数并设定边界条件
2
建立控制方程及定解条件
确定节点（区域离散化）
建立节点物理量的代数方程
设立迭代初值
求解代数方程组 否
收敛？ 是
解的分析
改进初场
3
1. 建立控制方程及定解条件
根据实际问题建立偏微分方程，同时给出边界条件。
2. 区域离散化
理论上可以通过任意的网格划分把求解区域划分成许多求解区域，以网格 线的交点作为需要确定的物理量的空间位置。实际应用中根据边界的形状采用 最简单、最有规律，和边界拟合程度最佳的方法来分割。
建立节点物理量的离散方程节点类型内节点边界节点泰勒级数展开法热平衡法泰勒级数展开法热平衡法热平衡法多运用于非均分网格划分下离散方程的建立其物理概念清晰推导过程简洁我们以二维稳态无内热源矩形均分下的温度场为例先用泰勒级数展开法对内节点由ab两个式子即可推出一阶导数和二阶导数的差分一般取中心差分更为精确一阶导数的中心差分
《高阶差分方程式》课件
04
高阶差分方程式的应用实例
金融领域的应用实例
股票价格预测
高阶差分方程式可以用于描述股 票价格的动态变化，通过历史数 据来预测未来的股票价格走势。
风险评估
在金融领域，高阶差分方程式可 以用于评估投资组合的风险，通 过分析资产价格的变动规律来预 测未来的市场波动。
期货价格建模
在期货市场中，高阶差分方程式 可以用于建立期货价格模型，以 预测未来期货价格的变化趋势。
数值分析法求解高阶差分方程式
数值分析法是一种基于数值计算的方法，通过将高阶差分方程式转化为数值求解问题，利用数值计算的方法得到近似解。常 用的数值分析法包括欧拉法、龙格-库塔法等。
数值分析法的优点是适用范围广，可以求解各种类型的高阶差分方程式。然而，数值分析法的精度和稳定性取决于所选择的 数值方法和计算步长，需要进行合理的选择和控制。
解释
高阶差分方程式在数学、物理、工程等领域有广泛应用，用于描述各种动态系统的行为。
高阶差分方程式的形式
线性形式
如果差分方程中不含有 (y) 的非线性项，则称为 线性差分方程。例如：(y_{n+1} - 2y_n + y_{n1} = 0)。
非线性形式
如果差分方程中含有 (y) 的非线性项，则称为非 线性差分方程。例如：(y_{n+1}^2 - y_n^2 = 0) 。
物理学中的应用
高阶差分方程式在物理学中的波动方程、热传导方程等领域有广泛应用，为解决实际问 题提供了有效工具。
工程学中的应用
在工程学中，高阶差分方程式被用于描述信号处理、图像处理等领域的问题，推动了相 关领域的科技进步。
高阶差分方程式未来的研究方向
01
高效算法研究
有限差分法
1.6.1 差分格式
上页
下页
返回
结束
如图，以二维场为例，在一由边界区域L界定的二维区域 界定的二维区域D 如图 ， 以二维场为例 ， 在一由边界区域 界定的二维区域 电位函数φ满足拉普拉斯方程且给定为第一类边界条件： 内，电位函数φ满足拉普拉斯方程且给定为第一类边界条件
∂ 2ϕ ∂ 2ϕ ∇ 2ϕ = + 2 ∂x ∂y 2
启动 赋边界节点已知电位值 赋予场域内各节点电位初始值 累计迭代次数 N=0 N=N+1 按超松弛法进行一 次迭代，求 ϕi(,N+1) j
上页
下页
返回
结束
Y N
所有内点 相邻二次迭代值的最大误差 是否小于 W 打印 N，ϕ(i, j ) 迭代解程序框图 停机
上页
下页
返回
结束
例 一长直接地金属槽截面如图。其侧壁与底面的电位均为零， 而顶盖电位ϕ4＝100。求槽内电位分布。 解： •二维场第一类边值问题。 •将二维场域划分成正方形 网格，步距h＝a/4。 •场域内任一点电位ϕ应满 足二维拉普拉斯方程的差 分计算格式。
α
4
k+ k) [ϕ i(− 1,1j ) + ϕ i(,kj+ 1 ) + ϕ i(+ 1, j + ϕ i(,kj )+ 1 − 4ϕ i(,kj ) ] −1
迭代收敛的速度与 α 有明显关系： 有明显关系： 迭代收敛的速度与 收敛的速度
收敛因子（α ） 迭代次数（ N） 1.0 >1000 1.7 269 1.8 1.83 1.85 1.87 1.9 174 143 122 133 171 2.0 发散
上页
下页
有限差分法基本原理PPT课件
uin1
uin
a
t x
(uin
un i 1
)
ui0 u (xi )
几种差分格式介绍
u a u 0 t x u(x,0) u(x)
FTFS格式（时间向前差分、空间向前差分）
uin1 uin uin1 uin 0
t
x
ui0 u (xi )
uin 1
uin
a
t x
(uin1
uin )
ui0 u (xi )
几种差分格式介绍
FTBS格式（时间向前差分、空间向后差分）
限差分方程的解是收敛T的(i。, n)
lim
x0,t
0
Ti
t
一般情况下，证明收敛性是非常难的，暂不予以证明。
3.稳定性 稳定性讨论的是差分解的误差在计算过程中的发展问题。
在 数值解中，引进误差是不可避免的，电子计算机也有舍入误差， 因此实际算得的有限差分方程的解是近似解。这种误差是要向其 他方向传播的，如果计算中引入的误差在以后逐层计算过程中影 响逐渐消失或者保持有界，则称差分方程是稳定的。否则就是不 稳定的。
Von Neumann稳定性分析方法简介
分析例题
T n1 i
Ti n
t x 2
(Ti
n 1
2Ti n
Ti
n 1
),
S
t x 2
Ti n1
STi n1
(1
2S )Tin
STi
n 1
上式T中i n 近似数值
有限差分法
有限差分法有限差分法finite difference method微分方程和积分微分方程数值解的方法。
基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。
然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。
有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。
此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。
对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。
另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。
此外，还有一个重要的概念必须考虑，即差分格式的稳定性。
因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。
前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。
只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。
关于差分格式的构造一般有以下3种方法。
最常用的方法是数值微分法，比如用差商代替微商等。
另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。
此外还可以用待定系数法构造一些精度较高的差分格式。
有限差分方法基础ppt课件
t
x
0
(x,0) (x)
这里 (x) 为某已知函数。同样，差分方程也必须有初始条件：
（2-7）
n1 i
n i
n i 1
n i 1
0
t
2x
0 i
(xi )
（2-8）
初始条件是一种定解条件。如果是初边值问题，定解条件中还应有适当的边界条件。差分方程和其定解条件一起， 称为相应微分方程定解问题的差分格式。
图1-3 均匀和非均匀网格实例2
22
第二节 差分方程、截断误差和相容性/差分方程（1/3）
差分相应于微分，差商相应于导数。差分和差商是用有限形式表 示的，而微分和导数则是以极限形式表示的。如果将微分方程中 的导数用相应的差商近似代替，就可得到有限形式的差分方程。 现以对流方程为例，列出对应的差分方程。
FTCS格式的截断误差为
Rin O(t, (x)2 )
FTFS和FTBS格式的截断误差为
Rin O(t, x)
3种格式对 t 都有一阶精度。
（2-12） （2-13）
30
第二节 差分方程、截断误差和相容性/相容性（1/3）
25
第二节 差分方程、截断误差和相容性/截断误差（1/6）
按照前面关于逼近误差的分析知道，用时间向前差商代替时间导数时的误差为 O(t) ,
用空间中心差商代替空间导数时的误差为 O((x)2 ) ，因而对流方程与对应的差分方程之间也存在一个误差，它是
Rin O(t) O((x)2 ) O(t, (x)2 )
表2
《有限差分法初步》课件
4. 三维有限差分法
1
三维有限差分法流程Fra bibliotek解释三维有限差分法的计算流程和步骤。
2
三维热传导方程的有限差分解法
探索如何利用有限差分法解决三维热传导方程。
3
显式法与隐式法
分析三维有限差分法中显式法和隐式法的性能和适用范围。
5. 有限差分法的误差和稳定性分析
1 截断误差
介绍有限差分法中的截断误差及其对计算结 果的影响。
显式法与隐式法
比较一维有限差分法中的显式法和隐式法的优缺点。
一维热传导方程的有限差分解法
探讨如何用有限差分法解决一维热传导方程。
3. 二维有限差分法
二维有限差分法流程
详细介绍二维有限差分法的计算 流程。
二维热传导方程的有限差 分解法
深入研究如何利用有限差分法解 决二维热传导方程。
显式法与隐式法
比较二维有限差分法中的显式法 和隐式法的应用场景。
2 稳定性分析
讨论有限差分法的稳定性问题及其对数值计 算的重要性。
6. 总结
有限差分法应用的优势和不足
总结有限差分法在实际应用中的优点和局限性。
发展前景
展望有限差分法在未来的应用前景和发展方向。
7. 参考文献
1. 引用的相关文献1 2. 引用的相关文献2
有限差分法初步
这是一份关于有限差分法的PPT课件，从引言到有限差分法的具体应用进行 了全面而深入的介绍，希望能够帮助大家更好地理解和应用这一方法。
1. 引言
有限差分法是一种计算数学方法，广泛应用于解决偏微分方程等问题，本节 介绍有限差分法的概述和使用场景。
2. 一维有限差分法
一维有限差分法流程
介绍一维有限差分法的整体流程和步骤。
有限差分方法
有限差分方法有限差分方法一种求偏微分(或常微分)方程和方程组定解问题的数值解的方法，简称差分方法。
微分方程的定解问题就是在满足某些定解条件下求微分方程的解。
在空间区域的边界上要满足的定解条件称为边值条件。
如果问题与时间有关，在初始时刻所要满足的定解条件，称为初值条件。
不含时间而只带边值条件的定解问题，称为边值问题。
与时间有关而只带初值条件的定解问题，称为初值问题。
同时带有两种定解条件的问题，称为初值边值混合问题。
定解问题往往不具有解析解，或者其解析解不易计算。
所以要采用可行的数值解法。
有限差分方法就是一种数值解法，它的基本思想是先把问题的定义域进行网格剖分，然后在网格点上，按适当的数值微分公式把定解问题中的微商换成差商，从而把原问题离散化为差分格式，进而求出数值解。
此外，还要研究差分格式的解的存在性和唯一性、解的求法、解法的数值稳定性、差分格式的解与原定解问题的真解的误差估计、差分格式的解当网格大小趋于零时是否趋于真解(即收敛性)，等等。
有限差分方法具有简单、灵活以及通用性强等特点，容易在计算机上实现。
偏微分方程初值问题的差分法许多物理现象随着时间而发生变化、如热传导过程、气体扩散过程和波的传播过程都与时间有关。
描述这些过程的偏微分方程具有这样的性质：若初始时刻t=t0的解已给定,则t t0时刻的解完全取决于初始条件和某些边界条件。
利用差分法解这类问题，就是从初始值出发，通过差分格式沿时间增加的方向，逐步求出微分方程的近似解。
双曲型方程的差分方法最简单的双曲型方程的初值问题是：式中嫓(x)为已知初值函数。
这初值问题的解是：u(x,t)=嫓(x-at)。
(2)由(2)可见，(1a)(1b)的解(2)当a>0时代表一个以有限的速度a沿特征线x-at=常数向右传播的波,而解u(x,t)在P(慜,惭)点的值完全由嫓(x)在x轴上的点A(慜-а惭,0)的值决定。
A点就是双曲型方程(1a)在P点的依赖域(图1)。