matlab报告用matlab研究抛体运动.docx

1、matlab报告用matlab研究抛体运动用matlab研究抛体运动2. 用matlab研究抛体运动2.1 引论MATLAB语言是一种集数值计算、符号运算、可视化建模、仿真和图形处理等多种功能的高级语言。使用MATLAB模拟物理现象为我们解决问题提供了一种新的方法，利用其方便的数值计算和作图功能，可以方便的模拟一些物理过程。对于处理非线性问题，既能进行数值求解，又能绘制有关曲线，方便实用，基于其功能强大，界面友善，语言自然，交互性强等优点，已成为教学和科研中最基础的软件之一，利用其解决复杂的数值计算问题，可以减少工作量，节约时间，图形绘制问题，真实直观，可以加深理解，提高工作效率将物体以一定的

2、初速度向空中抛出，仅在重力作用下物体所作的运动，它的初速度不为零，可分为平抛运动和斜抛运动。物理上提出的“抛体运动”是一种理想化的模型，即把物体看成质点，抛出后只考虑重力作用，忽略空气阻力。抛体运动加速度恒为重力加速度，相等的时间内速度变化量相等，并且速度变化的方向始终是竖直向下的。2.2抛体运动及应用2.2.1、 实验设计思路1、理论分析一般的处理方法是将其分解为水平方向和竖直方向，平抛运动水平方向是匀速直线运动，竖直方向是自由落体运动，斜抛运动水平方向是匀速直线运动，竖直方向是竖直上抛运动，在任意方向上分解有正交分解和非正交分解两种情加速度及位移等进行相应分析。无论怎样分解，都必须把运动的

3、独立性和独立作用原理结合进行系统分解，即将初速度、受力情、加速度及位移等进行相应分析。斜抛运动：水平方向速度 （1）竖直方向速度 （2）水平方向位移 （3）竖直方向位移 （4）平抛运动： 水平方向速度 （5） 竖直方向速度 (6)水平方向位移 (7)竖直方向位移 (8)合速度 （9） 合速度方向与水平夹角： (10)合位移 (11)位移方向与水平夹角： (12) 设某一抛射体的初速度为，抛射角为，将其运动在X,Y轴上进行正交分解，水平方向速度 (13)竖直方向 (14)质点的坐标是 (15) (16)从上两式消去，便得质点的轨迹运动方程 (17)抛射体能达到的最大高度为 (18)其到达最大高度

4、所需时间为 (19)空中飞行时间为 (20)抛射体的最大射程为 (21)它跟初速度和抛射角有关，在抛射角不变的情况下，射程与成正比，所以射程随初速度的增大而增大。在初速度不变的情况下，随着抛射角的增大，射程也增大，当度时，射程达到最大值，以后随着抛射角的增大，射程减小。 利用MATLAB的绘图功能，可以更直观的体现上述结论。（程序1）程序运行结果如图1所示。 图1 射程与抛射角、初速度的关系对于最大飞行路径所对应的抛射角问题（空气阻力忽略不计），X,Y坐标轴分别代表抛射体的射程与射高，在处，设在某一微小时段内抛射体的路径变量为，其对应的水平及竖直方向的变量为与，则 （22）设射程为R，则飞行路

