1、 物理模型建立

由课题问题可转化为实际数理问题,即典型弦振动问题,一根均匀弦两端分别在x=0以及x=L处固定,设初始速度为零,初始时刻弦的形状为一抛物线,抛物线的顶点为(L/2,h),讨论其弦振动的具体问题。<注:初值设定方式不收限制,但是不能全为零,否则无运动,在此以此初值是便于进一步的公式推导>
(L为弦长)

2、相关原理推导

2.1弦振动方程推导

由第一节所得出的物理实际问题可列出下列方程:

分离变量法:
u(x,t)=T(t)X(x)
代入原方程得:

求解可得:

将初始条件:代入得
且Cn等于:
因此:
至此,我们公式原理推导完毕,得出了弦振动的函数表达式:

2.2琴弦发音原理推导

经过多个文献的学习得出,小提琴声音是通过把琴弦的振动经过琴码传递到共鸣体,使得共鸣腔体产生共振,带动箱内外空气振动而产生的,而弦的振动则来源于琴弓与琴弦的相互作用。小提琴发声产生的振动大致分为弦的振动、琴码的振动和共鸣箱的振动三种,在此我们仅研究琴弦的振动,正是因此我们所仿真出来的发音效果会带来一定的偏差,请谅解。琴弦振动分为横振动、纵振动、扭转振动和倍频振动,根据课题要求与实际运算难度,在此我们只研究横振动。
横振动,将弦挑离其平衡位置再放掉,弦就开始做一个扁纱锭型的振动,它的振幅限制在两条明确的曲线之内。弦的横振动频率,可以用泰勒公式来表达,即
其中 a=√(F/ρ)
式中f为弦的振动频率,ρ、F、l,依次为弦的线密度、拉力和长度。且看出横振动频率与弦长成反比,与张力的平方根成正比,与弦线的密度的平方根成反比。

3、 相关参数的确立

有效琴弦长度L:《小提琴制作艺术》【英】克里斯.约翰逊,罗伊.考特诺尔著一书中讲到安装琴颈时琴颈站位与琴身站位比例为2比3,即F孔两个小缺口下口325毫米处。在此我们定弦长L=325mm=0.325m<注:此处我们选用的是4/4M尺寸的小提琴,即成人常用小提琴>
琴弦拉力F:由下表1可得H311型号E-Mi弦拉力为18.6磅。

表1琴弦参数
因此我们定为18.6磅= 8.436kg= 82.6728N
琴弦线密度:已知E-Mi钢弦的密度为7.8g/cm3=7800kg/m3,E弦直径为0.276mm则半径为0.138mm=0.000138m,由公式线密度=密度x琴弦的横截面积得,E弦的线密度是4.67x10^-4kg/m
至此,我们得到参数:L=0.325m,F=82.6728N,ρ线=4.67x10^-4kg/m

4、 运算检验与修正

首先判断基频,利用公式 将第三节的参数代入得

647Hz近似等于题干的标准基频660Hz,证明了我们数据的有效性与可靠性,在此我们做一下微小修正将拉力F从82.6728N改为86N可使基频调整至660Hz标准。
至此,修正后的参数为:L=0.325m,F=86N,ρ线=4.67x10^-4kg/m

5、Matlab仿真检验,习题解答

5.1 Matlab仿真琴弦发音

我们将利用matlab将振动函数与相关参数结合仿真出实际声音通过测量设备鉴别频率,过程如下见图5.1:

图5.1
测量结果如下见图5.2:

图5.2
因实际测量存在误差经多次测量去平均后得出matlab仿真发音频率确实为660Hz,在此需要关注图5.1代码31行与32行,由于2.2提到小提琴发声并不是单纯的琴弦振动所以仅仅模拟出琴弦函数声音偏差太大,在此我们用32行的典型振动函数并附上我们特定的参数仅用以进行发声判断,后期的图像处理我们仍将采用准确的弦振动方程得出的函数。
至此我们的理论阶段已基本完成可以进行更深度的探讨,下面开始习题的解答。

5.2 习题解答

习题1,当小提琴的发音为高八度A音(2x440Hz)时,弦长是多少?由题干可得多个变量要控制不变,仅仅允许弦长的变化,同时由公式

