fs=44100;?
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）;

mod4=（t4.^4）.*exp（-30*（t4.^0.5））;
mod4=mod4*（1/max（mod4））;
mod8=（t8.^4）.*exp（-50*（t8.^0.5））;
mod8=mod8*（1/max（mod8））;
mod16=（t16.^4）.*exp（-90*（t16.^0.5））;
mod16=mod16*（1/max（mod16））;

f0=2*146.8;?

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

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