1. 引言

二阶复合型非齐次线性微分方程在应用数学、理论物理、工程学、控制理论和最优化理论等领域广泛存在，其重要性体现在对各种物理、生物问题和工程问题的建模。在物理学中，这类微分方程被用于描述振动系统、热传导、电路等现象。在生物学研究中，这类微分方程可描述生物体内的动力学过程，例如药物分布和生态系统演变。在工程领域，它们则常被用于建模结构、材料性质和电路响应。因此，二阶复合型非齐次线性微分方程的求解方法吸引了众多学者研究。

Landahl和Hanson [1] 利用微分法解决了两点边值问题。一些作者 [2] [3] [4] 利用格林函数法求解两点边值问题。Coville和Dupaigne [5] 利用粒子群优化算法研究了线性常微分方程两点边值问题。利用格林函数方法研究微分方程边值问题是非常方便和容易的。基于文献 [3] 的研究，陆静 [6] 利用常数变易法研究了二阶常微分方程的求解问题，并给出了格林函数法在不同常微分方程两点边界条件下的求解方法、解的表达式及其性质。孙赫 [7] 利用常数变易法求得了一类二阶常微分方程周期边值问题的格林函数，并总结出一类高阶常微分方程边值问题的格林函数的求法。刘丽环 [8] 利用格林函数方法讨论了一类二阶常微分方程边值问题解的存在唯一性。该文献借助格林函数将边值问题等价于积分方程，并借助压缩映射不动点定理证明解的存在唯一性。刘慧 [9] 应用文献 [8] 的方法研究了一类二阶常微分方程在不同边值条件下的格林函数，从而研究积分方程解的存在性。田献珍 [10] 研究了二阶非齐次线性常微分方程在相同边界条件下和不同情况下的格林函数。利用格林函数给出了非齐次线性边值问题唯一解的精确表达式。李培超 [11] 利用通解和格林函数的方法给出了点源函数形式的二阶线性常微分方程(ODE)及其相应的非齐次项边值问题的解。这些方法可以应用于求解线性常微分方程以及相关的带有一般非齐次项的定解问题。

黄力 [12] 将二阶微分方程两点边值问题转化为弗雷德霍姆积分方程，利用泰勒矩阵方法得到了其近似解，给出了未知函数的n阶泰勒展开式，并通过矩阵运算得到了泰勒展开式的系数。马翠 [13] 基于变分原理提出了一个与二阶线性常微分方程问题等价的变分问题，并通过两点三次厄尔米特插值找到了一个近似函数。然后，将二阶线性常微分方程问题转化为多变量单目标优化问题，并采用粒子群算法进行求解。该解为二阶线性常微分方程的近似解。通过与瑞利–里兹法和有限元法的比较，该方法令人非常满意。

陈铁军 [14] 利用差分法研究边值问题，提出了三种边界条件，再根据边界条件对线性常微分方程进行求解。采用这种高精度差分法获得的结果具有较高精准度。刘杨 [15] 利用两种再生核方法求解线性常微分方程的初边值问题，并通过了一些数值实例比较了两种方法的精确度。黄佳玥 [16] 通过将一类边值问题转化为含单未知参数的初值问题，再利用初值问题的数值解法和有效的迭代法给出边值问题的数值解。叶康生 [17] 提出p型超收敛算法可应用于非线性常微分方程边值问题。泰勒展开可将非线性常微分方程边值问题转化为线性常微分方程边值问题，避免了迭代运算，提高了效率。李明英 [18] 针对二阶线性常微分方程两点边值问题，通过二阶导数的组合式近似构造了具有最高阶精度的等步长的五点差分格式，并将其对多个算例进行编程计算，再将得到的结果与精确解作比较，可发现所构成的差分格式达到了最高的八阶精度。该方法提高了计算精度。Riverafigueroa等人 [19] 给出了求解非齐次二阶常系数线性微分方程的一种直接方法。这种方法的优点是不要求初值问题解的唯一性和存在性定理。

