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