研究有限差分格式稳定性的其他方法 - 报告材料
2015 年秋季学期研究生课程考核
（读书报告、研究报告）
考核科目:偏微分方程数值解法
学生所在院（系）:理学院数学系
学生所在学科:数学
学生姓名:H i t e r
学号:1X S012000
学生类别:
考核结果阅卷人
研究有限差分格式稳定性的其他方法
摘要
偏微分方程的求解一直是大家比较关心的一个问题，而有限差分格式则是求解偏微分方程时常用并且有效的一个方法。因此，研究有限差分格式的性质就显得尤为重要。在课上我们已经跟着老师学习了运用Fourier方法研究有限差分格式的稳定性，但是在很多研究有限差分格式稳定性的问题中仅仅会用Fourier方法是不够的，所以在本篇论文中，将会介绍其他三种常用的研究有限差分格式稳定性的方法，分别是：Hirt启示型方法、直接方法（或称矩阵方法）和能量不等式方法。
关键字：偏微分方程；有限差分格式；稳定性
Abstract
The solution of partial differential equations has been more concerned with a problem, and the finite difference scheme is a common 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 commonly 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; stability
1 前言
微分方程的定解问题就是在满足某些定解条件下求微分方程的解。在空间区域的边界上要满足的定解条件称为边值条件。如果问题与时间有关，在初始时刻所要满足的定解条件，称为初值条件。不含时间而只带边值条件的定解问题，称为边值问题。与时间有关而只带初值条件的定解问题，称为初值问题。同时带有两种定解条件的问题，称为初值边值混合问题。定解问题往往不具有解析解，或者其解析解不易计算。所以要采用可行的数值解法。有限差分方法就是一种数值解法，它的基本思想是先把问题的定义域进行网格剖分，然后在网格点上，按适当的数值微分公式把定解问题中的微商换成差商，从而把原问题离散化为差分格式，进而求出数值解。此外，还要研究差分格式的解的存在性和唯一性、解的求法、解法的数值稳定性、差分格式的解与原定解问题的真解的误差估计、差分格式的解当网格大小趋于零时是否趋于真解（即收敛性），等等。有限差分方法具有简单、灵活以及通用性强等特点，容易在计算机上实现。在课上我们已经跟着老师学习了运用Fourier方法研究有限差分格式的稳定性，但是在很多研究有限差分格式稳定性的问题中仅仅会用Fourier方法是不够的，所以在本篇论文中，将会介绍其他三种常用的研究有限差分格式稳定性的方法，分别是：Hirt 启示型方法、直接方法和能量不等式方法。
2 Hirt启示性方法
2.1 方法概述
Hirt启示性方法是一种近似分析方法。主要是把差分格式在某确定点上作泰勒级数近似展开，把高阶误差略去，只留下最低阶的误差项。如果差分格式是相容的，那么这样得到的
新的微分方程（称之为第一微分近似或修正微分方程）与原来的微分方程相比只增加了一些含小参数的较高阶导数的附加项。Hirt 方法就是利用第一微分近似的适应性来研究差分格式的稳定性。Hirt 方法的判别准则是这样的：如果第一微分近似是适定的，那么原来微分方程的差分格式是稳定的，否则不稳定。其实所述的微分格式是原来微分方程问题的相容的差分格式，那么也可以看作第一微分近似问题的相容的差分格式。如果第一微分近似问题是不适定的，那么它的差分格式将不稳定[1]。
2.2 操作方法
先给出几个方程
0,,0,0>∈>=??+??t R x a x
u a t u （2.1） ΛΛ,2,1,0,,2,1,0,011=±±==-+-++n j h
u u a
u u n
j
n j n j
n j τ
（2.2）
01
1=-+--+h
u u a
u u n j n j n j
n j τ
（2.3）
考虑对流方程（2.1）的差分格式（2.3），在点）
，（n j t x 进行Taylor 技术展开，有 )(][2][)
,(),(2221h O x u h x u h
t u u t x u n
j n j n j n j +??-??=-- )(][2][)
,(),(2221ττ
O t
u h t u t u u t x u n
j n j n j n j +??-??=-+ 利用对流方程（2.1），有
22
222)(x
u a x u a t t u ??=??-??=?? 因此，在点）
，（n j t x 上，有差分方程（2.3）可以得到 )(2222222h O x
u a ah x u a t u ++??-=??+??ττ）（ 略去高阶误差项，得出第一微分方程近似
2
2222x u
a ah x u a t u ??-=??+??）（τ 要使上面的抛物型方程有意义，必须有
02
22>-τ
a ah 而上面的不等号改为等号，则就化为原来的对流方程。在这两种情况下，相应的问题是适定
的。即第一微分近似适定的条件是
02
22≥-τa ah
由此得出差分格式（2.3）的稳定性条件是1≤λa ，其中h
τ
λ=
。此结论与Fourier 方法分析
得到的结论是一致的。
下面我们再来分析逼近对流方程（2.1）（仍设0>a ）的差分格式（2.2）的稳定性。模仿上面的推导可以得到它的第一微分近似是
22222x
u a ah x u a t u ??+-=??+??）（τ 可以看出22x
u
??的系数小于0，因此第一微分近似是不适定的，从而推出差分格式（2.2）是
不稳定的。
3 直接方法
关于抛物型方程初值问题的差分格式的稳定性问题，可以用直接方法（或称矩阵方法）来研究。下面用具体例子来说明这个方法的基本思想及使用方法。
考虑常系数扩散方程的初值问题
???
???
?>==∈=>∈>??=??0,0),(),0(),0(),()0,(0
),,0(,0,022t t l u t u l x x u x u t l x a x u
a t
u （3.1） 采用显示差分格式来逼近，即
???
???
?>==-==-=>+-=--++0
,01,2,1),(1,2,1,0,20000
2
1
11n u u J j x u u J j n h u u u a u u J n j j n j n j n j n j n j ΛΛτ （3.2） 其中l Jh =。先把差分格式（3.2）写成
1,2,1
,)21(111-=+-+=-++J j u a u a u a u n j n j n j n j Λλλλ （3.3） 其中2
h τ
λ=
。可以把（3.3）写成向量形式，即
]00
[]][21212121[
][001
2
2
111
1212
1
1J
n
n
J
n
J n n
n J n J n n u u a u u u u a a a a a a a a a a u u u u M M O O
O
M λλλλλ
λλ
λλ
λ
λ+----=--+-+-++ （3.4）
如果令
T
n J n n n u u u u ),,,(1
21-=Λ 并考虑到00==n
J n u u ，则（3.4）式可以写成
n n Au u =+1 （3.5）
其中
]21212121[
λ
λλλλ
λλλλλa a a a a a a a a a A ----=O
O
O （3.6）
从显示格式出发，得到方程组（3.5）式，也可以理解为较为一般的形式，即对于逼近初值问题（3.2）的其他二层格式也可以化为（3.5）式的形式。当然此时A 不是（3.6）式所表示的形式。如果差分格式是二层隐式格式。则A 为C B 1
-这种形式。因此（3.5）式这种形式可理解为既包含二层显示格式又包含二层隐士格式的较为一般的形式。
引入误差向量~
n
n
n u u z -=，其中n
u 是差分方程（3.5）的精确值（理论值），~n
u 是差分方程（3.5）经数值求解得到的值（包括了舍入误差等）。显然，n
z 满足
n n Az z =+1 （3.7）
从而推出
0z A z n n = （3.8）
差分格式（3.5）的稳定性就要求
0,≥≤n K z n （3.9）
其中?为向量的2-数。由于
02
z A z n n ?≤
因此（3.9）式成立的充分必要条件为
M A n
≤2
（3.10）
上述采用2-数，当然也可以采用其他类型的数。对于稳定性条件（3.10），可以仿Fourier 方法中的推导，得到一些结论：
（1）谱半径条件
τρM A +≤1)( （3.11）
是差分格式稳定的一个必要条件，其中M 为常数。
（2）如果矩阵A 是一个正规矩阵，则（3.11）式也是格式稳定的一个充分条件。 下面讨论差分格式（3.5），（3.6）的稳定性。矩阵（3.6）是对称矩阵，所以只要使条件（3.11）成立即可。现在来计算A 的特征值。
令)1(-J 阶方阵
]01
101
10110[
O O O
=S
则A 可以表示为
S a I a A λλ+-=)21(
其中I 为)1(-J 阶单位矩阵。由此可知，关键是求出S 的特征值和特征向量。
设γ和T J w w w w ),,,(121-=Λ分别为S 的特征值和特征向量，
w Sw γ=
写成分量的形式有
?
?
?==-==+-++02
,,1,0,0021J j j j w w J j w w w Λγ （3.12） 先求出j w ，再求出S 的特征值γ。由于S 为对称矩阵，所以其特征值γ为实数。由Gerschgorin 定理知，
∑-≠≤-1
J k
j kj kk s s γ
其中kj s 为矩阵S 的元素。由此得到2≤γ。（3.12）式的第一式为常系数线性差分方程。设其解具有如下形式：
0,≠=μμj j w
将它代入（3.12）式的第一式，便得到关于μ的一元二次方程
012=+-λμμ
此方程称为（3.12）式的第一式的特征方程。由于2≤γ，所以其解为
2)2
(12γ
γ
μ-±=
i 其中1-=
i 。可以看到
1)2
(1)2
(222
=-+=γγμ
取2)2
(1sin ,2cos λ
?γ
?-==
，则?μi e ±=。因此差分方程（3.12）的解可以表示为
完整版有限差分方法概述.doc
有限差分法（ Finite Difference Method，简称FDM）是数值方法中最经典的方法，也是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较 早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分 为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上 述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后 差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 下面我们从有限差分方法的基本思想、技术要点、应用步骤三个方面来深入了解一下有限差分方法。 1.基本思想 有限差分算法的基本思想是把连续的定解区域用有限个离散点 构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。在采用数值计算方法求解偏微分方程时，再将每一处导数由有限差分近似公式替代，从而把求解偏微分方程的问题转换成求解代数方程的问题，即 所谓的有限差分法。 2.技术要点 如何根据问题的特点将定解区域作网格剖分；如何把原微分
差分放大电路仿真02605
苏州市职业大学实验报告姓名：学号：班级：
二、选好元器后，将所有元器件连接绘制成仿真电路（见图 1） R3 6.8k Q 三、仿真分析 1.静态工作点分析 1）调零。信号源先不接入回路中，将输入端对地短接，用万用表测量两个输出 节点，调节三极管的射极电位，使万用表的示数相同，即调整电路使左右完 全对称。测量电路及结果如图2所示 2）静态工作点调试。零点调好以后，可以用万用表测量 Q1、Q2管各电极电位, 结果如图 3 所示，测得 I B 1 15 A , I C 1 1.089mA , U CE 5.303V 。 2.测量差模放大倍数 将函数信号发生器XFG1的“ +”端接放大电路的R1输入端，“一”端接R2输入 端，COM 端接地。调节信号频率为1kHz ,输入电压10mV 调入双踪示波器，分别 接输入输出，如图4所示，观祭波形变化，示波器观祭到的差分放大电路输入、 输出波形如图5所示 R4 6.8k Q R1 ■ 酉 2 ?R6 >510 Q <3 ------- Q1 R8 12k Q 12 V 双端输入、 100Q Key=A 丄V2 -— 12 V 11 R5 5.1k 10 双端输出的长尾式差分放大电路 8 Q ■ 4 Q2 2N3903 R2 AAAr-| 2k Q 7 50% Rp1
4.607 V H-、4 -Q *: LR3 S : : ?6+BkQ : a ）： >R4 ：>G.?kn ............ R& '''' ---------- VA ---------- it：：12W5：：: 1 F ■! ■ I R1 .,,斗，- VA- :7W. . \ ■1 2M39G 3 :R2 ： : 2K1： 2N39G3 -” R6 5100 : ::5C% :10QQ ::Key=A 丄V2「::二12W TV '' 图2差分放大器电路调零
三轮DES差分分析实验报告-刘杰
DES 差分分析实验报告 四大队四队五班 刘杰 一、实验目的 差分密码分析是一种选择明文攻击，是现代分组密码分析的重要方法之一，也是理论分析密码算法和算法抗攻击测试的重要依据之一。本实验通过3轮DES 简化算法的差分分析来达到加深学员对差分分析方法原理的理解和利用该原理分析实际问题的操作能力。 二、实验内容 （1）3轮DES 简化算法的差分分析； （2）通过三组明密文对（每组两个相关明文和相应密文），利用差分原理提取密钥。 明 文 密 文 748502CD38451097 03C70306D8A09F10 3874756438451097 78560A0960E6D4CB 486911026ACDFF31 45FA285BE5ADC730 375BD31F6ACDFF31 134F7915AC253457 357418DA013FEC86 D8A31B2F28BBC5CF 12549847013FEC86 0F317AC2B23CB944 三、实验原理 设DES 两个明密文对:=00m L R ***=00m L R =33c L R *** =33c L R 计算过程: (,)(,)(,)(,)=⊕=⊕=⊕⊕322312300123R L f R k R f R k L f R k f R k
(,)(,)****=⊕⊕300123R L f R k f R k 令：*'=⊕000L L L (,)(,)(,)(* **''=⊕=⊕⊕⊕⊕333001012323R R R L f R k f R k f R k f R k 观察得：在本次实验原始数据中，明文对*=00R R ，即* '=⊕=00000000000R R R 则(,)(,)** ''=⊕=⊕⊕33302323R R R L f R k f R k 同时有：=00m L R ***=00m L R =23R L ** =23R L 则可计算出：*'=⊕000L L L *'=⊕333R R R (,)(,)* ''⊕=⊕232330f R k f R k R L 则可得出： S 盒输入差：(())(())()()* *⊕⊕⊕=⊕232333E R k E R k E L E L S 盒输出差：()*-''⊕=⊕13 0D D P R L 分析过程： 令：()()*⊕=3312345678E L E L B B B B B B B B ()-''⊕=13 012345678P R L C C C C C C C C ()=312345678E L A A A A A A A A =312345678 k J J J J J J J J ()⊕=3312345678E L k X X X X X X X X *()⊕=3312345678E L k Y Y Y Y Y Y Y Y 基本思路：（分别计算12345678J J J J J J J J ） {|,()()∈=⊕⊕=⊕=i i i i i i i J T e s t x A x y B S x S y C ,,,,,,,=12345678i 对于本次实验的3个具有明文差（*,0）的明密文对，则可构造上面的3个 Test 集合，显然 ()()( )∈12 i i i i J Test Test Test t ,,,,,,,=12345678i 一种确定Ji 的直接方法： 1.建立26=64长度的数组J[64]={0}; 2.对Testi(r)，r = 1,2,…,t ，若a ∈Testi(r)，则 J[a] = J[a] + 1。 3.若J[b] =3，则6比特串b 就是可能的密钥比特 Ji 。 四、实验环境 Microsoft visual c++ 五、实验步骤 （1）计算简化算法第3轮S 盒输入差
Poisson方程九点差分格式_米瑞琪
数值实验报告I 实验名称Poisson方程九点差分格式实验时间2016年 4 月 15 日姓名米瑞琪班级信息1303学号04成绩 一、实验目的，内容 1、理解Poisson方程九点差分格式的构造原理； 2、理解因为网格点的不同排序方式造成的系数矩阵格式的差异； 3、学会利用matlab的spdiags(),kron()函数生成系数矩阵； 二、算法描述 针对一个Poisson方程问题： 在Poisson方程五点差分格式的基础上，采用Taylor展开分析五点差分算子的截断误差，可以得到： 为了提高算子截断误差的精度，在(1)式中配凑出了差分算子的形式，将原Poisson方程代入（1）式有： 考虑，有：
将（3）代回（2）可得 得到Poisson方程的九点差分格式： 在计算机上实现（4）式，需要在五点差分格式 的基础上在等式两端分别增加一部分，将等式左侧新增的部分写成紧凑格式，有： 对于该矩阵，可以看成是两个矩阵的组合：
以及 则生成这两个矩阵可以采用Kroncker生成，方法类似于五点差分格式。 对于右端添加的关于f(x,y)的二阶导数，可以采用中心差分格式进行近似代替，即： 写成相应的紧凑格式有：
该式中的矩阵又可以分解为两个矩阵的和：
%计算误差 u_real=@(x,y)exp(pi*(x+y))*sin(pi*x).*sin(pi*y); for i=1:N1-1 u_m((i-1)*(N2-1)+1:i*(N2-1))=u_real(x(i),y); end u_v=u_m'; err_d=max(abs(u_d-u_v)); sol=reshape(u_d,N2-1,N1-1); mesh(X,Y,sol) 四. 数值结果 针对课本P93给出的问题，分别采用步长，将计算出的误差列表如下： 步长五点差分格式误差九点差分格式误差 可见采用九点差分格式可以进一步缩小误差，达到更高阶的精度。 五. 计算中出现的问题，解决方法及体会 在生成九点差分格式的时候，等号右端涉及到了对f的二阶偏导，我最初利用符号函数定义了f,随后求出其二阶偏导（仍然是符号函数）之后带入网格点，求f二阶偏导的精确解，但是代入过程相当繁琐，运行速度非常慢，最终我改变策略，选用f关于x,y的二阶中心差分格式替代精确值，最终得到了相对满意的结果。 教 师 评 语 指导教师：年月日
中心差分法的基本理论与程序设计
中心差分法的基本理论与程序设计 1程序设计的目的与意义 该程序通过用C语言（部分C++语言）编写了有限元中用于求解动力学问题的中心差分法，巩固和掌握了中心差分法的基本概念，提高了实际动手能力，并通过实际编程实现了中心差分法在求解某些动力学问题中的运用，加深了对该方法的理解和掌握。 2程序功能及特点 该程序采用C语言（部分C++语言）实现了用于求解动力学问题的中心差分法，可以求解得到运动方程的解答，包括位移，速度和加速度。计算简便且在算法稳定的条件下，精度较高。 3中心差分法的基本理论 在动力学问题中，系统的有限元求解方程（运动方程）如下所示： ()()()() Ma t Ca t Ka t Q t ++= 式中，() a t分别是系统的结点加速度向 a t是系统结点位移向量，() a t和() 量和结点速度向量，,, M C K和() Q t分别是系统的质量矩阵、阻尼矩阵、刚度矩阵和结点载荷向量，并分别由各自的单元矩阵和向量集成。 与静力学分析相比，在动力分析中，由于惯性力和阻尼力出现在平衡方程中，因此引入了质量矩阵和阻尼矩阵，最后得到的求解方程不是代数方程组，而是常微分方程组。常微分方程的求解方法可以分为两类，即直接积分法和振型叠加法。 中心差分法属于直接积分法，其对运动方程不进行方程形式的变换而直接进行逐步数值积分。通常的直接积分是基于两个概念，一是将在求解域0t T内的任何时刻t都应满足运动方程的要求，代之仅在一定条件下近似地满足运动方程，例如可以仅在相隔t?的离散的时间点满足运动方程；二是在一定数目的t?区域内，假设位移a、速度a、加速度a的函数形式。 中心差分法的基本思路是用有限差分代替位移对时间的求导，将运动方程中的速度和加速度用位移的某种组合表示，然后将常微分方程组的求解问题转换为
差动放大器实验报告
差动放大器实验报告 以下是为大家整理的差动放大器实验报告的相关范文，本文关键词为差动,放大器,实验,报告,篇一,实验,差动,放大器,南昌大学，您可以从右上方搜索框检索更多相关文章，如果您觉得有用，请继续关注我们并推荐给您的好友，您可以在工作报告中查看更多范文。 篇一：实验五差动放大器 南昌大学实验报告 实验五差动放大器 一、实验目的 1、加深对差动放大器性能及特点的理解 2、学习差动放大器主要性能指标的测试方法 二、实验原理 下图是差动放大器的基本结构。它由两个元件参数相同的基本共射放大电路组成。当开关K拨向左边时，构成典型的差动放大器。调零电位器Rp用来调节T1、T2管的静态工作点，使得输入信号ui＝0时，双端输出电压uo＝0。Re为两管共用的发射极电阻，它对差模信号无负反馈作用，因而不影响差模电压放大倍数，但对共模信号有较
强的负反馈作用，故可以有效地抑制零漂，稳定静态工作点。 图5-1差动放大器实验电路 1、静态工作点的估算典型电路Ic1=Ic2=1/2Ie恒流源电路Ic1=Ic2=1/2Ic3 2、差模电压放大倍数和共模电压放大倍数 双端输出:Re＝∞，Rp在中心位置时， Ad? 单端输出 △uoβRc ?? △ui Rb?rbe??β)Rp 2 Ad1? △uc11?Ad △ui2 Ad2? △uc21 ??Ad △ui2 当输入共模信号时，若为单端输出，则有 △uc1?βRcR
Ac1?Ac2????c △uiR?r?(1?β)(1R?2R)2Re bbepe 3、共模抑制比cmRR2 为了表征差动放大器对有用信号（差模信号）的放大作用和对共模信号的抑制能力，通常用一个综合指标来衡量，即共模抑制比AA cmRR?d或cmRR?20Logd?db? AcAc 三、实验设备与器材 1、函数信号发生器 2、示波器 3、交流毫伏表 4、万用表 5、实验箱 6、差动放大器集成块 四、实验内容 1、典型差动放大器性能测试 按图5-1连接实验电路，开关K拨向左边构成典型差动放大器。 1)测量静态工作点2)①调节放大器零点 信号源不接入。将放大器输入端A、b与地短接，接通±12V直流电源，用直流电压表测量输出电压uo，调节调零电位器Rp，使uo＝0。调节要仔细，力求准确。 ②测量静态工作点 零点调好以后，用直流电压表测量T1、T2管各电极电位及射极电阻Re两端电压uRe，记入表5－1。
有限差分法实验报告
工程电磁场 实验报告 ——有限差分法
用超松弛迭代法求解 接地金属槽内电位的分布 一、实验要求 按对称场差分格式求解电位的分布 已知： 给定边值：如图1-7示 图1-7接地金属槽内半场域的网格 给定初值)()(.1j 40 100 1j p 1 2j i -= --= ??? 误范围差: 510-=ε 计算：迭代次数N ，j i ,?，将计算结果保存到文件中 二、实验思想 有限差分法 有限差分法（Finite Differential Method ）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数?的泊松方程的问题转换为求解网格节点上? =?= V 100 ? 0 =?0 =?
的差分方程组的问题。 泊松方程的五点差分格式 )(4 1 4243210204321Fh Fh -+++=?=-+++?????????? 当场域中,0=ρ得到拉普拉斯方程的五点差分格式 )(4 1 044321004321??????????+++=?=-+++ 差分方程组的求解方法(1) 高斯——赛德尔迭代法 ][)(,)(,)(,)(,)(,2 k 1j i k j 1i 1k 1j i 1k j 1i 1k j i Fh 4 1 -+++=+++-+-+????? (1-14) 式中：??????=??????=,2,1,0,2,1,k j i ， ? 迭代顺序可按先行后列，或先列后行进行。 ? 迭代过程遇到边界节点时，代入边界值或边界差分 格式，直到所有节点电位满足ε??<-+)(,)(,k j i l k j i 为止。 （2）超松弛迭代法 ][) (,)(,)(,)(,)(,)(,)(,k j i 2k 1j i k j 1i 1k 1j i 1k j 1i k j i 1k j i 4Fh 4 ?????α??--++++=+++-+-+ (1-15) 式中：α——加速收敛因子)21(<<α 可见：迭代收敛的速度与α有明显关系 三、程序源代码 #include
差分法求解偏微分方程MAAB
南京理工大学 课程考核论文 课程名称：高等数值分析 论文题目：有限差分法求解偏微分方程 姓名：罗晨 学号： 成绩： 有限差分法求解偏微分方程 一、主要内容 1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程：具体求解的偏微分方程如下： 2.推导五种差分格式、截断误差并分析其稳定性； 3.编写MATLAB程序实现五种差分格式对偏微分方程的求解及误差分析；
4.结论及完成本次实验报告的感想。 二、推导几种差分格式的过程： 有限差分法（finite-differencemethods ）是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。 推导差分方程的过程中需要用到的泰勒展开公式如下： ()2100000000()()()()()()()......()(()) 1!2!! n n n f x f x f x f x f x x x x x x x o x x n +'''=+-+-++-+-(2-1) 求解区域的网格划分步长参数如下： 11k k k k t t x x h τ ++-=?? -=?(2-2) 2.1古典显格式 2.1.1古典显格式的推导 由泰勒展开公式将(,)u x t 对时间展开得 2,(,)(,)( )()(())i i k i k k k u u x t u x t t t o t t t ?=+-+-?(2-3) 当1k t t +=时有 21,112,(,)(,)( )()(())(,)()() i k i k i k k k k k i k i k u u x t u x t t t o t t t u u x t o t ττ+++?=+-+-??=+?+?(2-4) 得到对时间的一阶偏导数 1,(,)(,)()=()i k i k i k u x t u x t u o t ττ+-?+?(2-5) 由泰勒展开公式将(,)u x t 对位置展开得 223,,21(,)(,)()()()()(())2!k i k i k i i k i i u u u x t u x t x x x x o x x x x ??=+-+-+-??(2-6) 当11i i x x x x +-==和时，代入式(2-6)得
差分编译码实验报告
实验十三差分编译码实验 一、实验目的 掌握差分编码/译码原理 二、实验内容 1、学习差分编译码原理 2、用示波器观察差分编码结果和译码结果 三、基本原理 差分码是一种把符号‘0’和‘1’反映在相邻码元的相对变化上的波形。比如，若以相邻码元的电位改变表示符号‘1’，而以电位不改变表示符号‘0’，如图13－1所示。当然，上述规定也可以反过来。由图可见，这种码波形在形式上与单极性或双极性码波形相同，但它代表的信息符号与码元本身电位或极性无关，而仅与相邻码元的电位变化有关。差分波形也称相对码波形，而相应地称单极性或双极性波形为绝对码波形。差分码波形常在相位调制系统的码变换器中使用。 图13-1差分码波形 组成模块如下图所示： cclk d_out 端口说明： CCLK：编码时钟输入端 DIN：编码数据输入端 Diff－OUT：差分编码结果输出端 DCLK：译码时钟输入端
Diff－IN：差分译码数据输入端 DOUT：译码结果输出端 四、实验步骤 1、实验所用模块：数字编解码模块、数字时钟信号源模块。 实验连线： CCLK：从数字时钟信号源模块引入一高频时钟，如512K。 DIN：从数字时钟信号源模块引入一低频时钟，如16K。 DIFF－OUT与DIFF－IN短接。 DCLK与CCLK短接。 2、用示波器两探头同时观测DIN与DIFF－OUT端，分析差分编码规则。 3、用示波器两探头同时观测DIN与DOUT端，分析差分译码结果。 五、实验报告要求 设信息代码为1001101，码速率为128K，差分码的编码时钟为码速率的四倍，根据实验观察得到的规律，画出差分码波形。
实验四 两级放大电路实验报告
实验四 两级放大电路 一、实验目的 l 、掌握如何合理设置静态工作点。 2、学会放大器频率特性测试方法。 3、了解放大器的失真及消除方法。 二、实验原理 1、对于二极放大电路，习惯上规定第一级是从信号源到第二个晶体管BG2的基极，第二级是从第二个晶体管的基极到负载，这样两极放大器的电压总增益Av 为： 2V 1V 1 i 1 O 2i 2O 1i 2O ,i 2O S 2O V A A V V V V V V V V V V A ?=?==== 式中电压均为有效值，且2i 1O V V =，由此可见，两级放大器电压总增益是单级电压增益的乘积，由结论可推广到多级放大器。 当忽略信号源内阻R S 和偏流电阻R b 的影响，放大器的中频电压增益为： 1be 2 be 1C 1be 1L 11i 1O S 1O 1V r r //R 1 r R V V V V A β-='β-=== 2 be L 2C 2 2be 2L 21O 2O 1i 2O 2V r R //R r R V V V V A β-='β-=== 2 be L 2C 2 1be 2be 1C 12V 1V V r R //R r r //R A A A β?β=?= 必须要注意的是A V1、A V2都是考虑了下一级输入电阻（或负载）的影响，所以第一级的输出电压即为第二级的输入电压，而不是第一级的开路输出电压，当第一级增益已计入下级输入电阻的影响后，在计算第二级增益时，就不必再考虑前级的输出阻抗，否则计算就重复了。 2、在两极放大器中β和I E 的提高，必须全面考虑，是前后级相互影响的关系。 3、对两级电路参数相同的放大器其单级通频带相同，而总的通频带将变窄。 ) dB (A log 20G 式中G G G V u o 2u o 1u uo =+= 三、实验仪器 l 、双踪示波器。 2、数字万用表。 3、信号发生器。 4、毫伏表 5、分立元件放大电路模块 四、实验内容 1、实验电路见图4-1
差分方法实验报告
实验报告 课程名称：计算方法 院系：数学科学系 专业班级：数应1001 学号：1031110139 学生姓名：姚海保 指导教师：沈林 开课时间：2012至2013学年第一学期
一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。 学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求 1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。 2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。 3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。 4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩； 5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。
差动放大器实验报告
差动放大电路的分析与综合（计算与设计）实验报告 1、实验时间 10月31日(周五)17:50-21:00 2、实验地点 实验楼902 3、实验目的 1. 熟悉差动放大器的工作原理(熟练掌握差动放大器的静态、动态分析方法) 2. 加深对差动放大器性能及特点的理解 3. 学习差动放大电路静态工作点的测量 4. 学习差动放大器主要性能指标的测试方法 5. 熟悉恒流源的恒流特性 6. 通过对典型差动放大器的分析，锻炼根据实际要求独立设计基本电路的能力 7. 练习使用电路仿真软件，辅助分析设计实际应用电路 8. 培养实际工作中分析问题、解决问题的能力 4、实验仪器 数字示波器、数字万用表、模拟实验板、三极管、电容电阻若干、连接线 5、电路原理 1. 基本差动放大器 图是差动放大器的基本结构。它由两个元件参数相同的基本共射放大电路组成。 部分模拟图如下 1.直流分析数据 2.直流分析仿真数据 3.交流分析数据 4.交流分析仿真数据 具有平衡电位器的 差动放大器 图是差动放大器的结 构。它由两个元件参数相 近的基本共射放大电路组 成。 1.直流分析数据 2.直流分析仿真数据
3.交流分析数据 4.交流分析仿真数据 具有恒流源的差动放大器 图2-3是差动放大器的结构。它由两个元件参数相近的基本共射放大电路组成。 1.直流分析数据 2.直流分析仿真数据 3.交流分析数据 4.交流分析仿真数据 图3.1 差动放大器实验电路 当开关K 拨向右边时，构成具有恒流源的差动放大器。晶体管 T 3 与电阻3E R 共同组成镜象恒流源电路 , 为差动放大器提供恒定电流E I 。用晶体管恒流源代替发射极电阻 E R ，可以进一步提高差动 放大器抑制共模信号的能 力。 1、差动电路的输入输 出方式 根据输入信号和输出信号的不同方式可以有四种连接方式,即 : (l) 双端输入 -双端输出，将差模信号加在1s V 、2s V 两端 , 输出取自1o V 、2o V 两端。 (2) 双端输入 -单端输出，将差模信号加在1s V 、2s V 两端 , 输出取自1o V 或2o V 到地。 (3) 单端输入一双端输出，将差模信号加在1s V 上,2s V 接地 ( 或1s V 接地而信号加在2s V 上 ), 输出取自1o V 、2o V 两端。 (4) 单端输入 -单端输出 将差模信号加在1s V 上,2s V 接地 ( 或1s V 接地而信号加在2s V 上 ), 输出取自1o V 或2o V 到地。
加法器及差分放大器项目实验报告
加法器及差分放大器项目实验报告 一、项目内容和要求 （一）、加法器 1、任务目的： （1）掌握运算放大器线性电路的设计方法； （2）理解运算放大器的工作原理； （3）掌握应用仿真软件对运算放大器进行仿真分析的方法。 2、任务内容： 2.1 设计一个反相加法器电路，技术指标如下： （1）电路指标 运算关系:)25(21i i O U U U +-=。 输入阻抗Ω≥Ω≥K R K R i i 5,521。 （2）设计条件 电源电压Ec=±5V ； 负载阻抗Ω=K R L 1.5 （3）测试项目 A ：输入信号V U V U i i 5.0,5.021±=±=，测试4种组合下的输出电压； B ：输入信号V KHz U V U i i 1.0,1,5.021为正弦波±=信号，测试两种输入组合情况下的输出电 压波形。 C ：输入信号V U i 01=，改变2i U 的幅度，测量该加法器的动态范围。 D ：输入信号V U i 01=，V U i 1,2为正弦波，改变正弦波的频率，从1kHz 逐渐增加，步长为 2kHz ，测量该加法器的幅频特性。 2.2 设计一个同相加法器电路，技术指标如下： （1）电路指标 运算关系：21i i O U U U +=。 （2）设计条件 电源电压Ec=±5V ； 负载阻抗Ω=K R L 1.5 （3）测试项目 A ：输入信号V U V U i i 1,121±=±=，测试4种组合下的输出电压； B ：输入信号V KHz U V U i i 1,1,121为正弦波±=信号，测试两种输入组合情况下的输出电压 波形。 （二）、差分放大器 1、任务目的： （1）掌握运算放大器线性电路的设计方法； （2）理解运算放大器的工作原理； （3）掌握应用仿真软件对运算放大器进行仿真分析的方法。 2、任务内容 2.1 设计一个基本运放差分放大器电路，技术指标如下： （1）电路指标 运算关系:)(521i i O U U U --=。 输入阻抗Ω≥Ω≥K R K R i i 5,521。 （2）设计条件
一维波动方程的有限差分法
学生实验报告实验课程名称偏微分方程数值解 开课实验室数统学院 学院数统年级2013 专业班信计02班 学生姓名______________ 学号 开课时间2015 至2016 学年第 2 学期
数学与统计学院制 开课学院、实验室：数统学院实验时间：2016年6月20日
1、三层显格式建立 由于题中h 0.1, 0.1h，x 0,1 ,t 0,2，取N 10, M 200，故令网比r 0.1，h X j j h, j 0,1,2,L 10，t k k ,k O,1L ,200 ,在内网个点处，利用二阶中心差商得到如下格式: k 1 k U J 2U J 2- k 1 U j k k U j 1 2U j h2 k U j 1 o h2 略去误差项得到: k 1 U j 其中j 1,2丄9,k 对于初始条件 2 k r U J1 1,2,L ,199，局部截断误差为 U x,0 sin U J k U j k r U j 2 o k 1 U J h2。 (3) 对于初始条件-u x,0 t x，建立差分格式为： sin x j sin Jh , J 利用中心差商，建立差分格式为: 0,1,2,L 10 (4) 对于边界条件将差分格式延拓使综上（3 ）、 （4 ）、 k 1 u j 其中r山o.1 1 U J 2 1 U j 0,即U1二U j1, J 0,1,2,L 10 (5) 0,t 0,2 ，建立差分格式为： U N 0,k 0,1,L ,200 k 0为内点，代入（3）得到的式子再与（5）联立消去 1 1 2 0 ’ 2 0 1 5 r U, 1 1 r U, r J 2 J J 2 （7 ）得到三层显格式如下： U 0,t U 1,t k U0 (6 ) 、 2 k r U j 1 2 1 r2k 2 k U J r U J 1 k 1? U j , J U j （6） 1后整理得到: U j 1 (7) （局部截断误差为 1,2,L 9,k 1,2,L ,199 h2) 1 U j U J sin 1 2 0 2r U J 1 k U o X j k U N sin 2 0 r U j 0,k 0,1,2,L 10 Jh ,J 1 2r2u01, J 1,2,L 9 0,1L ,200 (8) 四?实验环境（所用软件、硬件等）及实验数据文件Matlab
(完整word版)差分放大器设计的实验报告
设计课题 设计一个具有恒流偏置的单端输入-单端输出差分放大器。 学校：延安大学
一： 已知条件 正负电源电压V V V V EE cc 12,12-=-+=+；负载Ω=k R L 20；输入差 模信号mV V id 20=。 二：性能指标要求 差模输入电阻Ω>k R id 10；差模电压增益15≥vd A ；共模抑制 比dB K CMR 50>。 三：方案设计及论证 方案一：
方案二
方案论证： 在放大电路中，任何元件参数的变化，都将产生输出电压的漂移，由温度变化所引起的半导体参数的变化是产生零点漂移的主要原因。采用特性相同的管子使它们产生的温漂相互抵消，故构成差分放大电路。差分放大电路的基本性能是放大差模信号，抑制共模信号好，采用恒流源代替稳流电阻，从而尽可能的提高共模抑制比。 论证方案一：用电阻R6来抑制温漂 ?优点：R6 越大抑制温漂的能力越强； ?缺点：<1>在集成电路中难以制作大电阻； <2> R6的增大也会导致Vee的增大（实际中Vee不
可能随意变化） 论证方案二 优点：（1）引入恒流源来代替R6,理想的恒流源内阻趋于无穷，直流压降不会太高，符合实际情况； （2）电路中恒流源部分增加了两个电位器，其中47R的用来调整电路对称性，10K的用来控制Ic的大小，从而调节静态工作点。 通过分析最终选择方案二。 四：实验工作原理及元器件参数确定 ?静态分析：当输入信号为0时， ?I EQ≈(Vee-U BEQ)/2Re ?I BQ= I EQ /(1+β) ?U CEQ=U CQ-U EQ≈Vcc-I CQ Rc+U BEQ 动态分析 ?已知：R1=R4,R2=R3
有限差分法
有限差分法有限差分法 finite difference method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛
武汉大学差动放大电路实验报告
武汉大学计算机学院教学实验报告 课题名称：电工实验专业：计算机科学与技术2013 年12 月14 日实验名称差动放大电路实验台号实验时数3小时姓名学号年级2013班3班 一、实验目的及实验内容 （本次实验所涉及并要求掌握的知识点；实验内容；必要的原理分析） 一、实验目的 1 、熟悉差动放大器工作原理 2、掌握差动放大器的基本测试方法 实验内容 1.计算下列差动放大器的静态工作点和电压放大 倍数电路图见5.1 信号源已替代 5.1 在图5.1的基础上画出单端输入时和共模输入时的电路图 二、实验环境及实验步骤 （本次实验所使用的器件、仪器设备等的情况；具体的实验步骤） 实验环境： 1.示波器 2.信号发生器 3.数字万用表 4.TPE-A3模拟电路实验箱 3、实验步骤： 1、将电路图5.1接线 2、测量静态工作点 3、测量差模电压放大倍数 4、测量共模电压放大倍数 5、在实验台上组成单端输入的差动电路进行下列实验
三、实验过程与分析 （详细记录实验过程中发生的故障和问题，进行故障分析，说明故障排除的过程和方法。根据具体实验，记录、整理相应的数据表格、绘制曲线、波形图等） 实验内容及数据记录 1、将电路图5.1接线 2、测量静态工作点 ①调零 将放大器输入端V11、V12接地，接通直流电源，调节调零电位器R P，使V O＝0。 ②测量静态工作点：测量V1，V2，V3各极各地电压， 并填入表5.1中。 5.1 对地 电压 Vc1 Vc2 Vc3 Vb1 Vb2 Vb3 Ve1 Ve2 Ve3 测量值 6.29 6.31 -0.74 0 0 - 7.77 -0.61 -0.61 - 8.39 3)测量差模电压放大倍数 在两个输入端各自加入直流电压信号，按有5.2要求测量并记录，由测量得到的数据计算出单端和输出的电压放大倍数。接入到V11t和V12，调节Dc信号源，使其输出为0.1和-0.1. (须调节直流电压源Ui1=0.1V ,Ui2=-0.1V) 4) 测量共模电压放大倍数 将输入端b1和b2 短接，接到信号源的输入端，信号源另一端接地。DC信号先后接OUT1和OUT2 测量有关数据后填入表5.32.，由测量得到的数据计算出单端和双端输出的电压放大倍数，并进一步计算出共模抑制比。 5.2 差模输入共模输入抑制 比测量值计算值测量值计算值计算 值Uc1 Uc2 Uo双Ad1 Ad2 Ad双Uc1 Uc2 Uco双Ac1 Ac2 Ac双CMRR +0.1V 10.08 2.55 7.46 -16. 8616.8 6-33. 71 6.29 6.31 -0.02 0.00 5 0.00 5 0 186.5 -0.1V 6.29 6.31 -0.02 0.00 50.00 5 0 186.5
实验八_差分放大器实验报告
差分放大电路 实验报告 姓名：黄宝玲 班级：计科1403 学号：201408010320 实验摘要（关键信息） 实验目的：由于差分放大器是运算放大器的输入级，清楚差分放大电路的工作原理，有助于理解运放的工作原理和方式。通过实验弄清差分放大器的工作方式和参数指标。这些概念有：差模输入和共模输入；差模电压增益Avd和共模电压增益Avc；共模抑制比Kcmr。 实验内容与规划： 1、选用实验箱上差分放大电路；输入信号为Vs=300mV，f=3KHz正弦波。 2、发射极先接有源负载，利用调零电位器使得输出端电压Vo=0。（Vo=Vc1-Vc2） 3、在双端输入和单端输入差模信号情况下，分别测量双端输出的输入输出波形，计算各自的差模放大倍数Avd。 4、在双端输入共模信号情况下，分别测量双端输出的输入输出波形，计算双端输出共模放大倍数Avc。 5、计算共模抑制比Kcm R 。 最好作好记录表格，因为要记录的数据较多。电路中两个三极管都为9013。 实验环境（仪器用品等） 1.仪器：示波器（DPO 2012B 100MHZ 1GS/s） 直流电源（IT6302 0~30V,3Ax2CH/0~5V,3A） 台式万用表（UT805A） 模拟电路实验箱（LTE-AC-03B）。 2、所用功能区：单管、多管、负反馈放大电路。 实验原理和实验电路 1、实验原理： 差分电路是具有这样一种功能的电路。该电路的输入端是两个信号的输入，这两个信号的差值，为电路有效输入信号，电路的输出是对这两个输入信号之差的放大。 概念梳理：
差模和共模是对于差动放大电路的两个输入端而言的。 A ）差模输入：差动放大电路的两管基极输入的信号幅度相等、极性相反，这样的信号称为差模信号，这样的输入称为差模输入。 差模信号Vid ：即差模输入的两个输入信号之差。 B ）共模输入：差动放大电路的两管基极输入的信号幅度相等、极性相同，这样的信号称为共模信号，这样的输入称为共模输入。 共模信号Vic ：即共模输入的两个输入信号的算数平均值。 C ）差模电压增益Avd ：指差动放大电路对差模输入信号的放大倍数。差模电压增益越大，放大电路的性能越好。 = D ）共模电压增益Avc ：指差动放大电路对共模输入信号的放大倍数。共模电压增益越小，放大电路的性能越好。 = E ）共模抑制比Kcmr ：指差模电压放大倍数与共模电压放大倍数之比，它表明差动放大电路对共模信号的抑制能力。 =20lg| |（dB ） =| | 2、实验电路： SW1 SW-SPDT Q1 NPN Q2 NPN Q3 NPN R1 510 R2 510 R3 10k R4 10k R5 10k R6 10k R7 10k R8 5.1K R9 68K R10 36K RV1 100 R9(1) R10(2) A B C D AM FM + -
Poisson方程九点差分格式_米瑞琪
数值实验报告I 、实 验目 的 ， 内 容 1、 理解Poisson 方程九点差分格式的构造原理； 2、 理解因为网格点的不同排序方式造成的系数矩阵格式的差异; 3、 学会利用matlab 的spdiags(),kron() 函数生成系数矩阵; 、算法描述 针对一个Poisson 方程问题: Au 二厂 f (x, y) 在Poisson 方程五点差分格式的基础上，采用 Taylor 展开分析五点差分算子的截断误差，可以 得到： AhU(x,> ￥j) 二 Au(xh yj) + —Jhi 厂 h2 a/ + h/ 12 亦 yj - J - 2u(x lF Yj - 1)+ u(xi -(, Yj - 1) hi^ufxn yj - J u(x jh yj) , 0u(xi,y 」)\ + hj 扩 \ h,十 Yj) 卜一 <2—癌L 时 0/ (1) 为了提高算子截断误差的精度，在⑴式中配凑出了差分算子的形式，将原Poisson 方程代入(1) 式有： yj) = ￥」) / 水 Uhi ，" 十 0(h 4) 昇舟1 ，有： /u(x I , yj) 1 i --^luxxCxj. Yj 十 1)- 2u xx (Xir yj + Uxxfxj, Vj - Jl + 0(h 『)=— hj h? u(Kj + 1, yj + 1)- 2u(x j t yj + 1)+ u(xj - v yj + i) h,护u(x lh yj 4 1) - ------- 2 考虑 h/ u(xj + ir yj) - 2u(xj. y 」)+ 12 i- Yj) 2hiVu(x b yj + 0(h 4)二 Au(x lT V J ) + — UlXi + 1