极客matlab弹奏清华大学校歌

清华百年校庆，牛人无数。有人用 matlab 演奏清华老校歌，还有板有眼。笔者在其基础上修改了代码，使其能够在播放完之后生成一个 wav 文件。对于不想自己去 matlab 里面试一下的同学，可以直接听我的这个版本。

http://rene1989.web-126.com/wp-content/uploads/2011/04/tsinghua.wav

如果想和原版对照，请点击：清华校歌欣赏

下面的代码中有两个可爱的：（形成的表情，请大家粘贴的时候注意哈！

%code

function playlxg()

%% play “老校歌”

rhythm1 = [1 1 3 5 5 6 1 6 5 5 3 3 5 3 1 6 1 3 5 6 6 6 1 5 3 2 3 2 1 2 5 5 14 5 6 6 7 6 5 1 1 6 1 5 6 5 6 6 5 3 2 2 3 5 1 1 1 3 2 3 2 1 6 6 5 3 2 2 3 1 1 1 6 6 5 5 6 5 2 3 5 1 1 6 6 5 5 6 5 2 3 2 1;

0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0;

1 0.5 0.5 1 1 1 0.5 0.5 1 1 1 1 0.5 0.5 1 1 0.5 0.5 2 1 1 0.5 0.5 1 1 1 0.5 0.5 1 1 1 0.5 0.5 1 1 1 0.5 0.5 4 1.5 0.5 1 1 1.5 0.5 2 1.5 0.5 1 1 1 0.5 0.5 2 1.5 0.5 1 1 1 0.5 0.5 2 1.5 0.5 1 1 1 0.5 0.5 2 1 3 1 3 1 1 1 1 1 1 2 1 3 1 3 1 1 1 1 1 0.5 0.5 2];

rhythm2 = [5 1 3 3 6 6 5 5 1 1 3 1 5 6 6 1 7 4 4 4 3 1 7 1 7 1 7 2 2 2 1 1 1 7 3 3 4 4 3 4 3 4 1 2 1 7 7 1 7 5 5 5 1 7 7 5 1 1 2 1 7 7 5 3 3 1 1 3 3 4 2 7 1 7 3 3 1 1 3 3 4 2 7 7 5;

-1 0 0 0 0 0 0 0 0 0 0 0 -1 -1 -1 0 -1 0 0 0 0 0 -1 0 -1 0 -1 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 0 0 0 -1 -1 0 -1 -1 -1 -1 0 -1 -1 -1 0 0 0 0 -1 -1 -1 0 0 0 0 0 0 0 0 -1 0 -1 0 0 0 0 0 0 0 0 -1 -1 -1;

1 1 1 1 1 1 1 1 1 1 0.5 0.5 1 1 0.5 0.5 2 1 1 1 1 1 1 0.5 0.5 1 1 1 1 1 1 1 1 4 1.5 0.5 1 1 1.5 0.5 2 1.5 0.5 1 1 1 0.5 0.5 2 1.5 0.5 1 1 1 1 2 1.5 0.5 1 1 1 1 2 1 3 1 3 1 1 1 1 1 1 2 1 3 1 3 1 1 1 1 1 1 2];

rhythm3 = [3 3 5 1 1 6 6 1 1 5 5 1 5 3 14 4 4 5 1 1 1 6 1 5 5 5 4 3 5 7 7 7 14 4 4 5 5 5 1 1 1 1 1 1 4 5 5 5 5 5 3 3 3 5 5 4 3 4 4 5 5 5 4 3 5 5 4 4 1 1 1 7 5 5 5 5 5 4 4 1 1 1 7 5 5 4 3;

0 0 0 1 1 0 0 1 1 0 0 1 0 0 0 0 0 0 1 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0;

1 0.5 0.5 1 1 1 1 1 1 1 1 0.5 0.5 1 1 0.5 0.5 2 1 1 0.5 0.5 1 1 1 0.5 0.5 1 1 1 1 1 1 1 1 4 1.5 0.5 1 1 1.5 0.5 2 1.5 0.5 1 1 1 1 2 1.5 0.5 1 1 1 1 2 1.5 0.5 1 1 1 1 2 1 3 1 3 1 1 1 1 1 1 2 1 3 1 3 1 1 1 1 1 0.5 0.5 2];

rhythm4 = [1 1 1 1 4 4 1 1 1 1 1 1 2 2 2 5 4 4 4 1 1 5 5 1 7 7 2 2 2 2 2 5 1 1 4 6 1 1 1 4 4 7 1 5 5 5 1 1 1 1 5 5 1 4 4 7 1 5 5 1 1 1 4 4 1 1 4 5 5 5 5 1 1 4 4 1 1 4 5 5 5 1;

0 0 0 0 -1 -1 0 0 0 0 0 0 0 0 0 -1 -1 -1 -1 0 0 -1 -1 0 -1 -1 0 0 0 0 0 -1 0 0 -1 -1 0 0 0 -1 -1 -1 0 -1 -1 -1 0 0 0 0 -1 -1 0 -1 -1 -1 0 -1 -1 0 0 0 -1 -1 0 0 -1 -1 -1 -1 -1 0 0 -1 -1 0 0 -1 -1 -1 -1 0;

1 1 1 1 1 1 1 1 1 1 1 1 1 0.5 0.5 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4 1.5 0.5 1 1 1.5 0.5 2 1.5 0.5 1 1 1 1 2 1.5 0.5 1 1 1 1 2 1.5 0.5 1 1 1 1 2 1 3 1 3 1 1 1 1 1 1 2 1 3 1 3 1 1 1 1 1 1 2];

y1 = generaterhythm(rhythm1, 3); % in F tune

y2 = generaterhythm(rhythm2, 3); % in F tune

y3 = generaterhythm(rhythm3, 3); % in F tune

y4 = generaterhythm(rhythm4, 3); % in F tune

%plot(y);

y=(y1+y2+y3/2+y4/2)/3;

sound(y); % play it!

wavwrite(y,’tsinghua.wav’);

%% generate a piece of rhythm

function y = generaterhythm(rhythm, basetune)

% GENERATERHYTHM function.

% rhythm is a 3 * len matrix

% for each column, the element on the first row

% shows the ‘sound name’ of one tune, the one

% on the second row shows the sound step of that tune,

% and the one on the third row shows how long

%in seconds it should lasts.

fs = 8000;

soundpos = [0 2 4 5 7 9 11 0:12];

y = zeros(1, sum(rhythm(3,:))*fs + 1); % initialize

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

%% generate a piece of sound with the frequence of freq

function y = generatetune (freq, time, fs)

t = 0:1/fs:(time - 1/fs);

y = 0.6 * sin(t*freq*2*pi) + 0.1 * sin(t*freq*4*pi) + 0.2 * sin(t*freq*6*pi);

for count1 = 1 : length(y)

y(1, count1) = y(1, count1) * amendment(count1/length(y));

end

%% amendment of the shape of the wave

function y = amendment(p) % 0 <= p <= 1

if p < 0.2

y = p * 5;

elseif p < 0.3

y = 1.8 - p * 4;

elseif p < 0.8

y = 0.6;

else

y = 3 - p * 3;

end

%end code