一类线性偏积分微分方程的数值解.docxVIP

一类线性偏积分微分方程的数值解 在这项工作中，我们研究了线性微分方程的相关分析。 此问题常出现在带有粘弹性流体模型及带有记忆功能的热传导物质,因此讨论这类偏积分微分方程的求解具有非常重要的意义.但这类方程中的大部分方程的数学精确解不易求出,所以有必要研究它们的数值解.研究数值解的传统方法有差分法、有限元法、谱方法等,本文将在时间t方向采用一种新的方法——Lubich的拉普拉斯变换数值逆. 1 构造所遵循的拉普拉斯逆变换 给出网格t=0,h,2h,…,Nh,卷积 可以离散为 其中wj(h)由幂级数 给出,这里F是f的拉普拉斯变换,δ(ξ)=∞∑j=0δjξj是生成线性多步法多项式的商数. 由罗朗定理有 其中Γρ:|ξ|=ρ,ρ0为常数. 设F(s)在区域|arg(s-c)|π-φ,φπ2,c∈R内解析,并且满足 |F(s)|≤Μ?|s|-μ,Μ∞,μ0,(1.5) δ(ξ)满足条件: (1)δ(ξ)在闭单位圆|ξ|≤1的某个邻域内解析且没有零点,除在ξ=1的一个零点外, (1.6a) (2)|argδ(ξ)|≤π-α,|ξ|1,αφ, (1.6b) (3)1hδ(e-h)1+o(hp),当p≥1. (1.6c) 当α=90°,90°,88°,73°,51°,18°时,相应的p=1,2,…,6,这时δ(ξ)由δ(ξ)=p∑i=1(1-ξ)i给出. 我们定义fj(h),由 F(δ(ξ)/h)=h∞∑j=0fj(h)ξj,(1.7) 得 fj(h)=wj(h)/h.(1.8) wj(h)由式(1.4)给出,可以证明fj(h)逼近F(s)的拉普拉斯逆变换f(t). 定理1 假定式(1.5),(1.6)成立,且设f(t)是F(s)的拉普拉斯逆变换.则 |fj(h)-f(jh)|≤C?tμ-1-p?hp,(t=jh),(1.9) 其中常数C与h∈(0,ˉh]和t∈[h,ˉt]无关,且ˉt+∞. 定理的证明见文献. 2 utx,tv1/2x-t 例1 解方程 {ut(x,t)-∫t0(t-s)-1/2uxx(x,s)ds=f(x,t),u(0,t)=u(1,t)=0,0≤t1,(Ι′) u(x,0)=w(x),0≤x≤1. 为了获得精确解 u(x,t)=x(1-x)(t3/2+1), 我们取 w(x)=x(1-x), 和 f(x,t)=(2+32x-32x2)t1/2, 确定试探函数空间Η01(0,1)={v(x)|v(x),v′(x)∈L2(0,1),v(0)=v(1)=0}. 将式(I′)的第一个方程两边乘以一个试探函数v(x),可得 ut(x,t)v(x)-∫0t(t-s)-1/2uxx(x,s)v(x)ds=f(x,t)v(x),(2.1) 两边取积分并整理得 ∫01ut(x,t)v(x)dx+∫01∫0t(t-s)-1/2ux(x,s)v′(x)dsdx=∫01f(x,t)v(x)dx.(2.2) 2.1 l-1维线性运动中各.1基函数的线性运动 在上引进等距节点0=x0x1…xL=1,记步长λ=1L,xn=nλ,0nL-1?[0?1]上的分段线性函数φλ在上连续且φλ(0)=φλ(1)=0,且在每个[xn-1,xn]上是线性函数(n=1,2,…,L-1),所有满足这些条件的函数构成集合Uλ,即Uλ={φλ|φλ∈C[0?1],且在每个[xn-1,xn]上是线性函数,φλ(0)=φλ(1)=0},容易验证Uλ是实数域上的一个L-1维的线性空间,Uλ的基函数为 φn(x)={0,0≤xn-1L,Lx-n+1,n-1L≤xnL,n+1-Lx,nL≤xn+1L,0,?n+1L≤x≤1,(n=1,2,?,L-1) 则Uλ中的任意函数都可以表示为基函数的线性组合.设 u(x,s)≈Uλ(x,s)=a1(s)φλ1(x)+a2(s)φλ2(x)+?+aL-1(s)φλL-1(x),(2.3)u(x,t)≈Uλ(x,t)=a1(t)φλ1(x)+a2(t)φλ2(x)+?+aL-1(t)φλL-1(x),(2.4) 于是 ux(x,s)=a1(s)(φλ1(x))′+a2(s)(φλ2(x))′+?+aL-1(s)(φλL-1(x))′,(2.5)ut(x,t)=a′1(t)φλ1(x)+a′2(t)φλ2(x)+?+a′L-1(t)φλL-1(x),(2.6) 将式(2.3)～(2.6)代入式(2.2)并取v(x)=φk(x)可得 ∫01∑n=1L-1a′n(t)φλn(x)φk(x)dx+∫0t[(t-s)-1/2∫01∑n=1L-1an(x)(φλn(x))′φ′k(x)dx]ds=∫01f(x,t)φk(x)dx,(k=1,2,?,L-1),(2.7) 其中 ∫01f(x,t)φk(x)dx=∫01(2+3