ux ( x , s) v∗( x ) d s d x =
1 0
f ( x , t ) v ( x ) dx .
( 2. 2)
2. 1
x 方向 : 线性有限元离散 在[ 0, 1] 上引进等距节点 0 = x 0 < x 1 < ∀ < x L = 1, 记步长 ∋=
1 , x n = n∋ , 0 < n < L - 1, [ 0, L 1] 上的分段线性函数 ∃ ( 0) = ∃ ( 1) = 0, 且在每个[ x n- 1 , x n ] 上是线性函数( n = ∋ 在 [ 0, 1] 上连续且 ∃ ∋ ∋
黎丽梅
( 湖南理工学院 数学系, 湖南 岳阳
t
414006 )
摘要 : 给出一种求一类线性偏积分微分 方程 u t ( x , t ) -
0
( t - s) u xx ( x , s) d s = f ( x , t ) 数值解的方法 , 空间 x
方向 采用线性有限元离散 , 时间 t 方向采用 Lubich 的拉普拉斯变换数值逆 , 得出数值解 的精度较高 , 计算也比较 简便 . 关键词 : 偏微分方程 ; 有限元 ; 拉普拉斯变换 ; 数 中图分类号 : O241. 8 文献标识码 : A
( 2. 7)
1 0
f ( x , t) ∃ k ( x ) dx =
2 1/ 2 (2+ 3 x - 3 x ) t ∃ k( x ) dx . 2 2
2. 2
t 方向: 取拉氏变换 记 an ( t ) 的拉氏变换为 (n ( s ) , 即 (n ( s) =
+ % 0
an ( t ) e
- st
( 2. 3) ( 2. 4)
u ( x , t ) + U ∋( x , t ) = a 1 ( t ) ∃∋( x ) + a 2 ( t ) ∃∋( x ) + ∀ + aL- 1 ( t ) ∃∋ ( x ) , 于是 ux ( x , s) = a 1 ( s ) ( ∃∋( x ) )∗ + a 2 ( s ) ( ∃∋( x ) ) ∗ + ∀ + aL- 1 ( s ) ( ∃∋ ( x ) ) ∗, ut ( x , t ) = a∗ ( t) 1
% -1 - p
j= 0 %
∃ ( 1-
p
) 给出 [ 1] .
i
∃ f j ( h)
j
,
( 1. 7)
( 1. 8)
&h ,
p
( t = j h) ,
( 1. 9)
其中常数 C 与 h ( ( 0, h] 和 t ( [ h , t ] 无关, 且 t < + % .
2 数值例子
t
例 1 解方程
ut ( x , t ) -
Lx - k + 1) d x +
2
( k + 1 - L x ) d x + x k+ 1 ( 1 + x k+ 1 )
( k + 1 - Lx ) ( Lx - k ) d x ,
则有方程组
204
北华大学学报 ( 自然科学版 )
第8卷
a d 0 ) 0 0
b a ∀ 0
0 d a ! 0 ∀
j= 0
∃ wj( h)
j
( 1. 3)
j= 0
∃
%
j j
是生成线性多步法多项式的商数 [ 1 3] .
∋
∀#
F( ( ) / h)
- j- 1
d ,
( 1. 4)
收稿日期 : 2007 01 15 基金项目 : 国家自然科学基金资助项目 ( 10271046) 作者简介 : 黎丽梅 ( 1974- ) , 女 , 讲师 , 硕士 , 主要从事计算数学研究 .
1 ∃∋( 1 2 L- 1
( 2. 5) ( 2. 6)
x ) + a∗ 2( t )
2 ∃∋(
x ) + ∀ + a∗ L- 1 ( t )
L- 1 ∃∋ (
x ),
将式 ( 2. 3) ~ ( 2. 6) 代入式 ( 2. 2) 并取 v ( x ) = ∃ k ( x ) 可得
1 L- 1
∃ 0
= d,
2 ( k + 1 - L x ) ( 2 + 3 x - 3 x ) dx + 2 2 k L
x k- 1 ( 1 - x k- 1 )
k+ 1 L k L
k- 1 ( L 2
k - L x ) ( Lx - k + 1) d x + x k ( 1 - x k )
k+ 1 L k L
k- 1 ( L
∀ 0 d ! d 0
0 ∀ 0 ! a d
0 0 ) d a
(1 ( s) (2 ( s) (3 ( s) ) (L- 2 (L- 1 =
g1 g2 g3 ) g L- 2 g L- 1
成立 , 从而可以解出 (1 ( s) , (2 ( s ) , ∀, (L- 1 ( s) .
3 取 ( ) = 1-
本文研究一类线性偏积分微分方程 :
t
ut ( x , t ) -
0
( t - s)
- 1/ 2
u xx ( x , s ) d s = f ( x , t ) , 0 x t 1. T, ( I)
u( 0 , t ) = u ( 1, t ) = 0, u( x , 0 ) = v ( x ) , 0
0
( t - s)
- 1/ 2
u xx ( x , s ) d s = f ( x , t ) , 0 t < 1, 0 x 1.
( I∗)
u( 0 , t ) = u ( 1, t ) = 0,
u ( x , 0) = w ( x ) , 为了获得精确解
u ( x , t) = x ( 1 - x ) ( t 我们取
t - 1/ 2
ut ( x , t ) v ( x ) 两边取积分并整理得
1 0
0
( t - s)
uxx ( x , s ) v ( x ) ds = f ( x , t ) v ( x ) ,
( 2. 1)
1 t
ut ( x , t ) v ( x ) dx +
0 0
( t - s)
- 1/ 2
∃ u n ( h)
j
%
j
,
( 3. 1) ( 3. 2)
u n ( h) =
w n ( h) , h
j j
u n ( jh ) = u n ( tj ) + u n = u n ( h) , 由式 ( 1. 4) 可得
j 1 w n ( h) = 2!& i
( 3. 3)
∋
∀
(n (
#
d t , 且 an ( 0) = w ( x n ) = x n ( 1 - x n ) , 将方
L- 1
程( 2 . 7) 对 t 取拉氏变换, 得
1 L- 1
n= 1
∃ [ s (n ( s ) 0
k = 1, 2 , ∀, L - 1. 令
a n( 0 ) ]
n ∃ ∋(
1 - 1/ 2 x ) ∃k ( x ) d x + ∀ ( ) s ∃ (n ( s) 2 n= 1
1 0
1 0
( ∃∋( x ) ) ∗ ∃ ∗ k( x ) d x =
n
s
- 3/ 2
(2+
3 3 2 xx )∃ k( x ) dx , 2 2 !Ls
k+ 1 L k L - 1/ 2
2s - 1/ 2 s + 2 !Ls = a, 3L 6L gk = s
- 3/ 2 2 ( Lx - k + 1) ( 2 + 3 x - 3 x ) d x + k- 1 2 2 L k L k L
1 Lubich 的拉普拉斯变换数值逆
给出网格 t = 0, h, 2 h, ∀, Nh, 卷积
t
f* g= 可以离散为
0
0
f ( s ) g( t - s) d s ,
( t # 0) ,
( 1. 1)
∃
jh
w j ( h) g ( t - jh ) ,
t %
( 1. 2)
其中 w j ( h) 由幂级数 F( ( ) / h) = 给出 , 这里 F 是 f 的拉普拉斯变换 , ( ) = 由罗朗定理 [ 4] 有 w j ( h) = 其中 ∀#: = # , #> 0 为常数. 1 2!& i
第3期
黎丽梅 : 一类偏积分微分方程的数值解法
203
1, 2, ∀, L - 1) , 所有满足这些条件的函数构成集合 U ∋, 即 U∋ = { ∃ ∋ ( C [ 0, 1] , 且在每个[ x n- 1 , x n ] ∋ ∃ 上是线性函数, ∃∋( 0) = ∃ ( 1) = 0} , 容易验证 U∋ 是实数域上的一个 L - 1 维的线性空间, U∋ 的基函数为 ∋ n- 1 0, 0 x < , L Lx - n + 1, ∃ n( x ) = n + 1 - Lx , 0, n- 1 L x < n, L ( n = 1, 2 , ∀, L - 1)
( 1. 5) ( 1. 6a) ( 1. 6b) ( 1. 6c)
1 的某个邻域内解析且没有零点 , 除在 < 1, &> ∃,