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

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

考虑空气阻力的抛射体的运动 ****(天水师范学院 物理与信息科学学院) 摘要：在力学中涉及许多复杂的计算问题，例如非线性问题，对其求解有时是很困难的 MATLAB是一种科学计算软件，主要适用于矩阵运算及控制和信息处理领域的分析设计， 它使用方便，输入简捷，运算高效，内容丰富，将MATLAB引入物理教学，利用可视化的 功能，可对物理学中抽象、不易理解的部分进行计算机模拟，从而使教学生动具体。以 提高学生的学习兴趣和积极性，得到明显的教学效果 关键词：MATLAB软件；力学；空气阻力；质点轨迹；抛射体 From the air resistance of the projectiles **** (College of Physics and Ination Science,Tianshui Normal University) Summary ： in mechanics of many complex calculations, for example, the solution of the nonlinear is sometimes very difficult MATLAB is a scientific terms, mainly applied to the matrix multiplication and control the field of analysis and ination processing, it is easy to use and operation of , efficient, inative, will be introduced in physics teaching, using MATLAB visual features, but for physics in the abstract, is not easy to understand a part of the computer models, so teaching is a specific. Keywords ： MATLAB software ； mechanics the particle ； air resistance ； track ； projectiles 引言： 在力学中涉及许多复杂的计算问题，例如非线性问题，对其求解有时是很困难的，MATLAB正是处理非线性问题的很好工具，既能进行数值计算，又能绘制有关曲线，非常 方便。经典力学原理基本上沿着两条路线进行。一条是基于牛顿运动定律，在静力分析中， 主要遵循力学平衡原理，加上组成结构材料的本构关系和应变，位移的几何协调关系可 以导出微分方程。另一条是基于功、能原理，它以能量原理(如最小势能原理，虚位移 原理等)为基础，可以导出需要求解的积分方程。 不管是解微分方程还是积分方程，均需求出函数y=f(x),使之满足方程并在边界上满足边 界条件。对于简单问题可以求得其解析解，但物理学中的实际问题是很复杂的，往往很 难求得其实用的解析解，因此，应用计算机得到其数值解成了可行的解决问题的途径。 常用的数值方法有差分法、有限元法、加权残值法、边界元法等，这些解法通常都有大 量的矩阵用算以及其他数值用算。MATLAB具有强大的科学计算功能，这使得人们可以用 它来代替FORTRAN等传统的编程语言。在计算要求相同的情况下，使用MATLAB编程，工 作量会大大减少。例如在物理学中应用FORTRAN求解振子的自由振动，需要调用jacobi 子程序求解矩阵的特征值及对应的特征向量，而且要求该矩阵必须为实对称矩阵，程序 繁杂，且对一般用户来说，要看懂程序算法实属不易，而采用MATLAB编制该自由振动的 子程序时只需要调用两个函数：求逆矩阵的inv函数及求特征值和特征向量的eig函数。 下面简单地介绍MATLAB在物理力学中的应用。 落体偏东和导弹右偏都是地球自转的非惯性效应. 关于落体偏东这个问 题的讨论一般都是在不计空气阻力影响的条件下进行, 而且已有许多种解释 方法; 以及落体偏东这个问题在考虑空气阻力影响下的讨论也已有多种解释. 其实, 空气阻力不仅对落体偏东有相当的影响, 同样对抛射体运动有不可忽 视的影响. 本文介绍抛射体在空气阻力两种影响下, 并且空气阻力与抛体速 度平方成正比时, 从牛顿动力学方程出发导出北半球抛射体的运动方程, 最后从抛射体的运动方程着手, 讨论了一些特例, 结果所得, 一般文献将均被 本文包含. 2、质点飞行的轨迹和距离考虑空气阻力是抛射体质心的飞行轨迹问题[1]。假设空气阻力的方向与速度向量相 反，大小与速度的平方成正比。计算质点飞行的轨迹和距离。0 x x v mg y v V 2 cv  图1抛射体的受力图 [2]mg cvv mg cv dt dv m cvv cv dt dv m v dt dy v dt dx y y x x y x               sin cos 2 2 其中c为空气阻力系数。用MATLAB来进行数值求解。设r= ，则原方程等价于下面的以及微分方程组；   y x v v y x , , ,                         g r r r m c r r r m c r r dt dr 4 2 4 2 3 3 2 4 2 3 4 3 3、MATLAB仿真与分析 程序如下： (1)函数文件 zf.m function rdot=zf(t,r)c=0.02;g=9.8;m=1; vm=sqrt(r(3)^2+r(4)^2); rdot=[r(3);r(4);-c*vm*r(3)/m;-c*vm*r(4)/m-g]; (2)主程序ex12＿6.m clear; y0=0;x0=0; %初始位置 v0=( 请输入初始输入速度(m/s): ); rho=( 请输入初始方向(度): ); tf=( 请输入飞行时间(s): ); vx0=vo*cos(rho*pi/180); %计算x、y方向的初始速度 vy0=v0*sing(rho*pi/180); [t,r]=ode45( zf ,[0,tf],[0;0;vx0;vy0;]); %解微分方程 H=max(r(:,2)) %求轨道的最高点 T=t(find(r(:,2)==H)); %到最高点的所需时间 L=min(r(find(r(:,2)<0),1)); %计算射程 plot([0,100],[0,0]);hold on %绘制X坐标线 xlabel( x );ylabel(y); plot(r(:,1),r(:,2)); %绘制运功轨迹 请输入初始输入速度(m/s):57 请输入初始方向(度):36 请输入飞行时间(s):6.5 H