当前位置：文档之家 › 差分格式稳定性及数值效应

差分格式稳定性及数值效应

  • 格式：docx
  • 大小：390.82 KB
  • 文档页数：5

下载文档原格式

  / 5

合集下载

研究报告有限差分格式稳定性的其他方法-报告

研究报告有限差分格式稳定性的其他方法-报告

研究有限差分格式稳定性的其他方法摘要偏微分方程的求解一直是大家比较关心的一个问题，而有限差分格式则是求解偏微分方程时常用并且有效的一个方法。

因此，研究有限差分格式的性质就显得尤为重要。

在课上我们已经跟着老师学习了运用Fourier方法研究有限差分格式的稳定性，但是在很多研究有限差分格式稳定性的问题中仅仅会用Fourier方法是不够的，所以在本篇论文中，将会介绍其他三种常用的研究有限差分格式稳定性的方法，分别是：Hirt启示型方法、直接方法（或称矩阵方法）和能量不等式方法。

关键字：偏微分方程；有限差分格式；稳定性AbstractThe solution of partial differential equations has been more concerned with a problem, and the finite difference scheme is a mon and effective method for solving partial differential equations. Therefore, it is very important to study the character of the finite difference scheme. We have followed the teacher to learn the use of Fourier method of finite difference scheme stability, but in a lot of research on the stability of finite difference scheme is only used Fourier method is not enough, so in this paper, will introduce the other three kinds of monly used in the study of finite difference scheme stability method, respectively is: Hirt enlightenment method, direct method (or matrix method) and energy inequality method.Key words: partial differential equation; finite difference scheme; stability1 前言微分方程的定解问题就是在满足某些定解条件下求微分方程的解。

无条件稳定 条件稳定 有限差分法

无条件稳定 条件稳定 有限差分法

无条件稳定条件稳定有限差分法
无条件稳定条件稳定有限差分法是一种常用的数值计算方法，广泛应用于科学和工程领域中的物理模拟、流体力学、电磁学等问题的求解。

该方法是基于有限差分逼近微分方程的思想，将求解区域离散化为有限个网格点，然后利用差分公式来逼近微分方程。

其中，稳定性是该方法最重要的性质之一。

无条件稳定的有限差分法是指在任何时刻、任何网格密度下都能保持数值解的稳定性。

这种方法的优点是计算方便，不需要对时间步长进行严格限制，但缺点是精度较低。

常见的无条件稳定的有限差分法有前向差分法和后向差分法。

条件稳定的有限差分法是指在一定条件下才能保持数值解的稳定性。

这种方法通常需要严格限制时间步长和网格密度，以保证数值解的精度和稳定性。

常见的条件稳定的有限差分法有CFL条件、稳定性限制条件等。

总之，无论是无条件稳定还是条件稳定的有限差分法，都在不同程度上应用广泛，是解决科学和工程问题的重要数值方法之一。

- 1 -。

数值天气预报涡度方程差分格式的守恒性和稳定性

数值天气预报涡度方程差分格式的守恒性和稳定性

数值天气预报涡度方程差分格式的守恒性和稳定性的报告，
600字
数值天气预报涡度方程差分格式是一种用来解决时间步长变化与风场变化之间的相互作用的数值解法，它的守恒性和稳定性是保证天气预报获得精确结果的关键。

本文将通过分析不同层次的涡度方程差分格式的守恒性和稳定性来讨论它在数值天气预报中的作用和意义。

首先，对于单层次的涡度方程差分格式，它的守恒性表明所有物理量在边界上完全满足守恒性定理，因此涡度方程差分格式能够保证模拟过程中物理量守恒。

其次，涡度方程差分格式的稳定性表明，所有物理量在空间和时间上的变化都是有限的，保证了模拟的精度不会随着时间的变化而变化，从而保证了模拟的准确性。

此外，多层次的涡度方程差分格式同样具有良好的守恒性和稳定性，即便不同层次的物理量之间存在交互作用，但守恒性和稳定性仍然得到保证。

