实验六常微分方程组求解.docVIP

常微分方程（组）求解 【实验目的】 1．求微分方程的符号解。 2．求微分方程的数值解。 【实验相关知识】 (一)微分方程的符号解. 微分方程的符号解也叫做解析解. 求微分方程（组）的符号解用命令dsolve. 命令格式如下： s=dsolve(方程1,方程2,…,初始条件1,初始条件2,…,自变量) 说明：用字符串表示方程，自变量缺省则默认为t. 导数用D表示, 2阶导数用D2表示,以此类推。返回值s是符号解. 例1. 求初值问题的解。 dsolve 命令格式是： dsolve(Dy = y*t/(t-5),y(0) = 2) ans = -2/3125*exp(t)*(t-5)^5 或者 s=dsolve(Dy = y*t/(t-5),y(0) = 2) s = -2/3125*exp(t)*(t-5)^5 二阶或更高阶方程的处理情况相类似。例如： , y(0) = -1, y(0) = 2 写成如下的形式输入MATLAB： dsolve(D2y - y = 0,y(0) = -1,Dy(0) = 2) ans = -3/2*exp(-t)+1/2*exp(t) 例2: 求满足初始条件的符号解. 写成如下的形式输入MATLAB： s=dsolve(D2y=sin(2*x)-y,y(0)=0,Dy(0)=1,x) simplify(s) %如果得到符号解比较复杂,可以试