5、径长度 （23） 根据前面的推论， (24)其中为抛射的初始速度，为抛射角，根据运动学原理，有 （25） （26）从（24）、(25）中消除，我们可得到该运动的抛物线方程： （27）从（24）中可知，为求解L,先得求出，因此在（4）式两边同时对求导，得： （28）将（27）代入式（24），等式两边同时积分，便得到了飞行路径长度与抛射角之间的关系： （29）根据式（28），为求得L的最大值，将（28）两边同时对求导 （30）令，可得到最大飞行路径所对应的抛射角的大小，但解此方程是比较困难的。为此，我们采用MATLAB的函数运算功能来解决这一问题。（程序2）程序如下，设其中的抛射初速度，。运行结果

6、如图2所示。图2 抛射角与飞行路径及其一阶导数曲线图图2给出了飞行路径随抛射角的变化曲线及飞行路径曲线的斜度，从图中可以得到，当（弧度）时，即度时，飞行路径最大，此时 （31）我们知道，在不考虑空气阻力的情况下，当抛射角度时，其射程最远，但此时其飞行路径并不是最远，而是当抛射角度时，其飞行路径最远，且其长度约为，实际上，由于空气阻力的存在，抛射体在空中是沿导弹曲线（弹头飞行时其重心所经过的路线）飞行的，它与抛物线不同，它的升弧与降弧不对称，在重力与空气阻力的共同影响下，弹道形成不均等的圆弧，升弧较长而直伸，降弧较短而弯曲.斜抛射出的炮弹的射程和射高都没有按抛体计算得到的值那么大，路线也不是理想

7、曲线。 物体在空气中受到的阻力，与物体运动速度大小有密切联系，速度越小，越接近理想情况，当物体速度低于200米每秒时，阻力与物体速度大小的平方成正比，速度介于400至600米每秒之间时，空气阻力与速度大小的三次方成正比，在速度很大的情况下，阻力与速度大小的高次方成正比。 将物体用一定的初速度沿水平方向抛出，不考虑空气的阻力，物体只在重力作用下所做的运动，叫做平抛运动。 竖直的重力与速度方向有夹角，做曲线运动；水平方向不受外力作用，是匀速运动，速度为Vo；竖直方向受重力作用，没有初速度，加速度为重力加速度g,是自由体运动。即做平抛运动的物体，在水平方向上由于不受力，将作匀速直线运动；在竖直方向上

8、的物体的初速度为0.且只受到重力作用，物体做自由落体运动，加速度为g。平抛运动的规律：（1）抛出t秒末的速度： 一抛出点为坐标原点，水平方向为x轴（正方向和初速度V0的方向相同），竖直方向为y轴，正方向向下，则： 水平分速度：Vx=Vo （32）竖直分速度：Vy=gt （33） 合速度：Vt= （34） tan= （35） (2)平抛运动的物体在任意时刻t的位置坐标：水平位移：x=Vot （36）竖直位移：y=g （37）合位移：s= （38） tan= （39）2、实验步骤(1)运用MATLAB编程得到平抛速度随时间的变化关系。(程序3)依据公示（32），（33），（34），(35) 图3

9、平抛运动速度随时间变化关系(2)运用MATLAB 编程到到平抛物体运动的曲线运用公式（32），（33），（34），(35),(37),(38),（39），我们可以求得物体在任意时刻的坐标并找到物体所在位置后，再用平滑曲线把这些点连起来，就得到平抛运动的轨迹。（程序4）运行结果如图4所示 图4 物体平抛轨迹曲线(3)利用matlab模拟物体斜抛运动通过该程序可以画出在任意位置以初始速度V和抛射角度抛出的轨迹。（程序5）按“run”运行时，弹出窗口点击图框中的“OK”，在“command window”中输出结果为：图5. 物体斜抛运动曲线(4)试计算抛射角为90度的特殊抛体运动任意时刻的位置和速

10、度一弹性小球，初始高度h=10m,向上初速度v0=15米每秒，与地面碰撞的速度衰减系数k=0.8，试计算任意时刻球的位置和速度。高度与时间的关系：， （40）速度与时间关系： （41）对等式两边积分，有， （42）， （43）由此可得数学方程：第一次落地前： （44） （45） （46）第二次落地前： （47） （48） (49) (50)第三次落地前： （51） （52） （53） （54）. . . . . .第n次落地前： （55） （56） （57) (58)如用手工进行计算，计算量极大，利用MATLAB编程（程序6）程序运行结果如图所示。 图6. 抛射角为90度的特殊抛体运动任意时刻

11、的位置和速度 (5)用matlab研究定点投篮命中率问题下图为过罚球点P和篮筐中心Q且垂直于地面的平面示意图。P,Q之间的距离L=4.60cm,Q点的高度H=3.05 cm，篮球直径d=24.6cm,篮筐直径D=45.0 cm.某人投篮的出手高度h位1.8-2.1m,出手速度v=8.09.0m/s。问题一：考虑球心对框心的点对点的投篮，求出手速度和出手方向的范围 问题二：若考虑球的大小和框的大小进行投篮，球入筐时可以偏离框心，求出手速度、角度及其最大偏值示意图1问题一不考虑篮球和篮筐的大小，不考虑空气阻力大小的影响，从未出手时的球心P为坐标原点，x轴为水平方向，y轴为竖直方向，篮球在t=0时以

12、出手速度v和出手角度投出，可视为质点的斜抛运动，其运动方程为： (59)其中g是重力加速度，由此可得球心的运动轨迹如下抛物线 (60)以x=L，y=H-h代入(60)式，就得到了球心命中框心的条件 （61）可以看出，给定出手速度v和出手高度h，就有两个满足条件，而（61）式有解的前提为： （62）可解得： （63）于是对于一定的出手高度h，使（63）式等号成立的v为最小的出手速度示意图2球入篮筐处的入射角度为，可从下式得到： （64）这里的导数由(60)式计算代入后得对应，有，设（65）求解（程序7）表1 对于不同的出手速度和出手高度的出手角度和入射角度速度 高度 最大出手角1 最小出手角 对

13、应入框角1 对应入框角28.0000 1.8000 62.4099 53.8763 42.7925 20.92138.0000 1.8500 62.7779 54.8876 41.8430 20.49158.0000 1.9000 63.1174 55.8206 40.9188 20.14318.0000 1.9500 63.4330 56.6865 40.0156 19.86498.0000 2.0000 63.7281 57.4941 39.1300 19.64788.0000 2.0500 64.0054 58.2505 38.2594 19.48478.0000 2.1000 64.2

14、670 58.9615 37.4017 19.36988.5000 1.8000 67.6975 62.1726 37.5049 12.62508.5000 1.8500 67.8663 62.6922 36.7546 12.68708.5000 1.9000 68.0288 63.1884 36.0075 12.77538.5000 1.9500 68.1854 63.6632 35.2632 12.88828.5000 2.0000 68.3367 64.1179 34.5214 13.02408.5000 2.0500 68.4829 64.5541 33.7819 13.18118.5000 2.1000 68.6244 64.9729 33.0444 13.35839.0000 1.8000 71.0697 67.1426 34.1327 7.65509.0000 1.8500 71.1736 67.4