得当f从660Hz变化到880Hz时,弦长L要变为原来的3/4,即弦长从0.325m变为0.325mx0.75=0.24375m。我们进行仿真验证结果如下见图5.3:
图5.3
频率测定结果如下见图5.4:

图5.4
第一问完成。
习题2,音准受哪些因素影响,影响多大?
又前几节的发音函数可知,理论上讲音准受弦长、拉力与线密度影响。<实际上还要考虑我们在matlab中的抽样频率,即时间序列,但现实生活中不用>该问题要用到我们在学到的误差传递公式即可轻松得出,见下图5.5所示:

图5.5
因此,
∆f/f=∆L/L+1/2∙∆F/F+1/2∙∆ρ/ρ
第二题解决。
习题3,检验自己在发生多大的设计误差是能听得出来。
该题,我们认为无论你改变哪个量最后都影响的基频f因此,通过改变f的大小来多次试验判断音准的变化。经过三位组员500hz-700hz,700hz-900hz,900hz-1000hz且每次变化5Hz我们得出以下数据见表1。
判断音准 频率界值(Hz)
Do 500-575
Ra 575-640
Mi 640-685
Fa 685-750
So 750-850
La 850-950
Xi 950-1000
表一:音频辨别表
习题4,时间序列取值间隔应注意哪些问题,为什么?
首先,时间序列的起点与终点即振动发音的起点与终点,不难理解我们最好要讲起点设为0s以后均可,因为理想情况下无阻力一直在振动,但是终点理论上讲要比起点大1/20K秒即可,因为人耳可听频率为20Hz-20KHz,但是由于matalab软件仿真不完全等同于现实问题,因此,经试验时间长度必须大于等于0.002秒才可以听出不同声音。<我们推测这一数值和电脑硬件性能也有很大关系,因为这影响着matlab其内部的算力>
至于采样频率,它当然越高越好,因为离散点越密集越像连续图像真实性越高,当然那样运算量也越大,经测验5kHz的采样率仅能达到人们讲话的声音质量。11kHz的采样率是播放小段声音的最低标准,是CD音质的四分之一。22kHz采样率的声音可以达到CD音质的一半,目前大多数网站都选用这样的采样率。
44kHz的采样率是标准的CD音质,可以达到很好的听觉效果。由于电脑运算量强大,在此我们完全可以把采样频率设置为100kHz来仿真。

6、拓展与延伸

6.1 调弦、按指与迷之神奇的a
每一根琴弦的空弦都对应各自的振动频率(决定音是否准),所谓的调弦,即在演奏前调整弦的拉力松紧和扭转也就是改变变量拉力F,当然也可以更换材料即改变变量线密度ρ。
按指可以认为是控制声音的发声与停止,(其实并未停止),当按下琴弦时,拉力瞬间增大,设想如果拉力无穷大那么就不存在横振动了,所以在按下琴弦时声音变小,阻力增大,如果及时恢复,声音也可近似恢复,但是按压时间过长就已经停止振动了,有必要说一下,按指在小提琴演奏中是常见的一种指法可以发出美妙的颤音。
迷之神奇的a a=√(F/ρ)
通过单位的量纲可得a的单位为m/s即速度单位,且通过查找文献得出,a为振动的传播速度。
6.2 图像与动画
我们用2.1节所推导的函数放入相关参数进行matlab仿真效果图如下所示:
点击下方图标即可播放,如无法播放根目录附有gif文件。

波形图与强度对比图如下见图6.1图6.2:

图6.1

图6.2
最后附上一首matlab制作的一闪一闪亮晶晶以此结束。(点击图标打开播放)

7、代码

最后,如果有问题请在评论去留言我会及时回复!谢谢

