%《MATLAB及其在理工課程中的應(yīng)用指南》(十一五規(guī)劃版)第六章例6－2－1程序
% 真空中的拋射體運(yùn)動
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 2007年5月
%
clear;
y0 = 0;  x0 = 0;         % 初始位置
vMag = input('輸入初始速度 (m/s):(書上為50) '); 		% 輸入初始速度
vDir = input(' 輸入初速方向(度):(書上為40或50) ');
yf = input('輸入目標(biāo)高度（米）:(書上為8)  '); 			% 輸入目標(biāo)高度
vx0 = vMag*cos(vDir* (pi/180));    % 計算x,y方向的初始速度
vy0 = vMag*sin(vDir* (pi/180));    % 
wy = -9.81; wx = 0;           % 重力加速度 (m/s^2)
tf=roots([wy/2,vy0,y0-yf]); % 解方程wy*t^2/2+vy0*t+y0=yh，計算落點(diǎn)tf
tf=max(tf); 			% 去除落點(diǎn)時間tf中的庸解
t=[0:0.1:tf,tf];		% 設(shè)定時間數(shù)組，因tf不大可能被0.1整除，必須加一個tf點(diǎn)
y = y0 + vy0*t + wy*t.^2/2;		 % 計算軌跡
x = x0 + vx0*t + wx*t.^2/2;
xf = max(x),plot(x,y),grid 		 % 計算射程，畫出軌跡
set(gcf,'color','w')