用Matlab件求常微分方程解(或通解).docx

用Matlab 件求常微分方程解(或通解) -＿... - ＿....-＿ ＿~－“．一· - ---- ———————————————————————————————— 作者： ———————————————————————————————— 日期： 2 吕梁学院《高等数学》实验报告 《高等数学》实验报告 实验人员:系(班): 学号: 姓名: 实验地点:电教楼五号机房 实验名称:Matlab 高等数学实验 实验时间:2014-6-3 16:30--18:30 实验名称:用Matlab 软件求常微分方程的解(或通解) 实验目的:熟练掌握 Matlab 软件求常微分方程的解(或通解) 实验内容:(给出实验程序与运行结果) 一、求微分方程的特解. ? d 2 y 1、? 1、 ? dx2 4 dy dx ? 3y ? 0 ?? y(0) ? 6, y (0) ? 10 程序: dsolve(D2y-4*Dy+3*y,y(0)=6,Dy(0)=10,x) ans = 4*exp(x)+2*exp(3*x) 3 ? d 2 y dy 2、?4 ? 4 ? y ? 0 ? dx2 dx ?? y(0) ? 2, y (0) ? 0 程序:dsolve(4*D2y+4*Dy+y,y(0)=2,Dy(0)=0,x) ans = ?3、?d 2 y ? 4 dy ? 29 y ? 0 ? 2*exp(-1/2*x)+exp(-1/2*x)*x ? dx2 dx ?? y(0) ? 0, y(0) ? 15 程序:dsolve(D2y+4*Dy+29*y=0,y(0)=9,Dy(0)=15,x) ans = 4、? d 2 y 4、 ? ? dx2  ? 4 dy ?13y ? 0 dx 33/5*exp(-2*x)*sin(5*x)+9*exp(-2*x)*cos(5*x) ?? y(0) ? 0, y(0) ? 3 程序:dsolve(D2y-4*dy+13*y=0,y(0)=0,Dy(0)=3,x) ans = 3/13*sin(13^(1/2)*x)*13^(1/2)-4/13*cos(13^(1/2)*x)*dy+4/13*dy 5、? d 2 y 5、 ? ? dx2 3 dy dx ? 4 y ? 0 ?? y(0) ? 0, y(0) ? ?5 程序:dsolve(D2y-3*Dy-4*y,y(0)=0,Dy(0)=-5,x) ans = exp(-x)-exp(4*x) 4 二、求齐次非线性微分方程的通解 1、 d 2 y dx2 2 dy dx ? 3y ? 3x ?1 程序:dsolve(D2y-2*Dy-3*y=3*x+1,x) ans = exp(-x)*C2+exp(3*x)*C1+1/3-x 2、 d 2 y dy? 5 ? 6 y ? 2、 d 2 y dy dx2 dx 程序:dsolve(D2y-5*Dy+6*y=x*exp(2*x),x) ans = exp(3*x)*C2+exp(2*x)*C1-1/2*x*exp(2*x)*(2+x) 3、 d 2 y dx2 4 y ? x cos x 程序:dsolve(D2y+4*y=x*cos(x),x) ans = sin(2*x)*C2+cos(2*x)*C1+2/9*sin(x)+1/3*x*cos(x) 4、 d 2 y dx2 y ? ex ? cos x 程序:dsolve(D2y+y=exp(x),x) ans = sin(x)*C2+cos(x)*C1+1/2*exp(x) 5 dsolve(D2y+y=cos(x),x) ans = sin(x)*C2+cos(x)*C1+1/2*cos(x)+1/2*sin(x)*x 则 原 式 = sin(x)*C2+cos(x)*C1+1/2*exp(x)+sin(x)*C2+cos(x)*C1+1/2*cos(x)+ 1/2*sin(x)*x 5、 d 2 y dx2 2 dy dx 5 y ? sin 2x 程序:dsolve(D2y-2*Dy+5*y=sin(2*x),x) ans = exp(x)*sin(2*x)*C2+exp(x)*cos(2*x)*C1+1/17*sin(2*x)+4/17*cos(2 *x) 三、微分方程实例 1、试求y ? x 的经过点 M（0,1）且在此点与直线y ? ? y ? x ?1由题意得? ? 1 ? y(0) ? 1, y(0) ? x ?1相切的积分曲线。 2 ? 2 程序:dsolve(D2y=x,y(0)=1,Dy(0)=1/2,x) ans = 1/6*x^3+1/2*x+1 实验心得: Matlab 是一个画图和