matlab国歌演奏源代码

clcclear all
fs = 44100; dt = 1/fs; f0 = 320; T16 = 0.125;
t16 = [0:dt:T16]; [temp m] = size（t16）;
t2 = linspace（08*T168*m）; [temp i] = size（t2）;
t4 = linspace（04*T164*m）; [temp j] = size（t4）;
t8 = linspace（02*T162*m）; [temp k] = size（t8）;
t12 = linspace（04/3*T164/3*m）; [temp l] = size（t12）;
% Modification functions
mod2 = sin（pi*t2/t2（end））;
mod4 = sin（pi*t4/t4（end））;
mod8 = sin（pi*t8/t8（end））;
mod12 = sin（pi*t12/t12（end））;
mod16 = sin（pi*t16/t16（end））;

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];

do0g = mod2.*cos（2*pi*ScaleTable（21）*f0*t2）;re0g = mod2.*cos（2*pi*ScaleTable（22）*f0*t2）;mi0g = mod2.*cos（2*pi*ScaleTable（23）*f0*t2）;fa0g = mod2.*cos（2*pi*ScaleTable（1）*f0*t2）;so0g = mod2.*cos（2*pi*ScaleTable（2）*f0*t2）;la0g = mod2.*cos（2*pi*ScaleTable（3）*f0*t2）;ti0g = mod2.*cos（2*pi*ScaleTable（4）*f0*t2）;
do1g = mod2.*cos（2*pi*ScaleTable（5）*f0*t2）;re1g = mod2.*cos（2*pi*ScaleTable（6）*f0*t2）;mi1g = mod2.*cos（2*pi*ScaleTable（7）*f0*t2）;fa1g = mod2.*cos（2*pi*ScaleTable（8）*f0*t2）;so1g = mod2.*cos（2*pi*ScaleTable（9）*f0*t2）;la1g = mod2.*cos（2*pi*ScaleTable（10）*f0*t2）;% tb1g = mod2.*cos（2*pi*ScaleTable（11）*f0*t2）;ti1g = mod2.*cos（2*pi*ScaleTable（12）*f0*t2）;
do2g = mod2.*cos（2*pi*ScaleTable（13）*f0*t2）;re2g = mod2.*cos（2*pi*ScaleTable（14）*f0*t2）;mi2g = mod2.*cos（2*pi*ScaleTable（15）*f0*t2）;fa2g = mod2.*cos（2*pi*ScaleTable（16）*f0*t2）;so2g = mod2.*cos（2*pi*ScaleTable（17）*f0*t2）;la2g = mod2.*cos（2*pi*ScaleTable（18）*f0*t2）;ti2g = mod2.*cos（2*pi*ScaleTable（19）*f0*t2）;
do3g = mod2.*cos（2*pi*ScaleTable（20）*f0*t2）;
I = zeros（1i）;

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）;
J = zeros（1j）;

do0e = mod8.*cos（2*pi*ScaleTable（21）*f0*t8）;re0e = mod8.*cos（2*pi*ScaleTable（22）*f0*t8）;mi0e = mod8.*cos（2*pi*ScaleTable（23）*f0*t8）;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*S