从弦振动方程说开去——解开音乐的密匙

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

本文就当作是课外读物，有一点高数基础，学过多元函数微分学和常微分方程的读者都可以读 。下面是发表时的原文：

终于做完了，这次绝对是做的最艰难的一次！更多内容请移步https://m.bilibili.com/video/BV1r34y1i7TH
      这次带来的是数理方程的应用专题，整体来讲，包括前面第0节，对比前几个专题，难度不算太大。
      本学期完成了数理方程的学习任务，但学下来一知半解，想来想去，或许是因为没有明确一个道路。我们从物理现象中抽象出了数学方程，而后研究它的解，却没有真正明白这些稀奇古怪的解的形式背后，究竟蕴含着什么——考试也鲜少涉及这方面内容。换句话说，从物理中来，却没有到物理中去，得到解之后没有过多思考解的物理意义是什么，没有理解这个解为什么长这个样子。专题中的第0节和第1节便解决了一维弦振动方程两种解的物理意义问题。
      同时本次专题侧重于应用方面，以数形结合的方法，从视听角度带来更直观的感受，理论与实践结合，易于理解波动方程问题之于音乐发展的重要性。也应了一些人的要求，把习题穿插在每一节当中，承上启下，符合思考的顺序，并且适当降低了题目难度。有一些难度的内容都在上述链接中详尽地做了讲解。本次参考教材惯例，把所有的定义出来的名词全部加粗。
      本想再加上第4节的东西：通过集合论与群论分析和弦，两个和弦的关系，调式，对位这些音乐中最常见的元素。这些是之前阅读从图书馆借来的《勋伯格对位法:对位艺术的探秘之钥》一书得到的启发，想借此发挥一番，只可惜各种条件所限，不能在专题中与大家见面了。
      音乐总归是给人听的，它是由各个历史时期的人们实验总结出来的。也就是说本次专题讨论的这些内容，只是用数学或者物理等自然科学方法来深入研究音乐当中蕴含的自然规律。我们做任何分析不能离开思考的出发点，即什么样的音乐是好听的，什么样的是不好听的；好听的为什么好听，不好听的为什么不好听；它们对应什么样的情感。我们可以用这种音乐的语言来表达什么，或者我们可以根据这些语言来看作者想要表达什么。
      至于为何是最艰难的一次，因为我在音乐方面属实一窍不通，但又挺感兴趣，算是又蔡又爱玩了。写着写着，发现越到后面越难下笔：做这次专题不同以往，大量篇幅都跳出了数学专业，需在搜集整理音乐资料后，结合数学思维，自己再去亲自补充证明，写出自己的理解，这对我真是一大难事。
      身为应用数学的学生，在同样学了一大堆东西之后，竟要为如何应用而抓破头皮。记得此前多次为其他同学解答数学疑惑，思之无果，只好无奈折衷。有一次小梁问我近代数学的问题，我仍需要去图书馆查阅诸多书籍，照搬前人经验。光是自己想出头绪，不知不觉也要花费一下午的工夫。这便是所谓解答疑惑，实在愧矣。
      说起这个，又想到文中探讨两个音的关系的时候，为了解决钢琴调音问题，引入了拍音的概念。思之令人发笑——我竟也用了一下午的时间去听那《beat_sample》中的四个采样有何不同。直到我利用Adobe Audition去解析时，才找到一些端倪。这令我豁然开朗，随即编写了一道习题并加上了配图，让工作得以顺利进行下去。这道习题也许有一定难度，读者需要耐心地听，找到细微的变化，方能找到正确答案。
      此次一来时间仓促，二来水平有限，三来篇幅有限，错误和漏洞在所难免，若发现错误，或想加以补充，还望各位不吝赐教！

MATLAB程序：
x1=0:.0001:2;
y1=sin(2*pi*440*x1)+sin(2*pi*442*x1);
figure(1);
plot(x1,y1);
xlabel('t轴');ylabel('y轴');
title('图2.1 频率为440Hz与442Hz叠加的振动图像');
text(0.4,1,'1个单位时间内有2个拍音')
figure(2);
x2=0:.001:0.5;
f=44;eps=4;
y2=sin(2*pi*f*x2);
y3=sin(2*pi*(f+eps)*x2);
plot(x2,y2+y3,'r');
xlabel('t轴');ylabel('y轴');
hold on;
y=2*sin(2*pi*(f+eps/2)*x2);
plot(x2,y,'b');
hold on;
z=cos(pi*eps*x2);
plot(x2,2*abs(z),'k');
axis([0,0.5,-3,3]);
hold off;
title('图2.2 命题2.2证明过程中用到的函数，取f=44Hz,ε=4Hz，0<=t<=0.5')