李顺初 [20] 在分析一类二阶复合线性齐次微分方程边值问题解的基础上，研究了相似核函数和解的相似结构，提出了求解该类边值问题的新方法。王俊超 [21] 通过简化二阶线性齐次微分方程线性非齐次边值问题的解析表达式，得到了解的形式相似性，并且利用简单边界条件，为构造具有复杂边界条件的方程解提供了一种新的方法。彭春 [22] 基于二阶线性齐次常微分方程第三类边值问题的相关理论，给出了求解二阶齐次常微分方程弹性边值问题的相似构造法，并将其应用于一些特殊方程的弹性边值问题中，得到了解的相似结构。为后续发展奠定了坚实的理论基础。但现有的研究在特定问题或边界条件下的适用性、精确性、计算效率等方面仍存在局限性。

基于上述研究成果，本文研究了二阶复合型非齐次线性微分方程的边值问题，给出了此类边值问题的求解方法——相似构造法 [20] ，总结了该方法的求解步骤。对于上述研究中的边值问题，传统的求解方法格林函数法、泰勒矩阵方法、变分原理、差分法等计算过程复杂，而相似构造法简化了计算过程，通过相对简单的数学运算将引解函数、相似核函数、内外边界和衔接条件系数组合即可得到边值问题的解。本文的研究将相似构造法应用于非齐次情况，并对该方法进行了扩展。这有助于读者更好地理解本文研究的独特性和贡献，同时扩充了微分方程边值问题的可解类。

2. 复合二阶线性非齐次常微分方程边值问题的求解

2.1. 复合二阶线性非齐次常微分方程边值问题

基于以上研究，本文针对复合二阶线性非齐次微分方程边值问题进行如下研究：

