Maple6-ch4.2-偏微分方程-不练.docVIP

§4.2 偏微分方程 pdsolve来处理方程（而不能是方程组）的解析解，能解的方程有限，解不出时就返回空解； PDEplot用图形描述解析解的形状，这个函数的功能是有限的，它只能对一阶的偏微分方程进行作图，它在库“PDEtools”中; pdetest 用来验解，使用形式与odetest一模一样. 4.2.1 偏微分方程的解析解 命令： pdsolve(PDE); pdsolve(PDE, f, HINT=…, INTEGRATE, build); 其中： PDE 是偏微分方程 f 是待求函数形式，这与dsolve中的参数是类似的 HINT=… 用来设置函数输出结果的形式，如设成两个独立的函数相乘的结果（如f(x)*g(x)）,f(x)+g(x)）HINT=strip(这种参数只在求一阶偏微分方程的解时有效，最后的结果将自变量和应变量用一个新的参数表示)等等，具体输出形式可根据具体方程和实际需要而定，这是一个可选参数； INTEGRATE 是一个可选参数，设置函数当待求的方程可用分离变量来求时，自动将函数得到的常微分方程集合求解综合起来。 Build 设置函数将最后返回一个显示解。 Ex.1 求一维波动方程(可用分离变量法求解)的解析解 pde1:=pde1: assume(a,positive); pde1:=diff(u(x,t),t,t)-a^2*diff(u(x,t),x,x)=0; pdsolve(pde1); Ex.2 求解一阶偏微分方程的解析解 解 pde2:=pde2: a:=a:b:=b: pde2:=diff(u(x,y),x)-a*diff(u(x,y),y)-b*u(x,y)=0; （1）先不设pdsolve的参数简单求解 sol1:=sol1: sol1:=pdsolve(pde2); pdetest(sol1,pde2); （2）对一阶可分离变量方程，可设参数HINT=f(x)*g(y)求解 ol2:=sol2: sol2:=pdsolve(pde2,u(x,y),HINT=f(x)*g(y)); 由于分离变量法给出的仍是将f(x)、g(y)含在微分方程组中，故需要用“build”函数给出原方程的解，这个函数既可以单独使用（如下），又可以作为“pdsolve”的参数（如sol4）。下面，用此函数求最后的解： sol3:=sol3: with(PDEtools): sol3:=build(sol2); pdetest(sol3,pde2); 验证上面所求的是解 sol4:=sol4: 将“build”作“pdsolve”的参数求解 sol4:=pdsolve(pde2,u(x,y),HINT=f(x)*g(y),build); 此解与sol3相同。 （3）由于这个一阶方程是线性齐次的，故其两个解的和也是解，因此它有f(x)+g(x)形式的解，故可取参数HINT=f(x)+g(y),build求解： sol5:=sol5: sol5:=pdsolve(pde2,u(x,y),HINT=f(x)+g(y),build); pdetest(sol5,pde2); 将此解与sol3或sol4比较可知，不同的参数选法，其解可能是不同的。 （4）对一阶方程，可以将解用图来显示其解的形状： with(PDEtools): a:=1:b:=1: PDEplot(pde2,u(x,y),[t,t,1/2*exp(3*t)],[t=0..1]); 注意：PDEplot这个函数的功能是有限的，它只能对一阶的偏微分方程进行作图。 它的一般格式为： PDEplot (PDE,inits, srange, options) 其中参数意义如下表： PDE 一阶偏微分方程，线性或非线性，但必须是只含一个n元未知函数 inits 一个n+1个表达式或方程组成的有序表，它们表示初值，用参数方式给出，包括n-1个参数，在n=3时，一般用s和t作为参数变量 srange 初值中参数的取值范围 options 辅助的参数，用“参数名=辅助参数”的形式给出 Ex.2` 再求解一阶偏微分方程的解析解，用单位圆作为初值曲线，绘积分曲面： with(PDEtools): PDE:=diff(u(x,y),x)*diff(u(x,y),y)-x*y+u(x,y)=0; PDEplot(PDE,[cos(t),sin(t