微分方程的maple求解.docVIP

数理方程与逻辑函数期中考试论文 PAGE PAGE 12 1、常用函数 1)求解常微分方程的命令dsolve. dsolve(常微分方程) dsolve(常微分方程，待解函数，选项) dsolve({常微分方程，初值}，待解函数，选项) dsolve({常微分方程组，初值}，{待解函数}，选项) 其中选项设置解得求解方法和解的表示方式。求解方法有type=formal_series(形式幂级数解)、type=formal_solution(形式解)、type=numeric(数值解)、type=series(级数解)、method=fourier(通过Fourier变换求解)、method=laplace(通过Laplace变换求解)等。解的表示方式有explicit(显式)、implicit(隐式)、parametric(参数式)。当方程比较复杂时，要想得到显式解通常十分困难，结果也会相当复杂。这时，方程的隐式解更为有用，一般也要简单得多。dsolve为标准库函数。 2)求解一阶线性常微分方程的命令linearsol. 在Maple中求解一阶线性方程既可以用dsolve函数求解，也可以用Detools函数包中的linearsol函数求解。linearsol是专门求解线性微分方程的命令，使用格式为: linearsol(线性方程，待解函数) linearsol的返回值为集合形式的解。 3)偏微分方程求解命令pdsolve. pdsolve(偏微分方程，待解变量，选项) pdsolve(偏微分方程，初值或边界条件，选项) pdsolve为标准库函数，可直接使用。 如果求解成功，将得到几种可能结果： 方程的通解； 拟通解(包含有任意函数，但不足以构造通解)； 一些常微分方程的集合； 2、方法 1)一阶常微分方程的解法 a分离变量法 = 1 \* ROMAN I直接分离变量法。如，方程右端是两个分别只含x或y的函数因式乘积，其通解为。 = 2 \* ROMAN II换元法之后再用分离变量法。对于以为中间变量的函数，如，令u=,则原方程变为，再用分离变量法可得。 b常数变易法 = 1 \* ROMAN I对于线性非齐次方程来说， 线性非齐次方程的通解=它所对应的齐次方程的通解+非齐次方程的一个特解。 如y+P(x)y=f(x)，若f(x)0,y+P(x)y=0为一阶线性齐次方程，其通解为，令代入非齐次方程，求出C(x),再的特解。 = 2 \* ROMAN II对于伯努利方程(非线性一阶)来说，先将其化为线性。 如,两端除以，得，令z=,则原方程可化为。 2)二阶线性常微分方程的解法 a二阶线性齐次方程,y+p(x)y+q(x)y=0 若与是方程的解，且(即线性无关)，则是通解，考虑常系数，即p.q都是常数，y+py+qy=0。其特征方程为。解为，。 = 1 \* ROMAN I0，两个不等实根，且常数时，。 = 2 \* ROMAN II0，一对共轭复根，，， ，常数，。 = 3 \* ROMAN III=0，两个相等实根，，，常数， 。 b二阶常系数线性非齐次微分方程，y+py+qy=r(x). 非齐次方程的通解=它所对应的齐次方程的通解+非齐次方程的一个特解。 利用常数变异法，令其特解为，则 ， 令=0…… = 1 \* GB3 ①，并求出，将并将它们都带入到原方程，得=r(x)…… = 2 \* GB3 ② 联立 = 1 \* GB3 ①， = 2 \* GB3 ②式得。 以上得出了特解，再将其与通解组合可得原方程的解。 c欧拉方程(变系数)，。 令，则，代入得 ，可以求解。 3、作图 1)常微分方程数值解作图命令odeplot 要作出常微分方程数值解的图像，要使用odeplot函数。odeplot在函数抱plots中，可通过with(plots)或plots[odeplot]调出。 odeplot(数值解，被绘函数，参数范围，选项) 2)偏微分方程作图命令PDEplot PDEplot(偏微分方程，初值，参数范围，选项) PDEplot位于PDEtools函数包中，使用前必须先调出PDEtools函数包。 三、各种方程的求解 第一部分:一阶常微分方程 1、可分离变量方程 例1: eq:=diff(y(x),x)=sin(x)/sin(y(x)); DEtools[odeadvisor](eq); dsolve(eq); dsolve(eq,implicit); dsolve({eq,y(0)=1}); dsolve({eq,y(0)=1},numeric,range=-2..2); plots[odeplot](%); 2、齐次方程 例2: eq:=D(y)(x)=y(x)/x+tan(