三区间复合Legendre方程边值问题解的相似结构
Similar Structure of Solutions of Boundary Value Problem of Three-Region Composite Legendre Equation
摘要: 针对三区复合Legendre方程的一类边值问题,对其解的表达形式进行了分析,发现其解式具有相似的结构;首先通过方程的两个线性无关解构造内、中、外区引解函数;再借助外边界条件、衔接性条件的系数以及引解函数构造内、中、外区相似核函数;最后将引解函数、相似核函数以及定解条件的系数进行组装得到其解,进而获得了求解此类边值问题的新方法——相似构造法,该方法简化了求解该类边值问题的复杂性,并提高了求解的效率。
Abstract: For a class of boundary value problems of three-region composite Legendre equation, the expres-sion of the solutions is analyzed, and it is found that its solution has a similar structure. Firstly, we can construct the functions of guide solution of the inner, middle and outer regions through the two linear independent solutions of the definite solution equation. Secondly, through the outer boundary condition, the connection condition coefficients and the functions of guide solution to con-struct the similar kernel functions in the inner, middle and outer regions. Finally, the solution is ob-tained by assembling the functions of guide solution, the similar kernel function and the coefficients of the definite solution conditions. Thus a new method for solving such boundary value problems is obtained, which is called the similarity construction method. This method simplifies the complexity of solving this kind of boundary value problem and improves the efficiency of solving them.
文章引用:付雪倩, 李顺初, 邵东凤, 刘盼, 范林. 三区间复合Legendre方程边值问题解的相似结构[J]. 理论数学, 2023, 13(2): 172-181. https://doi.org/10.12677/PM.2023.132020

1. 引言

工程中的许多偏微分方程的问题常常会转化为常微分方程问题来进行研究 [1] ,因此研究常微分方程的边值问题对实际生活具有重要意义。许多学者对微分方程边值问题进行了研究,张宇飞 [2] 等人研究了一类常微分方程组在定值右边界条件下的解的结构,发现其解的结构具有一定的内在联系。李顺初 [3] 研究了一类二阶齐次线性常微分方程的边值问题,通过分析整理得出该类方程的解具有相似结构。2010年,李顺初 [4] 通过对二阶齐次线性常微分方程(组)的边值问题和二阶齐次线性偏微分方程的混合问题进行了初探,提出了求解微分方程边值问题的方法——相似构造法。该方法运用广泛,董晓旭、郭丽洁、唐乙斌、李顺初、何签等人将该方法应用于三区复合Thomson方程 [5] 、Hermit方程 [6] 、Tschebyscheff方程 [7] 、三区间复合型第二种Weber方程 [8] 、三区间复合型第一种Weber方程 [9] 的边值问题的求解中。相似构造法还为解决工程问题带来了便利,比如:李顺初等人将该方法运用于均质油藏 [10] 、复合油藏 [11] 、多层油藏 [12] 等渗流模型的求解中,使其得到了不断的完善和发展。

在特殊函数中,Legendre方程具有重要的地位,他的解Legendre多项式在求解数学物理方程中具有广泛的运用 [13] [14] 。比如在求解带有轴对称性的带电体的静电势和球形域内电位分布常常要借助Legendre多项式来求解 [15] 。那么对Legendre方程的求解进行研究以及对其解式进行研究就显得十分有必要。但是一般的方法对其进行求解十分的繁琐,能否找到一种简便的方法进行求解呢?近年来,李顺初等人对Legendre方程 [16] ,连带Legendre方程 [17] ,复合Legendre方程 [18] 的解式进行了研究,并找到了其解式的相似结构,但是对于三区间复合及其以上的Legendre方程解式的研究还没有。基于以上的研究,下面将对三区复合Legendre方程的边值问题进行研究:

