最新微分方程作图PPT课件.ppt

《最新微分方程作图PPT课件.ppt》由会员分享，可在线阅读，更多相关《最新微分方程作图PPT课件.ppt（30页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、微分方程作图微分方程作图四川大学数学学院 徐小湛 May 2012with(DEtools):调用微分方程工具调用微分方程工具with(plots):调用绘图工具调用绘图工具de:=diff(y(x),x)=2*x*y(x);定义微分方程定义微分方程fxc:=DEplot(wffc,y(x),x=-2.2,y=-2.2):画斜率场画斜率场jfq:=contourplot(y/exp(x2),x=-2.2,y=-2.2):画积分曲线画积分曲线display(fxc,jfqx);用数学软件用数学软件Maple可以画出微分方程的的积可以画出微分方程的的积分曲线和方向场的图形。分曲线和方向场的图形。画

2、图的命令如下：画图的命令如下：2四川大学数学学院 徐小湛 May 2012四川大学数学学院 徐小湛 May 2012四川大学数学学院 徐小湛 May 2012四川大学数学学院 徐小湛 May 2012四川大学数学学院 徐小湛 May 2012四川大学数学学院 徐小湛 May 2012四川大学数学学院 徐小湛 May 2012微分方程：微分方程：作出微分方程的积分曲线的图形。作出微分方程的积分曲线的图形。原方程化为：原方程化为：9四川大学数学学院 徐小湛 May 2012with(DEtools):with(plots):wffc:=(x-y(x)3)*diff(y(x),x)+y(x)=0:d

3、solve(wffc);fangxiangcang:=DEplot(wffc,y(x),x=-2.2,y=-2.2,thickness=2):jifenquxian:=contourplot(x*y-y4/4,x=-2.2,y=-2.2,contours=20,color=blue,thickness=2):display(fangxiangcang,jifenquxian);10四川大学数学学院 徐小湛 May 2012with(DEtools):DEplot(x-y(x)3)*diff(y(x),x)+y(x)=0,y(x),x=-2.2,y=-2.2,y(0)=1,y(0)=0.3,y(

4、0)=1.5,y(0)=-0.5,y(0)=-1,y(0)=-1.5,linecolor=blue,black,gold,navy,green,maroon,color=violet,stepsize=0.01,scaling=constrained);11四川大学数学学院 徐小湛 May 2012with(DEtools):with(plots):wffc:=diff(y(x),x)=2*x*y(x):dsolve(wffc);fangxiangcang:=DEplot(wffc,y(x),x=-2.2,y=-2.2,thickness=2):jifenquxian:=contourplot

