姚老师最爱的两招：表格法与微分算子法
1 姚老师最爱的两招：表格法与微分算子法，因为效率高，所以喜欢，仅此而已！录入可是字字辛苦，希望大家珍惜哦！
一、 分部积分的表格法
分部积分主要针对被积函数为两类函数乘积的类型，主要可以归纳为反幂、对幂、幂三、幂指和三指五种，幂可以扩展为多项式函数，三主要指正弦和余弦两类三角函数，基本原则是把其中一类函数拿去凑微分，遵循“反对幂三指”、越往后越先凑微分的原则，前四种称为“终止模式”，最后一种称为“循环模式”。当涉及到幂函数（多项式函数）次数较高时，需多次用到分部积分，计算较繁且易出错，因此介绍一个推广公式：
定理：设(),()uuxvvx有1n阶连续导数，则
(1)()(1)(2)(3)1(1)''''''(1)nnnnnnnuvdxuvuvuvuvuvdx。（此定理及证明可略，仅告诉大家，我不是瞎编乱造，而是有理论依据的！）
【证：用数学归纳法。
当0n时，''uvdxuvuvdx。
设1nk时，(1)()(1)(2)(3)1(1)''''''(1)kkkkkkkuvdxuvuvuvuvuvdx
（*）
则当1nk时，(2)(1)(1)(1)'kkkkuvdxudvuvuvdx，
将上式的'u（*）式中的u，则有
(1)()(1)(2)1(2)'''''''(1)kkkkkkuvdxuvuvuvuvdx，
从而(2)(1)()(1)(2)2(2)''''''(1)kkkkkkkuvdxuvuvuvuvuvdx，得证。】
上述式子并不好记，它的一个直观表达就是表格法，如下表。
u的各阶导数 (1)''''''nuuuuu
(1)nv的各阶原函数 (1)()(1)(2)nnnnvvvvv
下面通过例子给予演示：
（1）“幂三”型
例1.1 52(325)cosxxxxdx
解： 2 52325xxx 4562xx 3206x 260x 120x 120 0
cosx sinx cosx sinx cosx sinx cosx
所以原式=
52432325sin562cos206sin60cos120sin120cosxxxxxxxxxxxxxxC
（2）“幂指”型
例1.2 432(21)xxxedx
解：
4321xx 3246xx 21212xx 2412x 24 0
2xe 22xe 24xe 28xe 216xe 232xe
所以原式=433222214612122412242481632xxxxxxxxeC
=4322123312xxxxxeC
（3）“反幂”型（尤其是反三角函数次数高于1时）
例1.3 2arcsinxxdx
解：令arcsinux，则sinxu，cosdxudu
所以原式=221sincossin22uuuduuudu
2u 2u 2 0
sin2u 1cos22u 1sin24u 1cos28u
从而原式=2111cos2sin22cos2488uuuuuC
2222111(arcsin)(12)1arcsin(12)428xxxxxxC。
（4）“对幂”型（尤其是对数函数次数高于1时）
例1.4 34(ln)xxdx
解：令lnxu，则uxe，udxedu， 3 故原式= 343444(ln)uuuxxdxeueduuedu（这是幂指类型了，用表格法自己解解看吧！）
= 43241333(lnlnlnln)44832xxxxxC。
（5）“三指”型（此为循环模式，想想与前面的有何不同？）
例1.5 2sinxxedx
解：
sinx cosx sinx
2xe 22xe 24xe
所以有222111sin(sincos)sin244xxxxedxxxexedx，
求解得22411sin(sincos)524xxxedxxxeC。
二、 求n阶常系数非齐次线性微分方程特解的微分算子法
n阶常系数非齐次线性微分方程为：
()(1)(2)121'()0nnnnnypypypypyfx （**）
求解非齐次方程（**）的特解*()yx常有三种方法：待定系数法、常数变易法和微分算子法。常数变易法在教材中一阶非齐次线性微分方程中已有介绍，待定系数法在二阶非齐次微分方程中着重讲解，因此，在此，主要讲微分算子法。
首先引进记号：222,,,nnndddDDDdxdxdx，
于是22()2','',,nnnndydydyyDyyDyyDydxdxdx，
（**）式变为111()()nnnnDpDpDpyfx。
记111()nnnnFDDpDpDp，
于是()()FDyfx，从而得特解*()1()()()fxyfxFDFD。 4 下面关键是要弄清楚1()FD这个算子是如何进行作用的呢？
通常()fx为幂、三、指、幂三、幂指、三指和幂三指几种类型，下面分别讨论（主要采用书上的例子）。
（1）幂：()fx为多项式函数，采用多项式除法进行计算，什么是多项式除法呢？
例2.1 2''3'25yyyx
解：原方程的一个特解为2*2532xyDD。
22222323137+2482313112231223932447344DDDDDDDDDDDDD
（后面可以继续写下去，但是想想，函数2()5fxx，还有必要吗？）
从而*22137+(5)248yDDx221371317(5)22248224xxxx。
（2）指：()xfxe，当不是特征方程的根时，将直接代入分母的D；当是特征方程的单根时，将分子乘以一个x，分母对D求导数，然后将代入分母；当是特征方程的二重根时，将分子乘以一个2x，分母对D求两阶导数，然后将代入分母。
例2.2.1 2''2'3xyyye
解：222*222322235xxxeeeyDD（因为2不是特征方程2230rr的根。） 5 例2.2.2 ''2'3xyyye
解：*223222124xxxxexexexeyDDD（因为1是特征方程2230rr的单根。）
例2.2.3 ''2'xyyye
解：2*221222xxxexexeyDDD（因为1是特征方程2210rr的二重根。）
（3）三：()fx为正弦函数或余弦函数时，由于欧拉公式连接了正、余弦函数，所以正、余弦函数可以转化为指数函数来求解，一般采用定理3.5进行求解。
例2.3 ''2'cosyyyx
解：cosRe()ixxe（实部）
所以先求解方程''2'ixyyye，
此时*1221(cossin)sincos21212222ixixixeeieiiyxixxxDDii，
故原方程的一个特解**11Re()sin2yyx。
（4）幂三、幂指、三指或幂三指，基本原则：三角函数看成复数域内指数函数的实部或虚部，从而转化成幂指类型，将指数函数提前，后面的算子中D换成D。
例2.4.1 3''6'95(1)xyyyxe（幂指）
解：3*332225(1)5(1)5(1)69(3)6(3)9xxxxexxyeeDDDDD（注：21D表示两次不定积分）
2133321215()552()62xxxxceexxcxcD
由于取的是一个特解，所以12,cc可以随便取，不妨取为0吧，从而得一个特解
*32355()62xyxxe。 6 例2.4.2 ''sin2xyyex（三指）
解：2(12)sin2Im()Im()xxixixexeee，
所以先求解(12)''ixyye，
解得(12)(12)*(12)1221(1)1(12)18ixixixeeyieDi
1(1)cos2sin2(cos2sin2)(cos2sin2)88xxeiexixxxixx
所以原方程的一个特解为**1Im()(cos2sin2)8xeyyxx。
例2.4.3 ''sin2yyxx（幂三）
解：2sin2Im()ixxxxe，所以先求解2''ixyyxe，
解得2*1222221(2)154ixixixxxeeDiiDDxeyD(这里用了多项式除法)
2144()(cos2sin2)()525525ixixieDxxix
44cos2sin2sin2cos2525525xxxxixx
所以原方程的一个特解为**14sin2coIm()s2525xxyxy。
例2.4.4 ''2'5sin2xyyyxex（幂三指）
解：2(12)sin2Im()Im()xxixixxexxeexe，
所以先求解(12)''2'5ixyyyxe，
解得
(12)*(12)(12)(12)122225(12)2(12)54(4)ixixixixxexxxyeeeDDDiDiDiDDDi(12)(12)(12)2(12)21111()111416416()()216216ixixixixDxxiiieeexxexxDDi（注：1D表示求不定积分，且用到了多项式除法哦！）
系数线性微分方程特解的微分算子法
第２８卷第４期
２００８年Ｂ月桂林电子科技大学学报
ＪｏｕｒｎＵ“ＣｕｍⅡＵⅡｊＶｅ噶ＩｔｙｏｆＥＩ睇ｔｍｎＩｃＴｅｃｈｎｏＩ蟾ｙＶ０１．２８。Ｎｏ，４
Ａｕｇ．２００８
二阶常系数线性微分方程特解的微分算子法
李绍刚，徐安农
（桂林电子科技大学数学与计算科学学院，广西桂林５４１００４）
摘要：微分算子法是求解常系数非齐次线性微分方程特解的有效方法，基于算子多项式的理论，针对二阶常系数
线性微分方程，论文给出了非线性项为指数函数、三角函数、幂函数及其混合函数的微分算子特解公式，实例表明
特解公式在解题中具有可应用性、有效性和简捷性．
关键词：线性微分方程，常系数；微分算子Ｉ特解
中图分类号：０１７５．１文献标识码：Ａ文章编号：１６７３—８０８Ｘ（２００８）０４一０３３０－０４
Ｄｉｆｆｅｒｅｎｔｉａｌｏｐｅｒａｔｏｒｍｅｔｈｏｄｆｏｒｐａｒｔｉｃｕｌａｒｓｏｌｕｔｉｏｎｆｏｒ
ｓｅｃｏｎｄ—ｏｒｄｅｒｃｏｎｓｔａｎｔｃｏｅｆｆｉｃｉｅｎｔｌｉｎｅａｒｄｉｆｆｅｒｅｎｔｉａＩｅｑｕａｔｉｏｎ
￡ｆＳ＾口ｏ—９４矗ｇ，ＸＵＡ忍一再。靠ｇ
（ｓｃｈ∞ＩｏｆＭ·ｔｈｅｍｆ西ｔｎｄｃｏｍｐｕｎｔｉｏ＂Ｉｓｃｉ呱ｅ．６ＩＩｉＩｉｎ
ＵｌＩｉ＂巧ｉｔｙｏｆＥＩ鳅∞ｌＩｉｃＴ鼬叫ｏｇｙ，白ｉ陆５４ｌ００４·Ｃｈ．ｍ）
Ａｂｓｔｒ址ｔ：ＤｉｆｆｅｒｅｎｔｉａｌｏｐｅＭｔｏｒｍｅｔｈｏｄ妇蛆ｅｆｆｅｃｔｉｖｅａｐｐｒ∞ｃｈｆｏｆ８０ｌｙｉｎｇｉＩＩｈｏｍｏｇｅｎ∞ｕｓｌｉｎｅａｒｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａ－
ｔｉｏｎｗｉｔｈｃｏｎｓｔａｎｔｃｏｅｆ“ｃｉｅｎｔｓ．１ｈｓｅｄｔｈｅｔｈｅｏｒｙＤｆｏｐｅｒａｔｏｒｐｏＩｙｎｏｍｉａｌａｎｄａｉｍｉｎｇａｔｓｅｃｏｎｄｏｒｄｅｒｉｒＩｌｌｏｍｏｇｅ—
ｌｉｎｅａｒｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｗｉｔｈｃｏｎｓｔａｎｔｃｏｅｆｆｉｃｉｅｎｔｓ，ｄｉｆｆｅｒｅｎｔｉａｌｏｐｅｒａｔｏｒｐａｒｔｉｃｕｌａｒｓｏｌｕｔｉｏｎｓｆｏｒｍｕｌａ
二阶常系数线性微分方程特解的微分算子法
第２８卷第４期 ２００８年８月 桂林电子科技大学学报 Ｊｏｕｒｎａｌ ｏｆ Ｇｕｉｌｉｎ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｅｌｅｃｔｒｏｎｉｃ Ｔｅｃｈｎｏｌｏｇｙ ＶｏＩ．２８，Ｎｏ．４ Ａｕｇ．２００８
二阶常系数线性微分方程特解的微分算子法
李绍刚，徐安农
（桂林电子科技大学数学与计算科学学院，广西桂林５４１００４）
摘要：微分算子法是求解常系数非齐次线性微分方程特解的有效方法，基于算子多项式的理论，针ｘ￣－－阶常系数 线性微分方程，论文给出了非线性项为指数函数、三角函数、幂函数及其混合函数的微分算子特解公式，实例表明 特解公式在解题中具有可应用性、有效性和简捷性。 关键词：线性微分方程；常系数；微分算子；特解 中图分类号：Ｏ１７５．１ 文献标识码：Ａ 文章编号：１６７３—８０８Ｘ（２００８）０４—０３３０—０４
Ｄｉｆｆｅｒｅｎｔｉａｌ ｏｐｅｒａｔｏｒ ｍｅｔｈｏｄ ｆｏｒ ｐａｒｔｉｃｕｌａｒ ｓｏｌｕｔｉｏｎ ｆｏｒ
ｓｅｃｏｎｄ－ｏｒｄｅｒ ｃｏｎｓｔａｎｔ ｃｏｅｆｆｉｃｉｅｎｔ ｌｉｎｅａｒ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｑｕａｔｉｏｎ
ＬＩ Ｓｈａｏ—ｇａｎｇ，ＸＵ Ａｎ—ｎｏｎｇ
（Ｓｃｈｏｏｌ ｏｆ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｃｏｍｐｕｔａｔｉｏｎａｌ Ｓｃｉｅｎｃｅ，Ｇｕｉｌｉｎ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｅｌｅｃｔｒｏｎｉｃ Ｔｅｃｈｎｏｌｏｇｙ・Ｇｕｉｌｉｎ ５４１００４，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｄｉｆｆｅｒｅｎｔｉａｌ ｏｐｅｒａｔｏｒ ｍｅｔｈｏｄ ｉｓ ａｎ ｅｆｆｅｃｔｉｖｅ ａｐｐｒｏａｃｈ ｆｏｒ ｓｏｌｖｉｎｇ ｉｎｈｏｍｏｇｅｎｅｏｕｓ ｌｉｎｅａｒ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｑｕａ— ｔｉｏｎ ｗｉｔｈ ｃｏｎｓｔａｎｔ ｃｏｅｆｆｉｃｉｅｎｔｓ．Ｂａｓｅｄ ｏｎ ｔｈｅ ｔｈｅｏｒｙ ｏｆ ｏｐｅｒａｔｏｒ ｐｏｌｙｎｏｍｉａｌ ａｎｄ ａｉｍｉｎｇ ａｔ ｓｅｃｏｎｄ ｏｒｄｅｒ ｉｎｈｏｍｏｇｅ— ｎｅｏｕ￥ｌｉｎｅａｒ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｑｕａｔｉｏｎ ｗｉｔｈ ｃｏｎｓｔａｎｔ ｃｏｅｆｆｉｃｉｅｎｔｓ，ｄｉｆｆｅｒｅｎｔｉａｌ ｏｐｅｒａｔｏｒ ｐａｒｔｉｃｕｌａｒ ｓｏｌｕｔｉｏｎｓ ｆｏｒｍｕｌａ ａｒｅ ｇｉｖｅｎ ｗｈｅｒｅ ｔｈｅ ｎｏｎｌｉｎｅａｒ ｉｔｅｍ ｈａｓ ｓｅｖｅｒａｌ ｔｙｐｅｓ ｏｆ ｆｕｎｃｔｉｏｎ ｓｕｃｈ ａｓ ｅｘｐｏｎｅｎｔｉａｌ，ｔｒｉｇｏｎｏｍｅｔｒｙ，ｐｏｗｅｒ，ｍｉｘｔｕｒｅ． Ｔｈｅ ｅｘａｍｐｌｅｓ ｐｒｏｖｅｄ ｔｈａｔ ｔｈｅ ｐａｒｔｉｃｕｌａｒ ｓｏｌｕｔｉｏｎ ｆｏｒｍｕｌａ ｈａｄ ｔｈｅ ｐｒｏｐｅｒｔｉｅｓ ｏｆ ａｐｐｌｉｃａｔｉｏｎ，ｖａｌｉｄｉｔｙ ａｎｄ ｃｏｎｃｉｓｅｎｅｓｓ ｉｎ ｓｏｌｖｉｎｇ ｐｒｏｂｌｅｍｓ． Ｋｅｙ ｗｏｒｄｓ：ｌｉｎｅａｒ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｑｕａｔｉｏｎ；ｃｏｎｓｔａｎｔ ｃｏｅｆｆｉｃｉｅｎｔ；ｄｉｆｆｅｒｅｎｔｉａｌ ｏｐｅｒａｔｏｒ；ｐａｒｔｉｃｕｌａｒ ｓｏｌｕｔｉｏｎ
算子法解微分方程
常系数非齐次线性微分方程的解法有很多，例如笔者的教材（《高
等数学 第六版》）所述的待定系数法和接下来给出的称之为“算
子法”以及另一种同样使用算子的方法。
1、首先介绍一种使用算子求解的方法：
考察二阶常系数非齐次线性微分方程
d2x/dt2＋a1dx/dt＋a0x＝b(t)
相应的齐次方程的通解是已知的，所以只须求出方程的一个特解
（由微分方程解的结构给出）。
设该方程的特征多项式q(λ)＝λ2＋a1λ＋a0分解为q(λ)＝(λ
－λ1) (λ－λ2)
则算子多项式q(D)也分解为q(D)＝(D－λ1) (D－λ2)
则原微分方程可写成 (D－λ1) (D－λ2)＝b(t)
依次解以下两个方程
(D－λ2) x1＝b(t)
(D－λ1) x＝x1 就可求得方程的特解。（其中x1看成是中间变量，只要通过求解
x1来求解x）
对于λ1和λ2是共轭虚数的情形，按上述步骤求得的方程特解
有可能是一个复值函数z(t)＝x(t)＋iy(t)。这时应有恒等式
d2z(t)/dt2＋a1dz(t)/dt＋a0z(t)＝b(t)
比较上式两边的实部，我们得到
d2x(t)/dt2＋a1dx(t)/dt＋a0x(t)＝b(t)
这样，不论λ1和λ2是实数或者是共轭虚数，我们都可能够求
出方程在实数范围内的特解，从而完全解决了这方程的求解问题。
给出教材上一个例子:求微分方程y``－5y`＋6y＝xe2x.（《高等
数学》P343）
解：该微分方程的算子多项式分解为 q(D)＝(D－2) (D－3)
设y1＝(D－2)y，代入知(D－3)y1＝xe2x（该式子是一阶常系数微
分方程），易求得
y1＝﹣(x＋1) e2x＋Ce3x（其中C为任意常数）.
所以 (D－2)y＝﹣(x＋1) e2x＋Ce3x.
得 y＝C1e2x＋C2e3x －(x2＋2x) e2x／2.
2、下面来说另一种更简便的方程，也就是“算子法”。不过在使
用算子法的时候，很多性质是必须了解的，在这里不作说明。“算
《有限差分法在微分方程中应用》课件论文
课
程
论
文
有限差分法在微分方程中的应用
本学期学习了《微分方程数值解》，本书中有限差分法给我留下的印象比较深刻，下边说说自己在方面的一点理解，请老师指正。
1.有限差分法的基本思想：
当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将原方程及边界条件中的微分用差分来近似，对于方程中的积分用求和或及机械求积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。
2.有限差分法求解偏微分方程的步骤：
区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格，这些离散点称作网格的节点；
近似替代，即采用有限差分公式替代每一个格点的导数。
逼近求解，换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程。
从原则上说，这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值进行插值计算来近似得到。理论上，当网格步长趋近于零时，差分方程组的解应该收敛于精确解，但由于机器字节的限制，网格步长不可能也没有必要取得无限小，那么差分法的收敛性或者说算法的稳定性就显得至关重要。因此，在运用有限差分法时，除了要保证精度外，还必须要保证其收敛性。
3.构造差分法的几种形式：
主要草用的是泰勒级数展开的方法。其基本差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。其中前两种形式为一阶计算精度，后一种为二阶计算精度。 4.有限差分法的应用:
4.1抛物线形的差分法中的一维常系数抛物线型方程
张宇微分算子法
张宇微分算子法
1. 引言
微分算子是数学中的一个重要概念，它在微积分和偏微分方程等领域有着广泛的应用。在这些领域中，我们经常需要对函数进行求导操作，而微分算子就是用来描述这种操作的工具。张宇微分算子法是一种基于张宇教授提出的方法来求解微分方程的技术。本文将详细介绍张宇微分算子法的原理、应用以及相关实例。
2. 原理
张宇微分算子法是一种基于特征方程和特征根的方法来求解线性常系数齐次线性微分方程的技术。对于给定的线性常系数齐次线性微分方程：
𝑎𝑛𝑦(𝑛)(𝑥)+𝑎𝑛−1𝑦(𝑛−1)(𝑥)+⋯+𝑎1𝑦′(𝑥)+𝑎0𝑦(𝑥)=0
其中，𝑦(𝑥)是未知函数，𝑎𝑖是常数，𝑦(𝑛)(𝑥)表示对函数𝑦(𝑥)求𝑛阶导数。
首先，我们可以将上述微分方程转化为一个特征方程：
𝑎𝑛𝑠𝑛+𝑎𝑛−1𝑠𝑛−1+⋯+𝑎1𝑠+𝑎0=0
其中，𝑠是特征根。解这个特征方程可以得到𝑛个互不相同的特征根𝑠1,𝑠2,…,𝑠𝑛。
然后，我们可以根据这些特征根来构造一个微分算子：
𝐷=(𝐷−𝑠1)(𝐷−𝑠2)…(𝐷−𝑠𝑛)
其中，𝐷表示对函数求导的微分算子。这样，原微分方程的通解可以表示为：
𝑦(𝑥)=𝐶𝑒𝑠1𝑥+𝐶𝑒𝑠2𝑥+⋯+𝐶𝑒𝑠𝑛𝑥
其中，𝐶是常数。
3. 应用
张宇微分算子法在求解线性常系数齐次线性微分方程时具有广泛的应用。它可以解决各种类型的微分方程，包括一阶、二阶、高阶等。
3.1 一阶线性常系数齐次微分方程
对于一阶线性常系数齐次微分方程：
𝑎1𝑦′(𝑥)+𝑎0𝑦(𝑥)=0
我们可以将其转化为特征方程：
𝑎1𝑠+𝑎0=0 解得特征根𝑠=−𝑎0𝑎1。然后构造微分算子：
𝐷=𝐷+𝑎0𝑎1
最后得到通解：
𝑦(𝑥)=𝐶𝑒−𝑎0𝑎1𝑥
3.2 二阶线性常系数齐次微分方程
褶积微分算子法数值模拟的边界问题研究
第３２卷第６期 物探化探计算技术 ２０１０年１１月
文章编号：ｌｏｏ１—１７４９（２０１０）Ｏ６－＿０５７ｌ－＿ｏ７
褶积微分算子法数值模拟的边界问题研究
贺同江 ，刘红艳 ，李小凡２
（１．天津市地震局，天津３００２０１；２．中国科学院地质与地球物理研究所，北京１０００２９）
摘要：褶积微分算子法是一种全新的数值模拟方法，已被广泛应用于复杂介质的地震波场数 值模拟，但是其边界反射问题一直没有解决。这里将最佳匹配层（ＰＭＬ）吸收边界条件引入到褶
积微分算子法中，此方法是在研究区域的边界上加入吸收层，使边界上传入吸收层的波，随传播
距离按指数规律衰减，不产生任何反射，以达到消除边界反射的目的。构造不同的模型，通过对
比分析证明，ＰＭＬ吸收边界条件能比较好地解决褶积微分算子法的边界问题，从而验证了完全
匹配层吸收效果的优越性。
关键词：褶积微分算子；数值模拟；ＰＭＬ吸收边界
中图分类号：Ｐ ６３１．４ 文献标识码：Ａ
０前言
褶积微分算子法作为一种全新的数值模拟方 法，已被广泛应用于复杂介质的地震波场数值模拟
中，该方法同时具有广义正交多项式方法的高精度
和短算子有限差分算法的高速度。通过对算子长 度的调节及算子系数的优化，可同时兼顾波场解的
全局信息与局部信息。该算法的计算速度快，计算
效率高，能够直观、高效地反映介质中波场的传播 规律 “Ｊ。但是在进行波动方程数值模拟时，考虑
到计算机的有限内存及有限计算时间，要对考虑问
题的无限区域进行截取，使进行中的数值模拟在有 限的区域内完成，即可以在有限的区域中得到较高
精度的无限空间解。为此需要引入人工边界来达 到此目的，但这样会在人工边界处产生人为反射，
如不消除或者压制这种虚假反射，就会影响数值模
拟的结果和精度。为此，边界条件问题是必须要解 决的，也是较难解决的问题。地震波数值模拟中， 在人工边界上使用吸收边界条件，应尽可能地降低
边界上所产生的反射波强度，使得反射波在边界上
微分方程几种求解方法
微分方程几种求解方法
微分方程是数学中的重要工具，用于描述自然界中关于变化的数学模型。微分方程的求解方法有多种，可以根据不同的特征和条件选择不同的方法。下面将介绍微分方程的几种常见求解方法。
1.可分离变量法
可分离变量法适用于形如 dy/dx = f(x)g(y) 的一阶微分方程。该方法的基本思路是将变量分离，即将方程写成 dx / f(x) = dy / g(y)，然后两边同时积分，从而得到方程的解。
2.齐次方程法
齐次方程指的是形如 dy/dx = f(x / y) 的一阶微分方程。齐次方程法的基本思路是变量替换，令 y = vx，然后将方程转化为关于 v 和 x
的一阶微分方程，再用可分离变量法求解。
3.线性方程法
线性方程是指形如 dy/dx + p(x)y = q(x) 的一阶微分方程。线性方程法的基本思路是找到一个积分因子，使得原方程变为恰当方程，然后进行积分求解。常见的积分因子有 e^(∫p(x)dx) 和 1 / (y^2)，选择合适的积分因子可以简化计算。
4.变量替换法
变量替换法适用于一些特殊形式的微分方程。通过合适的变量替换，可以将原方程转化为标准的微分方程形式，从而便于求解。常见的变量替换包括令 y = u(x) / v(x)，令 v = dy/dx等。
5.常数变易法 当已知一个特解时，可以利用常数变易法求解更一般的微分方程。该方法的基本思路是令y=u(x)y_0，其中y_0是已知的特解，然后将y代入原方程得到一阶线性非齐次方程，再用线性方程法进行求解。
6.欧拉法
欧拉法是一种数值求解微分方程的方法。它通过在函数的变化区间内分割小区间，并在每个小区间上用直线逼近函数的变化情况，从而得到微分方程的近似解。欧拉法的计算公式为y_(n+1)=y_n+h*f(x_n,y_n)，其中h为步长，f(x,y)为微分方程的右端。
7.泰勒级数法
泰勒级数法是一种近似求解微分方程的方法，利用函数的泰勒级数展开式进行计算。该方法的基本思路是将函数表示为一个无穷级数，并截断级数，保留一定的项数进行计算，从而得到微分方程的近似解。
微分方程解法的十种求法(非常经典)
微分方程解法的十种求法(非常经典)
本文将介绍微分方程的十种经典求解方法。微分方程是数学中重要的概念，广泛应用于物理学、工程学等领域。通过研究这十种求解方法，读者将更好地理解和应用微分方程。
1. 变量可分离法
变量可分离法是最常见和简单的微分方程求解方法之一。该方法适用于形如dy/dx=f(x)g(y)的微分方程，其中f(x)和g(y)是关于x和y的函数。通过将方程两边分离变量，即把f(x)和g(y)分别移到不同的方程一边，然后进行积分，最后得到y的表达式。
2. 齐次方程法
齐次方程法适用于形如dy/dx=F(y/x)的微分方程。通过令v=y/x，将微分方程转化为dv/dx=g(v)，其中g(v)=F(v)/v。然后再使用变量可分离法求解。
3. 线性微分方程法
线性微分方程法适用于形如dy/dx+a(x)y=b(x)的微分方程。通过乘以一个积分因子，将该方程转化为可以进行积分的形式。
4. 恰当微分方程法
恰当微分方程法适用于形如M(x,y)dx+N(x,y)dy=0的微分方程。通过判断M(x,y)和N(x,y)的偏导数关系，如果满足一定条件，则可以找到一个函数u(x,y)，使得u满足偏导数形式的方程，并且通过积分得到原方程的解。
5. 一阶线性常微分方程法
一阶线性常微分方程法适用于形如dy/dx+p(x)y=q(x)的微分方程。通过先求齐次线性方程的通解，然后再利用待定系数法找到特解，最后求得原方程的通解。
6. 二阶常系数齐次线性微分方程法
二阶常系数齐次线性微分方程法适用于形如d²y/dx²+a1dy/dx+a0y=0的微分方程。通过设y=e^(mx)，将微分方程转化为特征方程，然后求解特征方程得到特征根，利用特征根找到原方程的通解。
7. 二阶非齐次线性微分方程法
二阶非齐次线性微分方程法适用于形如d²y/dx²+a1dy/dx+a0y=F(x)的微分方程。通过先求齐次线性方程的通解，再利用待定系数法找到非齐次线性方程的特解，最后求得原方程的通解。
变系数线性微分方程算子解法的推广
２（）Ｉ１年２月 Ｆｅｂ．２Ｏｌｌ１ 汕头大学学报（自然科学版） Ｊｏｕｒｎａｌ ｏｆ Ｓｈａｎｔｏｕ Ｕｎｌｖｅｒｓｉｔｙ（Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅ） 第２６卷第１期 Ｖｔ＇１．２６ ＮＯ．１
文章编号：１００Ｉ一４２１７（ ｌ１ ）０１一ｏ０１ ８—０９
变系数线性微分方程算子解法的推广
方书盛
（汕头市达濠第二中学，广东汕头５Ｉ：ｓ０７１）
摘要：运用微分算子运算的基本原理，证明最高项系数含有函数的变系数线性微分算子可，． 分解为一阶微分算子及其幂的因式的乘积的充要条件，并得到常见的最高项系数含有函数的 二阶变系数线性微分算子可分解的充要条件；应用积分运算的基本方法，推导出由微分算子 分解式求出已知方程通解的方法和步骤，并得到了具体的计算公式；通过实例分析和比较说 明所得到解法的应用及其具有的优越性． 关键词：一般推广；变系数；线性微分算子；算子分解；算子解法 中图分类号：Ｏ Ｉ７５．１ 文献标识码：Ａ
０暑Ｉ言 。～ ， 《 ， 线性微分方程是最常见的微分方程，在科学技术上有着重要的应用．由微分方程的
理论可知，变系数线性微分方程在一般情况下是不可积的 ２］．为了满足理论研究和工
程实践的需要，近３０年来，人们不断探索变系数线性微分方程各种特殊的求解方法和
各种不同的可积类型，并已取得了一系列的成果１３￣９１．但这些解法都是在假设方程最高
项系数不含有函数的情形下得到的，而在一般情况下，这类方程的最高项系数是某一函
数，且通常不能直接化为最高项系数不含有函数的方程．因此， 研究最高项系数含有函 数的这类方程的解法是很有必要的．本文将变系数线性微分方程的算子解法推广到最高
项系数含有函数的一般情形，得到了这类方程的求解方法和可积的条件．
１ 主要结果
ｌ＿Ｊ 变系数线性微分算子的分解 变系数线性微分算子的分解是算子解法的一个重要环节．考虑最高项系数为某一函
数的ｎ阶变系数线性微分方程：
姚老师最爱的两招：表格法与微分算子法
1 / 1 姚老师最爱的两招：表格法与微分算子法，因为效率高，所以喜欢，仅此而已！录入可是字字辛苦，希望大家珍惜哦！
一、 分部积分的表格法
分部积分主要针对被积函数为两类函数乘积的类型，主要可以归纳为反幂、对幂、幂三、幂指和三指五种，幂可以扩展为多项式函数，三主要指正弦和余弦两类三角函数，基本原则是把其中一类函数拿去凑微分，遵循“反对幂三指”、越往后越先凑微分的原则，前四种称为“终止模式”，最后一种称为“循环模式”。当涉及到幂函数（多项式函数）次数较高时，需多次用到分部积分，计算较繁且易出错，因此介绍一个推广公式：
定理：设(),()uuxvvx有1n阶连续导数，则
(1)()(1)(2)(3)1(1)''''''(1)nnnnnnnuvdxuvuvuvuvuvdx。（此定理及证明可略，仅告诉大家，我不是瞎编乱造，而是有理论依据的！）
【证：用数学归纳法。
当0n时，''uvdxuvuvdx。
设1nk时，(1)()(1)(2)(3)1(1)''''''(1)kkkkkkkuvdxuvuvuvuvuvdx
（*）
则当1nk时，(2)(1)(1)(1)'kkkkuvdxudvuvuvdx，
将上式的'u（*）式中的u，则有
(1)()(1)(2)1(2)'''''''(1)kkkkkkuvdxuvuvuvuvdx，
从而(2)(1)()(1)(2)2(2)''''''(1)kkkkkkkuvdxuvuvuvuvuvdx，得证。】
上述式子并不好记，它的一个直观表达就是表格法，如下表。
u的各阶导数 (1)''''''nuuuuu
(1)nv的各阶原函数 (1)()(1)(2)nnnnvvvvv
下面通过例子给予演示：
（1）“幂三”型