{ ( 1 x 2 ) y 1 2 x y 1 + l 1 ( l 1 + 1 ) y 1 = 0 , ( 1 a x < b ) ( 1 x 2 ) y 2 2 x y 2 + l 2 ( l 2 + 1 ) y 2 = 0 , ( b < x < c ) ( 1 x 2 ) y 3 2 x y 3 + l 3 ( l 3 + 1 ) y 3 = 0 , ( c < x < d 1 ) [ E y 1 + ( 1 + E F ) y 1 ] | x = a = M y 1 | x = b = λ 1 y 2 | x = b , y 1 | x = b = λ 2 y 2 | x = b y 2 | x = c = μ 1 y 3 | x = c , y 2 | x = c = μ 2 y 3 | x = c [ G y 3 + H y 3 ] | x = d = 0 (1)

其中 a b c d E F M λ 1 λ 2 μ 1 μ 2 G H 均为实数, l i ( i = 1 , 2 , 3 ) 为正整数,且满足 M 0 G 2 + H 2 0

本文第一部分给出了研究所需要用到的基础知识,第二部分给出了本文的主要定理及其证明,第三部分通过对定理的证明归纳总结出相似构造法的求解步骤,第四部分将相似构造法的求解步骤进行举例应用,第五部分为本文的总结。

2. 预备知识

引理1 Legendre方程 ( 1 x 2 ) y i 2 x y i + μ i y i = 0 , ( 1 x 1 ) 它的本征值为 μ i = l i ( l i + 1 ) ( l = 0 , 1 , 2 , , i = 1 , 2 , 3 ) ,两个线性无关解为 P l i Q l i 则其方程的通解 [14] 为:

y i ( x ) = A i P l i ( x ) + B i Q l i ( x )

其中 A i B i 为任意常数。 P l i ( x ) 为Legendre多项式, Q l i ( x ) 为第二类Legendre函数, i = 1 , 2 , 3

引理2 Legendre函数的微分性质 [19]

P l i ( x ) = l i 1 x 2 [ P l i 1 ( x ) x P l i ( x ) ]

Q l i ( x ) = l i 1 x 2 [ Q l i 1 ( x ) x Q l i ( x ) ]

证明:见参考文献 [19]

引理3 由Legendre方程的两个线性无关解 P l i ( x ) Q l i ( x ) 构造二元函数

φ 0 , 0 i ( x , ξ ) = P l i ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i ( x ) (2)

则有

φ 1 , 0 i ( x , ξ ) = l i 1 x 2 { P l i 1 ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i 1 ( x ) x [ P l i ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i ( x ) ] } (3)

φ 0 , 1 i ( x , ξ ) = l i 1 ξ 2 { P l i ( x ) Q l i 1 ( ξ ) P l i 1 ( ξ ) Q l i ( x ) ξ [ P l i ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i ( x ) ] } (4)

φ 1 , 1 i ( x , ξ ) = l i 2 ( 1 x 2 ) ( 1 ξ 2 ) { P l i 1 ( x ) Q l i 1 ( ξ ) P l i 1 ( ξ ) Q l i 1 ( x ) x [ P l i ( x ) Q l i 1 ( ξ ) P l i 1 ( ξ ) Q l i ( x ) ] ξ [ P l i 1 ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i 1 ( x ) ] + x ξ [ P l i ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i ( x ) ] } (5)

其中 φ 0 , 0 i ( x , ξ ) , φ 1 , 0 i ( x , ξ ) , φ 0 , 1 i ( x , ξ ) , φ 1 , 1 i ( x , ξ ) 称作引解函数, i = 1 表示内区 1 a x < b i = 2 表示中区 b < x < c i = 3 表示外区 c < x < d 1

证明:根据引理2有

P l i ( x ) = l i 1 x 2 [ P l i 1 ( x ) x P l i ( x ) ]

Q l i ( x ) = l i 1 x 2 [ Q l i 1 ( x ) x Q l i ( x ) ]

φ 0 0 i ( x , ξ ) x 求偏导数有

φ 1 , 0 i ( x , ξ ) = x φ 0 , 0 i ( x , ξ ) = l i 1 x 2 [ P l i 1 ( x ) x P l i ( x ) ] Q l i ( ξ ) l i 1 x 2 [ Q l i 1 ( x ) x Q l i ( x ) ] P l i ( ξ ) = l i 1 x 2 { P l i 1 ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i 1 ( x ) x [ P l i ( x ) Q l i ( ξ ) P l i ( ξ ) Q l i ( x ) ] }

同理可证(4) (5)式。

3. 主要定理及其证明

定理1 若边值问题(1)有唯一解,那么它的内区、中区、外区解可以分别表示为:

y 1 ( x ) = M Φ 1 ( x ) E Φ 1 ( a ) + ( 1 + E F ) (6)

y 2 ( x ) = M 1 E Φ 1 ( a ) + ( 1 + E F ) φ 0 , 1 1 ( b , b ) λ 1 φ 1 , 1 1 ( a , b ) Φ 2 ( b ) λ 2 φ 1 , 0 1 ( a , b ) Φ 2 ( x ) (7)

y 3 ( x ) = M 1 E Φ 1 ( a ) + ( 1 + E F ) φ 0 , 1 1 ( b , b ) φ 0 , 1 2 ( c , c ) Φ 3 ( x ) [ λ 2 φ 1 , 0 1 ( a , b ) λ 1 φ 1 , 1 1 ( a , b ) Φ 2 ( b ) ] [ μ 2 φ 1 , 0 2 ( b , c ) μ 1 φ 1 , 1 2 ( b , c ) Φ 3 ( c ) ] (8)

其中 Φ 1 ( x ) 为内区相似核函数,表示为:

Φ 1 ( x ) = λ 2 φ 0 , 0 1 ( x , b ) λ 1 φ 0 , 1 1 ( x , b ) Φ 2 ( b ) λ 2 φ 1 , 0 1 ( a , b ) λ 1 φ 1 , 1 1 ( a , b ) Φ 2 ( b ) (9)

Φ 2 ( x ) 为中区相似核函数,表示为:

Φ 2 ( x ) = μ 2 φ 0 , 0 2 ( x , c ) μ 1 φ 0 , 1 2 ( x , c ) Φ 3 ( c ) μ 2 φ 1 , 0 2 ( b , c ) μ 1 φ 1 , 1 2 ( b , c ) Φ 3 ( c ) (10)

Φ 3 ( x ) 为外区相似核函数,表示为:

Φ 3 ( x ) = G φ 0 , 0 3 ( x , d ) + H φ 0 , 1 3 ( x , d ) G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) (11)

