含有阻尼项的弦振动方程及其仿真.docxVIP

含有阻尼项的弦振动方程及其仿真 内容提要： 本文通过对古典吉他的琴弦振动情况建立数学物理方程，得到一个含有阻尼项的双 曲型方程的初边值问题，对解用Matlab进行仿真。最后依据弦振动方程的结果，列举 了在这种情况下几种泛音的位置，并结合该方程，对右手给出指导。 关键词 数学物理方程，Matlab ,驻波。 引言： 在弦乐器表演屮常用到泛音这样的一个技巧，即左手虚按琴弦，滤掉一部分波在琴 弦上形成驻波。比如在弦的三分点进行滤波，则波长的三倍不能被弦长整除的波，将会 被滤掉。但是在拨弦乐器的教学中，关于泛音的位置一直是老师们口口相传。而月?某些 泛音准确位置并不在拨弦乐器的站（山口）上，所以缺乏理论指导。 在国内的研究领域屮，韩佩琪《弦乐器泛音的分析及应用》一文屮只是对弹拨乐器 的空弦状态下进行求解而且忽略了空气的阻力，而且并没有结合列出的解给出演奏技巧 上的指导。而邱桂明《阻尼作川卜?的弦振动研究》的初边值条件并不符合乐器的条件。 另外在周伟《古典吉他演奏教程》以及相关的-?些吉他教学视频中只是提及了左手虚按 的位置，关于右手的位置没有给出一个指导。综上来看，国内研究领域，对定弦振动泛 音的理论研究尚处于一个盲区。然血一维双曲型微分方程的理论已经比较完善给本文提 供了理论依据，给硏究带来了可行性。 一、模型建立: b b Fig_2」 如图所示：琴弦的初始状态: 其小〃是弹拨弦与初始位置间的距离，是弹拨点距离原点的距离，/表示弦的长度。 弦的两端是静止不动的，从而边值条件：为W(O,t)=W(/,t)=O 其中t表示振动吋间。 列出方程： 其中：错误！未找到引用源。，而T表示琴眩松弛时的张力，错误！未找到引用源。表 示琴眩线密度。 边值条件： u(Ort) = u(Ct) = 0 初值条件： lt=o = 0, ° v 二、问题的求解 从物理上知道，一个复朵的振动往往可以分解成许多简单的振动的叠加。如弦振动 所发出的声音可以分解成各种不同频率的单音叠加。相应于每种单音，弦振动时波形保 持不变，从而当时间变化是个点的振幅做同步的变化，所以可以冇如下形式： 带入到原方程会得到： 分离变量： 等式左右两边相等，左边仅是t的函数，右边仅是X的函数，左右两边要相等，只冇等于同一 个常数才可能。设此常数为错误！未找到引用源。。则得到两个常微分方程。 矿? + 為钿; 厂?十巩： 得到以下通解： f(x) = Cl * sin(^A $ 洱)+ C2 w g(r) = ca * 庐為司 + C4 * 声 因为阻尼系数很小，所以 — 4* a2* (、注 */*入 一“2 f) +(73他 4 * * * 入一 接下来带入初边值条件 “(0,0 = f(0)g( 得： C 讥心/W 得： /(%) =弼俘尤)我= 在代入: =(- =(- (C3cos Q + C4stn0) + ^4 * ?2 * A - p2(C4cos 0 一 C3 =f(x)(—耳C3 + V4* ?2 *A- 显然错误！未找到引用源。 所以错误！未找到引用源。 再代入： 用正弦级数來求:W- 1｝小‘—6^=1=仙(竽)-kbrccos(判儲-hlsin I + Icos Asr 用正弦级数來求: W - 1｝小‘ —6 ^=1 =仙(竽)- kbrccos (判儲 -hl sin I + Icos Asr 一 krrcos ( @ - i)k2TC2 旺3」 2 l4 * G* * 点A A知￡版字函加(* “齐住f -护寸+4加(卜宀(常 从解的表达式來看，和式的结果应当有界，而外面的指数函数可以使得波的振幅逐 渐减小，与含有阻尼的弹簧震动相类似。其中每一个单音的波长和与之对应的频率，分 别为％ 1/2错误！未找到引用源。 三、仿真模拟 由于求和结果不是慕本初等函数，我们将用MATLAB采取级数的部分和来模拟. 在这一部分，我们统一设弦长 由于求和结果不是慕本初等函数，我们将用MATLAB采取级数的部分和来模拟. 在这一部分，我们统一设弦长1二0.64m,拨弦位置b二0.55m,弦的张力T二24.5N,弦的线 密度().()()5kg/m3. 首先模拟基音：既k= 1,2,34,5 Fig_5_l (注：Fig_5_l中左图,[0,0.64]轴代表弦长,[0,0.03]轴代表时间，而纵轴在弦上一点, 随时间的前进的位移.由图中选取的是距离左端点0.55处，即拨弦处的样本，而琴真 正的发音是通过共鸣箱传出来的声音，所以右图不是声音真正的波形，但在表示频率方 面具冇代表性。) 以下我们将结合方程的解來研究泛音的性质。 2分点泛音就是滤掉在1/2处没有波节的音：即“2,4,6,8,10 其不动点为0.32m处， 即吉他屮的第12品位置，可以看岀在诸多的泛音屮，屮点