$\{\begin{array}{l}{y^{″}}_1+p_1\left(x\right){y^{\prime }}_1+q_1\left(x\right)y_1=f\left(x\right)\left(a\le x\le b\right)\\ {y^{″}}_2+p_2\left(x\right){y^{\prime }}_2+q_2\left(x\right)y_2=0\left(b\le x\le c\right)\\ {\left[E{y}_1+\left(1+EF\right){y^{\prime }}_1\right]}_{x=a}=D\\ {y_1|}_{x=b}={\lambda }_1{y_2|}_{x=b},{{y^{\prime }}_1|}_{x=b}={\lambda }_2{{y^{\prime }}_2|}_{x=b}\\ {\left[M{y}_2+N{y^{\prime }}_2\right]}_{x=c}=0\end{array}$(1)

其中： $D,E,F,M,N,{\lambda }_1,{\lambda }_2$是常数， ${\lambda }_1{\lambda }_2\ne 0$， $M^2+N^2\ne 0$， $p_i\left(x\right)\in C^1\left[a,b\right]$， $q_i\left(x\right)\in C^2\left[a,b\right]$ $\left(i=1,2\right)$。

2.2. 边值问题的求解

利用叠加原理，边值问题(1)可以分解为以下两个边值问题：

$\{\begin{array}{l}{y^{″}}_1+p_1\left(x\right){y^{\prime }}_1+q_1\left(x\right)y_1=0\left(a\le x\le b\right)\\ {y^{″}}_2+p_2\left(x\right){y^{\prime }}_2+q_2\left(x\right)y_2=0\left(b\le x\le c\right)\\ {\left[E{y}_1+\left(1+EF\right){y^{\prime }}_1\right]}_{x=a}=D\\ {y_1|}_{x=b}={\lambda }_1{y_2|}_{x=b},{{y^{\prime }}_1|}_{x=b}={\lambda }_2{{y^{\prime }}_2|}_{x=b}\\ {\left[M{y}_2+N{y^{\prime }}_2\right]}_{x=c}=0\end{array}$(2)

$\{\begin{array}{l}{y^{″}}_1+p_1\left(x\right){y^{\prime }}_1+q_1\left(x\right)y_1=f\left(x\right)\left(a\le x\le b\right)\\ {y^{″}}_2+p_2\left(x\right){y^{\prime }}_2+q_2\left(x\right)y_2=0\left(b\le x\le c\right)\\ {\left[E{y}_1+\left(1+EF\right){y^{\prime }}_1\right]}_{x=a}=0\\ {y_1|}_{x=b}={\lambda }_1{y_2|}_{x=b},{{y^{\prime }}_1|}_{x=b}={\lambda }_2{{y^{\prime }}_2|}_{x=b}\\ {\left[M{y}_2+N{y^{\prime }}_2\right]}_{x=c}=0\end{array}$(3)

定解方程 ${y^{″}}_1+p_1\left(x\right){y^{\prime }}_1+q_1\left(x\right)y_1=f\left(x\right)$的一般解为：

$y_1\left(x\right)=A_1{y}_{11}\left(x\right)+B_1{y}_{12}\left(x\right)+g\left(x\right)$(4)

定解方程 ${y^{″}}_2+p_2\left(x\right){y^{\prime }}_2+q_2\left(x\right)y_2=0$的一般解为：

$y_2\left(x\right)=A_2{y}_{21}\left(x\right)+B_2{y}_{22}\left(x\right)$(5)

定义引解函数如下：

${\phi }_{0,0}^i\left(x,\xi \right)=y_{i1}\left(x\right)y_{i2}\left(\xi \right)-y_{i2}\left(x\right)y_{i1}\left(\xi \right)\left(i=1,2\right)$(6)

${\phi }_{1,0}^i\left(x,\xi \right)=\frac{\partial }{\partial x}{\phi }_{0,0}^i\left(x,\xi \right)={y^{\prime }}_{i1}\left(x\right)y_{i2}\left(\xi \right)-{y^{\prime }}_{i2}\left(x\right)y_{i1}\left(\xi \right)\left(i=1,2\right)$(7)

${\phi }_{0,1}^i\left(x,\xi \right)=\frac{\partial }{\partial \xi }{\phi }_{0,0}^i\left(x,\xi \right)=y_{i1}\left(x\right){y^{\prime }}_{i2}\left(\xi \right)-y_{i2}\left(x\right){y^{\prime }}_{i1}\left(\xi \right)\left(i=1,2\right)$(8)

${\phi }_{1,1}^i\left(x,\xi \right)=\frac{{\partial }^2}{\partial x\partial \xi }{\phi }_{0,0}^i\left(x,\xi \right)={y^{\prime }}_{i1}\left(x\right){y^{\prime }}_{i2}\left(\xi \right)-{y^{\prime }}_{i2}\left(x\right){y^{\prime }}_{i1}\left(\xi \right)\left(i=1,2\right)$(9)

方程(7)~(9)通过 ${\phi }_{0,0}^i\left(x,\xi \right)$对x和 $\xi$求偏导得到，其中 $i=1$表示内区 $\left(a\le x\le b\right)$， $i=2$表示外区 $\left(b\le x\le c\right)$。

定理1：若边值问题(2)有唯一解，则内外区的解分别表示为 [10]

$y_1^{\left(1\right)}\left(x\right)=D\cdot \frac{1}{E{\Phi }_1\left(a\right)+1+EF}\cdot {\Phi }_1\left(x\right)\left(a\le x\le b\right)$(10)

$y_2^{\left(1\right)}\left(x\right)=D\cdot \frac{1}{E{\Phi }_1\left(a\right)+1+EF}\cdot \frac{{\phi }_{0,1}^1\left(b,b\right)}{{\lambda }_1{\Phi }_2\left(b\right){\phi }_{1,1}^1\left(a,b\right)-{\lambda }_2{\phi }_{1,0}^1\left(a,b\right)}\cdot {\Phi }_2\left(x\right)\left(b\le x\le c\right)$(11)

其中 ${\Phi }_2\left(x\right)$， ${\Phi }_1\left(x\right)$分别称为外区和内区的相似核函数，表示为：

${\Phi }_2\left(x\right)=\frac{M{\phi }_{0,0}^2\left(x,c\right)+N{\phi }_{0,1}^2\left(x,c\right)}{M{\phi }_{1,0}^2\left(b,c\right)+N{\phi }_{1,1}^2\left(b,c\right)}\left(b\le x\le c\right)$(12)

${\Phi }_1\left(x\right)=\frac{{\lambda }_2{\phi }_{0,0}^1\left(x,b\right)-{\lambda }_1{\Phi }_2\left(b\right){\phi }_{0,1}^1\left(x,b\right)}{{\lambda }_2{\phi }_{1,0}^1\left(a,b\right)-{\lambda }_1{\Phi }_2\left(b\right){\phi }_{1,1}^1\left(a,b\right)}\left(a\le x\le b\right)$(13)

定理2：若边值问题(3)有唯一解，则内外区的解分别表示为：

$\begin{array}{c}{y}_1^{\left(2\right)}\left(x\right)=\frac{1}{E{\Phi }_1\left(a\right)+1+EF}\times \left\{\left[-Eg\left(a\right)-\left(1+EF\right)g^{\prime }\left(a\right)\right]{\Phi }_1\left(x\right)+\left[{\lambda }_{2}g\left(b\right)-{\lambda }_1{g}^{\prime }\left(b\right){\Phi }_2\left(b\right)\right]{\Phi }_3\left(x\right)\right\}+g\left(x\right)\\ \left(a\le x\le b\right)\end{array}$(14)

$\begin{array}{c}{y}_2^{\left(2\right)}\left(x\right)=\frac{1}{E{\Phi }_1\left(a\right)+1+EF}\times \left\{-g^{\prime }\left(b\right){\Phi }_3\left(b\right)+g\left(b\right){\Phi }_4\left(b\right)-\left[Eg\left(a\right)+\left(1+EF\right)g^{\prime }\left(a\right)\right]{\phi }_{0,1}^1\left(b,b\right)\right\}{\Phi }_2\left(x\right)\\ \left(b\le x\le c\right)\end{array}$(15)

这里 ${\Phi }_3\left(x\right)$， ${\Phi }_4\left(x\right)$表示为：

${\Phi }_3\left(x\right)=\frac{E{\phi }_{0,0}^1\left(a,x\right)+\left(1+EF\right){\phi }_{1,0}^1\left(a,x\right)}{-{\lambda }_2{\phi }_{1,0}^1\left(a,b\right)+{\lambda }_1{\phi }_{1,1}^1\left(a,b\right){\Phi }_2\left(b\right)}\left(a\le x\le b\right)$(16)

${\Phi }_4\left(x\right)=\frac{E{\phi }_{0,1}^1\left(a,x\right)+\left(1+EF\right){\phi }_{1,1}^1\left(a,x\right)}{-{\lambda }_2{\phi }_{1,0}^1\left(a,b\right)+{\lambda }_1{\phi }_{1,1}^1\left(a,b\right){\Phi }_2\left(b\right)}\left(a\le x\le b\right)$(17)

证明：将(4)~(5)式代入内外区的边界条件和两个收敛条件，得：

$\left[E{y}_{11}\left(a\right)+\left(1+EF\right){y^{\prime }}_{11}\left(a\right)\right]A_1+\left[E{y}_{12}\left(a\right)+\left(1+EF\right){y^{\prime }}_{12}\left(a\right)\right]B_1=-Eg\left(a\right)-\left(1+EF\right)g^{\prime }\left(a\right)$(18)

$y_{11}\left(b\right)A_1+y_{12}\left(b\right)B_1-{\lambda }_1{y}_{21}\left(b\right)A_2-{\lambda }_1{y}_{22}\left(b\right)B_2=-g\left(b\right)$(19)

${y^{\prime }}_{11}\left(b\right)A_1+{y^{\prime }}_{12}\left(b\right)B_1-{\lambda }_2{y^{\prime }}_{21}\left(b\right)A_2-{\lambda }_2{y^{\prime }}_{22}\left(b\right)B_2=-g^{\prime }\left(b\right)$(20)

$\left[M{y}_{21}\left(c\right)+N{y^{\prime }}_{21}\left(c\right)\right]A_2+\left[M{y}_{22}\left(c\right)+N{y^{\prime }}_{22}\left(c\right)\right]B_2=0$(21)

根据边值问题(3)解的存在唯一性，可以发现线性方程的系数行列式(方程(18)~(21))关于待定系数不等于零，且

$\begin{array}{c}\Delta =\left|\begin{array}{cccc}E{y}_{11}\left(a\right)+\left(1+EF\right){y^{\prime }}_{11}\left(a\right)& E{y}_{12}\left(a\right)+\left(1+EF\right){y^{\prime }}_{12}\left(a\right)& 0& 0\\ y_{11}\left(b\right)& y_{12}\left(b\right)& -{\lambda }_1{y}_{21}\left(b\right)& -{\lambda }_1{y}_{22}\left(b\right)\\ {y^{\prime }}_{11}\left(b\right)& {y^{\prime }}_{12}\left(b\right)& -{\lambda }_2{y^{\prime }}_{21}\left(b\right)& -{\lambda }_2{y^{\prime }}_{22}\left(b\right)\\ 0& 0& M{y}_{21}\left(c\right)+N{y^{\prime }}_{21}\left(c\right)& M{y}_{22}\left(c\right)+N{y^{\prime }}_{22}\left(c\right)\end{array}\right|\\ =E\left\{{\phi }_{0,0}^1\left(a,b\right)\left[-{\lambda }_{2}M{\phi }_{1,0}^2\left(b,c\right)-{\lambda }_{2}N{\phi }_{1,1}^2\left(b,c\right)\right]-{\phi }_{0,1}^1\left(a,b\right)\left[-{\lambda }_{1}M{\phi }_{0,0}^2\left(b,c\right)-{\lambda }_{1}N{\phi }_{0,1}^2\left(b,c\right)\right]\right\}\\ +\left(1+EF\right)\left\{{\phi }_{1,0}^1\left(a,b\right)\left[-{\lambda }_{2}M{\phi }_{1,0}^2\left(b,c\right)-{\lambda }_{2}N{\phi }_{1,1}^2\left(b,c\right)\right]-{\phi }_{1,1}^1\left(a,b\right)\left[-{\lambda }_{1}M{\phi }_{0,0}^2\left(b,c\right)-{\lambda }_{1}N{\phi }_{0,1}^2\left(b,c\right)\right]\right\}\\ =\left[E{\Phi }_1\left(a\right)+1+EF\right]\\ \times \{{\phi }_{1,0}^1\left(a,b\right)\left[-{\lambda }_{2}M{\phi }_{1,0}^2\left(b,c\right)-{\lambda }_{2}N{\phi }_{1,1}^2\left(b,c\right)\right]-{\phi }_{1,1}^1\left(a,b\right)\left[-{\lambda }_{1}M{\phi }_{0,0}^2\left(b,c\right)-{\lambda }_{1}N{\phi }_{0,1}^2\left(b,c\right)\right]\}\end{array}$

由克莱姆法则得到 $A_1,B_1,A_2,B_2$的值如下：

$\begin{array}{l}{A}_1=\frac{1}{\Delta }\{\left[-Eg\left(a\right)-\left(1+EF\right)g^{\prime }\left(a\right)\right]\{y_{12}\left(b\right)\left[-{\lambda }_{2}M{\phi }_{1,0}^2\left(b,c\right)-{\lambda }_{2}N{\phi }_{1,1}^2\left(b,c\right)\right]\\ -{y^{\prime }}_{12}\left(b\right)\left[-{\lambda }_{1}M{\phi }_{0,0}^2\left(b,c\right)-{\lambda }_{1}N{\phi }_{0,1}^2\left(b,c\right)\right]\}-\left[E{y}_{12}\left(a\right)+\left(1+EF\right){y^{\prime }}_{12}\left(a\right)\right]\\ \times \left\{-g\left(b\right)\left[-{\lambda }_{2}M{\phi }_{1,0}^2\left(b,c\right)-{\lambda }_{2}N{\phi }_{1,1}^2\left(b,c\right)\right]+g^{\prime }\left(b\right)\left[-{\lambda }_{1}M{\phi }_{0,0}^2\left(b,c\right)-{\lambda }_{1}N{\phi }_{0,1}^2\left(b,c\right)\right]\right\}\end{array}$(22)

$\begin{array}{l}{B}_1=\frac{1}{\Delta }\{\left[E{y}_{11}\left(a\right)+\left(1+EF\right){y^{\prime }}_{11}\left(a\right)\right]\{-g\left(b\right)\left[-{\lambda }_{2}M{\phi }_{1,0}^2\left(b,c\right)-{\lambda }_{2}N{\phi }_{1,1}^2\left(b,c\right)\right]\\ +g^{\prime }\left(b\right)\left[-{\lambda }_{1}M{\phi }_{0,0}^2\left(b,c\right)-{\lambda }_{1}N{\phi }_{0,1}^2\left(b,c\right)\right]\}+\left[Eg\left(a\right)+\left(1+EF\right)g^{\prime }\left(a\right)\right]\\ \times \left\{y_{11}\left(b\right)\left[-{\lambda }_{2}M{\phi }_{1,0}^2\left(b,c\right)-{\lambda }_{2}N{\phi }_{1,1}^2\left(b,c\right)\right]-{y^{\prime }}_{11}\left(b\right)\left[-{\lambda }_{1}M{\phi }_{0,0}^2\left(b,c\right)-{\lambda }_{1}N{\phi }_{0,1}^2\left(b,c\right)\right]\right\}\end{array}$(23)

$\begin{array}{l}{A}_2=\frac{1}{\Delta }\left[M{y}_{22}\left(c\right)+N{y^{\prime }}_{22}\left(c\right)\right]\{\left[E{y}_{11}\left(a\right)+\left(1+EF\right){y^{\prime }}_{11}\left(a\right)\right]\left[-y_{12}\left(b\right)g^{\prime }\left(b\right)+{y^{\prime }}_{12}\left(b\right)g\left(b\right)\right]\\ -\left[E{y}_{12}\left(a\right)+\left(1+EF\right){y^{\prime }}_{12}\left(a\right)\right]\left[-y_{11}\left(b\right)g^{\prime }\left(b\right)+{y^{\prime }}_{11}\left(b\right)g\left(b\right)\right]\\ -\left[Eg\left(a\right)+\left(1+EF\right)g^{\prime }\left(a\right)\right]{\phi }_{0,1}^1\left(b,b\right)\}\end{array}$(24)

$\begin{array}{l}{B}_2=\frac{-1}{\Delta }\left[M{y}_{21}\left(c\right)+N{y^{\prime }}_{21}\left(c\right)\right]\{-E{\phi }_{0,0}^1\left(a,b\right)g^{\prime }\left(b\right)+E{\phi }_{0,1}^1\left(a,b\right)g\left(b\right)-\left(1+EF\right){\phi }_{1,0}^1\left(a,b\right)\\ \times g^{\prime }\left(b\right)+\left(1+EF\right){\phi }_{1,1}^1\left(a,b\right)g\left(b\right)-\left[Eg\left(a\right)+\left(1+EF\right)g^{\prime }\left(a\right)\right]{\phi }_{0,1}^1\left(b,b\right)\}\end{array}$(25)

将 $A_1,B_1,A_2,B_2$的值代入方程(4)~(5)，并且利用相似核函数 ${\Phi }_1\left(x\right)$， ${\Phi }_2\left(x\right)$， ${\Phi }_3\left(x\right)$和 ${\Phi }_1\left(a\right)$， ${\Phi }_2\left(b\right)$， ${\Phi }_3\left(b\right)$， ${\Phi }_4\left(b\right)$的值分别得到边值问题(3)的内区和外区的解，即方程(14)~(15)。

定理3：利用叠加原理，边值问题(1)的内外区的解为

$\begin{array}{l}{y}_1\left(x\right)=y_1^{\left(1\right)}\left(x\right)+y_1^{\left(2\right)}\left(x\right)\\ =\frac{1}{E{\Phi }_1\left(a\right)+1+EF}\\ \times \{\left[D-Eg\left(a\right)-\left(1+EF\right)g^{\prime }\left(a\right)\right]{\Phi }_1\left(x\right)+\left[{\lambda }_{2}g\left(b\right)-{\lambda }_1{g}^{\prime }\left(b\right){\Phi }_2\left(b\right)\right]{\Phi }_3\left(x\right)\}+g\left(x\right)\\ \left(a\le x\le b\right)\end{array}$(26)

$\begin{array}{l}{y}_2\left(x\right)=y_2^{\left(1\right)}\left(x\right)+y_2^{\left(2\right)}\left(x\right)\\ =\frac{1}{E{\Phi }_1\left(a\right)+1+EF}\cdot \{\{\frac{D}{{\lambda }_1{\Phi }_2\left(b\right){\phi }_{1,1}^1\left(a,b\right)-{\lambda }_2{\phi }_{1,0}^1\left(a,b\right)}\\ \stackrel{}{}-\left[Eg\left(a\right)+\left(1+EF\right)g^{\prime }\left(a\right)\right]\}{\phi }_{0,1}^1\left(b,b\right)-g^{\prime }\left(b\right){\Phi }_3\left(b\right)+g\left(b\right){\Phi }_4\left(b\right)\}{\Phi }_2\left(x\right)\\ \left(b\le x\le c\right)\end{array}$(27)

3. 常微分方程边值问题的求解步骤

根据上述求解(1)、(2)、(3)边值问题的过程，得到求解复合二阶线性非齐次微分方程边值问题的步骤。具体步骤如下：

Step 1. 求解定解方程

得到边值问题(1)的每一个定解方程的一般解为(4)~(5)。

Step 2. 构造引解函数

利用 $y_{i1}\left(x\right),y_{i2}\left(x\right)\left(i=1,2\right)$分别构造了内外区的引解函数 ${\phi }_{0,0}^i\left(x,\xi \right)\left(i=1,2\right)$，其它引解函数可以通过 ${\phi }_{0,0}^i\left(x,\xi \right)$对x和 $\xi$求偏导得到。

Step 3. 构建内、外区域相似核函数

根据公式(12)、(13)、(16)和(17)，分别利用外区和内区的引解函数、齐次外边界条件的系数M，N、非齐次内边界条件的系数E，F、内外区域两个衔接条件的系数 ${\lambda }_1$， ${\lambda }_2$构造相似核函数。

Step 4. 求边值问题的解

根据公式(26)和(27)，边界问题(1)的内区和外区的解分别由非齐次内边界条件的系数D，E，F、相似核函数 ${\Phi }_1\left(x\right)$， ${\Phi }_2\left(x\right)$， ${\Phi }_3\left(x\right)$、内外区引解函数的值 ${\Phi }_1\left(a\right)$， ${\Phi }_2\left(b\right)$， ${\Phi }_3\left(b\right)$， ${\Phi }_4\left(b\right)$和内外区域两个衔接条件的系数 ${\lambda }_1$， ${\lambda }_2$组合得到。

4. 结论

1) 本文解决了复合二阶线性非齐次常微分方程的边值问题。利用叠加原理，得到边值问题解的相似结构。

2) 根据上述方法的步骤，可以看出该方法是求解一般复合二阶线性微分方程边值问题的一种方便、有效、有创新性的方法。

基金项目

西华大学人才引进项目(Z202068)资助。

参考文献