第四章-偏微分方程的有限差分法课件

《第四章-偏微分方程的有限差分法课件》由会员分享，可在线阅读，更多相关《第四章-偏微分方程的有限差分法课件（80页珍藏版）》请在装配图网上搜索。

1、计算物理学 Harbin Institute of Technology Yangkun 1/75第四章第四章 偏微分方程的有限差分法偏微分方程的有限差分法4.1 有限差分法原理有限差分法原理4.2 热传导方程的差分解法热传导方程的差分解法4.3 波动方程的波动方程的差分解法差分解法计算物理学 Harbin Institute of Technology Yangkun 2/754.1 有限差分法原理有限差分法原理抛物线形抛物线形双曲型双曲型椭圆形椭圆形不可逆过程不可逆过程可逆过程可逆过程平衡过程平衡过程热传导方程热传导方程波动方程波动方程位势方程位势方程 物理学中许多物理规律都用偏微分方程描

2、述，物理学中许多物理规律都用偏微分方程描述，偏微分方程主要分为以下三类：偏微分方程主要分为以下三类：fauuctudfauuctud22fauuc上式中a,c,f以及未知函数u为定义在求解区域上的实（复）函数计算物理学 Harbin Institute of Technology Yangkun 3/754.1 有限差分法原理有限差分法原理2222 fdfxdxfd fxdx有限差分有限差分解法解法差分近似代替微分，差商近似代替微商差分近似代替微分，差商近似代替微商 这样就把求解区域内连续分布函数离散化成这样就把求解区域内连续分布函数离散化成求网络节点上的分立函数值，从而把所需求解的求网络节点

3、上的分立函数值，从而把所需求解的微分方程微分方程变为一组相应的变为一组相应的差分方程差分方程，进一步可，进一步可以求解离散节点上的函数值。以求解离散节点上的函数值。数学基础数学基础泰勒（泰勒（Taylor）展开）展开计算物理学 Harbin Institute of Technology Yangkun 4/754.1 有限差分法原理有限差分法原理差商公式的构造差商公式的构造223323223323232323232323()()2!3!()()2!3!22(2)()22!3!22(2)()22!3!dfh d fh d ff xhf xhdxdxdxdfh d fh d ff xhf xhd

4、xdxdxhhdfd fd ff xhf xhdxdxdxhhdfd fd ff xhf xhdxdxdx利用泰勒级数展开定义差商利用泰勒级数展开定义差商计算物理学 Harbin Institute of Technology Yangkun 5/754.1 有限差分法原理有限差分法原理()()ff xhf xxh误差为误差为O(h)差商公式：差商公式：一阶向前差商：一阶向前差商：()()ff xf xhxh一阶向后差商：一阶向后差商：223323223323()()2!3!()()2!3!dfh d fh d ff xhf xhdxdxdxdfh d fh d ff xhf xhdxdxdx

5、计算物理学 Harbin Institute of Technology Yangkun 6/754.1 有限差分法原理有限差分法原理二阶二阶向前差商：向前差商：式式(2)-式式(1)X2222(2)2()()ff xhf xhf xxh223323232323()()(1)2!3!22(2)()2 (2)2!3!dfh d fh d ff xhf xhdxdxdxhhdfd fd ff xhf xhdxdxdx误差为误差为O(h2)差商公式：差商公式：计算物理学 Harbin Institute of Technology Yangkun 7/754.1 有限差分法原理有限差分法原理二阶二阶

6、向后差商向后差商:式式(2)-式式(1)X2222()2()(2)ff xf xhf xhxh223323232323()()(1)2!3!22(2)()2(2)2!3!dfh d fh d ff xhf xhdxdxdxhhdfd fd ff xhf xhdxdxdx计算物理学 Harbin Institute of Technology Yangkun 8/754.1 有限差分法原理有限差分法原理(2)4()3()2ff xhf xhf xxh一阶向前差商：一阶向前差商：222223323(2)2()()()()2!3!ff xhf xhf xxhdfh d fh d ff xhf xhd

