资源描述:
《考虑空气阻力的抛射体运动mtlab仿真》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、专业资料考虑空气阻力的抛射体的运动****(天水师范学院物理与信息科学学院)摘要:在力学中涉及许多复杂的计算问题,例如非线性问题,对其求解有时是很困难的MATLAB是一种科学计算软件,主要适用于矩阵运算及控制和信息处理领域的分析设计,它使用方便,输入简捷,运算高效,内容丰富,将MATLAB引入物理教学,利用可视化的功能,可对物理学中抽象、不易理解的部分进行计算机模拟,从而使教学生动具体。以提高学生的学习兴趣和积极性,得到明显的教学效果关键词:MATLAB软件;力学;空气阻力;质点轨迹;抛射体Fromtheairresistanceoftheproject
2、iles****(CollegeofPhysicsandInformationScience,TianshuiNormalUniversity)Summary:inmechanicsofmanycomplexcalculations,forexample,thesolutionofthenonlinearissometimesverydifficultMATLABisascientificterms,mainlyappliedtothematrixmultiplicationandcontrolthefieldofanalysisandinformati
3、onprocessing,itiseasytouseandoperationofinput,efficient,informative,willbeintroducedinphysicsteaching,usingMATLABvisualfeatures,butforphysicsintheabstract,isnoteasytounderstandapartofthecomputermodels,soteachingisaspecific.Keywords:MATLABsoftware;mechanicstheparticle;airresistanc
4、e;track;projectiles引言:word完美格式专业资料在力学中涉及许多复杂的计算问题,例如非线性问题,对其求解有时是很困难的,MATLAB正是处理非线性问题的很好工具,既能进行数值计算,又能绘制有关曲线,非常方便。经典力学原理基本上沿着两条路线进行。一条是基于牛顿运动定律,在静力分析中,主要遵循力学平衡原理,加上组成结构材料的本构关系和应变,位移的几何协调关系可以导出微分方程。另一条是基于功、能原理,它以能量原理(如最小势能原理,虚位移原理等)为基础,可以导出需要求解的积分方程。不管是解微分方程还是积分方程,均需求出函数y=f(x),使之满
5、足方程并在边界上满足边界条件。对于简单问题可以求得其解析解,但物理学中的实际问题是很复杂的,往往很难求得其实用的解析解,因此,应用计算机得到其数值解成了可行的解决问题的途径。常用的数值方法有差分法、有限元法、加权残值法、边界元法等,这些解法通常都有大量的矩阵用算以及其他数值用算。MATLAB具有强大的科学计算功能,这使得人们可以用它来代替FORTRAN等传统的编程语言。在计算要求相同的情况下,使用MATLAB编程,工作量会大大减少。例如在物理学中应用FORTRAN求解振子的自由振动,需要调用jacobi子程序求解矩阵的特征值及对应的特征向量,而且要求该矩
6、阵必须为实对称矩阵,程序繁杂,且对一般用户来说,要看懂程序算法实属不易,而采用MATLAB编制该自由振动的子程序时只需要调用两个函数:求逆矩阵的inv函数及求特征值和特征向量的eig函数。下面简单地介绍MATLAB在物理力学中的应用。落体偏东和导弹右偏都是地球自转的非惯性效应.关于落体偏东这个问题的讨论一般都是在不计空气阻力影响的条件下进行,而且已有许多种解释方法;以及落体偏东这个问题在考虑空气阻力影响下的讨论也已有多种解释.其实,空气阻力不仅对落体偏东有相当的影响,同样对抛射体运动有不可忽视的影响.本文介绍抛射体在空气阻力两种影响下,并且空气阻力与抛体
7、速度平方成正比时,从牛顿动力学方程出发导出北半球抛射体的运动方程,word完美格式专业资料最后从抛射体的运动方程着手,讨论了一些特例,结果所得,一般文献将均被本文包含.2、质点飞行的轨迹和距离考虑空气阻力是抛射体质心的飞行轨迹问题[1]。假设空气阻力的方向与速度向量相反,大小与速度的平方成正比。计算质点飞行的轨迹和距离。图1抛射体的受力图[2]0xmgVword完美格式专业资料其中c为空气阻力系数。用MATLAB来进行数值求解。设r=,则原方程等价于下面的以及微分方程组;3、MATLAB仿真与分析程序如下:(1)函数文件zf.mfunctionrdot=
8、zf(t,r)c=0.02;g=9.8;m=1;word完美格式专业资料vm=s