例6-4.m

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

clear;
%初始位置
x0=0;
y0=0;  
%输入初始速度的大小                                             
vMag=input('初始速度(m/s)：');	
%输入初始速度的方向					
vDir=input('初速方向(度):');	
%输入目标高度yh						
yh=input('目标高度(m):');
%计算x,y方向的初始速度							
vx0=vMag*cos(vDir*(pi/180));                           
vy0=vMag*sin(vDir*(pi/180));
%重力加速度
wy=-9.81;                                            
wx=0;
%解二次线性方程，计算落点时间th
th=roots([wy/2,vy0,y0-yh]);  
%去除th两个解中的庸解                            
th=max(th);                                           
t=0:0.1:th;
%计算轨迹
y=y0+vy0*t+wy*t.^2/2;                                
x=x0+vx0*t+wx*t.^2/2;
%计算射程，画出轨迹
xh=max(x),plot(x,y);