数理方程关于某振动方程地分析报告matlab

《数理方程关于某振动方程地分析报告matlab》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案数理方程基于MATLAB的问题分析报告一、问题的提出、背景、意义文档大全实用标准文案振动是指物体经过它的平衡位置所作的往复运动或某一物理量在其平衡值附近的来回变动。而波动则是一种能量传播的方式。虽然形式不同，但是两者的联系十分紧密，振动是波动的根源，波动是振动的传播形式。因此在分析问题乃至实际操作中，往往是把两者放在一起分析的，首先讨论振动的各方面特性，这样就相当于已知了波动一点上的相应特性，再对波动进行分析时，就只用讨论距离的影响了。一般来说，振动只受时间影响，加上距离的参数，最终波动就只受两个变量影响，而且也

2、知道了它们是无关的，就可以使用分离变量法进行求解。弦振动是波动的一类特殊形式，它在音乐物理学、材料学、地理学、物质分析学等许多领域都得到了应用，而弦振动所属声学又是力学的一个非常独立的分支，因此它在各领域的作用几乎是不可取代的。由于近年来的各方面硬件设施和软件的发展，曾经停止发展很长一段时间的对弦振动的分析又开始体现出它独特的优势。在产生音乐的过程中，琴弦的振动是很常见的一种方式，本文就将对琴弦振动进行一定的研究，通过对弦振动方程的理解，给出不同初始条件，并分析出琴弦不同地方产生波的特性，再用MATLAB做好程序，画出相应的

3、图像，经比较后得到琴弦的拨发与产生声音的联系。二、问题分析思路2.1建立偏微分方程分析一根琴弦的振动问题，通过针对具体要分析的问题，可以列出弦振动方程以及初始条件（L为弦的长度，因为是两端固定的弦，初始条件一定有），用分离变量法很容易求得它相应的解，即弦振动的函数。2.2对琴弦参数的求解已知常量T=128N，普通钢琴弦密度，根据琴弦传播速度公式，可以求得速度v。文档大全实用标准文案2.3求解对象由弦振动的函数可以得到弦上不同点的振动情况。随机选取几个点，得到它们的振动情况，并比较。2.4作图方法通过MATLAB仿真出不同点的

4、图像，比较图像的幅值周期等参数。（开始考虑到有两种方式，一种直接通过上一个步骤求出的解使用简单的MATLAB命令画出图，另一种则是通过MATLAB解方程后再画出相应的图像，事实上第一种MATLAB是做不到的，于是用第二种）2.5仿真结果仿真出弦振动的频谱图，即以频率和振幅为横纵坐标的图，得到不同频率与振幅的关系，对图可以进行一系列的分析，得到相应的结果。2.6方程解的现实意义由于琴弦振动实际意义，我们将弦振动的实际音效也用MATLAB做出来了，这样更能直观的体会到琴弦振动条件不同带来的影响。但是发出的声音不如实际生活那么和谐

5、美妙（缺少腔体等音乐元件）。三、具体求解步骤3.1标准齐次弦振动的求解如前文所提，对于这样一个标准的齐次弦振动问题，分离变量法是我们主要所采取的解题方法。设方程具有的解的形式为：（3-1）将变量t与变量x分离开后，代入原方程，得到：文档大全实用标准文案（3-2）（3-3）令：（3-4）此时，得到两个常微分方程：（3-5），（3-6）代入边界条件，得到：（3-7）由于不是我们需要的解，对T（t）不能恒为0，所以对于X（x），我们可以得到：（3-8）这样一来，我们可以得到常微分方程满足边界条件的平凡解。当时，原方程的边值问题就只

6、有零解。当时，原方程的通解为：（3-9）代入边界条件，得：（3-10）（3-11）解得的结果为,A=0,。为了使X（x）不恒为0，应有，亦即，则：文档大全实用标准文案（3-12）相应的特征函数为，其中Bn为任意非零常数，对应每一个特征值方程的解是（3-13）其中，Cn,Dn为任意常数。我们得到原方程一系列特解为：（3-14）为了求出满足的解，我们将作傅立叶拓展，把每一项全部叠加起来，则：（3-15）为了确定系数Cn，Dn，将方程代入初始条件，得（3-16）之后即可解出Cn，Dn：（3-17）（3-18）3.2实际弦振动的求解

7、对于第二节一开始提出的一维实际琴弦振荡问题，我们将实际参数代入公式中。这里，取考虑到弦乐器的常见技法就是拨弦，拨弦即用手指把琴弦拨离平衡位置。使其振动发声。这相当于在X=a处把弦拉高到高度h，然后松开，使其自由振动，即弦振动的初始位移不为零而初速度为零。文档大全实用标准文案(1)假设在琴弦的正中间拨弦，则a=L/2，取值为434mm，拨弦高度h为4mm。可以得到：那么，此时的波动方程表达式为：（3-19）X为坐标，t为时间，如果我们取弦上三个具有代表性的点，根据琴弦的对称性，就可以大致了解整个弦的振动情况；为此，我们不妨选取

8、x=L/6,x=L/3,x=L/2三个点作为特征点。此时，分别令x=L/6,x=L/3,x=L/2，代入（3-19），有：（3-20）（3-21）（3-22）（2）若不是在琴弦的正中间拨弦，是在a=L/3处拨弦，则此时的波动方程为：（3-23）仅仅是最前端系数发生变化，不影响我们对问题的研