7、xdxdx2222fd fxdx计算物理学 Harbin Institute of Technology Yangkun 9/754.1 有限差分法原理有限差分法原理(2)4()3()2ff xhf xhf xxh一阶向后差商：一阶向后差商：222223323()2()(2)()()2!3!ff xf xhf xhxhdfh d fh d ff xhf xhdxdxdx2222fd fxdx计算物理学 Harbin Institute of Technology Yangkun 10/754.1 有限差分法原理有限差分法原理()()2off xhf xhxh一阶中心差商：一阶中心差商：2233

8、23223323()()2!3!()()2!3!dfh d fh d ff xhf xhdxdxdxdfh d fh d ff xhf xhdxdxdx计算物理学 Harbin Institute of Technology Yangkun 11/754.1 有限差分法原理有限差分法原理222()2()()off xhf xf xhxh二阶中心差商：二阶中心差商：223323223323()()2!3!()()2!3!dfh d fh d ff xhf xhdxdxdxdfh d fh d ff xhf xhdxdxdx计算物理学 Harbin Institute of Technology

9、Yangkun 12/754.1 有限差分法原理有限差分法原理差分格式的收敛性和稳定性差分格式的收敛性和稳定性收敛性：收敛性：稳定性：稳定性：当步长当步长h0时，差分方程的解趋向于微分方时，差分方程的解趋向于微分方程的解。程的解。误差在运算过程中不会失控，即累计误差不误差在运算过程中不会失控，即累计误差不会无限增加。会无限增加。计算物理学 Harbin Institute of Technology Yangkun 13/754.1 有限差分法原理有限差分法原理 从数学上讲，没有限制的微分方程会有从数学上讲，没有限制的微分方程会有无穷多个解，不能构成一个定解问题。无穷多个解，不能构成一个定解问

10、题。从物理上讲，描述物理问题的微分方程仅从物理上讲，描述物理问题的微分方程仅适用于描述在一个连续体或物理场的内部发生适用于描述在一个连续体或物理场的内部发生的物理过程，仅靠这些微分方程不足以确定物的物理过程，仅靠这些微分方程不足以确定物理过程的理过程的具体具体特征。特征。因此，要想解决实际的物理问题，必须因此，要想解决实际的物理问题，必须知道一个连续体或物理场的初始状态和边界知道一个连续体或物理场的初始状态和边界受到的外界影响。受到的外界影响。计算物理学 Harbin Institute of Technology Yangkun 14/754.1 有限差分法原理有限差分法原理初始条件：初始条

