用Matlab解微分方.pptVIP

* 用 Matlab 解微分方程 一、微分方程的解析解 求微分方程（组）的解析解用函数dsolve。 dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’) 结果：u =tan(t+C1) 输入命令: y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x) 结 果 : y =3*exp(-2*x)*sin(5*x) 解 解 输入命令 ： [x,y,z]=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z, t)； 结果为： x =exp(2*t)*C1+C2*exp(-t)-C2*exp(2*t)+exp(2*t)*C3-3*exp(-t) y =-C1*exp(-2*t)+exp(2*t)*C1+C2*exp(-2*t)+C2*exp(-t)-C2*exp(2*t)+exp(2*t)*C3-C3*exp(-t) z =-C1*exp(-2*t)+exp(2*t)*C1-C2*exp(2*t)+C2*exp(-2*t)+exp(2*t)*C3 二、微分方程的数值解 （一）常微分方程数值解的定义 在生产和科研中所处理的微分方程往往很复杂且大多得不出一般解。而在实际上对初值问题，一般是要求得到解在若干个点上满足规定精确度的近似值，或者得到一