工程应用数学11-4.ppt

§11.4 MATLAB与微分方程求解 微分方程可以通过函数dsolve求解， 该函数的调用格式为： r = dsolve(eq1,eq2,…,cond1,cond2,…,v) r = dsolve(eq1,eq2,…,cond1,cond2,…,v) 输入参数： eq1,eq2,…表示微分方程， v为独立变量， cond1,cond2,…表示边界条件和/或初始条件 * §11.4 MATLAB与微分方程求解 默认的独立变量是t，用户也可以使用别 的变量来代替t，只要把别的变量放在输入变 量的最后即可． 字母D代表微分算子，即d/dt，字母D后 面所跟的数字代表几阶微分，如D=d/dx， D2=d2/dx2，…．跟在微分算子后面的字母 是被微分的变量，如D3y代表对y(t)的三阶微 分． * §11.4 MATLAB与微分方程求解 初始和边界条件由字符串表示，如 y(a)=b，Dy(c)=d，D2y(e)=f分别表示 例如 f1 = dsolve(D2y - Dy =exp(x)) f2 = dsolve((Dy)^2 + y^2 = 1,s) * §11.4 MATLAB与微分方程求解 % 带一个定解条件 D3 = dsolve(Dy = a*y, y(0) = b) % 带两个定解条件 D4 = dsolve(D2y = -a^2*y, y(0) = 1, Dy(pi/a) = 0) % 求解线性微分方程组 [x,y] = dsolve(Dx = y, Dy = -x) [u,v] = dsolve(Du=u+v,Dv=u-v) * §11.4 MATLAB与微分方程求解 例11-22 求解微分方程： 的通解． 解：在MATLAB中输入以下命令： dsolve(Dy=1+y^2,x) 例11-23 求通解 解：在MATLAB中输入以下命令： y=dsolve(Dy=(sin(x)-cos(x))*sqrt(1-y^2),x) * §11.4 MATLAB与微分方程求解 例11-24 求方程 满足条件 的特解. 解：在MATLAB中输入以下命令： y = dsolve(Dy=2*x*y/(x^2+1),y(0)=1,x) * §11.4 MATLAB与微分方程求解 例11-25 求微分方程 的通解. 解：在MATLAB中输入以下命令： y = dsolve(Dy - 2*y = exp(-x),x) * §11.4 MATLAB与微分方程求解 例11-26 求微分方程 满足初始条件 的特解 解：在MATLAB中输入以下命令： y=dsolve(x*Dy+y=cos(x),y(pi)=1,x) * §11.4 MATLAB与微分方程求解 例11-27 求方程 满足初始条件 的特解 解：在MATLAB中输入以下命令： y=dsolve(D2y+4*y=0,y(0)=1,Dy(0)=0,x) *