matlab神曲卡农

% Cripple Pachebel‘s Canon on Matlab
% Have fun

fs = 44100; % sample rate
dt = 1/fs;

T16 = 0.125;

t16 = [0:dt:T16];
[temp k] = size（t16）;

t4 = linspace（04*T164*k）;
t8 = linspace（02*T162*k）;

[temp i] = size（t4）;
[temp j] = size（t8）;

% Modification functions
mod4 = sin（pi*t4/t4（end））;
mod8 = sin（pi*t8/t8（end））;
mod16 = sin（pi*t16/t16（end））;

f0 = 2*146.8; % reference frequency

ScaleTable = [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];
% 1/4 notes
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 = mod4.*cos（2*pi*ScaleTable（11）*f0*t4）;
ti1f = mod4.*cos（2*pi*ScaleTable（12）*f0*t4）;
do2f = mod4.*cos（2*pi*ScaleTable（13）*f0*t4）;
re2f = mod4.*cos（2*pi*ScaleTable（14）*f0*t4）;
mi2f = mod4.*cos（2*pi*ScaleTable（15）*f0*t4）;
fa2f = mod4.*cos（2*pi*ScaleTable（16）*f0*t4）;
so2f = mod4.*cos（2*pi*ScaleTable（17）*f0*t4）;
la2f = mod4.*cos（2*pi*ScaleTable（18）*f0*t4）;
ti2f = mod4.*cos（2*pi*ScaleTable（19）*f0*t4）;
do3f = mod4.*cos（2*pi*ScaleTable（20）*f0*t4）;
blkf = zeros（1i）;

% 1/8 notes
fa0e = mod8.*cos（2*pi*ScaleTable（1）*f0*t8）;
so0e = mod8.*cos（2*pi*ScaleTable（2）*f0*t8）;
la0e = mod8.*cos（2*pi*ScaleTable（3）*f0*t8）;
ti0e = mod8.*cos（2*pi*ScaleTable（4）*f0*t8）;
do1e = mod8.*cos（2*pi*ScaleTable（5）*f0*t8）;
re1e = mod8.*cos（2*pi*ScaleTable（6）*f0*t8）;
mi1e = mod8.*cos（2*pi*ScaleTable（7）*f0*t8）;
fa1e = mod8.*cos（2*pi*ScaleTable（8）*f0*t8）;
so1e = mod8.*cos（2*pi*ScaleTable（9）*f0*t8）;
la1e = mod8.*cos（2*pi*ScaleTable（10）*f0*t8）;
tb1e = mod8.*cos（2*pi*ScaleTable（11）*f0*t8）;
ti1e = mod8.*cos（2*pi*ScaleTable（12）*f0*t8）;
do2e = mod8.*cos（2*pi*ScaleTable（13）*f0*t8）;
re2e = mod8.*cos（2*pi*ScaleTable（14）*f0*t8）;
mi2e = mod8.*cos（2*pi*ScaleTable（15）*f0*t8）;
fa2e = mod8.*cos（2*pi*ScaleTable（16）*f0*t8）;
so2e = mod8.*cos（2*pi*ScaleTable（17）*f0*t8）;
la2e = mod8.*cos（2*pi*ScaleTable（18）*f0*t8）;
ti2e = mod8.*cos（2*pi*ScaleTable（19）*f0*t8）;
do3e = mod8.*cos（2*pi*ScaleTable（20）*f0*t8）;
blke = zeros（1j）;

% 1/16 notes
fa0s = mod16.*cos（2*pi*ScaleTable（1）*f0*t16）;
so0s = mod16.*cos（2*pi*ScaleTable（2）*f0*t16）;
la0s = mod16.*cos（2*pi*ScaleTable（3）*f0*t16）;
ti0s = mod16.*cos（2*pi*ScaleTable（4）*f0*t16）;
do1s = mod16.*cos（2*pi*ScaleTable（5）*f0*t16）;
re1s = mod16.*cos（2*pi*ScaleTable（6）*f0*t16）;
mi1s = mod16.*cos（2*pi*ScaleTable（7）*f0*t16）;
fa1s = mod16.*cos（2*pi*ScaleTable（8）*f