例6-7(2).m

这是一个MATLAB7.0基础与提高例题的所有源码

clear;
y0=0;
x0=0;                 						%初始位置
vMag=input('输入初始速度(m/s)：');        		%输入初始速度
vDir=input('输入初始方向(度)');
tf=input('输入飞行时间(秒)');            			%输入飞行时间
vx0=vMag*cos(vDir*(pi/180));           			%计算x，y方向的初始速度
vy0=vMag*sin(vDir*(pi/180));
r0=[0;0;vx0;vy0];
[t,r]=ode45('e31a',[0,tf],r0),           			%数值积分（调用函数程序e31a.m）
plot(r(:,1),r(:,2)),hold on                   		%计算轨迹
%ode45规定返回的结果中：t是列向量，各时刻的r为4列向量
%注意下一个语句的意义：找y<0的下标所对应的x的最小值，以粗略计算射程
xmax=min(r(find(r(:,2)<0),1))
plot([0,150],[0,0])                    			%画出x坐标线