二维热传导方程隐式差分格式系数矩阵特征值的求取.docVIP

大庆石油学院学报 JOURNALOFDAQINGPETROLEUMINSTITUTE第29卷Vol.29第3期No.32005年6月Jun.2005 二维热传导方程隐式差分格式 系数矩阵特征值的求取 冯凤萍1,周瑞芬2 (1.大庆石油学院数学系,黑龙江大庆　163318;　2.大庆石油学院机械科学与工程学院,黑龙江大庆　163318)　　摘　要:给出了二维热传导方程隐式差分格式.对该差分格式的未知量和右端项以矩阵形式表示,的矩阵表示新形式.,该类线性方程组时的收敛速度分析奠定了基础. 关　键　词:隐式差分格式;矩阵特征值;稳定性中图分类号:O177　　　:10000　引言 .对高维方程差分格式系数矩阵特征值的求取是比较困难的.文献[1,2]给出了求矩阵特征值的方法,但其结果只求得矩阵特征值的近似值.文献[3]给出求Poisson方程第一边值问题差分格式矩阵的特征值的求取方法.基于文献[3],笔者给出了二维热传导方程第一边值混合问题隐式差分格式的矩阵特征值的求取方法. 1　二维热传导方程及其隐式差分格式 1.1　第一边值混合问题 222+=a,(x,y)∈(0,1)×(0,1),t0,x2y2t u(0,y,t)=1(y,t),u(1,y,t)=2(y,t), u(x,0,t)=ψ1(x,t),u(x,1,t)=ψ2(x,t),(1) u(x,y,0)=g(x,y). 式中:u(x,y,t)为待求函数;1(y,t),2(y,t),ψ1(x,t),ψ2(x,t),g(x,y)均为已知函数;a为已知常数.该问题是利用已知边界条件和初始条件求u(x,y,t). 1.2　隐式差分格式的构造 把正方形的每边n等份,设Δx,Δy,Δt分别为X-Y方向网格步长和时间网格步长,并约定记号mΔx,jΔy,mΔt),Δ2xuim,j=uim+1,j-2uim,j+uim-1,j,Δ2yuim,j=uim,j+1-2uim,j+uim,j-1.ui,j=u(i Δx,jΔy ,m+Δt)处用差分逼近微分得首先使式(1)的第1个方程在点(i2 m+1m22Δ2Δ2Δ2Δ2mm+1mm+1(2)=+.2ui,j+2ui,j2ui,j+2ui,j(Δx)(Δy)Δt2(Δx)2(Δy) 此差分格式局部截断误差为O((Δx)2+(Δy)2+(Δt)2),且为隐格式. 收稿日期审稿人:张长海;编辑:关开澄 　　基金项目:黑龙江省自然科学基金项目(A0019) 　　作者简介:冯凤萍(1965-),女,讲师,主要从事应用数学方面的研究. ?88? ? 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 第3期　　　　　　　　　　　　　冯凤萍等:二维热传导方程隐式差分格式系数矩阵特征值的求取 取Δx=Δy,令r=,将式(2)化简得 2(Δx)2 1 (1-rΔ2x-rΔ2y)uim+Δ2x+rΔ2y)uim,j.=(1+r,j (3) mmmmmm 　　将边界条件和初始条件离散后,并令um=(u1,1,u1,2,…,u1,n-1,…,un-1,1,un-1,2,…,un-1,n-1),则差分格式的求解成为求解一个线性代数方程组 Au m+1 =Bu m +f,(4) 式中:A为系数矩阵;矩阵B、向量f和向量um均为已知值.通过线性代数方程组(4)即可利用第m层的待求函数值求得第m+1层的函数值. 2　隐式差分格式系数矩阵特征值的求取 由于是二维差分格式,系数矩阵A,常困难.令U=[ui,j],F=[fi,j],则线性代数方程组(4)(5 )Tn-1U+U其中, -2 -r2r+-r (6)Tn-1=. ωωω ωω-r -r2r+ 此时,系数矩阵 Tn-1+2In-1-In-1 A= -In-1 ωω ωω -In-1 -In-1Tn-1+ 2In-. 　　Tn-1是特殊的三对角矩阵,关于它们的特征值已有结果[2],即Tn-1的特征值 λ+0.5,j=1,2,…,n-1;j=21-cos n 　　相应的特征向量 zj= sin,sin,…,sin nnN T ,j=1,2,…,n-1. 　　利用Tn-1的特征值和特征向量可求出矩阵A的特征值和特征向量. T 定理1　矩阵A的特征值为ξpq=λzq按列拉直的列向量,其中p=1.p+λq,特征向量vpq是矩阵zp?2,…,n-1;q=1,2,…,n-1. 证明　因为 TTT λzq+λzq=ξzq,pzp?qzp?pqzp? (λzq+zp?(λzq,pzp)?qzq)=ξpqzp?(Tn-1zp)?zq+zp?(λqzq) T T T T T