(数理方程)Matlab模拟琴弦振动发声并显示振动图像(matlab)_人工智能相关推荐

  1. matlab模拟薄膜振动,(数理方程)Matlab模拟琴弦振动发声并显示振动图像

    1.物理模型建立 由课题问题可转化为实际数理问题,即典型弦振动问题,一根均匀弦两端分别在x=0以及x=L处固定,设初始速度为零,初始时刻弦的形状为一抛物线,抛物线的顶点为(L/2,h),讨论其弦振动的 ...

  2. matlab 模拟风,三维空间相关风场的计算机模拟及Matlab程序实现

    第 14 卷第 2 期2008 年 6 月 空 间 结 构 SPATIAL STRUCTURES Vol. 14 No. 2 Jun. 2008 收稿日期:200610208. 作者简介:刘文洋(19 ...

  3. 【光学】基于matlab模拟涡旋光与球面波的干涉含Matlab源码

    1 简介 基于涡旋光与球面波的干涉原理,提出一种物体微位移的光学测量方法.改进马赫泽德干涉光路,其中一束光照射至空间光调制器产生涡旋光束作为参考光,另一束光经透镜变为球面波后照射至物体上,两束光干涉后 ...

  4. matlab模拟巴特沃斯滤波器设计,巴特沃斯滤波器matlab实现

    描述 巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦,没有起伏,而在阻频带则逐渐下降为零. 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大 ...

  5. matlab 对mnist手写数字数据集进行判决分析_人工智能TensorFlow(十四)MINIST手写数字识别...

    MNIST是一个简单的视觉计算数据集,它是像下面这样手写的数字图片: MNIST 每张图片还额外有一个标签记录了图片上数字是几,例如上面几张图的标签就是:5.0.4.1. MINIST数据 MINIS ...

  6. matlab振动方程,数理方程关于振动方程的分析matlab.doc

    数理方程关于振动方程的分析matlab.doc 1数理方程基于MATLAB的问题分析报告一.问题的提出.背景.意义振动是指物体经过它的平衡位置所作的往复运动或某一物理量在其平衡值附近的来回变动.而波动 ...

  7. 数理方程及MATLAB解算学习笔记

    数理方程及MATLAB解算学习笔记 文章目录 数理方程及MATLAB解算学习笔记 第一章 MATLAB基础知识 1.class查询数值类型 2.永久性数值变量 3.创建特殊矩阵的专用指令 4.基本初等 ...

  8. 【数理方程】定解问题

    数理方程 数学物理方程是研究物理.力学.工程技术及其他自然科学时,经一些简化得到的一些偏微分方程,反映客观世界物理量之间关系.这些偏微分方程研究物理量在空间的分布规律和在时间的变化规律. 课程内容:3 ...

  9. 加加速度matlab,关于加加速度的若干机械运动分析及MATLAB模拟

    收稿日期:2004 - 04 - 06 基金项目:遵义师范学院科研基金资助项目(2003012) 作者简介:董水金(1946 -) ,男 ,江苏无锡人 ,遵义师范学院物理系副教授 ,主要从事电子技术和 ...

最新文章

  1. oracle sql MERGE 关键字
  2. 算法(第4版)Algorithms, Fourth Edition
  3. 新的 FreeBSD 核心团队选举结果出炉
  4. 如何設定 VS2008 讓 HTML 在格式化排版時不要破壞版面
  5. Omi框架学习之旅 - 生命周期 及原理说明
  6. uml 时序图_UML学习-14种UML图
  7. python 批量下载图片_Python 批量下载图片示例
  8. PS教程第二十三课:套索工具
  9. Linux用户、用户组、文件权限设置
  10. 编程三角形面积公式_利用边长计算三角形面积 — 海伦公式推导
  11. flask response对象
  12. iOS核心动画详解swift版----基础动画
  13. 异步请求---Get
  14. php 5.5 preg replace,解决ecshop在php5.5以上环境preg_replace报错方法
  15. 单片机c51语言显示help,单片机入门-C51语言LCD12864液晶显示器的基本使用介绍
  16. 我的世界神秘时代研究笔记计算机,我的世界神秘时代5教程全研究笔记连线一览...
  17. 谷歌统计Google Analytics使用入门
  18. IOS11更新教程大分享(附带降级教程)
  19. 给大忙人看的码农翻身记
  20. 抓取微博热搜榜数据并保存在Excel中

热门文章

  1. Linux 之 FIFO
  2. simics 使用指南
  3. 嵌入式工程师学习路线
  4. JetsonNano 网络配置记录
  5. P7599 [APIO2021] 雨林跳跃
  6. n代表二进制数字的长度,m代表二进制数字中1的个数,返回所有符合条件的数
  7. ubuntu 无线网卡启动命令
  8. 安安猜价格聪明机器人_3.1.3 安安猜价格行为编辑
  9. 李彦宏:一个人的百度
  10. USB2.0 HUB眼图调试经验总结