11、件：与时间相联系与时间相联系边界条件：边界条件：边界边界受到外界的影响受到外界的影响偏微分方程的定解条件偏微分方程的定解条件常见的物理问题可以归结为三大类边界条件常见的物理问题可以归结为三大类边界条件),(),(2010zyxftuzyxfutt计算物理学 Harbin Institute of Technology Yangkun 15/754.1 有限差分法原理有限差分法原理2 第二类边界条件第二类边界条件（诺依曼（诺依曼Neumann）1 第一类边界条件第一类边界条件（狄利克雷（狄利克雷Dirichlet）),(0truu热传导热传导问题：问题：边界边界上上温度分布已知温度分布已知热传导

12、热传导问题：问题：通过通过边界边界单位面积上的热流量已知单位面积上的热流量已知),(0trqnun表示表示的外法线的外法线q0定义在定义在上的已知函数上的已知函数uuuu ni nj nnxy 计算物理学 Harbin Institute of Technology Yangkun 16/754.1 有限差分法原理有限差分法原理由热力学傅立叶定律得：由热力学傅立叶定律得：QukStn 热流量：热流量：QuqkS tn 单位面积上的单位面积上的热流量：热流量：K:热传导系数热传导系数单位时间内通过给定截面的热量，正比例于垂直于单位时间内通过给定截面的热量，正比例于垂直于该界面方向上的温度变化率和

13、截面面积，而热量传该界面方向上的温度变化率和截面面积，而热量传递的方向则与温度升高的方向相反。递的方向则与温度升高的方向相反。计算物理学 Harbin Institute of Technology Yangkun 17/754.1 有限差分法原理有限差分法原理3 第三类边界条件第三类边界条件（洛平（洛平Robin）热传导热传导问题：问题：边界表面边界表面与外界之间的热量交换已知与外界之间的热量交换已知),(000trcnubuaa0,b0.c0定义在定义在上的已知上的已知函数函数外界温度为外界温度为u0，热交换规律遵循，热交换规律遵循热传导实验定律热传导实验定律：单位时间内，从边界单位面积传

14、递给周围的热流量单位时间内，从边界单位面积传递给周围的热流量正比于边界表面和外界的温度差。正比于边界表面和外界的温度差。计算物理学 Harbin Institute of Technology Yangkun 18/754.1 有限差分法原理有限差分法原理 对于实际物理问题，边界条件往往是很复杂的，对于实际物理问题，边界条件往往是很复杂的，可能是一种或不同边界区域几种边界条件的组合，可能是一种或不同边界区域几种边界条件的组合，甚至不能用这三类边界条件描述。甚至不能用这三类边界条件描述。0quu：热交换系数：热交换系数 u：边界温度：边界温度QuqkS tn 单位面积上的单位面积上的热流量：热流

15、量：0ukuun0uukun计算物理学 Harbin Institute of Technology Yangkun 19/754.2 热传导方程的差分解法热传导方程的差分解法 物理学中对热传导、热辐射以及气体扩散现物理学中对热传导、热辐射以及气体扩散现象的描述，常可以归结为同一类型的抛物线型方象的描述，常可以归结为同一类型的抛物线型方程，通常采用二阶偏微分方程描述，这类方程统程，通常采用二阶偏微分方程描述，这类方程统称为热传导方程。称为热传导方程。fauuctud计算物理学 Harbin Institute of Technology Yangkun 20/754.2 热传导方程的差分解法热

16、传导方程的差分解法一维各向同性、均匀介质，且无热源的热传导方程：一维各向同性、均匀介质，且无热源的热传导方程：220 0tTuuxltx 为了求解为了求解u(x,t)，还必须利用边界条件和初，还必须利用边界条件和初始条件。始条件。定解条件定解条件：边界条件边界条件和和初始条件初始条件。定解问题定解问题：解存在、唯一并且连续依赖初始条件。：解存在、唯一并且连续依赖初始条件。4.2.1一维热传导方程的差分解法一维热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 21/754.2 热传导方程的差分解法热传导方程的差分解法 22210,0ut

17、gtu l tuutxxtugx对于一维热传导问题对于一维热传导问题（第一类边界条件）（第一类边界条件）00 xtlT 数值解就是在求解区域数值解就是在求解区域:0 0tGxTl 中某些中某些离散点（离散点（xi,ti）上求出上求出u(xi,ti)足够近似的解。足够近似的解。计算物理学 Harbin Institute of Technology Yangkun 22/754.2 热传导方程的差分解法热传导方程的差分解法1 把求解区域离散化（确定离散点）把求解区域离散化（确定离散点）00 i=0,1,Nk=0,1,MkilhxxihihTkNttkM,(,),ikiki kx ti ku x

18、tuxtxitkTl,i ku计算物理学 Harbin Institute of Technology Yangkun 23/754.2 热传导方程的差分解法热传导方程的差分解法2 推导差分递推公式推导差分递推公式222,2,(,)u xh tu x tu xh tu x txh在节点（在节点（xi,tk）上）上2221,1,2,2,(,)2ikikikikx xt tiki kiku xh tu x tu xh tu x txhuuuh二阶二阶向前差商向前差商O(h2)计算物理学 Harbin Institute of Technology Yangkun 24/754.2 热传导方程的差分

19、解法热传导方程的差分解法同样，在节点（同样，在节点（xi,tk）上）上,1,(,)ikikikx xt ti ki ku x tu x tx ttuu一一阶阶向前差商向前差商O(h)计算物理学 Harbin Institute of Technology Yangkun 25/754.2 热传导方程的差分解法热传导方程的差分解法一维热传导方程可以近似为一维热传导方程可以近似为,1,1,1,22i ki kiki kikuuuuuh令令2h,11,1,(1 2)i kiki kikuuuuO(h)计算物理学 Harbin Institute of Technology Yangkun 26/75

20、4.2 热传导方程的差分解法热传导方程的差分解法,11,1,(1 2)i kiki kikuuuuxtxitk,1i ku初始条件初始条件边界条件边界条件计算物理学 Harbin Institute of Technology Yangkun 27/754.2 热传导方程的差分解法热传导方程的差分解法一维热传导方程显示差分递推公式为一维热传导方程显示差分递推公式为,11,1,00,1,2(12)()()()i kiki kikikl kuuuuuihug kugki=0,1,Nk=0,1,M计算物理学 Harbin Institute of Technology Yangkun 28/754.

21、2 热传导方程的差分解法热传导方程的差分解法显示差分递推公式的稳定性：显示差分递推公式的稳定性：212h对于一维热传导方程，差分格式为稳定差对于一维热传导方程，差分格式为稳定差分格式的充分条件是分格式的充分条件是即即22hki,kkikikiuu,计算物理学 Harbin Institute of Technology Yangkun 29/754.2 热传导方程的差分解法热传导方程的差分解法为了提高数值解的精度，必须减小为了提高数值解的精度，必须减小22hhxdx 相应就要变小，这必然增加计算量。这就相应就要变小，这必然增加计算量。这就是显示差分格式的缺点，但它的优点是计算是显示差分格式的缺

22、点，但它的优点是计算简单。简单。计算物理学 Harbin Institute of Technology Yangkun 30/754.2 热传导方程的差分解法热传导方程的差分解法差分格式计算步骤：差分格式计算步骤：给定给定 计算计算 计算初始值：计算初始值：计算边界值：计算边界值：用差分格式计算用差分格式计算,l T h2;2h,0()iuih0,1,2(),()kl kug kugk,1i ku2;h,11,1,00,1,2(1 2)()()()i kiki kikikl kuuuuuihug kugk22 00uutt Txx l 计算物理学 Harbin Institute of Te

23、chnology Yangkun 31/754.2 热传导方程的差分解法热传导方程的差分解法,11,1,(1 2)i kiki kikuuuui=0,1,Nk=0,1,M1,12,1,0,2,13,2,1,3,14,3,2,N 1,1,N 1,2,(1 2)(1 2)(1 2)(1 2)kkkkkkkkkkkkkN kkNkuuuuuuuuuuuuuuuu计算物理学 Harbin Institute of Technology Yangkun 32/754.2 热传导方程的差分解法热传导方程的差分解法1,10,1,2,12,3,13,1,1N,1,1 21 21 21 2kkkkkkkNkkN

24、kuuuuuuuuuu1,12,1,0,2,13,2,1,3,14,3,2,N 1,1,N 1,2,(1 2)(1 2)(1 2)(1 2)kkkkkkkkkkkkkN kkNkuuuuuuuuuuuuuuuu计算物理学 Harbin Institute of Technology Yangkun 33/754.2 热传导方程的差分解法热传导方程的差分解法X=diag(v,k)若若v为为n个元素向量，返回一个阶数为个元素向量，返回一个阶数为n+abs(k)的方阵的方阵X，将，将v作为作为X的第的第k个对角个对角元，元，k=0代表主对角元，代表主对角元，k0表示在主对角表示在主对角元之上，元之上

25、，k0表示在主对角元以下。表示在主对角元以下。v=ones(1,5);X1=diag(v)X2=diag(v,1)X3=diag(v,-1)计算物理学 Harbin Institute of Technology Yangkun 34/754.2 热传导方程的差分解法热传导方程的差分解法（1-2*）*diag(ones(1,N-1)+*（diag(ones(1,N-2),1)+diag(ones(1,N-2),-1)）1,10,1,2,12,3,13,1,1N,1,1 21 21 21 2kkkkkkkNkkNkuuuuuuuuuu计算物理学 Harbin Institute of Techn

26、ology Yangkun 35/754.2 热传导方程的差分解法热传导方程的差分解法例例 求热传导方程混合问题求热传导方程混合问题22,01,0(,0)4(1),01(0,)0,(1,)0,0uuxttxu xxxxututt计算物理学 Harbin Institute of Technology Yangkun 36/754.2 热传导方程的差分解法热传导方程的差分解法221 0.05 4hhh,11,1,(1 2)i kiki kikuuuu1,11,00,1,2,.1 1,2,.04(1)ikm kihiuuimhknu22h计算物理学 Harbin Institute of Tech

27、nology Yangkun 37/754.2 热传导方程的差分解法热传导方程的差分解法,14(1)iihuih,11,1,1,1,1(1 2)(1 2)i kiki kikkik ik ik iuuuuuuuu计算物理学 Harbin Institute of Technology Yangkun 38/759.2 4.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 39/754.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangk

28、un 40/754.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 41/754.2 热传导方程的差分解法热传导方程的差分解法二维各向同性、均匀介质，且无热源的热传导方程：二维各向同性、均匀介质，且无热源的热传导方程：2222uuutxy000 xlystT,0,u x yx y初始条件：初始条件：边界条件：？边界条件：？4.2.2 二二维热传导方程的差分解法维热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 42/754.2 热传导方程的差分解法热

29、传导方程的差分解法同一维类似，把求解区域离散化同一维类似，把求解区域离散化 ijklhNShxihyjtkMhxyxiyj,ijki j ku x y tu,i j kuktt 计算物理学 Harbin Institute of Technology Yangkun 43/754.2 热传导方程的差分解法热传导方程的差分解法在节点（在节点（xi,yj,tk）上）上2221,1,2222,1,1,2,2,(,)2,2,(,)2ikikijkijkijkx xt tij ki j kij kijkijkijkx xt ti jki j ki jku xh y tu x y tu xh y tu x

30、 y txhuuuhu x yh tu x y tu x yh tu x y tyhuuuh计算物理学 Harbin Institute of Technology Yangkun 44/754.2 热传导方程的差分解法热传导方程的差分解法在节点（在节点（xi,yj,tk）上）上,1,(,)ikijkijkx xt ti j ki j ku x y tu x y tx y ttuu计算物理学 Harbin Institute of Technology Yangkun 45/754.2 热传导方程的差分解法热传导方程的差分解法二维热传导方程可以近似为二维热传导方程可以近似为,1,1,1,1,1

31、,222i j ki j kij ki kij ki jki ki jkuuuuuuuuh令令2h计算物理学 Harbin Institute of Technology Yangkun 46/75 t k+1 j+1 j k j-1 i-1 i i+1图(9.3)二 维 热 传 导 差 分 解 法4.2 热传导方程的差分解法热传导方程的差分解法,1,1,1,1,1,(1 4)i j ki j kij kij ki jki jkuuuuuu差分递推公式为差分递推公式为,1i j ku1,ij ku1,ij ku,i j ku,1,i jku,1,i jku计算物理学 Harbin Instit

32、ute of Technology Yangkun 47/754.2 热传导方程的差分解法热传导方程的差分解法边界条件：边界条件：一、右图中阴影部分一、右图中阴影部分为绝热壁，此绝热壁为绝热壁，此绝热壁可以用第二类边界条可以用第二类边界条件描述，即热流量为件描述，即热流量为零。零。第二类边界条件：第二类边界条件：通过通过边界表面单位面积上的热流量已知边界表面单位面积上的热流量已知),(0trqnuuuuu ni nj nnxy 计算物理学 Harbin Institute of Technology Yangkun 48/754.2 热传导方程的差分解法热传导方程的差分解法差分近似为差分近似为

33、1,0,1,00j kj kN j kNj kuuhuuh12,0,11,2,1,1,10 0 1,2,1i j kii j ki NuxuxjMMMjM计算物理学 Harbin Institute of Technology Yangkun 49/754.2 热传导方程的差分解法热传导方程的差分解法递推公式为：递推公式为：120,1,1,1,2,1,1,11,2,1 j kj kN j kNj kjMMMuuuujM计算物理学 Harbin Institute of Technology Yangkun 50/754.2 热传导方程的差分解法热传导方程的差分解法二、二、i=0边界，边界，M1

34、jM2 区区域为与高温恒温热源相连域为与高温恒温热源相连接的口，温度可取归一化接的口，温度可取归一化值值1。j=0和和j=M边界与低温恒温边界与低温恒温热源相连，温度始终为热源相连，温度始终为0。计算物理学 Harbin Institute of Technology Yangkun 51/754.2 热传导方程的差分解法热传导方程的差分解法二维热传导方程显示差分递推公式为二维热传导方程显示差分递推公式为计算物理学 Harbin Institute of Technology Yangkun 52/754.2 热传导方程的差分解法热传导方程的差分解法214h稳定差分格式的充分条件是稳定差分格式

35、的充分条件是即即24h计算物理学 Harbin Institute of Technology Yangkun 53/754.2 热传导方程的差分解法热传导方程的差分解法,1,1,1,1,1,(1 4)i j ki j kij kij ki jki jkuuuuuu计算物理学 Harbin Institute of Technology Yangkun 54/754.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 55/754.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Instit

36、ute of Technology Yangkun 56/754.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 57/754.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 58/754.2 热传导方程的差分解法热传导方程的差分解法 view(az,el)az:方位角方位角el:仰角仰角view(0,90)view(-37.5,30）计算物理学 Harbin Institute of Technology Yangkun

37、59/754.2 热传导方程的差分解法热传导方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 60/75 4.3 波动方程的差分解法波动方程的差分解法一维均匀弦线的自由振动方程：一维均匀弦线的自由振动方程：22222vtx00 xltT 波动方程的差分解法也利用构造网格节点的方法波动方程的差分解法也利用构造网格节点的方法NlihNiihxx,1,0,MTkMkktt,1,0,(,)i kikx t计算物理学 Harbin Institute of Technology Yangkun 61/75 4.3 波动方程的差分解法波动方程的差分解法

38、用二阶中心差分近似方法得：用二阶中心差分近似方法得：221,1,22222,1,122222iki kiki ki ki kxxhtt计算物理学 Harbin Institute of Technology Yangkun 62/75 4.3 波动方程的差分解法波动方程的差分解法代入波动方程得：代入波动方程得：,1,11,1,22222i ki ki kiki kikvh令：令：vh222,1,1,1,12(1)i ki kikiki k 计算物理学 Harbin Institute of Technology Yangkun 63/75 4.3 波动方程的差分解法波动方程的差分解法222,1

39、,1,1,12(1)i ki kikiki k t M k+1 k k-1 x i-1 i i+1 N图(10.1)一 维 波 动 方 程 的 差 分 解 法边界条件边界条件初始条件初始条件计算物理学 Harbin Institute of Technology Yangkun 64/75 4.3 波动方程的差分解法波动方程的差分解法a 初始条件：初始条件：12(,0)()(,0)()xxxxt对于初始时刻速度，也须用差分格式给出：对于初始时刻速度，也须用差分格式给出：0 xl计算物理学 Harbin Institute of Technology Yangkun 65/75 4.3 波动方程

40、的差分解法波动方程的差分解法a1 向前差分：向前差分：,0,0,1,02,1,0212()()()()iiiiiiiiiixttxxx,01,112()()()iiiiixxx误差：误差：()O h计算物理学 Harbin Institute of Technology Yangkun 66/75 4.3 波动方程的差分解法波动方程的差分解法a2 中心差分：中心差分：,0,0,1,1,12,12()22()iiiiiiiixttx由由 222,1,1,1,12(1)i ki kikiki k 得（得（k=0）222,1,01,01,0,12(1)iiiii +计算物理学 Harbin Inst

41、itute of Technology Yangkun 67/75 4.3 波动方程的差分解法波动方程的差分解法22,1,01,01,0222,1111112(1)()()2(1)()()()()2iiiiiiiiiixxxxx,0122,1111112()(1)()()()()2iiiiiiixxxxx误差：误差：2()O h整理得：整理得：计算物理学 Harbin Institute of Technology Yangkun 68/75 4.3 波动方程的差分解法波动方程的差分解法b 边界条件：边界条件：12(0,)()(,)()tg tl tg t一维波动方程的差分格式有如下两种一维波

42、动方程的差分格式有如下两种0tT 计算物理学 Harbin Institute of Technology Yangkun 69/75 4.3 波动方程的差分解法波动方程的差分解法第一种：第一种：22,1,1,1,1,01,1120,1,22(1)()()()()()i ki kikiki kiikN kihihihg kg k误差：误差：()O h计算物理学 Harbin Institute of Technology Yangkun 70/75 4.3 波动方程的差分解法波动方程的差分解法第二种：第二种：22,1,1,1,1,0122,111120,1,22(1)()(1)()(1)(1)

43、()2()()i ki kikiki kiikN kihihihihihg kgk误差：误差：2()O h计算物理学 Harbin Institute of Technology Yangkun 71/75 4.3 波动方程的差分解法波动方程的差分解法 第二种差分格式精度要高于第一种，是经常第二种差分格式精度要高于第一种，是经常采用的方法。采用的方法。1vh 理论上可以证明，两种差分格式稳定条件是：理论上可以证明，两种差分格式稳定条件是：计算物理学 Harbin Institute of Technology Yangkun 72/75 4.3 波动方程的差分解法波动方程的差分解法,vhlTh

44、v,ikxihtk,1i k波动方程差分格式的计算步骤如下：波动方程差分格式的计算步骤如下：;计算计算;计算计算;计算初值和边值计算初值和边值;给定给定计算计算计算物理学 Harbin Institute of Technology Yangkun 73/75 4.3 波动方程的差分解法波动方程的差分解法2222,01,0(,0)(,0)sin,(1),01(0,)(1,)0,0yyxttxy xy xxxxxtytytt例：计算下列一维波动方程例：计算下列一维波动方程计算物理学 Harbin Institute of Technology Yangkun 74/75 4.3 波动方程的差分解

45、法波动方程的差分解法22,1,1,1,11,1,22(1)()00sin(1)sin(1)(1)(1(1)i ki kikiki kkm kiiyyyyyyyyihyihihih11,0.05,2hvhv计算物理学 Harbin Institute of Technology Yangkun 75/75 4.3 波动方程的差分解法波动方程的差分解法22,1,1,1,11,1,22(1)()00sin(1)sin(1)(1)(1(1)i ki kikiki kkm kiiyyyyyyyyihyihihih计算物理学 Harbin Institute of Technology Yangkun 7

46、6/75 4.3 波动方程的差分解法波动方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 77/75 4.3 波动方程的差分解法波动方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 78/75 4.3 波动方程的差分解法波动方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 79/75 4.3 波动方程的差分解法波动方程的差分解法计算物理学 Harbin Institute of Technology Yangkun 80/75 上机上机 422010.1 01 010(,0)(,)0(,)0 xxxuuxttxu xeuu x tnuu x tn 编程计算下列一维扩散方程的解编程计算下列一维扩散方程的解要求：要求：1 推导差分递推公式推导差分递推公式2 图形显示计算结果，并与解析解比较。图形显示计算结果，并与解析解比较。a x-t-u 三维曲线三维曲线 b t=0 0.5 1 时刻时刻 x-u曲线曲线 c x=0.5处处t-u曲线曲线exp(0.1)uxt