实验九微分方程及应用.docVIP

实验九 微分方程及应用 【】【】MATLAB求常微分方程的解的方法。 【实验内容】 1．熟悉各种简单常微分方程及解法； 2．利用MATLAB求解常见常微分方程； 【】一、实验的基本理论与方法 二、实验使用的MATLAB函数 equ) ◆dsolve(equ, var) 上述命令调用格式中，equ为待求解的常微分方程，第一种调用格式视变量为自变量进行求解；第二种调用格式中var为指定变量，dsolve将以var为自变量进行常微分方程的求解。 dsolve(equ, condition1,condition2,,conditionm, var) 或 dsolve(equ, condition1, condition2,, conditionm, var)：求解有初始条件的常微分方程。 以上两种调用格式所得结果完全相同，其中：equ为常微分方程；condition1, condition2,,conditionm为初始条件；var为指定变量。 注：平时以形式出现的常微分方程，在Matlab中需重新改写。Matlab中用表示对变量求导数，表示对求一阶导数，表示对求阶导数，因此，这一常微分方程在Matlab中需描述为： 三、实验指导 。 可以执行以下命令：y=dsolve(x*Dy-y*log(y)=0,x) 则执行结果为： y = exp(exp(C1)*x) 其中exp(C1)表示所求出的解为通解。若一时粗心，把执行命令写为： y=dsolve(x*Dy-y*log(y)=0) 则所得到的执行结果为： y = exp(exp((t+C1)/x)) 显然结果相差甚远。 例2、求解微分方程。 可以执行以下命令：y=dsolve(Dy-2*y/(x+1)-(x+1)^(5/2)=0,x) 则执行结果为： y = 2/3*(x+1)^(3/2)*x^2+4/3*(x+1)^(3/2)*x+2/3*(x+1)^(3/2)+C1*x^2+2*C1*x+C1 例3、求解微分方程。 可以执行以下命令：y=dsolve(D3y=exp(2*x)-cos(x),x) 则执行结果为： y = 1/8*exp(2*x)+sin(x)+1/2*C1*x^2+C2*x+C3 例4、求解微分方程。 可以执行以下命令：y=dsolve(y*D2y-Dy^2=0,x) 则执行结果为： y = [ exp((x+C2)/C1)] 例5、求解微分方程。 可以执行以下命令：y=dsolve(D4y-2*D3y+5*D2y=0,x) 则执行结果为：y = C1+C2*x+C3*exp(x)*cos(2*x)+C4*exp(x)*sin(2*x) 例6、求解微分方程。 可以执行以下命令：y= dsolve(D2y-5*Dy+6*y=x*exp(2*x),x) 则执行结果为： y= -1/2*exp(2*x)*(x^2+2*x+2)+C1*exp(2*x)+C2*exp(3*x) 例7、求解微分方程。 可以执行以下命令：y=dsolve(D2y+y=cos(2*x),x) 则执行结果为： y = (1/2*sin(x)+1/6*sin(3*x))*sin(x)+(1/6*cos(3*x)-1/2*cos(x))*cos(x)+C1*sin(x)+C2*cos(x) 利用pretty命令将结果化为手写格式： pretty(y) 以上命令的执行结果为： (1/2 sin(x) + 1/6 sin(3 x)) sin(x) + (1/6 cos(3 x) - 1/2 cos(x)) cos(x) + C1 sin(x) + C2 cos(x) 例8、求解微分方程。 可以执行以下命令：y=dsolve(x*Dy+2*y=0,y(2)=1,x) 则执行结果为： y = 4/x^2 例9、求解常微分方程满足的特解。 可以执行以下命令：w=dsolve(D3w=-w,w(0)=1,Dw(0)=0,D2w(0)=0) 则执行结果为： w = 1/3*exp(-t)+2/3*exp(1/2*t)*cos(1/2*3^(1/2)*t) 若想使得到的结果为自变量的表达式，只需执行以下命令： w=dsolve(D3w=-w,w(0)=1,Dw(0)=0,D2w(0)=0,x) 以上命令的执行结果为： w = 1/3*exp(-x)+2/3*exp(1/2*x)*cos(1/2*3^(1/2)*x) 【实验练习】 1、求解下列微分方程。 （1）；（2）； （3）；（4）； （5）； （6）求解微分方程满足的特解。 （7），满足当时，的特解。 高等数学实验指导书 80 N阶导数表示方法。 方程 自变量 只有方程没有自变量