沈阳理工大学装备工程学院综合课程设计说明书
q0=70*pi/180; tao0=45*pi/180;
faiv=k*q0-tao0-(k-1)*x(2); f(1)=vt*cos(x(2))-vt*x(3)*cos(faiv); f(2)=(vt*x(3)*sin(faiv)-vt*sin(x(2)))/x(1); f(3)=0; f=f(:);
3 分析不同弹道参数的程序
clear; clc; k=4; tf=60;
x0=[4700,1.202,k]; [t,x]=ode45(@kfun,[0 tf],x0); zt=4700*sin(70*pi/180);
xt=4700*cos(70*pi/180)+12.*t(1:60); xm=xt-x(1:60,1).*cos(x(1:60,2)); zm=zt-x(1:60,1).*sin(x(1:60,2)); r=x(:,1); q=x(:,2)*57.3; figure (1) plot(xt,zt,xm,zm,'--'); xlabel('X/m'); ylabel('Z/m');
title('比例k=4对应的弹道曲线');
figure (2) plot(t,r); xlabel('t/s'); ylabel('r/m');
31
title('');
figure (3) plot(t,q); xlabel('t/s'); ylabel('q/m'); title('');
沈阳理工大学装备工程学院综合课程设计说明书
32