证明:由引理1知,边值问题(1)的三个定解方程的通解为:

y i ( x ) = A i P l i ( x ) + B i Q l i ( x ) , i = 1 , 2 , 3 (12)

其中 A i B i 为任意常数。

由引理2知

y i ( x ) = l i 1 x 2 { A i [ P l i 1 ( x ) x P l i ( x ) ] + B i [ Q l i 1 ( x ) x Q l i ( x ) ] } (13)

y ( x ) , y ( x ) 代入边值问题(1)的内边值条件、衔接性条件和外边值条件可以得到

A 1 { E P l i ( a ) + ( 1 + E F ) l 1 1 a 2 [ P l 1 1 ( a ) a P l 1 ( a ) ] } + B 1 { E Q l 1 ( a ) + ( 1 + E F ) l 1 1 a 2 [ Q l 1 1 ( a ) a Q l 1 ( a ) ] } = M (14)

A 1 P l 1 ( b ) + B 1 Q l 1 ( b ) A 2 λ 1 P l 2 ( b ) B 2 λ 1 Q l 2 ( b ) = 0 (15)

A 1 l 1 1 b 2 [ P l 1 1 ( b ) b P l 1 ( b ) ] + B 1 l 1 1 b 2 [ Q l 1 1 ( b ) b Q l 1 ( b ) ] A 2 l 2 λ 2 1 b 2 [ P l 2 1 ( b ) b P l 2 ( b ) ] B 2 l 2 λ 2 1 b 2 [ Q l 2 1 ( b ) b Q l 2 ( b ) ] = 0 (16)

