常微分方程的数值解与解析解.docVIP

常微分方程的解析解 常微分方程的解析解也就是常微分方程的精确解，也称为常微分方程的符号解；一般可理解为求微分方程的通解或者特解的解析式或表达式；但只有少数的微分方程存在解析解。 在MATLAB中，由函数dsolve()求解常微分方程（组）的解析解，其具体格式如下： X=dsolve(‘方程1’，‘方程2’，…‘方程n’，‘初始条件’，‘自变量’) 函数dsolve用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解。 例1：求解常微分方程的MATLAB程序为：dsolve(Dy=1/(x+y),x)，注意，系统缺省的自变量为t，因此这里要把自变量写明。 结果为：-lambertw(-C1*exp(-x-1))-x-1 其中：Y=lambertw(X)表示函数关系Y*exp(Y)=X。 例2：求解常微分方程的MATLAB程序为：Y2=dsolve(y*D2y-Dy^2=0’,’x’) 结果为： Y2 =[ exp((x+C2)/C1)] [ C2] 我们看到有两个解，其中一个是常数。 例3：求常微分方程组通解的MATLAB程序为： [X,Y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=exp(2*t),t) 例4：求常微分方程组通解的MATLAB程序为： [X,Y]=dsolve(Dx+2*x-Dy=10*c