《matlab报告——用matlab研究抛体运动.doc》由会员分享，可在线阅读，更多相关《matlab报告——用matlab研究抛体运动.doc（21页珍藏版）》请在金锄头文库上搜索。
1、用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)其到达最大高度所需时间为 (19)空中飞行时间为 (20
4、)抛射体的最大射程为 (21)它跟初速度和抛射角有关，在抛射角不变的情况下，射程与成正比，所以射程随初速度的增大而增大。在初速度不变的情况下，随着抛射角的增大，射程也增大，当度时，射程达到最大值，以后随着抛射角的增大，射程减小。 利用MATLAB的绘图功能，可以更直观的体现上述结论。（程序1）程序运行结果如图1所示。 图1 射程与抛射角、初速度的关系对于最大飞行路径所对应的抛射角问题（空气阻力忽略不计），X,Y坐标轴分别代表抛射体的射程与射高，在处，设在某一微小时段内抛射体的路径变量为，其对应的水平及竖直方向的变量为与，则 （22）设射程为R，则飞行路径长度 （23） 根据前面的推论， (24
5、)其中为抛射的初始速度，为抛射角，根据运动学原理，有 （25） （26）从（24）、(25）中消除，我们可得到该运动的抛物线方程： （27）从（24）中可知，为求解L,先得求出，因此在（4）式两边同时对求导，得： （28）将（27）代入式（24），等式两边同时积分，便得到了飞行路径长度与抛射角之间的关系： （29）根据式（28），为求得L的最大值，将（28）两边同时对求导 （30）令，可得到最大飞行路径所对应的抛射角的大小，但解此方程是比较困难的。为此，我们采用MATLAB的函数运算功能来解决这一问题。（程序2）程序如下，设其中的抛射初速度，。运行结果如图2所示。图2 抛射角与飞行路径及其一阶
6、导数曲线图图2给出了飞行路径随抛射角的变化曲线及飞行路径曲线的斜度，从图中可以得到，当（弧度）时，即度时，飞行路径最大，此时 （31）我们知道，在不考虑空气阻力的情况下，当抛射角度时，其射程最远，但此时其飞行路径并不是最远，而是当抛射角度时，其飞行路径最远，且其长度约为，实际上，由于空气阻力的存在，抛射体在空中是沿导弹曲线（弹头飞行时其重心所经过的路线）飞行的，它与抛物线不同，它的升弧与降弧不对称，在重力与空气阻力的共同影响下，弹道形成不均等的圆弧，升弧较长而直伸，降弧较短而弯曲.斜抛射出的炮弹的射程和射高都没有按抛体计算得到的值那么大，路线也不是理想曲线。 物体在空气中受到的阻力，与物体运动
7、速度大小有密切联系，速度越小，越接近理想情况，当物体速度低于200米每秒时，阻力与物体速度大小的平方成正比，速度介于400至600米每秒之间时，空气阻力与速度大小的三次方成正比，在速度很大的情况下，阻力与速度大小的高次方成正比。 将物体用一定的初速度沿水平方向抛出，不考虑空气的阻力，物体只在重力作用下所做的运动，叫做平抛运动。 竖直的重力与速度方向有夹角，做曲线运动；水平方向不受外力作用，是匀速运动，速度为Vo；竖直方向受重力作用，没有初速度，加速度为重力加速度g,是自由体运动。即做平抛运动的物体，在水平方向上由于不受力，将作匀速直线运动；在竖直方向上的物体的初速度为0.且只受到重力作用，物体
8、做自由落体运动，加速度为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 平抛运动速度随时间变化关系(2)运用MAT
9、LAB 编程到到平抛物体运动的曲线运用公式（32），（33），（34），(35),(37),(38),（39），我们可以求得物体在任意时刻的坐标并找到物体所在位置后，再用平滑曲线把这些点连起来，就得到平抛运动的轨迹。（程序4）运行结果如图4所示 图4 物体平抛轨迹曲线(3)利用matlab模拟物体斜抛运动通过该程序可以画出在任意位置以初始速度V和抛射角度抛出的轨迹。（程序5）按“run”运行时，弹出窗口将图框中的相关数据更改为：点击图框中的“OK”，在“command window”中输出结果为：图5. 物体斜抛运动曲线(4)试计算抛射角为90度的特殊抛体运动任意时刻的位置和速度一弹性小球，初始高度h=10m,向上初速度v0=15米每秒，与地面碰撞的速度衰减系数k=0.8，试计算任意时刻球的位置和速度。高度与时间的关系：， （40）速度与时间关系：