此外，由于多层次的涡度方程差分格式本身具有较低的计算复杂度，因此可以实现更加精细的模拟。

综上所述，涡度方程差分格式是一种具有良好守恒性和稳定性的数值解法，它能够有效地模拟大气中的涡度，用于数值天气预报时能够发挥重要作用。

因此，在利用涡度方程差分格式进行数值天气预报时，必须首先评估它的守恒性和稳定性，以保证模拟结果的准确性和可靠性。

非定常流体力学中间差分格式稳定性分析研究

非定常流体力学中间差分格式稳定性分析研究

非定常流体力学中间差分格式稳定性分析研究随着计算机技术的发展，数值模拟已经成为研究非定常流体力学的重要手段。

其中差分法是最常用的一种计算方法。

而中心差分法是差分法中最为常用的方法之一。

在数值计算中，稳定性是非常重要的一个问题。

本文将从非定常流体力学的角度出发，分析中心差分格式的稳定性问题。

一、中心差分法中心差分法是一种最为常用的差分法，其具体计算过程是将计算点的函数值表示为它自身与周围计算点值的线性组合，其中，每个计算点的函数值均采用相同的线性组合模式。

这个模式就是中心差分法的核心。

中心差分法可以用于求解一些常见的偏微分方程，例如泊松方程、热传导方程、对流扩散方程，以及非定常流体力学中的纳维－斯托克斯方程等。

二、非定常流体力学的求解非定常流体力学是流体运动学和动力学的研究，其中：研究的是在时间和空间上变化的流场。

在非定常流体力学中，求解纳维－斯托克斯方程是相当难的。

要解决这一问题，可以采用数值模拟的方法。

由于非定常流体力学的求解过程涉及到高维空间和复杂的数学模型，因此需要具有高性能的计算机和优秀的数值方法。

中心差分法作为一种常见的数值方法，可以用于求解非定常流体力学。

不过，如果不考虑其稳定性问题，这种方法也是会出现一些问题的。

三、中心差分格式的稳定问题在数值计算中，稳定性问题是非常重要的一个问题。

稳定性是指对精度的要求。

一种数值计算方法，如果该方法对初始误差非常敏感，或者计算过程中误差放大得太快，那么这种方法就是不稳定的。

因此，中心差分格式的稳定性问题需要引起我们的关注。

中心差分格式的稳定性取决于流场的不稳定性，并且与形式构成的方程相关。

由于中心差分格式本身是一种稳定的方法，但它的稳定性却取决于数值格式和解的一些特性，如模型方程、网格尺寸等因素。

为了解决中心差分格式的稳定性问题，我们可以采用标量稳定性分析和矩阵稳定性分析两种方法。

通过这两种方法的研究和分析，我们可以更好地了解中心差分格式的稳定性问题，并实现更为精准的求解。

对差分数学模型进行评价和分析

对差分数学模型进行评价和分析

对差分数学模型进行评价和分析
差分数学模型是一种常用的数理模型，用于描述某一自变量在不同时间或空间下的变化情况。

其优点在于简单易懂、计算机处理方便、适用范围广泛。

但是，对于差分数学模型的评价需要考虑以下几个方面：
1. 模型的拟合效果：通过模型的拟合效果来评价模型的可靠性和精度。

通常使用残差和拟合均方误差来评估模型的拟合效果，如果残差较小、拟合均方误差较小，说明模型拟合效果较好。

2. 模型的稳定性：模型的稳定性是指当输入变量的值发生微小变化时，输出结果是否会有大幅度的变化。

稳定性好的模型会使得结果更加可靠。

3. 模型的可解释性：模型的可解释性是指模型是否能够清晰、简洁地解释自变量与因变量之间的关系。

如果模型的解释能够使人理解和接受，那么这个模型就会显得更加可靠。

4. 模型的适用性：模型的适用性是指模型可用于解决实际问题的情况。

模型应该能够适应不同的数据和变化趋势，并且能够进行预测和决策。

