Advances in Applied Mathematics 应用数学进展, 2012, 1, 41-48 http://dx.doi.org/10.12677/aam.2012.12006 Published Online November 2012 (http://www.hanspub.org/journal/aam.html) Similar Constructive Method of Solution for a Class of Boundary Value Problem of Thomson Equation Congyin Fan, Shunchu Li, Dongxu Xu Institute of Applied Mathematics, Xihua University, Chengdu Email: 1147138817@qq.com Received: Oct. 8th, 2012; revised: Oct. 10th, 2012; accepted: Oct. 23rd, 2012 Abstract: Aiming at the solution for a class of boundary value problem of Thomson equation, we found that its solution can be obtained by structuring similar kernel function with the coefficients of non-homogeneous left boundary condition. While the similar kernel function is constructed by both the function of guide solu- tion and the coefficients of homogeneous right boundary condition. Thus, a new method of solving such kind of boundary value problem is obtained—Similar Constructive Method. Keywords: Boundary Value Problem; Thomson Equation; Similar Constructive Method; Similar Kernel Function; Function of Guide Solution 一类 Thomson 方程边值问题解 的相似构造法 范聪银,李顺初,许东旭 西华大学应用数学研究所,成都 Email: 1147138817@qq.com 收稿日期:2012 年10 月8日;修回日期:2012 年10 月10 日;录用日期:2012 年10 月23日 摘 要:针对 Thomson方程的一类边值问题的解式进行研究,发现其解式可由相似核函数和非齐次左 边值条件中的系数进行组装得到,而相似核函数由引解函数以及齐次右边值条件中的系数来构造,由 此获得了求解该类边值问题的一种新方法——相似构造法。 关键词:边值问题;Thomson 方程;相似构造法;相似核函数;引解函数 1. 引言 在数论中,任何一个实数都能 用一个连 分式 表示出来[1]。而一个微 分方 程对于某 个固 定类型的 边值 问题而 言,其解曲线簇中的曲线都具有几何上相似或放射相似性;那么其解式能否用一个连分式或连分式的乘积表示 出来呢?自 2004 年李顺初等人分别对一些二阶齐次线性常微分方程、二阶齐次线性常微分方程组、二阶齐次线 性偏微分方程、二阶齐次线性偏微分方程组以及一些油气藏工程中的渗流方程的解进行了研究,不仅对此做出 了肯定的回答,还逐渐形成了微分方程解的相似结构理论以及解的相似构造法[2-17]。 Thomson 方程被誉为表面化学四大方程之一[18],其边值问题尤其在化学工业领域中得到广泛应用,如:石 油化工,感光化学等。为使有关 Thomson 方程的边值问题能在实际工程中的应用变得更加方便,本文就 Thomson 方程边值问题[19] Copyright © 2012 Hanspub 41 范聪银 等 一类 Thomson 方程边值问题解的相似构造法 222 i0 1 0 xa xb xy xyxy EyEF yD Gy Hy (1) (DEFGHab 、、、、、、、为已知常数且 22 2 0,0, 0, i1DGH ab )进行研究。首先对 Thomson 方程 进行变量替换找到该方程的两个线性无关解并用这两个解来构造引解函数;其次求边值问题(1)在非齐次左边值 条件一种特殊情形下的解,得到相似核函数(用引解函数和齐次右边值条件中的系数表示出来);再次求出边值问 题(1)的解(用相似核函数和非齐次左边值条件中的系数组装得到)。最后给出相似构造法的具体步骤,并举例说 明该方法及其优越性。本文的研究不仅有利于 Thomson 方程边值问题在实际工程中的应用,还能进一步完善微 分方程解的相似结构理论以及扩宽相似构造法的应用范围。 2. 预备知识 引理 1:Thomson 方程可经过变量替换[19] i x ,化为: 2 222 2 dd 0 d d yy y 证明:若对Thomson 方程的自变量 x 作替换 i x ,因 22 22 dd dd i, i dd dd yy y yy xx y 则Thomson 方程可化为: 2 222 2 dd 0 d d yy y (2) (2)式是以自变量为 的变型Bessel 方程。 引理 2:Thomson 方程的通解可以写为: iyAIx BKx i (3) 其中: A B、是任意实数, 分别为 IK 、 阶的第一和第二类变型Bessel 函数。 证明:因为变型Bessel 方程(2)的通解可以表示为[19]: yAI BK 则由引理 1知,令 i x ,即可得到Thomson 方程的通解为: iiyAIx BKx 为能方便简洁地把边值问题(1)的解表示出来,下面引入引解函数并给出其微分性质。 引理 3:作二元函数: 1 ,,, 1 mn mnm nmn x y tKtxItyItxKty (4) 分别令二元函数(4)式中的参数 i,1;,tm n ; 1 ,即可得到引解函数: ,,,iiiii , x yKxIy IxKy (5) 1, 11 ,,iiiii , x yKxIyIxKy (6) ,11 1 ,,iiiii , x yKxIy IxKy (7) Copyright © 2012 Hanspub 42 范聪银 等 一类 Thomson 方程边值问题解的相似构造法 1, 11111 ,,iiiii , x yKxI xI xK x (8) 且它们具有如下的微分性质: ,, 1, ,,i,,i i,,i,xy xyxy xx (9) ,,,1 ,,i,,i i,,i,xy xyxy yy (10) ,1 ,1,1 ,,i,,i i,,i,xy xyxy xx (11) 1,,1, 1 ,,i,,i i,,i,xy xyxy yy (12) 证明:对(4)中的 进行赋值易得(5)~(8)式。 mnt、、 下面证明引解函数的微分性质。 根据变型Bessel 函数的微分性质[20]: 1 d d I tItI tt t; 1 d d K tKtKt tt 对(5)式两边的自变量 x 求导得: ,1 1 ,1, ,,ii iiii iii ,,i i,,i x yKxKxIyI xIxKy xx x xy xy x 即(9)式得证。 同理可证(10)~(12)式。 3. 主要定理及其证明 先求边值问题(1)在非齐次左边值条件的一种特殊情形下的解,它将是边值问题(1)的解的相似结构中的相似 核函数。 定理 1:若边值问题 222 i0 1 0 xa xb xy xyxy y Gy Hy (13) (其中参数的限制与(1)中的相同)有解,则其解是唯一的且可以表示为: ,,,1 2 ,1,,1,, 11, 1 ,,i,,ii,,i ,,ii ,,i,,ii ,,ii ,,ii,,i yx GxyHxyxy b Gxyxy Hxyxyxyxy aabba (14) Copyright © 2012 Hanspub 43 范聪银 等 一类 Thomson 方程边值问题解的相似构造法 证明:由引理 2知Thomson 方程的通解为(3)式,故 11 diiiiiii d yAIxBKxAIx IxBKx Kx xxx i (15) 由1 xa y 可得: 11 ii iiiiiiAI aBKaAIa IaBKa Ka aa =1 (16) 再由边值问题(13)中的齐次右边界条件 0 xb Gy Hy 得: 1 1 1 1 ii i iiii iiii i 0 A GIbHIbBGKbHKb AGIb HbIbbIb BGKbHbKb bKb (17) (16)、(17)式是关于待定系数 , A B的线性方程组,其系数矩阵为: ii ii i Ia Ka GIbHI bGK bHK b R i 简记为 ,不难得到其增广矩阵 11 12 21 22 RR RR R11 12 21 22 1 0 RR RR R。由边值问题(13)的解的存在性可得 12rr或RR 。 首先,证明边值问题(13)的解是唯一的。 假设边值问题(13)的解不唯一,那么有 1rr RR。现分情况进行讨论。 ①当矩阵 R 有零行时,由变型 Bessel函数的零点存在定理[20]可知:不可能同时为 0,故只讨论 的情况。由 11 12 RR、 2122 0RR 2122 0RR ,可得到如下关于 的线性方程组: GH、 ii ii GI bHI b GK bHK b 0 0 (18) 容易得到方程组(18)的系数行列式为: 1 ii ii I bIb K bKb 根据变型Bessel 函数 Wronski 公式[20]知: 1 10 ib ,则齐次线性方程组(18)只有零解,即,这与 相矛盾。 0GH 22 0GH ②当都不为 0时,因 11 122122 RRRR、、、 1rR ,故有矩阵 R 的两行或两列对应成比例。 当两行对应成比例时,不妨设比例系数为 0 ,那么经过初等行变换后增广矩阵R可变为: 1112 1 00 RR R 显然此时 rrRR,这与 1rrRR矛盾。 当两列对应成比例时,由于矩阵 R 的两列对应成比例,那么 1r R, 2r R,这与 1rrRR矛盾。 Copyright © 2012 Hanspub 44 范聪银 等 一类 Thomson 方程边值问题解的相似构造法 ③如果矩阵 R 出现零列,那么增广矩阵R只出现下面两种情况: 12 11 22 21 01 0 01 0 RR RR 或RR 1 0 再由①中的讨论可知,这与边值问题(13)的解的存在相矛盾,所以矩阵 R 不会出现零列。 由①、②、③可知 2rr RR,即边值问题(13)的解是唯一的。 其次,证明边值问题(13)的解可以表示为(14)式。 由边值问题(13)有唯一解知,关于待定系数, A B的线性方程组(16)、(17)的系数行列式,再由引理 3, 并经化简得: 0 ,1, 2 ,,11,1,1 ,, ii,, i ,, ii,, i,,ii,, i Gab ab a Habababab ab ab (19) 利用 Cramer 法则求得: 1 1iii Hi A GK bK bHKb b (20) 1 1iii H BGIbIbHIb b i (21) 现在把(20)、(21)式代入(3)式,可得边值问题(13)的解: 1 1 ii iiii iii xI xK x HH i y GKbKbHKbGIbIbH Ib bb 再利用(19)式以及引理 3,即可得到(14)式。 定理 2:若边值问题(1)有解,则其解是唯一的且可以表示为: 11 1 yD x Fa EFa (22) 证明:与定理 1同理可证边值问题(1)的解是唯一的。此处略。 下面证明边值问题(1)的解可以表示为(22)式。 由引理 2知Thomson 方程的通解为(3)式,由边值问题(1)的非齐次左边在值条件 可 得: 1xa EyEF yD 1 1 i1i iii1i ii A EIaEFIaIaB EKaEFKaKaD aa (23) 由边值问题(1)的解是唯一的知,关于待定系数 , A B的线性方程组(17)、(23)的系数行列式 ,再利用引 理3,并经化简得: 0 ,,,1 ,1, 2 ,1,,11,1 ,,i,,i i,,i1,,ii,,i ,, ii,, ii,, ii,, i E GabHababEFGabab ba Hababab ab ab ba (24) Copyright © 2012 Hanspub 45 范聪银 等 一类 Thomson 方程边值问题解的相似构造法 利用 Cramer 法则求得: 1 1iiii A GKbKbH Kb b (25) 1 1iii H BGIb IbHIb b i (26) 把(25)和(26)式代入(3)式可得边值问题(1)的解: 1 1 i iii i iiii xI xH yGKbKbHK b Kx H GI bI bHIb b ib 再利用(24)式以及引理 3,并经整理后得到(22)式。 由定理 1和定理2,可得到如下在实际应用中十分有用几个的推论。 推论 1:在边值问题(1)或(13)中,若齐次右边界条件为 0yb (即0, 0HG ),则相应的相似核函数: , ,1, ,, i ,, ii,, i xb x ab ab a (27) 推论 2:在边值问题(1)或(13)中,若齐次右边界条件为 0yb (即0, 0GH ),则相应的相似核函数: ,,1 2 ,1,,11,1 ,, i,, i ,, ii,, ii,, ii,, i xbi xb b x ab ab abab ab ba (28) 推论 3:边值问题(1)的解式(22) 的结构中的第一个连分式有如下性质: 1 xa D yxFy x E F a (29) (29)式反映了解在非齐次左边值处本质性的特征,在实际应用中起着十分重要的作用。 4. 相似构造法步骤 经过对引理 2、定理 1和定理 2的分析知,相似核函数实质上是边值问题(1)在一种特殊的非齐次左边值条 件 1 xa y 下的解且该解可由引解函数 ,,,i, 1;,1 mn xy mn , 和齐次右边值条件中的系数表 示出来,而边值问题(1)的解又可由相似核函数以及非齐次左边值条件中的系数进行组装得到。 根据以上分析可建立求解Thomson 方程边值问题的一种新方法(相似构造法),其具体步骤如下: 第一步:根据(5)~(8)式,由 Thomson 方程的两个线性无关的解 i i I xK x 、,作引解函数: ,,,i, 1;, 1 mn xy mn 。 第二步:根据(14)式,由引解函数及齐次右边值条件 0 xb Gy Hy 的系数 生成相似核函数:GH、 x , 并计算 。 a 第三步:根据(22)式,由非齐次左边界条件 1xa EyEF yD 的系数 DE 与相似核函数进行组装 而得到边值问题(1)的解。 F、、 实例 求解边值问题: Copyright © 2012 Hanspub 46 范聪银 等 一类 Thomson 方程边值问题解的相似构造法 22 1 2 i1 1 20 x x xy xyxy yx yx 0 (30) 解:第一步由 11 i i I xK x、是边值问题(30)定解方程的两个线性无关解,根据(5)~(8)式作引解函数: 1,11 111 2,12 121 1,21 212 2,22 222 ,,iiiii , ,,iii+ii , ,,iii+ii, ,,iiiii , xyKxIy IxKy xyKxIy IxKy xyKxIy IxKy x yKxIyIxKy 第二步由(14)式构造相似核函数: 111 11212 111 12121 1 12122222 i2ii2i 2ii2i 2ii2i i 2ii2iii2ii2i 2ii2ii2i2ii 2ii2i xKxI IxKKxIIxK KI IKKxI IxK KI IKKI IK 111 11212 111 12121 1 121 22222 1i2ii2i 2ii2i 2ii2i i 2ii2iii2ii2i +2ii 2ii2i2ii 2ii2i KI IKKIIK KI IKKxI IxK KI IKKI IK 第三步由(22)式可得边值问题(30)的解: 111 11212 111 12121 1 121 22222 i2ii2i 2ii2i 2ii2i i 2ii2iii2ii2i +2ii 2ii2i2ii 2ii2i yK xIIxKK xII xK KI IKKxIIxK KI IKKI IK 5. 结论 根据相似构造法的步骤知:只需要知道边值问题的定解方程的两个线性无关解,首先构造二元函数并由二 元函数生成引解函数,其次由引解函数以及非齐次左边界条件中的系数构造相似核函数,最后由相似核函数以 及其次右边界条件中的系数进行组装得到边值问题的解。应用相似构造法解该类边值问题可以把复杂的问题简 单化,到达事半功倍的效果。 一方面相似构造法能够避免繁琐的求解过过程,另一方面该方法还有利于计算机程序员编程;同时它不仅 是微分方程边值问题解的理论创新,更是一种思想创新。另外,从最终构造出来的解的结构式看很容易发现边 值问题的解与边界条件中的系数之间的关系,这对于该类边值问题在工程中的应用有着非常重要的实际意义。 参考文献 (References) [1] 闵嗣鹤, 严士健. 初等数论(第三版)[M]. 北京: 高等教育出版社, 2003. [2] 李顺初, 伊良忠, 郑鹏社. 微分方程定解问题解的相似结构[J]. 四川大学学报(自然科学版), 2006, 43(4): 933-934. Copyright © 2012 Hanspub 47 范聪银 等 一类 Thomson 方程边值问题解的相似构造法 Copyright © 2012 Hanspub 48 [3] 贾闽惠, 李顺初. Bessel方程的边值问题解的相似结构[J]. 大学数学, 2005, 21(5): 37-39. [4] 李顺初. 二阶齐次线性微分方程的边值问题的解的相似结构[J]. 西华大学学报(自然科学版), 2009, 28(5): 40-41, 90. [5] 严娟, 李顺初, 邢承林. 二阶欧拉方程的一类边值问题的相似结构[J]. 西华大学学报(自然科学版), 2009, 28(6): 86-88. [6] 李全勇, 李顺初, 迟颖. 常系数齐次线性微分方程组边值问题解的相似结构[J]. 四川兵工学报, 2010, 31(4): 126-129. [7] 李顺初. 微分方程解的相似结构的初探与展望[J]. 西华大学学报(自然科学版), 2010, 29(2): 223-226, 238. [8] 李顺初, 郑鹏社, 张宇飞. 复合油藏试井分析解的相似结构[J]. 数学的实践与认识, 2008, 38(3): 23-28. [9] 唐乙斌, 李顺初, 严娟, 李全勇. Tschebyscheff 方程边值问题解的相似结构[J]. 四川兵工学报, 2011, 32(1): 155-156. [10] J. D. Tian, S. C. Li. The formal similarity of solutions in the Laplace space on the class of quasilinear partial differential equation. Mathemati- cal Theory and Applications, 2004, 24(2): 66-73. [11] M. H. Jia, S. C. Li. The formal similarity of solutions in the Laplace space on the class of fluid flow differential equation. Journal of Electronic Science and Technology of China, 2005, 3(2): 172-174. [12] J.-P. Su, S.-C. Li and C.-J. Li. The similar of solutions in the Laplace space of composite parabolic partial differential equation. Journal of Zaozhuang University, 2009, 26(2): 6-11. [13] 李顺初. 一类偏微分方程组的 Laplace 空间解的形式相似性[J]. 西华大学学报(自然科学版), 2007, 26(4): 83-86. [14] 李顺初, 郑鹏社, 张宇飞. 复合油藏试井分析解的相似结构[J]. 数学的实践与认识, 2008, 38(3): 23-28. [15] 李顺初, 郑鹏社, 张宇飞. 合采油藏试井分析解的相似结构[J]. 高校应用数学学报, 2009, 24(2): 234-238. [16] 李顺初, 张建军. 分形双孔介质油藏试井分析解的相似结构[J]. 西华大学学报(自然科学版), 2006, 25(1): 40-43. [17] W. Adamson. A textbook of physical chemistry. New York: Academic Press, 1973: 308-315. [18] 刘式适, 刘式达. 特殊函数(第二版)[M]. 北京: 气象出版社, 2002. [19] 李顺初, 黄柄光. Lpalace 变换与 Bessel 函数及试井分析理论基础[M]. 北京: 石油工业出版社, 2000. [20] 郭玉翠. 数学物理方法学习指导[M]. 北京: 清华大学出版社, 2006. |