A 2 P l 2 ( c ) + B 2 Q l 2 ( c ) A 3 μ 1 P l 3 ( c ) B 3 μ 1 Q l 3 ( c ) = 0 (17)

A 2 l 1 1 c 2 [ P l 2 1 ( c ) c P l 2 ( c ) ] + B 2 l 2 1 c 2 [ Q l 2 1 ( c ) c Q l 2 ( c ) ] A 3 l 3 μ 2 1 c 2 [ P l 3 1 ( c ) c P l 2 ( c ) ] B 3 l 2 μ 2 1 c 2 [ Q l 3 1 ( c ) c Q l 3 ( c ) ] = 0 (18)

A 3 { G P l 3 ( d ) + H l 3 1 d 2 [ P l 3 1 ( d ) d P l 2 ( d ) ] } + B 3 { G Q l 3 ( d ) + H l 3 1 d 2 [ Q l 3 1 ( d ) d Q l 3 ( d ) ] } = 0 (19)

从而由(14)~(19)式可以得到关于待定系数 A 1 , B 1 , A 2 , B 2 , A 3 , B 3 的线性方程组。结合(2)~(5)式得到其系数行列式为

Δ = μ 2 λ 2 φ 1 , 0 2 ( b , c ) [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] [ E φ 0 , 0 1 ( a , b ) + ( 1 + E F ) φ 1 , 0 1 ( a , b ) ] μ 1 λ 2 φ 1 , 1 2 ( b , c ) [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] [ E φ 0 , 0 1 ( a , b ) + ( 1 + E F ) φ 1 , 0 1 ( a , b ) ] μ 2 λ 1 φ 0 , 0 2 ( b , c ) [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] [ E φ 0 , 1 1 ( a , b ) + ( 1 + E F ) φ 1 , 1 1 ( a , b ) ] + μ 1 λ 1 φ 0 , 1 2 ( b , c ) [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] [ E φ 0 , 1 1 ( a , b ) + ( 1 + E F ) φ 1 , 1 1 ( a , b ) ] (20)

由于边值问题(1)有唯一解 [16] ,即 Δ 0 。根据Cramer法则可得

A 1 = M Δ { μ 2 λ 2 φ 1 , 0 2 ( b , c ) Q l 1 ( b ) [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] μ 1 λ 2 φ 1 , 1 2 ( b , c ) Q l 1 ( b ) [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] μ 2 λ 1 φ 0 , 0 2 ( b , c ) l 1 1 b 2 [ Q l 1 1 ( b ) b Q l 1 ( b ) ] [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] + μ 1 λ 1 φ 0 , 1 2 ( b , c ) l 1 1 b 2 [ Q l 1 1 ( b ) b Q l 1 ( b ) ] [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] } (21)

B 1 = M Δ { μ 2 λ 2 φ 1 , 0 2 ( b , c ) P l 1 ( b ) [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] μ 1 λ 2 φ 1 , 1 2 ( b , c ) P l 1 ( b ) [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] μ 2 λ 1 φ 0 , 0 2 ( b , c ) l 1 1 b 2 [ P l 1 1 ( b ) b P l 1 ( b ) ] [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] + μ 1 λ 1 φ 0 , 1 2 ( b , c ) l 1 1 b 2 [ P l 1 1 ( b ) b P l 1 ( b ) ] [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] } (22)

A 2 = M Δ { μ 2 φ 0 , 1 1 ( b , b ) Q l 2 ( c ) [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] + μ 1 φ 0 , 1 1 ( b , b ) l 2 1 c 2 [ Q l 2 1 ( c ) c Q l 2 ( c ) ] [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] } (23)

B 2 = M Δ { μ 2 φ 0 , 1 1 ( b , b ) P l 2 ( c ) [ G φ 1 , 0 3 ( c , d ) + H φ 1 , 1 3 ( c , d ) ] + μ 1 φ 0 , 1 1 ( b , b ) l 2 1 c 2 [ P l 2 1 ( c ) c P l 2 ( c ) ] [ G φ 0 , 0 3 ( c , d ) + H φ 0 , 1 3 ( c , d ) ] } (24)

