差热传导方程的分解法c语言,热传导方程差分解法的最佳网格_李先枝

第34卷第3期2004年9月河南大学学报(自然科学版)

JournalofHenanUniversity(NaturalScience)Vol.34 No.3Sep.2004

热传导方程差分解法的最佳网格

李先枝

(郑州师范高等专科学校,河南郑州450044)

摘 要:论述了用有限差分方法求解一维热传导方程带权六点格式的最优网格选择问题,即在固定的运算次数下,对给定的权系数一,如何选择空间步长vx与时间步长vt,使得在此网格下,六点格式有更高阶的截断误差,并且保持差分解的正性。

关键词:数值解法;正性;热情导方程中图分类号:O175.7 文献标识码:A

文章编号:1003-4978(2004)03-0016-03

TheBestMeshoftheDifferenceSolution

oftheHea-tconductionEquation

LIXian-zhi

(ZhengzhouTeacher'sCollege,HenanZhengzhou450044,China)

Abstract:Thisarticlediscussesthemethodoffinitedifferencetosolvetheselectionproblemoftheoptimalmeshaboutone-dimensionalheat-conductionequation,weightwithsixformats,thatishowtochoosethespacestepvxandtimestepvtunderattheconditionofthefixedoperationdegreeandthegivenweightcoefficientHinordertomakethesixformatsownthehighertruncationerrorandkeepthedifferencesolutionpositivepropertyunderthismesh.Keywords:numericalsolution;positiveproperty;heat-conductionequation

考虑热传导方程的初边值问题

2=,xI(0,1),t>0,9x2

u(0,t)=u(1,t)=0,t\0,u(x,0)=u0(x),xI[0,1].

其中u0(x)是充分光滑的1

为了建立差分格式,引入空间步长vx=

8(h,S)=

n+1

(1)

=h,时间步长vt=SIR+的差分网格:M+1(xj,tn)Bxj=jh,tn=nS,j=0,1,,,M+1,nIN.

n

n

n

n

已知微分方程(1)的古典显格式(a)与古典隐格式(b)1

ujj

-ujuj+1-2uj+uj-1=,vx

n

n+1

n+1

n

(a)(b)

jj+1jj-1

=.vtvx2

设HI[0,1],作运算:/(a)@(1-H)+(b)@H0,便得到带权六点格式(c):

1n+11nnn1

un++un+uj+un+-ujnj+1-2ujj-11-2uj+uj-1j

=H+(1-H).vxvxn+1

(c)

收稿日期:2004-02-18

作者简介:李先枝(1967-),女,河南登封人,陕西师范大学教育硕士,郑州师范高等专科学校数学系讲师.