matlab演奏天空之城

%天空の城

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

T16 = 0.125; 

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

t1 = linspace（016*T1616*k）;
t34 = linspace（0 12*T16 12*k）;
t2 = linspace（08*T168*k）;
t38 = linspace（06*T166*k）;
t4 = linspace（04*T164*k）; 
t8 = linspace（02*T162*k）; 

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

% Modification functions 
mod1 = sin（pi*t1/t1（end））;
mod34 = sin（pi*t34/t34（end））;
mod38 = sin（pi*t38/t38（end））;
mod2 = sin（pi*t2/t2（end））; 
mod4 = sin（pi*t4/t4（end））; 
mod8 = sin（pi*t8/t8（end））; 

f0 = 2*146.8; % reference frequency

ScaleTable = [1 9/8 5/4 4/3 3/2 5/3 15/8 ...
2 9/4 5/2 8/3 3 10/3 15/4];

%1 notes
la04 = mod1.*cos（2*pi*ScaleTable（6）*f0*t1）; 
mi14 = mod1.*cos（2*pi*ScaleTable（10）*f0*t1）; 

%1/2 notes
do0d = mod2.*cos（2*pi*ScaleTable（1）*f0*t2）; 
re0d = mod2.*cos（2*pi*ScaleTable（2）*f0*t2）; 
mi0d = mod2.*cos（2*pi*ScaleTable（3）*f0*t2）; 
fa0d = mod2.*cos（2*pi*ScaleTable（4）*f0*t2）; 
so0d = mod2.*cos（2*pi*ScaleTable（5）*f0*t2）; 
la0d = mod2.*cos（2*pi*ScaleTable（6）*f0*t2）; 
si0d = mod2.*cos（2*pi*ScaleTable（7）*f0*t2）; 
do1d = mod2.*cos（2*pi*ScaleTable（8）*f0*t2）; 
re1d = mod2.*cos（2*pi*ScaleTable（9）*f0*t2）; 
mi1d = mod2.*cos（2*pi*ScaleTable（10）*f0*t2）; 
fa1d = mod2.*cos（2*pi*ScaleTable（11）*f0*t2）; 
so1d = mod2.*cos（2*pi*ScaleTable（12）*f0*t2）; 
la1d = mod2.*cos（2*pi*ScaleTable（13）*f0*t2）; 
si1d = mod2.*cos（2*pi*ScaleTable（14）*f0*t2）; 

%3/8 notes
do0t = mod38.*cos（2*pi*ScaleTable（1）*f0*t38）; 
re0t = mod38.*cos（2*pi*ScaleTable（2）*f0*t38）; 
mi0t = mod38.*cos（2*pi*ScaleTable（3）*f0*t38）; 
fa0t = mod38.*cos（2*pi*ScaleTable（4）*f0*t38）; 
so0t = mod38.*cos（2*pi*ScaleTable（5）*f0*t38）; 
la0t = mod38.*cos（2*pi*ScaleTable（6）*f0*t38）; 
si0t = mod38.*cos（2*pi*ScaleTable（7）*f0*t38）; 
do1t = mod38.*cos（2*pi*ScaleTable（8）*f0*t38）; 
re1t = mod38.*cos（2*pi*ScaleTable（9）*f0*t38）; 
mi1t = mod38.*cos（2*pi*ScaleTable（10）*f0*t38）; 
fa1t = mod38.*cos（2*pi*ScaleTable（11）*f0*t38）; 
so1t = mod38.*cos（2*pi*ScaleTable（12）*f0*t38）; 
la1t = mod38.*cos（2*pi*ScaleTable（13）*f0*t38）; 
si1t = mod38.*cos（2*pi*ScaleTable（14）*f0*t38）; 

%3/4 notes
so0q = mod34.*cos（2*pi*ScaleTable（5）*f0*t34）; 
si0q = mod34.*cos（2*pi*ScaleTable（7）*f0*t34）; 

%1/4 notes
do0f = mod4.*cos（2*pi*ScaleTable（1）*f0*t4）; 
re0f = mod4.*cos（2*pi*ScaleTable（2）*f0*t4）; 
mi0f = mod4.*cos（2*pi*ScaleTable（3）*f0*t4）;
fa0f = mod4.*cos（2*pi*ScaleTable（4）*f0*t4）; 
so0f = mod4.*cos（2*pi*ScaleTable（5）*f0*t4）;
la0f = mod4.*cos（2*pi*ScaleTable（6）*f0*t4）; 
si0f = mod4.*cos（2*pi*ScaleTable（7）*f0*