MATLAB计算阻尼振动

弹簧振子的阻尼振动

[问题]一弹簧振子的质量为m,倔强系数为k。振子还受到与速度大小成正比、方向相反的阻力,比例系数为γ。当振子从静止开始运动时,初位移为A。物体的运动规律是什么?不同的阻尼下的运动曲线和速度曲线有什么特点?

[数学模型]

根据牛顿运动定律,物体运动的微分方程为

m

2

dxdt

2

2

kx

dxdt

, (6.2.1)

取k/m = ω0,γ/m = 2β,ω0就是无阻尼时物体的固有角频率,β是阻尼因子。物体的运动方程可表示为

dxdt

rt

22

2

dxdt

0x 0。 (6.2.2)

2

设微分方程的解为x = e,代入上式可得特征方程

r2 - 2βr + ω02 = 0。 (6.2.3)

特征方程的解为

r

, (6.2.4)

α可以是实数和零以及虚数,则r1 = -β + α,r2 = -β – α,r1和r2可以是

实数或复数。微分方程的解为

x C1e1 C2e

rt

r2t

e

t

(C1e

t

C2e

t

), (6.2.5)

其中C1和C2是由初始条件决定的常数。物体的速度为

v

dxdt

C1r1e1 C2r2e

rt

r2t

e

t

[C1( )e

t

C2( )e

t

]。 (6.2.6)

当t = 0时,x = A,v = 0,因此可得

A = C1 + C2,0 = C1(-β + α) + C2(-β - α), (6.2.7)

如果β ≠ ω0,即α ≠ 0,解得两个常数分别为

C1 A,C2 A。

2 2

因此物体的位移为

x

A2

e

t

[( )e

t

( )e

t

]。 (6.2.8)

[讨论]①当β > ω0时,即α > 0,上式就是过阻尼的情况。

②当β → ω0时,即α → 0,不论用罗必塔法则还是用公式eαt → 1 + αt和e-αt → 1 - αt,都可得

x A(1 0t)e

0t

。 (6.2.9)

这是临界阻尼的情况。

③当β < ω0

时,设

α = iω

,为虚数单位,利用欧拉公式

1

matlab弹簧振子的阻尼振动,MATLAB计算方法和技巧6_2阻尼振动相关推荐

  1. matlab动画_弹簧振子振动的matlab动画演示

    用Matlab阐述物理中的胡克定律,为弹簧振子震动的matlab动画示意图,另有一个为不同质量不同弹簧系数的比较.spring.m,compare.m rectangle('position',[12 ...

  2. MATLAB设计课题推荐及选择老师技巧

    MATLAB设计课题推荐及选择老师技巧 前言:时光稍纵即逝,一转眼又到了大四开学季,目前即将面临的是bs选题,以作为毕业前的最后一道工程项目.兴许你会说选题是下半年的事情,还早呢,下半年再说.其实你错 ...

  3. matlab求解焦点,Matlab焦点光场分布数值计算方法研究+源程序(2)

    Matlab焦点光场分布数值计算方法研究+源程序(2) 时间:2018-03-08 10:05来源:毕业论文 6.2 MATLAB程序运行结果分析29 结论 32 致谢 33 参考尔献34 1 引言 ...

  4. Matlab学习一本通,matlab基础教程

    链接:https://pan.baidu.com/s/1uTCbiRfIxcrt6lmiy6_QlQ  提取码:f2dn  Matlab学习一本通,matlab基础教程 <MATLAB R201 ...

  5. matlab dsp实验报告,matlab实验报告14.pdf

    matlab实验报告14 Matlab实验报告 院系名称 :信息科学与工程学院 专业班级 :通信工程 1303 指导老师 :陈科文 ,支国明 ,张金焕 ,周扬 学生姓名 : 学号 : 目录 实验一 熟 ...

  6. matlab com组件 发布,matlab生成com组件

    (中国科学院声学研究所北京100084) 摘要:介绍了采用MaflabCOMBuilder制作COM组件的方法,通过支持COM方式的编译工具调用生成的COM组件,能简 单方便地实现Matlab利其他. ...

  7. 面向对象:MATLAB的自定义类 [MATLAB]

    https://www.cnblogs.com/gentle-min-601/p/9785812.html 面向对象:MATLAB的自定义类 [MATLAB] 这几天刚刚开始学习MATLAB的面向对象 ...

  8. matlab接口编程配置,MATLAB与外部程序接口编程

    基本信息 书名:MATLAB与外部程序接口编程 定价:42.00元 作者: 出版社:机械工业出版社 出版日期:2009-01-00 ISBN:9787111257066 字数: 页码: 版次: 装帧: ...

  9. matlab中插入希腊,matlab中希腊字母打法

    如何在matlab中输入希腊字母_数学_自然科学_专业资料.matlab数据可视化处理,绘图技巧 如何在 matlab 中输入希腊字母 matlab 中用转义符来输入希腊字母的方法 上标用 ..... ...

最新文章

  1. (深入理解)model.eval() 、model.train()以及torch.no_grad() 的区别
  2. 【Java】什么是CAS、synchronized升级概述、偏向锁/轻量级锁详解 - 笔记
  3. 写个函数用来对二维数组排序
  4. Ext Designer Preview生成代码在html上显示步骤
  5. 3.15 晚会—「饿了么」之殇
  6. js禁用按钮_探索js让你的网页“自己开口说话”
  7. JAVA day06 酒店管理系统
  8. Android:读取assets目录下文件(二)
  9. java 文件分割器_Java分割器
  10. Selenium爬虫 -- WebDriver多标签页创建与切换
  11. 2020 智慧城市解决方案(智慧城市系统及相关技术)
  12. 20190818 思维导图简易教程
  13. 什么是3D建模?用到哪些软件?
  14. java搜索引擎框架_搜索引擎框架介绍
  15. 采蘑菇电脑c语言,英菲尼迪终于升级英菲尼迪Q50L,内行人告诉你怎么选还配备主动降噪、胎压显示!凯美瑞都比不上它! 早买早享受...
  16. html之行间样式和外部样式以及内部样式
  17. 有关/mnt/asec /mnt/secure文件夹及app2sd原理
  18. java中length,length()和size()各自含义
  19. 三元一次方程组计算机解法,math 三元一次方程组的解法
  20. FTP工作原理及内网用端口映射方式建FTP的注意事项

热门文章

  1. 机器学习西瓜书第一章总结
  2. 【flutter】打渠道包-多渠道脚本打包
  3. 树莓派HDMI无输出
  4. word表格分页时怎样能自动生成表头
  5. DANN:利用神经网络算法评估变异位点的有害程度
  6. ps保存html和图像格式不显示,photoshop保存web格式不能显示该怎样解决
  7. 在for循环中运行setTimeout的三种情况
  8. python运维自动化老男孩_老男孩第十四期 python 自动化运维第二周
  9. IE的各种设置列表及禁止修改IE主页
  10. Maximal Information Coefficient (MIC)最大互信息系数