考虑空气阻力的抛射体运动mtlab仿真.doc

考虑空气阻力的抛射体的运动(天水师范学院物理与信息科学学院)摘要在力学中涉及许多复杂的计算问题，例如非线性问题，对其求解有时是很困难的MATLAB是一种科学计算软件，主要适用于矩阵运算及控制和信息处理领域的分析设计，它使用方便，输入简捷，运算高效，内容丰富，将MATLAB引入物理教学，利用可视化的功能，可对物理学中抽象、不易理解的部分进行计算机模拟，从而使教学生动具体。以提高学生的学习兴趣和积极性，得到明显的教学效果关键词MATLAB软件；力学；空气阻力；质点轨迹；抛射体FROMTHEAIRRESISTANCEOFTHEPROJECTILES(COLLEGEOFPHYSICSANDINATIONSCIENCE,TIANSHUINORMALUNIVERSITY)SUMMARYINMECHANICSOFMANYCOMPLEXCALCULATIONS,FOREXAMPLE,THESOLUTIONOFTHENONLINEARISSOMETIMESVERYDIFFICULTMATLABISASCIENTIFICTERMS,MAINLYAPPLIEDTOTHEMATRIXMULTIPLICATIONANDCONTROLTHEFIELDOFANALYSISANDINATIONPROCESSING,ITISEASYTOUSEANDOPERATIONOF,EFFICIENT,INATIVE,WILLBEINTRODUCEDINPHYSICSTEACHING,USINGMATLABVISUALFEATURES,BUTFORPHYSICSINTHEABSTRACT,ISNOTEASYTOUNDERSTANDAPARTOFTHECOMPUTERMODELS,SOTEACHINGISASPECIFICKEYWORDSMATLABSOFTWARE；MECHANICSTHEPARTICLE；AIRRESISTANCE；TRACK；PROJECTILES引言在力学中涉及许多复杂的计算问题，例如非线性问题，对其求解有时是很困难的，MATLAB正是处理非线性问题的很好工具，既能进行数值计算，又能绘制有关曲线，非常方便。经典力学原理基本上沿着两条路线进行。一条是基于牛顿运动定律，在静力分析中，主要遵循力学平衡原理，加上组成结构材料的本构关系和应变，位移的几何协调关系可以导出微分方程。另一条是基于功、能原理，它以能量原理(如最小势能原理，虚位移原理等)为基础，可以导出需要求解的积分方程。不管是解微分方程还是积分方程，均需求出函数YFX,使之满足方程并在边界上满足边界条件。对于简单问题可以求得其解析解，但物理学中的实际问题是很复杂的，往往很难求得其实用的解析解，因此，应用计算机得到其数值解成了可行的解决问题的途径。常用的数值方法有差分法、有限元法、加权残值法、边界元法等，这些解法通常都有大量的矩阵用算以及其他数值用算。MATLAB具有强大的科学计算功能，这使得人们可以用它来代替FORTRAN等传统的编程语言。在计算要求相同的情况下，使用MATLAB编程，工作量会大大减少。例如在物理学中应用FORTRAN求解振子的自由振动，需要调用JACOBI子程序求解矩阵的特征值及对应的特征向量，而且要求该矩阵必须为实对称矩阵，程序繁杂，且对一般用户来说，要看懂程序算法实属不易，而采用MATLAB编制该自由振动的子程序时只需要调用两个函数求逆矩阵的INV函数及求特征值和特征向量的EIG函数。下面简单地介绍MATLAB在物理力学中的应用。落体偏东和导弹右偏都是地球自转的非惯性效应关于落体偏东这个问题的讨论一般都是在不计空气阻力影响的条件下进行,而且已有许多种解释方法以及落体偏东这个问题在考虑空气阻力影响下的讨论也已有多种解释其实,空气阻力不仅对落体偏东有相当的影响,同样对抛射体运动有不可忽视的影响本文介绍抛射体在空气阻力两种影响下,并且空气阻力与抛体速度平方成正比时,从牛顿动力学方程出发导出北半球抛射体的运动方程,最后从抛射体的运动方程着手,讨论了一些特例,结果所得,一般文献将均被本文包含2、质点飞行的轨迹和距离考虑空气阻力是抛射体质心的飞行轨迹问题1。假设空气阻力的方向与速度向量相反，大小与速度的平方成正比。计算质点飞行的轨迹和距离。0XXVMGYVV2CV图1抛射体的受力图2MGCVMGCVDTVMCVCVTDVTYVDTXYYXXYXSINOS22其中C为空气阻力系数。用MATLAB来进行数值求解。设R，则原方程等价于下面的以及微分方程组；YXV,GRMCRDT4233433、MATLAB仿真与分析程序如下1函数文件ZFMFUNCTIONRDOTZFT,RC002G98M1VMSQRTR32R42RDOTR3R4CVMR3/MCVMR4/MG2主程序EX12＿6MCLEARY00X00初始位置V0 请输入初始输入速度M/S RHO 请输入初始方向度 TF 请输入飞行时间S VX0VOCOSRHOPI/180计算X、Y方向的初始速度VY0V0SINGRHOPI/180T,RODE45 ZF ,0,TF,00VX0VY0解微分方程HMAXR,2求轨道的最高点TTFINDR,2H到最高点的所需时间LMINRFINDR,20,1计算射程PLOT0,100,0,0HOLDON绘制X坐标线XLABEL X YLABELYPLOTR,1,R,2绘制运功轨迹请输入初始输入速度M/S57请输入初始方向度36请输入飞行时间S65H218766T17619L739735其运动轨迹如图010203040506070809010050403020100102030体1体体体体体体体体体体体体体体体Y4、结论当抛射体运动考虑空气阻力的影响时，它的落地点将会发生偏移，与只考虑重力时的将会有所不同，所以在导弹的发射中为了进一步提高准确度都会考虑空气阻力对它的影响3，当然在我们的一般的研究中不会考虑空气阻力的影响参考文献1楼荣训,陈静源关于空气阻力对落体偏东影响的研究J浙江师大学报,1994,1741～452周衍柏理论力学教程M北京人民教育出版社,1986256～2573作者朱戎文章题名［空气阻力和地球自转对抛射体影响的系统研究］期刊名陕西师范大学学报(自然科学版)，2003年，03期课程论文题目考虑空气阻力的抛射体的运动学院物理与信息科学学院姓名张新专业物理学学号271040148研究类型研究综述指导教师提交日期