第9章：清华大学谷源涛MATLAB音乐合成.pdfVIP

大作业一 音乐合成 第一部分 1.1 合成东方红 合成源程序如下，程序参考谷老师上课内容(自己编完代码后跟老师的差不多，而且实在是 不知道还能做什么改进了。。。) function k= play() rhythm = [5 5 6 2 1 1 6 2;0 0 0 0 0 0 -1 0 ; 0.5 0.25 0.25 1 0.5 0.25 0.25 1;]; y=generaterhythm(rhythm, 8); sound(y); end function y = generaterhythm(rhythm, basetune) fs = 8000; soundpos = [0 2 4 5 7 9 11 0:12]; y = zeros(1, sum(rhythm(3,:))*fs + 1); curpos = 1; for count = 1 : length(rhythm) cursoundname = rhythm(1, count); cursoundpos = soundpos(cursoundname); curfreq = 220 * 2 .^ ((cursoundpos + basetune + 3) / 12 + rhythm(2, count)); cursound = generatetune(curfreq, rhythm(3, count), fs); y(1,curpos:(curpos + length(cursound) - 1)) = cursound; curpos = curpos + length(cursound); end end function y = generatetune (freq, time, fs) t = 0:1/fs:(time - 1/fs); y = zeros(1, length(t)); y=sin(t*freq*2*pi); end 1.2 用包络修正每个乐音，消除噪音 在1.1中听到有“啪”的杂声，下面通过加包络来消噪音。 使“啪”音效出的最简便的方法是，为每个音加上外形相似的包络，使不同音的连接处 的响度为零即可，题目中提示使用e指数保罗所以，我们添加包络程序如下 function y = amendment(p, l) % p is start position of a syllable and l is length of asyllable if p/l0.7 y=p/l*exp(-7*p/l); else y=0.7/0.3*exp(-4.9)*(1-p/l); end end 程序运行后，可明显发现没有“啪”音的出现，但由于不管音的长短，最终都要衰减到 零，这使音乐的真实感大大降低。 加过包络后的波形图如图1.2.1所示 图1.2.1 在两个不同音的交界处放大图如1.2.2所示，可见由于交界处响度为零，所以不会出现 大的高频分量 图1.2.2 在试验程序中发现，当交界处的响度非常小而非零时，“啪”音也不会出现，例如当交 界处响度为最大的10%时，“啪”也会消除，波形如图1.2.3所示 图1.2.3 1.3 用最简单的方法将1.2中的音乐分别升高和降低一个八度，并用resample 函数将上述音 乐升高半个音阶 根据乐理知识，每一个整八度的频率差为2倍，也就是说，升高一个八度就是把频率改 为原来的两倍，降低一个八度是将频率改为原来的二分之一，所以在频率标记上乘以2或除 以2即可。 例如 y=sin(t*2*freq*2*pi); %twice frequency y=sin(t*freq*pi); %half frequency 若要上升半个音阶，则频率上升为原来的 倍即可，所以利用resample函数 对原数据进行采样率改变，即 y=resample(y,100,106);