%%%%?SingASong
%%%%?by?圓
%%%%?樂理基礎
%%%%?音高
%%%%?C=do?D=re?E=mi?F=fa?G=sol?A=la?B=si
%%%%?鋼琴上一共有88個音高，將這些音高進行分組可分為
%%%%?小字一組，小字二組，小字三組，小字四組，小字五組
%%%%?小字組，大字組，大字一組，大字二組
%%%%?音樂中將中央C上放的A音定為標準音高，其頻率為f=440
%%%%?中央C是指小字一組c即a1的頻率為440
%%%%?音值
%%%%?全音符、二分音符、四分音符、八分音符、十六分音符、三十二分音符
%%%%?附點二分音符、附點四分音符、附點八分音符、附點十六分音符
clc
clear?all

t0=3;??????????????????????%%一個全音符的時長

fs=44100;??????????????????%%采樣率該值為MP3的標準采樣率
dt=1/fs;

t16=0:dt:t0/16;????????????%%十六分音符
N=length（t16）;
t8=linspace（0t0*2/162*N）;???%%八分音符
t4=linspace（0t0*4/164*N）;???%%四分音符
t2=linspace（0t0*8/168*N）;???%%二分音符
t1=linspace（0t0*16/1616*N）;??%%全音符

t8d=linspace（0t0*3/163*N）;???%%附點八分音符
t4d=linspace（0t0*6/166*N）;????%%附點四分音符
t2d=linspace（0t0*12/1612*N）;???%%附點二分音符
%%%%?音高與頻率的關系
p=1:88;
f=440*2.^（（p-49）/12）;???????%%?p為音高
%%%%?Modification?函數
mod1=sin（pi*t1/t1（end））;
mod2=sin（pi*t2/t2（end））;
mod4=sin（pi*t4/t4（end））;
mod8=sin（pi*t8/t8（end））;
mod2d=sin（pi*t2d/t2d（end））;
mod4d=sin（pi*t4d/t4d（end））;
mod8d=sin（pi*t8d/t8d（end））;
mod16=sin（pi*t16/t16（end））;
%?mod16d=sin（pi*t16d/t16d（end））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?全音符
%%大字二組
A2_1?=mod1.*cos（2*pi*f（1）*t1）;
A2s_1=mod1.*cos（2*pi*f（2）*t1）;
B2_1?=mod1.*cos（2*pi*f（3）*t1）;
%%大字一組
C1_1?=mod1.*cos（2*pi*f（4）*t1）;
C1s_1=mod1.*cos（2*pi*f（5）*t1）;
D1_1?=mod1.*cos（2*pi*f（6）*t1）;
D1s_1=mod1.*cos（2*pi*f（7）*t1）;
E1_1?=mod1.*cos（2*pi*f（8）*t1）;
F1_1?=mod1.*cos（2*pi*f（9）*t1）;
F1s_1=mod1.*cos（2*pi*f（10）*t1）;
G1_1?=mod1.*cos（2*pi*f（11）*t1）;
G1s_1=mod1.*cos（2*pi*f（12）*t1）;
A1_1?=mod1.*cos（2*pi*f（13）*t1）;
A1s_1=mod1.*cos（2*pi*f（14）*t1）;
B1_1?=mod1.*cos（2*pi*f（15）*t1）;
%%大字組
C_1?=mod1.*cos（2*pi*f（16）*t1）;
Cs_1=mod1.*cos（2*pi*f（17）*t1）;
D_1?=mod1.*cos（2*pi*f（18）*t1）;
Ds_1=mod1.*cos（2*pi*f（19）*t1）;
E_1?=mod1.*cos（2*pi*f（20）*t1）;
F_1?=mod1.*cos（2*pi*f（21）*t1）;
Fs_1=mod1.*cos（2*pi*f（22）*t1）;
G_1?=mod1.*cos（2*pi*f（23）*t1）;
Gs_1=mod1.*cos（2*pi*f（24）*t1）;
A_1?=mod1.*cos（2*pi*f（25）*t1）;
As_1=mod1.*cos（2*pi*f（26）*t1）;
B_1?=mod1.*cos（2*pi*f（27）*t1）;
%%小字組
c_1?=mod1.*cos（2*pi*f（28）*t1）;
cs_1=mod1.*cos（2*pi*f（29）*t1）;
d_1?=mod1.*cos（2*pi*f（30）*t1）;
ds_1=mod1.*cos（2*pi*f（31）*t1）;
e_1?=mod1.*cos（2*pi*f（32）*t1）;
f_1?=mod1.*cos（2*pi*f（33）*t1）;
fs_1=mod1.*cos（2*pi*f（34）*t1）;
g_1?=mod1.*cos（2*pi*f（35）*t1）;
gs_1=mod1.*cos（2*pi*f（36）*t1）;
a_1?=mod1.*cos（2*pi*f（37）*t1）;
as_1=mod1.*cos（2*pi*f（38）*t1）;
b_1?=mod1.*cos（2*pi*f（39）*t1）;
%%小字一組
c1_1?=mod1.*cos（2*pi*f（40）*t1）;
c1s_1=mod1.*cos（2*pi*f（41）*t1）;
d1_1?=mod1.*cos（2*pi*f（42）*t1）;
d1s_1=mod1.*cos（2*pi*f（43）*t1）;
e1_1?=mod1.*cos（2*pi*f（44）*t1）;
f1_1?=mod1.*cos（2*pi*f（45）*t1）;
f1s_1=mod1.*cos（2*pi*f（46）*t1）;
g1_1?=mod1.*cos（2*pi*f（47）*t1）;
g1s_1=mod1.*cos（2*pi*f（48）*t1）;
a1_1?=mod1.*cos（2*pi*f（49）*t1）;
a1s_1=mod1.*cos（2*pi*f（50）*t1）;
b1_1?=mod1.*cos（2*pi*f（51）*t1）;
%%小字二組
c2_1?=mod1.*cos（2*pi*f（52）*t1）;
c2s_1=mod1.*cos（2*pi*f（53）*t1）;
d2_1?=mod1.*cos（2*pi*f（54）*t1）;
d2s_1=mod1.*cos（2*pi*f（55）*t1）;
e2_1?=mod1.*cos（2*pi*f（56）*t1）;
f2_1?=mod1.*cos（2*pi*f（57）*t1）;
f2s_1=mod1.*cos（2*pi*f（58）*t1）;
g2_1?=mod1.*c