第二讲Maple初步和应用续.pptVIP

第二部分 Maple 的简单应用 * * 第二讲 Maple 简介和应用 单种群模型的解析求解和稳定性分析 单种群模型的解析求解和稳定性分析 dsolve({diff(x(t),t)=r*x(t)*(1-x(t)/K),x(0)=x0}); dsolve({diff(x(t),t)=x(t)*(r-ln(x(t)/K)),x(0)=x0}); 程序 Logisticsol.mws 单种群模型的解析求解和稳定性分析 solve(r*x*(1-x/K)=0,{x}); 求平衡态： 程序 Logisticequ.mws 强调求解功能强大 第二功能性反应函数与捕食被捕食系统 内部平衡态的稳定性分析 程序HollingLV1.mws 程序pattern.mws DEtools[phaseportrait] #画向量场及积分曲线； ([diff(y(x),x)=-y(x)],y(x), #定义微分方程y’=-y; x=-2..2, #指出x的范围 [[y(-2)=2],[y(-2)=1],[y(-2)=-2]], #给出3个初始值 dirgrid=[17,17], #定义网格点密度 arrows=LINE, #定义线段类型 axes=NORMAL); #定义坐标系类型 DEtools[dfieldplot] ([diff(y(x),x)=x^2-y(x)],y(x), x=-2..2,y=-2..2, dirgrid=[9,9], arrows=LINE, axes=NORMAL); 画向量场及积分曲线 with(DEtools): phaseportrait(cos(x)*diff(y(x),x$3)-diff(y(x),x$2)+Pi*diff(y(x),x)=y(x)-x,y(x),x=-2.5..1.4,[[y(0)=1,D(y)(0)=2,(D@@2)(y)(0)=1]],y=-4..5,stepsize=.05); 寻找帮助命令 DEtools phaseportrait([D(x)(t)=y(t)-z(t),D(y)(t)=z(t)-x(t),D(z)(t)=x(t)-y(t)*2], [x(t),y(t),z(t)],t=-2..2,[[x(0)=1,y(0)=0,z(0)=2]],stepsize=.05, scene=[z(t),x(t)],linecolour=sin(t*Pi/2),method=classical[foreuler]); phaseportrait(D(y)(x)=-y(x)-x^2,y(x),x=-1..2.5,[[y(0)=0],[y(0)=1],[y(0)=-1]],title=`Asymptotic solution`,colour=magenta,linecolor=[gold,yellow,wheat]);