matlab微分方程组解法.pptVIP

数学实验;实验目的;求微分方程的数值解;1、目标跟踪问题一：导弹追踪问题;微分方程的解析解; 解 输入命令: y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,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=simple(x) % 将x化简 y=simple(y) z=simple(z);微分方程的数值解;（二）建立数值解法的一些途径;2、使用数值积分;3、使用泰勒公式;; 1、在解n个未知函数的方程组时，x0和x均为n维向量，m-文件中的待解方程组应以x的分量形式写成.;解: 令 y1=x，y2=y1’;解 1、建立m-文件rigid.m如下： function dy=rigid(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-0.51*y(1)*y(2);;;由(1),(2)消去t整理得模型:;解法二(数值解);解法三(建立参数方程求数值解);4. 解导弹运动轨迹的参数方程;5. 结??见图1;慢跑者与狗;2