某大神用matlab写了D调《卡农》，之后出现了《最炫民族风》、《突然的自我》等。今天闲着没事写了一个《青花瓷》的，纯粹改的《卡农》，只是针对这首歌在程序里加了附点音符和二分音符，并且改成A调，只可惜后半部分转降B没做出来。不废话了，上程序！大家轻点拍砖。。。

fs = 44100; % sample rate

dt = 1/fs;

T16 = 0.25;

t16 = [0:dt:T16];

[temp k] = size(t16);

t8=linspace(0,2*T16,2*k);

t8fd=linspace(0,3*T16,3*k);

t4=linspace(0,4*T16,4*k);

t4fd=linspace(0,5*T16,5*k);

t2fd=linspace(0,7*T16,7*k);

t2=linspace(0,6*T16,6*k);

[temp i] = size(t4);

[temp j] = size(t8);

% Modification functions

mod2 = sin(pi*t2/t2(end));

mod2fd=sin(pi*t2fd/t2fd(end));

mod4 = sin(pi*t4/t4(end));

mod4fd=sin(pi*t4fd/t4fd(end));

mod8 = sin(pi*t8/t8(end));

mod8fd=sin(pi*t8fd/t8fd(end));

mod16 = sin(pi*t16/t16(end));

f0 = 2*146.8; % reference frequency

ScaleTable = 3*[2/3 3/4 5/6 15/16 ...

1 9/8 5/4 4/3 3/2 5/3 9/5 15/8 ...

2 9/4 5/2 8/3 3 10/3 15/4 4 ...

1/2 9/16 5/8]/((2^(1/12)))^5;

% half note

do1h = mod2.*cos(2*pi*ScaleTable(5)*f0*t2);

% dotted minim

so1hfd = mod2fd.*cos(2*pi*ScaleTable(9)*f0*t2fd);

% quarter note

do0f = mod4.*cos(2*pi*ScaleTable(21)*f0*t4);

re0f = mod4.*cos(2*pi*ScaleTable(22)*f0*t4);

mi0f = mod4.*cos(2*pi*ScaleTable(23)*f0*t4);

fa0f = mod4.*cos(2*pi*ScaleTable(1)*f0*t4);

so0f = mod4.*cos(2*pi*ScaleTable(2)*f0*t4);

la0f = mod4.*cos(2*pi*ScaleTable(3)*f0*t4);

ti0f = mod4.*cos(2*pi*ScaleTable(4)*f0*t4);

do1f = mod4.*cos(2*pi*ScaleTable(5)*f0*t4);

re1f = mod4.*cos(2*pi*ScaleTable(6)*f0*t4);

mi1f = mod4.*cos(2*pi*ScaleTable(7)*f0*t4);

fa1f = mod4.*cos(2*pi*ScaleTable(8)*f0*t4);

so1f = mod4.*cos(2*pi*ScaleTable(9)*f0*t4);

la1f = mod4.*cos(2*pi*ScaleTable(10)*f0*t4);

tb1f &#