A 3 = M Δ { G Q l 3 ( d ) + H l 3 1 d 2 [ Q l 3 1 ( d ) d Q l 2 ( d ) ] } φ 0 , 1 1 ( b , b ) φ 0 , 1 2 ( c , c ) (25)

B 3 = M Δ { G P l 3 ( d ) + H l 3 1 d 2 [ P l 3 1 ( d ) d P l 2 ( d ) ] } φ 0 , 1 1 ( b , b ) φ 0 , 1 2 ( c , c ) (26)

将(21)~(26)式带入到Legendre方程的通解(12)式中,并将结果进行整理,即得到边值问题(1)的解式(6)~(8)式。根据定理1容易得到如下推论。

推论1如果Legendre方程边值问题(1)的内边界条件为 y 1 | x = a = 1 ,即 E = 0 , F 0 ,那么边值问题(1)的内区解为

y 1 ( x ) = Φ 1 ( x ) , ( 1 a x < b )

推论2如果Legendre方程边值问题(1)的外边界条件为 y 3 | x = d = 0 ,即 G 0 H = 0 ,那么边值问题(1)的外区相似核函数为

Φ 3 ( x ) = φ 0 , 0 3 ( x , d ) φ 1 , 0 3 ( c , d ) , ( c < x < d )

推论3如果Legendre方程边值问题(1)的外边界条件为 y 3 | x = d = 0 ,即 G = 0 H 0 ,那么外区相似核函数为

Φ 3 ( x ) = φ 0 , 1 3 ( x , d ) φ 1 , 1 3 ( c , d ) , ( c < x < d )

推论4在Legendre方程边值问题(1)中有

[ y 1 + F y 1 ] | x = a = M E + 1 F + Φ 1 ( a )

4. 相似构造法的具体步骤

根据上述定理1的证明过程,我们不难归纳总结出求解三区复合Legendre方程边值问题的相似构造法的具体步骤如下:

第一步,构造内区、中区、外区引解函数。根据Legendre方程的两个线性无关解 P l i ( x ) Q l i ( x ) 构造内区、中区、外区引解函数 φ m , n i ( x , ξ ) ( i = 1 , 2 , 3 , m , n = 0 , 1 ) ,如(5)~(8)式。

第二步,构造内区、中区、外区的相似核函数。根据外边值条件的系数G,H以及外区引解函数 φ m , n 3 ( x , ξ ) ( m , n = 0 , 1 ) 构造外区相似核函数,如(11)式所示,并求出 Φ 3 ( c ) ;根据衔接性条件的系数 μ 1 μ 2 以及中区引解函数 φ m , n 2 ( x , ξ ) ( m , n = 0 , 1 ) 构造中区相似核函数,如(10)所示,并求出 Φ 2 ( b ) ;根据衔接性条件的系数 λ 1 λ 2 以及内区引解函数 φ m , n 1 ( x , ξ ) ( m , n = 0 , 1 ) 构造内区相似核函数,如(9)式所示,并求出 Φ 1 ( a )

第三步,获得三区复合Legendre边值问题(1)的解。将内边值条件的系数M,E,F与内区相似核函数 Φ 1 ( x ) 以及 Φ 1 ( a ) 进行组装得到内区解,即(6)式;将内边值条件的系数M,E,F、衔接性条件的系数 λ 1 λ 2 、内区引解函数、中区相似核函数 Φ 2 ( x ) 以及 Φ 1 ( a ) Φ 2 ( b ) 进行组装得到中区解,即(7)式;将内边值条件的系数M,E,F、衔接性条件的系数 λ 1 λ 2 μ 1 μ 2 、内区引解函数,中区引解函数、外区相似核函数 Φ 3 ( x ) 以及 Φ 1 ( a ) Φ 2 ( b ) Φ 3 ( c ) 进行组装得到外区解,即(8)式。