5. 模型的复杂度：模型的复杂度是指不同变量之间的关系是否能够用较为简单的数学公式进行解释。

如果模型过于复杂，就会使得结果难以理解和接受。

综上所述，对差分数学模型进行评价和分析应综合考虑以上几点，以得出一个全面、客观、准确的结论。

讨论对流程的差分格式的精度及稳定性的认识

讨论对流程的差分格式的精度及稳定性的认识

讨论对流程的差分格式的精度及稳定性的认识下载提示：该文档是本店铺精心编制而成的，希望大家下载后，能够帮助大家解决实际问题。

文档下载后可定制修改，请根据实际需要进行调整和使用，谢谢！本店铺为大家提供各种类型的实用资料，如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等，想了解不同资料格式和写法，敬请关注！Download tips: This document is carefully compiled by this editor. I hope that after you download it, it can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you! In addition, this shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!讨论对流程的差分格式的精度及稳定性的认识1. 引言对流程是流体力学中的重要现象，研究对流过程的数值模拟需要考虑数值方法的精度和稳定性。

差分格式的稳定性与收敛性1

差分格式的稳定性与收敛性1
得到的，而在此边值问题的解是
(11)
v( x)
M ( x a)(b x) . 2
因为 v( x) 是 x 的二次函数，它的四阶导数为零，从式(2)、 (3)看到 v( x) 在点 xm 的二阶中心差商与 v '' ( xm ) 相等，因此差分方 程(10)的解等于边值问题(11)的解，即
Lu u ''' q( x)u f ( x), a x b, u (a) , u (b) ,
其中 q ( x) 、 f ( x) C a, b , q( x) 0. 将区间 a, b 分成 N 等份，记分点为 xm a mh, m 0,1, , N , 这里步长 h b a .利用泰勒公式，得
和
(1) (2) 的解分别记为 u m 和 um ，其中差分算子 Lh 由式(5)定义，则方程 组(5)的解 u m 为 (1) (2) um um um
(8) (9)
由极值原理可知
(1) um max , , m 1, 2, , N 1 .
3
(2) 接下来再估计 um ，考虑差分方程
Rm (u) Lu( xm ) Lhu( xm ) ，称 Rm (u) 是用差分算子 Lh 代替微 记 分算子 L 所产生的截断误差.由式(2)，二阶中心差商代替二阶 导数所产生的截断误差 Rm ，从式(4)和式(5)可以得出
Rm Lh (u ( xm ) um ) ， Rm 称为差分方程（5）的截断误差.
N
(1)
1 [(u ( xm1 ) 2u ( xm ) u ( xm1 )] u '' ( xm ) Rm 2 h Nhomakorabea其中

差分格式的定性分析

差分格式的定性分析
c 2 (1 − r )∆t ∂ 2u 耗散性和稳定性： 耗散项： 2r ∂x 2
∆t 若截断误差项含 O ， ∆x ∆t 则须 → 0才相容 ∆x ∆x
若c<0则恒为负，相应差分格式为逆耗散格式，按 Hirt论断，一切逆耗散格式均不稳； 该耗散项为正的条件是 r ≤ 1 或 色散性： 色散项
9
1 Neumann-Richtmyer人为耗散法
(1950，JAP, 21, 232)
在动量方程中加入人为粘性压力项 q N ，其人为粘性系数比例于速度梯度
x （若激波在x向传播）
2 2 ∂u ρ(b∆x) qN = ∂x 0
∂u 当 <0 ∂x 当 ∂u ≥0 ∂x
方案2
U U
max i i min
= max( U id− 1 , U id , U id+ 1 ; U in− 1 , U in , U in+ 1 ) = min( U id− 1 , U id , U id+ 1 ; U in− 1 , U in , U in+ 1 )
max i
此法提高 U
减小 U
式中ρ为密度，u为速度，b为可调参数，一般取b=1.5-2 最后激波区的厚度大致为 δ x = π 即4－5个格距
2 b∆ x ≅ 2 .5b∆x γ +1
10
2 Lapidus人为耗散法
对各分量方程同时加入人为耗散项
(1967, JCP, 2, 154)
∂U ∂F ∂G + + + T x + Ty = S ∂t ∂x ∂y
强耗散格式（或引入人为耗散） 弱耗散格式 激波波头的振荡

差分格式的稳定性讨论

差分格式的稳定性讨论

差分格式的稳定性讨论
关于差分格式的稳定性讨论
本文在差分格式稳定性的概念的基础上,按照稳定性的定义来验证某个差分格式,验证差分格式是否稳定,往往比较复杂.讨论稳定性有很多方法,如矩阵方法,离散扰动方法.Hitt呲启示性方法和Fouder方法等等.每个方法各有优劣,本文讨论应用最广泛,也较为方便的Fourier方法.
作者：李卓识作者单位：吉林农业大学,吉林,长春,130118 刊名：网络财富英文刊名：INTEMET FORTUNE 年，卷(期)：2009 ""(23) 分类号：G64 关键词：差分格式稳定性 Lax等价定理 Foutier方法。

【计算机应用】_差分格式的稳定性_期刊发文热词逐年推荐_20140727

【计算机应用】_差分格式的稳定性_期刊发文热词逐年推荐_20140727

2012年 序号
科研热词 推荐指数 1 解的存在唯一性 3 2 紧致差分格式 3 3 收敛性 3 4 1 5 非线性schroedinger方程 1 6 非线性schr(o)dinger方程 1 7 非定常 1 8 涡量-速度法 1 9 格子波尔兹曼方法 1 10 有限差分法 1 11 广义分布函数 1 12 对角占优 1 13 不可压navier-stokes方程组 1 14 一维burgers守恒型方程 1 15 weno差分格式 1 16 unique solvability 1 17 nonlinear schrsdinger equation1 18 linearized compact difference 1 scheme 19 convergence 1
2013年 序号 1 2 3 4 5
2013年 科研热词 紧致差分格式 稳定性 收敛性 唯一性 richardson外推 5
科研热词 推荐指数 非线性leland方程 1 稳定性 1 数值试验 1 并行计算 1 交替分段crank-nicolson(asc-n)格式 1
2008年 序号 1 2 3 4
科研热词 稳定性 子域精细积分隐格式 三次样条函数 一维抛物型方程
推荐指数 1 1 1 1
2009年 序号
科研热词 1 紧致格式 2 污染传播 3 数值方法
推荐指数 1 1 1
2010年 序号 1 2 3 4
科研热词 推荐指数 稳定性 1 收敛性 1 延迟抛物微分方程 1 crank-nicolson差分格式 1

一类半线性抛物方程的差分格式及收敛性和稳定性

一类半线性抛物方程的差分格式及收敛性和稳定性

维普资讯
５２ ５
工
程
数
学
学
报
第 ２ 卷 ５
（ ） ， 二阶可导，且存在常数 Ｃ ＞０ Ｋ ０ Ｈ２ （） １ 及 ，使得当ｌ Ｃ ＋ 时有 ｓ ｉ ｏ
ｉ￣ Ｉｓ，， ｓ，ｆ（Ｉ Ｃ Ｉ ｛ （ｌＩ（ Ｉｉ ｓ） １ ｉ ｆ ） ，） ） ａ
引理１ 】 设ｕ Ｗ∈ｖ ，则有 【 。 ， ｈ
ｍ — —１
一
ｈ
：
２ （ 仇） ｔ （ １ ）ＳＷ一／） 叫 ：ｈ ： ｕ一／ （ｘ ｉ１２
ｉ １ ：
（） ９ （０ １） （１ １）
ｉ １
训。 。Biblioteka ｍ 一１ 一－ Ｉ Ｉ ｉ 训
１ －
维普资讯
第２卷 第３ ５ 期
２ ０ 年 ０ 月 ０８ ６
工
程
数
学
学
报
Ｖ １２ ｏ ３ ｏ．５Ｎ ．
Ｊ ｎ ０ ８ ｕ ｅ２ ０
ＣＨＩ ＮＥＳＥ Ｊ ＯＵＲＮＡＬ ＯＦ ＥＮＧＩ ＮＥＥＲＩ ＮＧ ＡＴＨＥＭ ＡＴＩ Ｍ ＣＳ
造 了一个计算 简便 的线性 化二层格 式，证明了该差分格式解 的存在唯一性、收敛 性、和 无条件
稳定 性，并给出 了在离散 模和 。 ｏ模意义下 收敛 阶数为Ｏ（ ＋丁 ） ｈ 。最后给 出了数值例子验
证 了分 析 结 果 。
由非线性抛物型方程理论［知道 ，当， ｕ ， （） 】 （） 及边界函数满足适 当条件 时，（） ３存在唯 １一 ） （
这 里 ｄ＞ ０ 为常 数 ，Ｑ Ｔ＝ ｆ，１ ０Ｔ１ ０ １ ×ｆ， 。
该模型在化学反应和扩 散、神经 网络传 导、生物竞争 ［３ １】 － 等方面有许 多应用 。研究该 问题 的

差分格式的定性分析

差分格式的定性分析

通常只关心稳定条件本身，则直接用下述方法：
ε
m
( x
j
, t
n
) = e
α t
n
e
ik
m
x
j
= e
α n ∆ t
e
ik
m
j∆ x
代入差分方程，解出放大因子
e
α ∆ t
稳定性条件是放大因子的模小于等于1。
7
10. 人为耗散和通量改正
即使满足稳定性条件，当计算进行到一定时间步数时，仍 出现不稳定，此源于非线性效应。即使格式自身存在的数 , 值耗散（或叫隐耗散）, 在激波区一般不足以克服波头振 荡等计算不稳定现象。故需要在格式中人为加入耗散项。
强耗散格式（或引入人为耗散） 弱耗散格式 激波波头的振荡
过分平滑
为兼顾稳定性和精度，有各种办法，其中广泛应用的是 通量改正法 （Flux-Corrected Transport） 1. 正耗散：从守恒型方程出发，采用强耗散格式，保证无振荡 2. 逆耗散：对差分解在非激波区进行逆耗散，以抵消差分解的耗散误差
A 1 = F H1 − F L 1 第二步：用逆耗散通量 i ± i± i±
2
2
消除耗散
2
U
n +1 i
=U
d i
∆t − (A 1 − A 1 ) i− ∆x i + 2 2
13
U
弱耗散格式
强耗散格式
x 过度的改正将出现波头振荡，故须对逆耗散通量
A
i± 1 2
加以限制
14
AC 1 = C
i± 2
9
1 Neumann-Richtmyer人为耗散法
(1950，JAP, 21, 232)

差分格式稳定性分析和基于Matlab编程对抛物线方程的数值计算

差分格式稳定性分析和基于Matlab编程对抛物线方程的数值计算
其解析解：
u( x, t ) exp( 2t )sin( x)
解：当取 h=0.1，t=0.01 和 h=0.1，t=0.001 计算时，网比 r 分别为 1 和 0.1，古典显式格式要求网比小于等于 1/2 才稳定，而古典隐式格 式无条件稳定，所以本题采用古典隐式格式计算。 古典隐式差分格式： (1 2r )u j ,k r (u j 1,k u j 1,k ) u j ,k 1 （ 1）首先取 h=0.1，t=0.01，即 r=1，并将结果与解析解比较。 计算代码：
k 1 k k k 1 0 u j 0 u j 1 0 u j 1 0 1 u j 0 k 1 k k k 1 0 0 0 0 1 0 j j 1 j 1 j k u j k Wj k j 1 0 k 1 0 k 0 k 0 1 k 0 W j 0 0 W j 1 0 0 W j 1 1 Wj 1 0
W jk V k ei x 1 0 k 1 (ei h e i h ) 1 k V 0 V 1 1 0
则传播因子为：
2 cos( h) 1 G ( , ) 1 1 1 0
解：采用古典隐式格式计算，取 h=0.02， t=0.0001，即 r=0.25， 画出 物体表面温度在 0.5 秒时间内的分布图，计算代码如下：
clc;clear all; h=0.02;t=0.0001;n=1; r=t/h^2; s=0:0.02:1;s1=0:0.02:0.5;s2=0.5:0.02:1; k=0:0.0001:0.5; c=1:51;d=1:5001; u=zeros(5001,51); B=zeros(48,1);

差分方法稳定性介绍

差分方法稳定性介绍

03
多尺度问题的求解
多尺度问题广泛存在于科学和工程领 域，对差分方法的稳定性提出更高要 求。未来研究中，将更加注重多尺度 问题的求解方法和技术研究。
THANKS
感谢观看
差分方法稳定性介绍
• 引言 • 差分方法的基本原理 • 差分方法的稳定性分析 • 差分方法的误差分析 • 提高差分方法稳定性的措施 • 差分方法稳定性的应用案例 • 总结与展望
01
引言
差分方法的概念
差分方法
差分方法是一种数值计算方法，用于求解微分方程的近似解。它通过构造差分 格式来逼近微分方程的导数，从而将微分方程转化为代数方程进行求解。
差分方法的稳定性分析
稳定性定义
数值稳定性
差分方法在数值计算过程中，对于初 始条件和边界条件的小扰动，解的变 化能够保持有界，即不会因计算步数 的增加而无限放大。
渐近稳定性
当计算步数趋于无穷时，差分方法的 解能够收敛到真实解，即误差能够逐 渐减小并趋于零。
稳定性判据
要点一
Lax-Richtmyer稳定性判据
对于线性偏微分方程，如果差分格式能够保持离散能量不 增长，则该格式是稳定的。该判据提供了判断差分格式稳 定性的一个充分条件。
要点二
Courant-Friedrichs-Lewy （CFL…
对于显式差分格式，为了保证计算的稳定性，时间步长与 空间步长之间需要满足一定的关系，即CFL条件。该条件 给出了时间步长的上限。
边界条件的处理
Dirichlet边界条件
直接给出边界上的函数值，处理简单。
Neumann边界条件
给出边界上的法向导数值，需要通过差分 近似进行处理。
Robin边界条件
周期边界条件

差分格式稳定性及数值效应

差分格式稳定性及数值效应
1、在编程过程中，是否只有先选取一个较大的区域，之后缩小来寻找比较清晰的区பைடு நூலகம்。
2、其次，在不同的区间中是一个周期函数，选取不同的区间是否有影响？
3、在修正迎风格式中，当a较大时，区间的极大的产生了偏离，如何能够快速的找到我们需要的区间进行分析，还是只能选取较大区域后一步步缩小范围？
(3)a=4时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式的震荡更加强烈。
由上得出，稳定性对差分格式求解偏微分方程有重大意义。一个差分格式是否好，是否可用，首先要判定它是否稳定并找到稳定性条件。修正迎风格式强大的稳定性在解决一阶线性双曲线方程中有着很强的实用价值。
另外在实验中也有一些其他的问题：
取a=1,2,4, h=0.1,τ=0.08,t=4,对不同的差分格式（迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式，修正迎风格式）及不同的a值进行迭代计算，来讨论分析差分格式的稳定性。
实验原理：
1.迎风格式：
运算格式：
x-Friedrichs格式：
运算格式：
x-Wendroff格式：
(1)a=1时，迎风格式，Lax-Friedrichs格式，修正迎风格式的计算结果与解析解近似情况较好，而Lax-Wendroff格式则在间断点处形成双波现象，这符合Lax-Wendroff格式为二阶迭代格式的性质。
(2)a=2时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式都出现了比较强烈的震荡。
3在修正迎风格式中当较大时区间的极大的产生了偏离如何能够快速的找到我们需要的区间进行分析还是只能选取较大区域后一步步缩小范围
差分格式稳定性及数值效应
F1107102 5112409006陆逢源实验目的：

第五讲——显式差分和隐式差分(5)

第五讲——显式差分和隐式差分(5)

( n 1)
Bc
( n)
T1n 1 n 1 T2 T3n 1 n 1 T4 T n 1 5n 1 T 6
0 0 0 0 1 s 0 0 s 2 2s 0 s 2 2s s 0 B 0 s 2 2s s 0 0 0 0 s 2 2s 0 0 0 0 0 0 0 0 0 s 1
0 0 0 0 0 1 s 0 0 0 s 2 2s 0 s 2 2s s 0 0 A 0 s 2 2s s 0 0 0 0 0 s 2 2s s 0 0 0 0 0 1
内部节点：
边界节点：
0 0 0 0 1 s 0 0 s 2 2s 0 s 2 2s s 0 B 0 s 2 2s s 0 0 0 0 s 2 2s 0 0 0 0 0 0 0 0 0 s 1
G(Y (t ), Y (t t )) 0
例子：
1. 显式差分格式：
左端：n+1时刻的值； 右端：n时刻的值。
特点：结构简洁，直接求解，求解速度快。
但是，时间步长需满足：
显式差分格式才能得到稳定的数值解，否则，数值解将会不稳定而振荡。
显示差分格式示意图
2. 隐式差分格式：
时间一阶精度 空间二阶精度
回顾
1. 有限差分法基础 2. 差分格式 3. 差分方程 4. 边界条件的处理
5. 相容性、稳定性和收敛性
回顾
1. 有限差分法的相容性、稳定性和收敛性 相容性：针对差分格式而言，在时间步长和空间步长趋近于零的情况下， 如果差分格式的截断误差（差分格式与原有偏微分方程之差）的模趋近于零， 则该差分格式与原偏微分方程是相容的，或称该差分方程与原偏微分方程 具有相容性。
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
另外在实验中也有一些其他的问题：
1、在编程过程中，是否只有先选取一个较大的区域，之后缩小来寻找比较清晰的区域。
2、其次，在不同的区间中是一个周期函数，选取不同的区间是否有影响？
3、在修正迎风格式中，当a较大时，区间的极大的产生了偏离，如何能够快速的找到我们需要的区间进行分析，还是只能选取较大区域后一步步缩小范围？
差分格式稳定性及数值效应
F11.以对流方程为例，分析4种差分格式的误差。
2.了解4种差分格式的稳定性
实验问题：
对于一阶线性双曲型方程：
取a=1,2,4, h=0.1,τ=0.08,t=4,对不同的差分格式（迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式，修正迎风格式）及不同的a值进行迭代计算，来讨论分析差分格式的稳定性。
实验原理：
1.迎风格式：
运算格式：
x-Friedrichs格式：
运算格式：
x-Wendroff格式：
这种格式构造是采用Taylor级数展开和微分方程本身得到，运算格式：
4.修正迎风格式：
其中 是 取整数部分， = 。根据之后的理论分析可以得到这是一个无条件稳定结构。
五实验结果：分别按顺序为迎风格式，Friedrichs,Wendroff,修正迎风
(2)a=2时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式都出现了比较强烈的震荡。
(3)a=4时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式的震荡更加强烈。
由上得出，稳定性对差分格式求解偏微分方程有重大意义。一个差分格式是否好，是否可用，首先要判定它是否稳定并找到稳定性条件。修正迎风格式强大的稳定性在解决一阶线性双曲线方程中有着很强的实用价值。
a=1：
a=2:
a=4:
六实验总结：
本次实验，通过4种差分格式求解T=4时的解并与解析解画图比较，可以看出：
(1)a=1时，迎风格式，Lax-Friedrichs格式，修正迎风格式的计算结果与解析解近似情况较好，而Lax-Wendroff格式则在间断点处形成双波现象，这符合Lax-Wendroff格式为二阶迭代格式的性质。