5、(y/exp(x2),x=-2.2,y=-2.2,contours=20,color=blue,thickness=2):display(fangxiangcang,jifenquxian);方向场与积分曲线方向场与积分曲线12四川大学数学学院 徐小湛 May 2012with(DEtools):DEplot(diff(y(x),x)=2*x*y(x),y(x),x=-2.2,y(0)=-1,y(0)=-0.5,y(0)=0,y(0)=0.5,y(0)=1,y(0)=1.5,y=-4.4,linecolor=gold,black,blue,red,brown,green,color=grey,

6、stepsize=0.01,scaling=constrained);方向场与积分曲线方向场与积分曲线13四川大学数学学院 徐小湛 May 2012wffc:=3*x*y(x)2*diff(y(x),x)=x3+y(x)3:dsolve(wffc,implicit);fangxiangcang:=DEplot(wffc,y(x),x=-2.2,y=-2.2,thickness=2):jifenquxian:=contourplot(x2-2*y3/x,x=-2.2,y=-2.2,contours=20,color=blue,thickness=2):display(fangxiangcang,

7、jifenquxian);的方向场及积分曲线的方向场及积分曲线14四川大学数学学院 徐小湛 May 2012方程的通解：方程的通解：wffc:=diff(y(x),x$2)-diff(y(x),x)-2*y(x)=0:tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-1.1),_C2=-1.1):plot(toplot,x=-1.1,y=-10.10,thickness=3,color=red);通解中的部分曲线通解中的部分曲线15四川大学数学学院 徐小湛 May 2012wffc:=diff(y(x),x$2)-diff(

8、y(x),x)-2*y(x)=0:tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-2.2),_C2=-2.2):plot(toplot,x=-1.1,y=-10.10,thickness=3,color=red);更多的曲线更多的曲线方程的通解：方程的通解：16四川大学数学学院 徐小湛 May 2012求特解：求特解：方程的通解：方程的通解：特解特解:17四川大学数学学院 徐小湛 May 2012通解中的部分曲线通解中的部分曲线和特解曲线和特解曲线wffc:=diff(y(x),x$2)+2*diff(y(x),x)+y(

9、x)=0:tongjie:=dsolve(wffc,y(x):tejie:=dsolve(wffc,y(0)=4,D(y)(0)=-2,y(x):tongjiequxian:=seq(seq(rhs(tongjie),_C1=3.5),_C2=0.4):p1:=plot(tongjiequxian,x=-2.2,y=0.8,thickness=1,color=red):p2:=plot(rhs(tejie),x=-2.2,y=0.8,thickness=5,color=blue):display(p1,p2,scaling=constrained);18四川大学数学学院 徐小湛 May 201

10、2wffc:=diff(y(x),x$2)+2*diff(y(x),x)+y(x)=0:tongjie:=dsolve(wffc,y(x):tejie:=dsolve(wffc,y(0)=4,D(y)(0)=-2,y(x):tongjiequxian:=seq(seq(rhs(tongjie),_C1=2.6),_C2=0.8):p1:=plot(tongjiequxian,x=-2.2,y=0.8,color=blue):p2:=plot(rhs(tejie),x=-2.2,y=0.8,thickness=5,color=red):display(p1,p2);19四川大学数学学院 徐小湛

11、May 2012方程的通解：方程的通解：20四川大学数学学院 徐小湛 May 2012通解中的部分曲线通解中的部分曲线wffc:=diff(y(x),x$2)-4*diff(y(x),x)+13*y(x)=0:tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-1.1),_C2=-1.1):plot(toplot,x=-1.1,y=-10.10,thickness=3,color=red);21四川大学数学学院 徐小湛 May 2012通解中的部分曲线通解中的部分曲线wffc:=diff(y(x),x$2)-4*diff(y(

12、x),x)+13*y(x)=0:tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-2.2),_C2=-2.2):plot(toplot,x=-2.1.5,y=-20.20,thickness=3,color=red);22四川大学数学学院 徐小湛 May 2012wffc:=diff(y(x),x$2)-4*diff(y(x),x)+13*y(x)=0:tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-3.3),_C2=-3.3):plot(toplot,

13、x=-2.1,y=-20.20,thickness=2,color=blue);23四川大学数学学院 徐小湛 May 2012方程：方程：方程的通解：方程的通解：24四川大学数学学院 徐小湛 May 2012wffc:=diff(y(x),x$2)+y(x)=2*x2-3:tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-3.3),_C2=-3.3):plot(toplot,x=-2.2,y=-12.6,thickness=2,color=blue);25四川大学数学学院 徐小湛 May 2012方程：方程：通解：通解：26

14、四川大学数学学院 徐小湛 May 2012wffc:=diff(y(x),x$2)-2*diff(y(x),x)-3*y(x)=exp(-x):tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-3.3),_C2=-3.3):plot(toplot,x=-2.2,y=-12.12,color=blue,thickness=2);27四川大学数学学院 徐小湛 May 2012方程：方程：通解：通解：28四川大学数学学院 徐小湛 May 2012wffc:=diff(y(x),x$2)-2*diff(y(x),x)+y(x)=(1+x)*exp(x):tongjie:=dsolve(wffc,y(x):toplot:=seq(seq(rhs(tongjie),_C1=-3.3),_C2=-3.3):plot(toplot,x=-4.3,y=-10.10,color=blue);29结束语结束语谢谢大家聆听！谢谢大家聆听！30