根据以上的求解步骤可以绘出如图1的求解流程图。

Figure 1. Solution flow chart

图1. 求解流程图

5. 举例

求解如下的边值问题:

{ ( 1 x 2 ) y 1 2 x y 1 + 2 y 1 = 0 , ( 1 < x < 0 ) ( 1 x 2 ) y 2 2 x y 2 + 6 y 2 = 0 , ( 0 < x < 1 2 ) ( 1 x 2 ) y 3 2 x y 3 + 12 y 3 = 0 , ( 1 2 < x < 1 ) [ y 1 + 2 y 1 ] | x = 1 = 1 y 1 | x = 0 = y 2 | x = 0 , y 1 | x = 0 = 2 y 2 | x = 0 y 2 | x = 1 2 = y 3 | x = 1 2 , y 2 | x = 1 2 = 2 y 3 | x = 1 2 [ y 3 + 2 y 3 ] | x = 1 = 0 (27)

通过对比边值问题(1)和边值问题(27)可以得到 l 1 = 1 l 2 = 2 l 3 = 3 a = 1 b = 0 c = 1 2 d = 1 E = 1 F = 1 M = 1 λ 1 = 1 λ 2 = 2 μ 1 = 1 μ 2 = 2 G = 1 ,利用相似构造法对边值问题(27)进行求解。

内区定解方程 ( 1 x 2 ) y 1 2 x y 1 + 2 y 1 = 0 的两个线性无关解为 P 1 ( x ) Q 1 ( x ) ,中区定解方程 ( 1 x 2 ) y 2 2 x y 2 + 6 y 2 = 0 的两个线性无关解为 P 2 ( x ) Q 2 ( x ) ,外区定解方程 ( 1 x 2 ) y 3 2 x y 3 + 12 y 3 = 0 的两个线性无关解为 P 3 ( x ) Q 3 ( x )

第一步,根据边值问题(27)的内区、中区、外区定解方程的两个线性无关解分别构造相应的引解函数。即内区引解函数为:

φ 0 , 0 1 ( x , ξ ) = P 1 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 1 ( x )

φ 1 , 0 1 ( x , ξ ) = 1 1 x 2 { P 0 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 0 ( x ) x [ P 1 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 1 ( x ) ] }

φ 0 , 1 1 ( x , ξ ) = 1 1 ξ 2 { P 1 ( x ) Q 0 ( ξ ) P 0 ( ξ ) Q 1 ( x ) ξ [ P 1 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 1 ( x ) ] }

φ 1 , 1 1 ( x , ξ ) = 1 ( 1 x 2 ) ( 1 ξ 2 ) { P 0 ( x ) Q 0 ( ξ ) P 0 ( ξ ) Q 0 ( x ) x [ P 1 ( x ) Q 0 ( ξ ) P 0 ( ξ ) Q 1 ( x ) ] ξ [ P 0 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 0 ( x ) ] + x ξ [ P 1 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 1 ( x ) ] }

中区引解函数为:

φ 0 , 0 2 ( x , ξ ) = P 2 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 2 ( x )

φ 1 , 0 2 ( x , ξ ) = 2 1 x 2 { P 1 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 1 ( x ) x [ P 2 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 2 ( x ) ] }

φ 0 , 1 2 ( x , ξ ) = 1 1 ξ 2 { P 2 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 2 ( x ) ξ [ P 2 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 2 ( x ) ] }

φ 1 , 1 2 ( x , ξ ) = 4 ( 1 x 2 ) ( 1 ξ 2 ) { P 1 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 1 ( x ) x [ P 2 ( x ) Q 1 ( ξ ) P 1 ( ξ ) Q 2 ( x ) ] ξ [ P 1 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 1 ( x ) ] + x ξ [ P 2 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 2 ( x ) ] }

外区引解函数为:

φ 0 , 0 3 ( x , ξ ) = P 3 ( x ) Q 3 ( ξ ) P 3 ( ξ ) Q 3 ( x )

φ 1 , 0 3 ( x , ξ ) = 3 1 x 2 { P 2 ( x ) Q 3 ( ξ ) P 3 ( ξ ) Q 2 ( x ) x [ P 3 ( x ) Q 3 ( ξ ) P 3 ( ξ ) Q 3 ( x ) ] }

φ 0 , 1 3 ( x , ξ ) = 3 1 ξ 2 { P 3 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 3 ( x ) ξ [ P 3 ( x ) Q 3 ( ξ ) P 3 ( ξ ) Q 3 ( x ) ] }

φ 1 , 1 1 ( x , ξ ) = 9 ( 1 x 2 ) ( 1 ξ 2 ) { P 2 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 2 ( x ) x [ P 3 ( x ) Q 2 ( ξ ) P 2 ( ξ ) Q 3 ( x ) ] ξ [ P 2 ( x ) Q 3 ( ξ ) P 3 ( ξ ) Q 2 ( x ) ] + x ξ [ P 3 ( x ) Q 3 ( ξ ) P 3 ( ξ ) Q 3 ( x ) ] }

第二步,构造相似核函数。由于 G = 1 H = 1 结合外区引解函数构造外区相似核函数为:

Φ 3 ( x ) = φ 0 , 0 3 ( x , 1 ) + φ 0 , 1 3 ( x , 1 ) φ 1 , 0 3 ( 1 2 , 1 ) + φ 1 , 1 3 ( 1 2 , 1 )

由于 μ 1 = 1 μ 2 = 2 并结合中区引解函数和 Φ 3 ( 1 2 ) 构造中区相似核函数为:

Φ 2 ( x ) = 2 φ 0 , 0 2 ( x , 1 2 ) φ 0 , 1 2 ( x , 1 2 ) Φ 3 ( 1 2 ) 2 φ 1 , 0 2 ( 0 , 1 2 ) φ 1 , 1 2 ( 0 , 1 2 ) Φ 3 ( 1 2 )

由于 λ 1 = 1 λ 2 = 2 并结合内区引解函数和 Φ 2 ( 0 ) 构造内区相似核函数为:

Φ 1 ( x ) = 2 φ 0 , 0 1 ( x , 0 ) φ 0 , 1 1 ( x , 0 ) Φ 2 ( 0 ) 2 φ 1 , 0 1 ( 1 , 0 ) φ 1 , 1 1 ( 1 , 0 ) Φ 2 ( 0 )

第三步,根据 E = 1 F = 1 M = 1 并将各部分进行组装,得到

内区解为:

y 1 = Φ 1 ( x ) Φ 1 ( 1 ) + 2

中区解为:

y 2 = 1 Φ 1 ( 1 ) + 2 φ 0 , 1 1 ( 0 , 0 ) φ 1 , 1 1 ( 1 , 0 ) Φ 2 ( 0 ) 2 φ 1 , 0 1 ( 1 , 0 ) Φ 2 ( x )

外区解为:

y 3 = 1 Φ 1 ( 1 ) + 2 φ 0 , 1 1 ( 0 , 0 ) φ 0 , 1 2 ( 1 2 , 1 2 ) Φ 3 ( x ) [ 2 φ 1 , 0 1 ( 1 , 0 ) φ 1 , 1 1 ( 1 , 0 ) Φ 2 ( 0 ) ] [ 2 φ 1 , 0 2 ( 0 , 1 2 ) φ 1 , 1 2 ( 0 , 1 2 ) Φ 3 ( 1 2 ) ]

6. 结论

1) 对于边值问题(1)只要知道其两个线性无关解,就可以构造出内、中、外区引解函数,进而构造出对应的相似核函数;

