MATLAB计算平抛阻尼运动

小球在空气中平抛运动的轨迹

[问题]一小球在空气中作平抛运动,初速度为v0,所受的阻力与速率成正比:f = -kv,k称为阻力系数。画出小球运动的轨迹。

[数学模型]

如B7.1图所示,小球受到重力mg,方向向下;空气阻力f,方向与速度方向相反。根据牛顿第二定律可列出直角坐标方程 d2xdxd2ydy,m2 mg k。 (7.1.1) m2 kdtdtdtdt

由于vx = dx/dt,vy = dy/dt,上式可化为 B7.1图 dvdvxkk vx,y g vy。 (7.1.2) dtmdtm

分离变量得

dvydvxk dt, dt, vxmg kvy/m

积分得

lnvx lnCx kvkmt, [ln(g y) lnCy] t, mkm

kmgkt),vy [1 exp( t)]。 (7.1.3) mkm当t = 0时,vx = v0,vy = 0,可得Cx = -v0,Cy = -g,因此 vx v0exp( 当t = 0时,x = 0,y = 0,积分上式可得

mv0kmgm2gkx [1 exp( t)],y t 2[exp( t) 1]。 (7.1.4) kmkkm

这是小球的运动方程,也是以时间t为参数的轨道方程。

[算法]方法一:用解析式。取t0 = m/k为时间单位,取s0 = m2g/k2为坐标单位,坐标方程可表示为

x = s0v0*[1 – exp(-t*)],y = s0[t* + exp(-t*) - 1]。 (7.1.4*) 其中,t* = t/t0,为约化时间;v* = kv0/mg,是小球的约化初速度。

[程序]zqy6_3ode.m的第一部分如下。

%阻力与速度成正比的平抛运动的轨迹(用解析解)

clear %清除变量

v0=input('请输入水平初速度kv0/mg:'); %键盘输入水平初速度

t=0:0.1:4; %时间向量

x=v0*(1-exp(-t)); %有空气阻力的横坐标

y=t+exp(-t)-1; %有空气阻力的纵坐标

fs=16; %字体大小

figure %开创图形窗口

matlab 平抛运动,MATLAB计算平抛阻尼运动相关推荐

  1. c语言平抛编程试题,平抛运动练习题及答案

    平抛运动规律巩固练习 一.选择题(不定项): 1.关于平抛运动,下列说法正确的是 ( ) A.不论抛出位置多高,抛出速度越大的物体,其水平位移一定越大 B.不论抛出位置多高,抛出速度越大的物体,其飞行 ...

  2. matlab卷积矩阵绝对值,MATLAB矩阵分析和计算

    MATLAB矩阵分析和计算 编辑 锁定 讨论 上传视频 本词条缺少概述图,补充相关内容使词条更完整,还能快速升级,赶紧来编辑吧! <MATLAB矩阵分析和计算>是清华大学出版社出版的一本图 ...

  3. MATLAB调用refprop计算物性参数详解

    MATLAB调用refprop计算物性参数详解 欢迎使用Markdown编辑器 欢迎使用Markdown编辑器 REFPROP(REference Fluid PROPerties)是一款国际权威工质 ...

  4. 基于MATLAB的动态前景目标特征提取与运动跟踪

    基于MATLAB的动态前景目标特征提取与运动跟踪 摘 要:动态前景目标识别和提取是计算机视觉领域的重要内容.对动态图像进行前景目标提取与运动跟踪,通过改进高斯混合背景模型,提出一种基于自适应特征加权的 ...

  5. 悬置位移matlab,matlab悬置解耦计算工具

    matlab悬置解耦计算工具从名字就可以知道这款软件的功能就是帮助用户在软件分析车辆悬置解耦数据,你可以在软件将数据添加分析,选择需要分析的数值范围,从而立即在软件完成解耦计算,支持对数据查看,立即分 ...

  6. [MATLAB学习笔记]基于MATLAB的坐标系变换及飞行器姿态运动显示

    描述三维空间物体的运动通常是在指定的坐标系下进行,在不同的坐标系下物体运动的轨迹和姿态不尽相同.求解物体在不同坐标系下的位置和姿态的关键在于求解不同坐标系之间的变换矩阵. 以从地心坐标系到体坐标系为例 ...

  7. matlab闰年,基于MATLAB的万年历计算

    基于MATLAB的万年历计算 黄小川(200431610102) 测绘工程专业2004级9班 摘要:运用MATLAB编程实现1936-2032年农历与公历间的相互转换,并求当日是星期几. 关键词:万年 ...

  8. 在matlab中可以计算式子的最大值吗,matlab求最大值

    matlab遗传算法求最大值问题例题完整代码_IT认证_资格考试/认证_教育专区.matlab遗传算法求最大值问题例题完整代码 本文件的目的是减少您打字的烦恼 %主程序:用...... 11MATLA ...

  9. matlab圆形器件,计算围绕点+ Matlab的圆形箱

    我的问题与此链接stackoverflow问题有关 本质上,重复在此绘制的图形..我在图像中有一个中心点(x,y),我必须围绕它绘制四个1-4单位半径的圆,它们之间有8个角度. 在此图中,有12个角仓 ...

最新文章

  1. 5亿整数的大文件,怎么排?
  2. java临时的api,JAVA API-day03
  3. python闭包的原理_web前端:js 闭包原理
  4. 【MongoDB学习笔记21】MongoDB的复合索引
  5. Latex 中插入超链接 插入网址
  6. HappyLeetcode64:Sqrt(x)
  7. 题外:分类篇(音乐风格分类)基于BP神经网络
  8. GDKOI2018游记
  9. UnityShader22:序列帧动画
  10. oracle 九阴真精,《九阴真经》真的很强吗?其实它杂而不精,顶级高手都不愿意用...
  11. Unity游戏接入Steam成就
  12. C陷阱与缺陷(二)语义“陷阱”、连接
  13. 文献(2): 综述_癌症相关成纤维细胞(CAF)的异质性【建议收藏】
  14. Unity小游戏教程系列 | 创建小型太空射击游戏(一)
  15. php生成word,并下载
  16. 用完加速器国内的网址打不开了?
  17. 关于Qt 缺少xcb 解决办法
  18. win7计算机管理中看不到新加的硬盘,win7系统看不到第二块硬盘的解决方法.
  19. PHP爆绝对路径方法总结帖
  20. 宜品优享抽奖系统小程序,多方式真是展示抽奖功能

热门文章

  1. 黄毅然的网络安全学习
  2. win10网络上只能看见小部分电脑
  3. 10+31=100团队项目简介
  4. python文件第三关
  5. 方博机器人_贝尔机器人的学员这么优秀 难道你真的不知道吗
  6. html5输入框自动放大镜,HTML5实现放大镜效果
  7. 3D【1】人脸重建:Face Alignment Across Large Poses: A 3D Solution中侧脸数据生成实验
  8. iOS仿照Yahoo天气:油条天气
  9. 九月十月百度,迅雷,华为,阿里巴巴最新校招笔试面试三十题(10.18)
  10. ie8兼容性视图对ajax,ie8的兼容性视图设置在哪