资源描述:
《试论《有限差分法在微分方程中的应用》课程论文》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、试论《有限差分法在微分方程中的应用》课程论文试论《有限差分法在微分方程中的应用》课程论文_论文范文导读:课程论文有限差分法在微分方程中的应用本学期学习了《微分方程数值解》,本书中有限差分法给我留下的印象比较深刻,下边说说自己在方面的一点理解,请老师指正。1.有限差分法的基本思想:当系统的数学模型建立后,我们面对的主要问题就是微分积分方程的求解。基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件,并把差分方程组的解作为微分方程定解问题的近似解。将原方程及边界条件中的微分用差分来近似,对于方程
2、中的积分用求和或及机械求积公式来近似代替,从而把原微分积分方程和边界条件转化成差分方程组。2.有限差分法求解偏微分方程的步骤:区域离散,即把所给偏微分方程的求解区域细分成由有限个格点组成的X格,这些离散点称作X格的节点;近似替代,即采用有限差分公式替代每一个格点的导数。逼近求解,换而言之,这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程。从原则上说,这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格,或者通过离散点上的函数值进行插值计算来近似得到。理论上,当X格
3、步长趋近于零时,差分方程组的解应该收敛于精确解,但由于机器字节的限制,X格步长不可能也没有必要取得无限小,那么差分法的收敛性或者说算法的稳定性就显得至关重要。因此,在运用有限差分法时,除了要保证精度外,还必须要保证其收敛性。3.构造差分法的几种形式:主要草用的是泰勒级数展开的方法。其基本差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。其中前两种形式为一阶计算精度,后一种为二阶计算精度。4.有限差分法的应用:4.1抛物线形的差分法中的一维常系数抛物线型方程考虑最简单的以为常系数抛物线型方程u
4、uLua2f(x)tx2(x,t)其中是(x.t)平面内的给定区域,可以是有节区域或无解区域;a>0是常数,L是微分算子。根据定解条件的不同,可以将上述方程分为两类:1.初值问题在区域{(x,t)
5、x,t0}上求解方程满足初始条件u(x,0)(x),x的解。2,初边值问题(混合问题)在区域{(x,t)
6、0xl,0xT}内求方程满足初始条件u(x,0)(x),0xl和下列边界条件之一的解。第一边届条件u(0,t)1,u(l,t)2第二边界条件0tTux(0,t)1,ux(l,t)2第三边界条件0tT(u1(t)ux)
7、x
8、0r1(t)(u2(t试论《有限差分法在微分方程中的应用》课程论文_论文范文(2)导读:f1(a12hb1)u0h2f2g)ux)
9、x0r2(t)i0,i1,2,0tT用适当的差商代替方程中相应的偏导数,可得到以下几种最简差分格式:古典显示格式:Lu(1)hkj1ukukjjakkukj12ujuj1h2fjk古典隐式格式:k1L(2)huj1ukukjja1k11ukukj12ujj1h2fjk1加权六点隐式格式:1ukukjj1k11kkukukukj12ujj1j12ujuj1k1kaa(1)fj(1)fj22hhLu
10、(3)hkj4.2椭圆型方程边值问题的差分法考虑如下两点边值问题aubucufu(0)u0,u(1)u1x(0,1)其中u0,u1为常数,系数a=a(x),b=b(x),c=c(x),f=f(x)为一致的充分光滑函数,且满足a(x)>0,c(x)>0.首先将区间离散化,我们采用剖分部分,取正整数M,将区间M等分,的M+1个节点:0x0x1xM1其中xjjh(j0,1,M),h1。设U为定义在节点xjjh(j0,1,M)上的X格函M数并用Ui近似u(xj)。下边可以得到两点边值问题的有限差分:AhUjajUj12U
11、jUj1h2bjUj1Ujh1cjUjfjU0u0,U1u1j0,1,M上述方程组成为差分方程,它的解就是两点边值问题的差分解,当j0,1,M1时忧伤出差分方程可得:11(ajhbj)Uj1(2ajh2cj)Uj(ajhbj)Uj1h2fj22由此得出线性代数方程:AUg(M1)(M1)其中U(U1,U2,UM1,),g(g1,g2,gM1),因此ARTT122ahb(ahb1)1112(a1hb)2a2hb2(a2hb2)2222A(aij)1(a3hb3)2a3h2b322aM1h2bM1gRM1定义为12hf1(a12
12、hb1)u0h2f2g试论《有限差分法在微分方程中的应用》课程论文_论文范文(3)导读:只列出结果。1。Courant-Isaacson-Rees格式(迎风格式)1ukukjjkukjuj1kukj1ujahah10其中a(
13、a
14、a)22.Lax-Friedrichs格式。11kkuk(u