2) 观察边值问题(1)的解式可以发现,其内、中、外区三个解的表达形式具有相似性,其解式的系数只与内边界条件、衔接性条件、外边界条件的系数有关,并且其解式是由引解函数、相似核函数、边界条件的系数进行组装而成;

3) 在边值问题(1)中,随着边界条件系数的改变,其解的表示形式不会发生改变,改变的是相应的相似核函数,从而对边值问题的求解提供了一种简单、快捷的方法。

参考文献

[1] Dong, X.X., Liu, Z. and Li, S.C. (2019) Similar Constructing Method for Solving Nonlinear Spherical Seepage Model with Quadratic Pressure Gradient of Three-Region Composite Fractal Reservoir. Computational and Applied Mathemat-ics, 38, 1-27.
https://doi.org/10.1007/s40314-019-0847-z
[2] 张宇飞, 李顺初, 郑鹏社. 一类常微分方程组在定值右边界条件下的解的结构[J]. 四川工业学院学报, 2004(S1): 31-33+36.
[3] 李顺初. 二阶齐次线性微分方程的边值问题的解的相似结构[J]. 西华大学学报(自然科学版), 2009, 28(5): 40-41+90.
[4] 李顺初. 微分方程解的相似结构初探与展望[J]. 西华大学学报(自然科学版), 2010, 29(2): 223-226+238.
[5] Wang, Y., Dong, X.X. and Li, S.C. (2016) Study on How to Resolve the Boundary Value Problem of Three-Region Composite Thomson Equation. Journal of Computational Methods in Sciences & Engineering, 16, 111-124.
https://doi.org/10.3233/JCM-160606
[6] 郭丽洁, 李顺初, 严娟. Hermit方程边值问题解的相似结构[J]. 济宁学院学报, 2010, 31(3): 9-10.
[7] 唐乙斌, 李顺初, 严娟, 李全勇. Tschebyscheff方程边值问题解的相似结构[J]. 四川兵工学报, 2011, 32(1): 155-156.
[8] 李顺初, 何签, 夏星, 桂钦民. 三区间复合型第二种Weber方程边值问题的相似构造法[J]. 徐州工程学院学报(自然科学版), 2021, 36(4): 1-9.
[9] 何签, 李顺初, 董晓旭, 等. 三区间复合型第一种Weber方程边值问题求解的新方法[J]. 内蒙古师范大学学报: 自然科学版, 2022, 51(1): 9.
[10] L, S.C., Zhao, C.C., Zheng, P.S., et al. (2019) Analysis of Oil and Gas Flow Characteristics in the Reservoir with the Elastic Outer Boundary. Journal of Petroleum Science and Engineering, 175, 280-285.
https://doi.org/10.1016/j.petrol.2018.12.042
[11] He, Q., Dong, X.X., et al. (2022) The Elasticity of the Outer Boundary and the Solution of Two-Region Composite Reservoir Seepage Model. Petroleum Science and Technology, 40, 2773-2791.
[12] 王强. 多层复合油藏渗流模型新解法的研究[D]: [硕士学位论文]. 成都: 西华大学, 2016.
[13] 杨守文, 王海军. 勒让德方程本征值的确定[J]. 大学数学, 2021, 37(2): 33-36.
[14] 曹佳薇. 基于移位Legendre多项式对三种粘弹性材料本构模型的动力学分析[D]: [硕士学位论文]. 秦皇岛: 燕山大学, 2021.
[15] 余海洋, 方世跃. 关于勒让德多项式递推公式的研究[J]. 四川理工学院学报: 自然科学版, 2008, 21(2): 27-2932.
[16] 夏文文, 李顺初, 桂东冬. Legendre方程边值问题的解的相似结构[J]. 嘉应学院学报, 2015, 33(2): 5-9.
[17] 罗梅, 李顺初. 连带Legendre微分方程边值问题解的相似结构[J]. 重庆工商大学学报(自然科学版), 2015, 32(11): 34-37.
[18] 夏星, 李顺初, 邵东凤, 等. 复合Legendre方程边值问题解的相似结构法[J]. 内江师范学院学报, 2022, 37(2): 41-45+52.
[19] 刘适式, 刘适达. 特殊函数[M]. 北京: 气象出版社, 2002.