二阶常微分方程边值问题求解的常数变易法.docx

二阶常微分方程边值问题求解的常数变易法.docx

《二阶常微分方程边值问题求解的常数变易法.docx》由会员分享，可在线阅读，更多相关《二阶常微分方程边值问题求解的常数变易法.docx（3页珍藏版）》请在冰点文库上搜索。

二阶常微分方程边值问题求解的常数变易法

二阶常微分方程边值问题求解的常数变易法

数学物理方程与特殊函数复习资料二阶常微分方程边值问题求解的常数变易法20XX年-8-31数理方程所解决的问题与高等数学（微积分）教科书中的常微分方程有很大区别，其中最显著的特点是多数微分方程的条件是边值问题，即知道未知函数在自变量变化区域的边界上的取值。

这就是所谓的边值问题。

最简单的是二阶常微分方程的两点边值问题。

二阶常微分方程的解是一个一元函数，关于这个一元函数的信息，知道的不多，除了微分方程本身提供的之外，还有未知函数在一个区间的两个端点处的值。

微积分所教给我们的技巧是先求出常微分方程的通解，再根据两个条件确定通解中的两个任意常数。

进入这门课之初，先回顾初值问题，再思考边值问题。

在边值问题中，数理方程课程内容中出现了一个历史上非常著名的函数，即格林函数。

对力的分析中普遍使用一个方程：

F=ma。

这是著名的牛顿第二定律，其中，F表示力，m表示物体的质量，而a表示物体运动的加速度。

由于加速度的物理意义可解释为物体运动时位移变量对时间的二阶导数，再结合使用虎克定律，就得出简单的振动所满足的二阶常微分方程y2y0如果考虑外力作用，该方程化为更一般的情况y2yf（x）y（0）,y（0）两个初始条件可解释为已经知道初始位移和初始速度。

求解上面方程需要用常数变易法。

先回顾一阶常微分方程求解的方法，然后再讨论二阶常微分方程的常数变易法。

一、一阶常微分方程初值问题的常数变易法一阶常微分方程常数变易法，用于解源函数不为零的常微分方程问题y（x）ry（x）f（x）,x0y（0）先求解简化的（源函数为零）的方程：

y（x）ry（x）0由分离变量：

dydyrdxry，ydx积分：

lnyrxc，y（x）Cexp（rx）应用常数变易法，假设简化前的方程的解具有与简化后方程的解有相同形式，将常数替换为待定的函数，即y（x）u（x）exp（rx）求导数，得y（x）u（x）exp（rx）ru（x）exp（rx）u（x）exp（rx）ry（x）数学物理方程与特殊函数复习资料将其代入化简前的方程，得等式u（x）exp（rx）f（x），u（x）exp（rx）f（x）积分，得u（x）xexp（r）f（）dC代入表达式y（x）u（x）exp（rx），得y（x）[exp（r）f（）dC]exp（rx）x应用初始条件，得解函数y（x）exp（rx）exp[r（x）]f（）dx从两部分解读解函数的意义。

第一部分利用了初始条件的信息，第二部分利用了微分方程右端项的信息。

它们分别是两个子问题的解yry0yryf（x），y（0）y（0）0二、二阶常微分方程初值问题的常数变易法应用常数变易法解二阶常微分方程初值问题y（x）2y（x）f（x）,x0y（0）,y（0）先考虑化简后的方程：

yy0。

其辅助方程为：

2m220，mi,x）C2sin（x）通解为：

y（x）C1cos（假设化简前的方程的解具有相同形式，将常数替换为待定的函数，即y（x）u（x）cos（x）v（x）sin（x）在上式中，有两个未知函数待定。

如果直接代入微分方程，可产生一个等式，由一个等式不能唯一确定两个函数。

如果人为增加一个等式，就可以构造出二元线性方程组，朗斯基行列式方法是成功的确定两个待定函数的方法，方法如下，对假设的函数求一阶导数，得y[ucos（x）vsin（x）][usin（x）vcos（x）]在上面表达式中，令第一个方栝号为零，产生一个等式ucos（x）vsin（x）0同时，由yusin（x）vcos（x）继续求导数，得y[usin（x）vcos（x）][2ucos（x）2vsin（x）]数学物理方程与特殊函数复习资料代入方程，产生第二个等式usin（x）vcos（x）f将两个等式联立，得线性代数方程组ucos（x）vsin（x）0usin（x）vcos（x）f或写成矩阵形式sin（x）u0cos（x）sin（x）cos（x）vf上式中系数矩阵行列式称为朗斯基行列式，由于利用克莱姆法则解方程组，有cos（x）sin（x）sin（x）cos（x）10fsin（x）cos（x）f（x）sin（x），2cos（x）sin（x）0ff（x）cos（x）u1/1x1f（x）sin（x），v2/1x1f（x）cos（x）积分，得两个待定函数表达式u（x）f（）sin（）dC1，v（x）f（）cos（）dgxC2代入常数变易法假设的函数中，得yC1cos（x）C2sin（x）1cos（x）sin（）f（）dxx1sin（x）cos（）f（）dx利用初始条件确定任意常数C1和C2，显然C1，C2/代入并利用三角函数和差化积公式，得1xy（x）cos（x）sin（x）sin[（x）]f（）d0三、二阶常微分方程两点边值的常数变易法与初值问题相似的二阶常微分方程问题如下y（x）y（x）f（x）,x（0,/2）y（0）,y（/2）与前面问题相比，只有1，且第二初值条件改为了右端边界条件。

利用常数变易法得xxyC1cosxC2sinxcosxsinf（）dxsinxcosf（）d数学物理方程与特殊函数复习资料利用边值条件确定任意常数，得C1，C2代入，整理得x/2cosf（）d/2ycosxsinxcosxsinf（）dsinxxcosf（）d所以，有两点边值问题的解ycosxsinx其中，/2G（x,）f（）dcosxsin,0xG（x,）sinxcos,x/2这一函数被称为格林函数。

另一个简单的两点边值问题如下yf（x）,x（0,1）y（0）,y

（1）这一问题物理背景可解释为：

在高速公路上狂奔的汽车，通过雷达测量得到了从起点到终点每一时刻的加速度f（x），需计算汽车在各个时间的位置。

简化后的方程为：

y0，通解为：

yC1（1x）C2x。

应用常数变易法，令化简前的方程解为yu（x）（1x）v（x）x由于朗斯基行列式为而x1x111所以0fx1xf（x），2x10f（1x）f（x）uxf（x），v（1x）f（x）积分并代入解，得yC1（1x）C2x（1x）f（）dx

（1）f（）dxx利用边界条件，得C1，C2

（1）f（）d1代入整理得数学物理方程与特殊函数复习资料y（1x）x（1x）f（）dx

（1）f（）dxx1所以y（1x）xG（x,）f（）d1其中格林函数为（1x）,0xG（x,）x

（1）,x1如下两点边值问题中的边界条件为零值，即0，0。

则方程的解为yG（x,）f（）d1格林函数的数学意义在于，对于任意的源函数f（x），两点边值问题的解可以通过格林函数和源函数的乘积的积分形式表示出来。

使这一问题的解的存在唯一性变得显然。