（转载）
% 天空之城完美钢琴版——Chfs = 44100;
%% sample rate
dt = 1/fs;
T16 = 0.125;
t16 = 0:dt:T16;
[~, k] = size(t16);
t1 = linspace(0,12*T16,12*k);%3拍
t2 = linspace(0,8*T16,8*k);%2拍
t4 = linspace(0,4*T16,4*k);%1拍
t8 = linspace(0,2*T16,2*k);%半拍
t17 = linspace(0,6*T16,6*k);%1拍半
t18 = linspace(0,(8/3)*T16,(8/3)*k);%1拍3连音其中1音
t19 = linspace(0,24*T16,24*k);%6拍
t20 = linspace(0,10*T16,10*k);%2拍半
[~, q] = size(t1);
[~, l] = size(t2);
[~, i] = size(t4);
[~, j] = size(t8);
[~, p] = size(t17);
[~, m] = size(t18);
[~, n] = size(t19);
[temp u] = size(t20);
mod1=(t1.^4).*exp(-30*(t1.^0.25));
mod1=mod1/max(mod1);
mod2=(t2.^4).*exp(-30*(t2.^0.3));
mod2=mod2/max(mod2);
mod4=(t4.^4).*exp(-30*(t4.^0.5));
mod4=mod4/max(mod4);
mod8=(t8.^4).*exp(-50*(t8.^0.75));
mod8=mod8/max(mod8);
mod16=(t16.^4).*exp(-110*(t16.^1));
mod16=mod16/max(mod16);
mod17=(t17.^4).*exp(-25*(t17.^0.5));
mod17=mod17/max(mod17);
mod18=(t18.^4).*exp(-35*(t18.^0.6));
mod18=mod18/max(mod18);
mod19=(t19.^4).*exp(-30*(t19.^0.2));
mod19=mod19/max(mod19);
mod20=(t20.^4).*exp(-30*(t20.^0.275));
mod20=mod20/max(mod20);
f0 = 2*146.8;
%% reference frequency（参考频率）
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 20/3 10/2 ...
17/12 19/12 17/6 ...%28
1/4 9/32 5/16 1/3 3/8 5/12 15/32 17/24 ...%36
19/48 19/24 19/16 19/8 17/8];
%% 1/2音符+1/8音符
so2u = mod20.*cos(2*pi*ScaleTable(17)*f0*t20);
%% 全音符+1/2音符
mi2n = mod19.*cos(2*pi*ScaleTable(15)*f0*t19);
%% 1拍3连音
la1m = mod18.*cos(2*pi*ScaleTable(10)*f0*t18);
ti1m = mod18.*cos(2*pi*ScaleTable(12)*f0*t18);
do2m = mod18.*cos(2*pi*ScaleTable(13)*f0*t18);
mi2m = mod18.*cos(2*pi*ScaleTable(15)*f0*t18);
blkm = zeros(1,m);
%% 3/4 音符
la1q = mod1.*cos(2*pi*ScaleTable(10)*f0*t1);
mi2q = mod1.*cos(2*pi*ScaleTable(15)*f0*t1);
so2q = mod1.*cos(2*pi*ScaleTable(17)*f0*t1);
ti2q = mod1.*cos(2*pi*ScaleTable(19)*f0*t1);
%% 1/2 音符
ss0g = mod2.*cos(2*pi*ScaleTable(38)*f0*t2);
do0g = mod2.*cos(2*pi*ScaleTable(21)*f0*t2);
la0g = mod2.*cos(2*pi*ScaleTable(3)*f0*t2);
do1g = mod2.*cos(2*pi*ScaleTable(5)*f0*t2);
mi1g = mod2.*cos(2*pi*ScaleTable(7)*f0*t2);
so1g = mod2.*cos(2*pi*ScaleTable(9)*f0*t2);
la1g = mod2.*cos(2*pi*ScaleTable(10)*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);
blkg = zeros(1,l);
%% 1/4 音符
so4f = mod4.*cos(2*pi*ScaleTable(33)*f0*t4);
ss0f = mod4.*cos(2*pi*ScaleTable(38)*f0*t4);
do0f = mod4.*cos(2*pi*ScaleTable(21)*f0*t4);
mi0f = mod4.*cos(2*pi*ScaleTable(23)*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);
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);
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);
ti2f = mod4.*cos(2*pi*ScaleTable(19)*f0*t4);
do3f = mod4.*cos(2*pi*ScaleTable(20)*f0*t4);
la3f = mod4.*cos(2*pi*ScaleTable(24)*f0*t4);
sf1f = mod4.*cos(2*pi*ScaleTable(26)*f0*t4);
sf2f = mod4.*cos(2*pi*ScaleTable(28)*f0*t4);
ss1f = mod4.*cos(2*pi*ScaleTable(27)*f0*t4);
blkf = zeros(1,i);
%% 1/8 音符
do4e = mod8.*cos(2*pi*ScaleTable(29)*f0*t8);
re4e = mod8.*cos(2*pi*ScaleTable(30)*f0*t8);
mi4e = mod8.*cos(2*pi*ScaleTable(31)*f0*t8);
fa4e = mod8.*cos(2*pi*ScaleTable(32)*f0*t8);
so4e = mod8.*cos(2*pi*ScaleTable(33)*f0*t8);
la4e = mod8.*cos(2*pi*ScaleTable(34)*f0*t8);
ti4e = mod8.*cos(2*pi*ScaleTable(35)*f0*t8);
ss4e = mod8.*cos(2*pi*ScaleTable(37)*f0*t8);
ss0e = mod8.*cos(2*pi*ScaleTable(38)*f0*t8);
sr1e = mod8.*cos(2*pi*ScaleTable(39)*f0*t8);
sr2e = mod8.*cos(2*pi*ScaleTable(40)*f0*t8);
sd2e = mod8.*cos(2*pi*ScaleTable(41)*f0*t8);
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*ScaleTable(5)*f0*t8);
re1e = mod8.*cos(2*pi*ScaleTable(6)*f0*t8);
mi1e = mod8.*cos(2*pi*ScaleTable(7)*f0*t8);
fa1e = mod8.*cos(2*pi*ScaleTable(8)*f0*t8);
so1e = mod8.*cos(2*pi*ScaleTable(9)*f0*t8);
la1e = mod8.*cos(2*pi*ScaleTable(10)*f0*t8);
ti1e = mod8.*cos(2*pi*ScaleTable(12)*f0*t8);
do2e = mod8.*cos(2*pi*ScaleTable(13)*f0*t8);
re2e = mod8.*cos(2*pi*ScaleTable(14)*f0*t8);
mi2e = mod8.*cos(2*pi*ScaleTable(15)*f0*t8);
so2e = mod8.*cos(2*pi*ScaleTable(17)*f0*t8);
la2e = mod8.*cos(2*pi*ScaleTable(18)*f0*t8);
ti2e = mod8.*cos(2*pi*ScaleTable(19)*f0*t8);
do3e = mod8.*cos(2*pi*ScaleTable(20)*f0*t8);
sf1e = mod8.*cos(2*pi*ScaleTable(26)*f0*t8);
sf2e = mod8.*cos(2*pi*ScaleTable(28)*f0*t8);
ss1e = mod8.*cos(2*pi*ScaleTable(27)*f0*t8);
blke = zeros(1,j);
%% 1/16 音符
do4s = mod16.*cos(2*pi*ScaleTable(29)*f0*t16);
fa4s = mod16.*cos(2*pi*ScaleTable(32)*f0*t16);
so4s = mod16.*cos(2*pi*ScaleTable(33)*f0*t16);
la4s = mod16.*cos(2*pi*ScaleTable(34)*f0*t16);
do0s = mod16.*cos(2*pi*ScaleTable(21)*f0*t16);
re0s = mod16.*cos(2*pi*ScaleTable(22)*f0*t16);
mi0s = mod16.*cos(2*pi*ScaleTable(23)*f0*t16);
fa0s = mod16.*cos(2*pi*ScaleTable(1)*f0*t16);
so0s = mod16.*cos(2*pi*ScaleTable(2)*f0*t16);
la0s = mod16.*cos(2*pi*ScaleTable(3)*f0*t16);
ti0s = mod16.*cos(2*pi*ScaleTable(4)*f0*t16);
do1s = mod16.*cos(2*pi*ScaleTable(5)*f0*t16);
re1s = mod16.*cos(2*pi*ScaleTable(6)*f0*t16);
mi1s = mod16.*cos(2*pi*ScaleTable(7)*f0*t16);
so1s = mod16.*cos(2*pi*ScaleTable(9)*f0*t16);
ti1s = mod16.*cos(2*pi*ScaleTable(12)*f0*t16);
do2s = mod16.*cos(2*pi*ScaleTable(13)*f0*t16);
re2s = mod16.*cos(2*pi*ScaleTable(14)*f0*t16);
mi2s = mod16.*cos(2*pi*ScaleTable(15)*f0*t16);
blks = zeros(1,k);
%% 3/8 音符
fa1p = mod17.*cos(2*pi*ScaleTable(8)*f0*t17);
la1p = mod17.*cos(2*pi*ScaleTable(10)*f0*t17);
ti1p = mod17.*cos(2*pi*ScaleTable(12)*f0*t17);
do2p = mod17.*cos(2*pi*ScaleTable(13)*f0*t17);
re2p = mod17.*cos(2*pi*ScaleTable(14)*f0*t17);
mi2p = mod17.*cos(2*pi*ScaleTable(15)*f0*t17);
fa2p = mod17.*cos(2*pi*ScaleTable(16)*f0*t17);
so2p = mod17.*cos(2*pi*ScaleTable(17)*f0*t17);
la2p = mod17.*cos(2*pi*ScaleTable(18)*f0*t17);
ti2p = mod17.*cos(2*pi*ScaleTable(19)*f0*t17);
do3p = mod17.*cos(2*pi*ScaleTable(20)*f0*t17);
mi3p = mod17.*cos(2*pi*ScaleTable(25)*f0*t17);
blkp = zeros(1,p);
%% 不奏响
blkblock = [blkg blkg blkg blkg...
blkg blkg blkg blkg];
blkblock1 = [blkm blkm blkm blkm blkm blkm...
blkm blkm blkm blkm blkm blkm...
blkblock...
blkm blkm blkm blkm blkm blkm...
blkm blkm blkm blkm blkm blkm];
%% 主旋律
piano1 = [la1e ti1e do2p ti1e do2e mi2p...
ti1g blkf mi1f la1p so1e la1e do2p...
so1g blkf mi1f fa1p mi1e fa1f do2f...
mi1g blke do2e do2e do2e ti1p sf1e sf1f ti1f...
ti1g blkf la1e ti1e do2p ti1e do2f mi2f...
ti1g blkf mi1f la1p so1e la1f do2f...
so1g blkp mi1e fa1f do2e ti1p do2f...
re2e re2f mi2e do2g do2e ti1e la1f ti1f ss1f...
la1g blkf do2e re2e mi2p re2e mi2f so2f...
re2g blkf so1e so1e do2p ti1e do2f mi2f...
mi2q blkf la1e ti1e do2f ti1e do2e re2e re2e...
do2p so1e so1g fa2f mi2f re2f do2f...
mi2n blkf mi2f la2p blke so2p so2e...
mi2e re2e do2g do2f re2f do2e re2p so2f...
mi2q mi2f la2p la2e so2p so2e...
mi2e re2e do2g blke do2e re2f do2e re2p ti1f...
la1q la2e ti2e do3p ti2e do3e mi3p...
ti2q mi2e mi2e la2p so2e la2e do3p...
so2q mi2e mi2e fa2p mi2e fa2f do3f...
mi2g blke do3e do3e do3e ti2p sf2e sf2f ti2f...
ti2q la1e ti1e do2p ti1e do2f mi2f...
ti1g blkf mi1f la1p so1e la1f do2f...
so1g blkp mi1e fa1f do2e ti1p do2f...
re2e re2f mi2e do2g do2e ti1e la1f ti1f ss1f...
la1q la1e ti1e do2e ti1e la1f ti1f ss1f...
la1q la3f];
%% 高音伴奏
piano7 = [blkp do2e ti1e la1e ti1f so1e la1e so1e do2f so2u...
blke do2e ti1e la1e ti1f fa1e mi1e so1e do2f so2u...
blke do2e ti1e la1e blke do2e blkf mi1e do2f so2u...
blke sf1e mi1e sr1e blke fa1e blkf blke ss1e sf1e mi1e ss1e mi1e fa1e ss1e...
blke do2e ti1e la1e ti1f so1e la1e so1e do2f so2u...
blke do2e ti1e la1e ti1f fa1e mi1e so1e do2f so2u...
blke do2e ti1e la1e blke do2e blkf blke do2e ti1e la1e blke do2e blkf...
blke fa1e blkf mi1f blkf blke do2e ti1e la1e blke do2e blkf...
blke so1e blke mi1e blke so1e blkf blke ti1e blke so1e blke ti1e la1e so1e...
blke la1e blke mi1e blke la1e blkf blke so1e blke mi1e blke so1s so1s blks so1s blke...
blke fa1e fa1f blke fa1e blke fa1e blke mi1e mi1f blke mi1s mi1s blks mi1s blke...
blke fa1e blke fa1e blke fa1e blke fa1e blke sf1e do2e ti1e sf1e do2e ti1e sf1e...
blke mi1e do2e ti1e mi1e do2e ti1e mi1e blke do2e ti1e la1e blke ti1e so1e la1e...
blke do2e ti1e do2e la1e ti1e so1e la1e blke ti1e la1e so1e blke ti1e fa1e so1e...
blke do2e ti1e do2e la1e ti1e ss1e la1e blke do2e ti1e la1e blke ti1e so1e la1e...
blke do2e ti1e do2e la1e ti1e so1e la1e blke ti1e la1e ti1e so1f blkf...
blke do2e ti1e do2e mi1e re1e do1e mi1e blke do2e ti1e la1e ti1f so1e la1e...
so1e do2f so2u blke do2e ti1e la1e ti1f fa1e mi1e...
so1e do2f so2u blke do2e ti1e la1e blke do2e blkf...
mi1e do2f so2u blke sf1e mi1e sr1e blke fa1e blkf...
blke ss1e sf1e mi1e ss1e mi1e fa1e ss1e blke do2e ti1e la1e ti1f so1e la1e...
so1e do2f so2u blke do2e ti1e la1e ti1f fa1e mi1e...
so1e do2f so2u blke do2e ti1e la1e blke do2e blkf...
blke do2e ti1e la1e blke do2e blkf blke fa1e blkf mi1f blkf...
blke do2e ti1e la1e blke do2e blkf blke fa1e blkf mi1f blkf...
blke do2e ti1e la1e blke do2e blkf];
piano8 = [blkp mi2e re2e do2e re2f ti1e do2e ti1e blkg blkf blke...
blke mi2e re2e do2e re2f la1e so1e ti1e blkg blkf blke...
blke mi2e re2e do2e blke mi2e blkf ti1e blkg blkf blke...
blke sr2e sd2e ti1e blke sr2e blkf blke mi2e re2e do2e mi2e ti1e do2e re2e...
blke mi2e re2e do2e re2f ti1e do2e ti1e blkg blkf blke...
blke mi2e re2e do2e re2f la1e so1e ti1e blkg blkf blke...
blke mi2e re2e do2e blke mi2e blkf blke mi2e re2e do2e blke mi2e blkf...
blke la1e blkf ss1f blkf blke mi2e re2e do2e blke mi2e blkf...
blke mi2e blke do2e blke mi2e blkf blke re2e blke ti1e blke re2e do2e ti1e...
blke do2e blke la1e blke do2e blkf blke ti1e blkp ti1s ti1s blks ti1s blke...
blke la1e la1f blke ti1e blke ti1e blke do2e do2f blke do2s do2s blks do2s blke...
blke re2e blke re2e blke re2e blke re2e blke la1e blkf la1e blkf la1e...
blke ss1e blkf ss1e blkf ss1e blke mi2e re2e do2e blke re2e ti1e do2e...
blke mi2e re2e mi2e do2e re2e ti1e do2e blke re2e do2e ti1e blke re2e la1e ti1e...
blke mi2e re2e mi2e do2e re2e ti1e do2e blke mi2e re2e do2e blke re2e ti1e do2e...
blke mi2e re2e mi2e do2e re2e ti1e do2e blke re2e do2e re2e ti1f blkf...
blke mi2e re2e mi2e do2e ti1e la1e blke blke mi2e re2e do2e re2f ti1e do2e...
ti1e blkg blkf blke blke mi2e re2e do2e re2f la1e so1e...
ti1e blkg blkf blke blke mi2e re2e do2e blke mi2e blkf...
ti1e blkg blkf blke blke sr2e sd2e ti1e blke sr2e blkf...
blke mi2e re2e do2e mi2e ti1e do2e re2e blke mi2e re2e do2e re2f ti1e do2e...
ti1e blkg blkf blke blke mi2e re2e do2e re2f la1e so1e...
ti1e blkg blkf blke blke mi2e re2e do2e blke mi2e blkf...
blke mi2e re2e do2e blke mi2e blkf blke la1e blkf ss1f blkf...
blke mi2e re2e do2e blke mi2e blkf blke la1e blkf ss1f blkf...
blke mi2e re2e do2e blke mi2e blkf];
piano9 = [blkp blke blke blke blkf blke blke blke blkg blkf blke...
blke blke blke blke blkf blke blke blke blkg blkf blke...
blke blke blke blke blke blke blkf blke blkg blkf blke...
blke blke blke blke blke blke blkf blke blke blke blke blke blke blke blke...
blke blke blke blke blkf blke blke blke blkg blkf blke...
blke blke blke blke blkf blke blke blke blkg blkf blke...
blke blke blke blke blke blke blkf blke blke blke blke blke blke blkf...
blke mi2e blkf mi2f blkf blke blke blke blke blke blke blkf...
blke blke blke blke blke blke blkf blke blke blke blke blke blke blke blke...
blke blke blke blke blke blke blkf blke blke blkp blks blks blks blks blke...
blke blke blkf blke blke blke blke blke blke blkf blke blks blks blks blks blke...
blke re2e blke re2e blke re2e blke re2e blke re2e blkf re2e blkf re2e...
blke re2e blkf re2e blkf re2e blke blke blke blke blke blke blke blke...
blke blke blke blke blke blke blke blke blke blke blke blke blke blke blke blke...
blke blke blke blke blke blke blke blke blke blke blke blke blke blke blke blke...
blke blke blke blke blke blke blke blke blke blke blke blke blkf blkf...
blke blke blke blke blke blke blke blke blke blke blke blke blkf blke blke...
blke blkg blkf blke blke blke blke blke blkf blke blke...
blke blkg blkf blke blke blke blke blke blke blke blkf...
blke blkg blkf blke blke blke blke blke blke blke blkf...
blke blke blke blke blke blke blke blke blke blke blke blke blkf blke blke...
blke blkg blkf blke blke blke blke blke blkf blke blke...
blke blkg blkf blke blke blke blke blke blke blke blkf...
blke blke blke blke blke blke blkf blke mi2e blkf mi2f blkf...
blke blke blke blke blke blke blkf blke mi2e blkf mi2f blkf...
blke blke blke blke blke blke blkf];
%% 低音伴奏
piano2 = [blke blke la4e mi0e la0e ti0e do1g...
mi4e mi0e la0e mi0e ti1g fa4e do0e fa0e so0e do0g...
do0e so0e do1e re1e mi1g fa4e do0e fa0e so0e la0g...
do4e so4e do0e re0e mi0f la0f mi4e ss4e ti4e mi0e ss0g...
mi4e mi0e ss0e ti0e mi1g la4e mi0e la0e ti0e do1f la0f...
mi4e mi0e la0e mi0e ti1f so1f fa4e do0e fa0e so0e la0f do0f...
do4e so4e do0e re0e mi0f do0f re4e la4e re0e fa0e do1f la0f...
la4e mi0e la0e ti0e do1f la0f fa4e do0e fa0e la0e so4e re0e so0e ti0e...
la4e mi0e la0e ti0e do1f la0f do0s so0s do1s re1s mi1e re1e do1e so0e mi0e do0e...
so4s re0s so0s la0s ti0e so0e re0e ti4e so4e mi4e la4s mi0s la0s ti0s do1e la0e mi0e do0e ti4e la4e...
do4s so4s do0s re0s mi0e so0e do1e so0e mi0e do0e la4e blkp ti4e blkp...
do0e so0e do1e re1e mi1f do1f fa4e do0e fa0e la0e ti0f la0f...
la4e mi0e la0e mi0e ss0f blkf la4e mi0e la0e mi0e ss0f blkf fa4s do0s fa0s la0s do1e fa1e re1e ti0e so0e re0e...
la4s mi0s la0s ti0s do1e mi1e do1e la0e mi0e do0e so4s re0s so0s ti0s re1e ti0e so0e re0e ti4e so4e...
do0e so0e do1e mi1e ss0f blkf fa4s do0s fa0s la0s do1e fa1e re1e ti0e so0e re0e...
la4s mi0s la0s ti0s do1e mi1e do1e la0e mi0e do0e so4s re0s so0s ti0s re1e ti0e so0e re0e so4f...
la4e mi0e la0e ti0e do1f la0f blkf blkf blkf blkf...
blkf blkf blkf blkf blkf blkf blkf blkf...
blkf blkf blkf blkf blkf blkf blkf blkf...
blkf blkf blkf blkf blkf blkf blkf blkf...
blkf blkf blkf blkf la4e mi0e la0e ti0e do1f la0f...
mi4e mi0e la0e mi0e ti1f so1f fa4e do0e fa0e so0e la0f do0f...
do4e so4e do0e re0e mi0f do0f re4e la4e re0e fa0e do1f la0f...
la4e mi0e la0e ti0e do1f la0f fa4e do0e fa0e la0e so4e re0e so0e ti0e...
la4e mi0e la0e ti0e do1f la0f fa4e do0e fa0e la0e so4e re0e so0e ti0e...
la4e mi0e la0e ti0e do1f la0f];
%% 前奏部分
piano4 = [ti1m ti1m la1m ti1m ti1m la1m...
ti1m ti1m la1m ti1m ti1m la1m...
mi1e blkp mi1e blkp mi1e blkp mi1e blkp...
mi1e blkp mi1e blkp mi1e do2e ti1s do2s la1e la1e ti1e blke ti1e...
ti1m ti1m la1m ti1m ti1m la1m...
ti1m ti1m la1m ti1m ti1m la1m];
piano5 = [do2m do2m do2m do2m do2m do2m...
do2m do2m do2m do2m do2m do2m...
la1e do2e ti1s do2s la1e la1e do2e blke la1e la1e do2e ti1s do2s ti1e la1e do2e blke la1e...
la1e do2e ti1s do2s la1e la1e do2e blke la1e la1e mi2e re2s mi2s blkf re2e blke re2e...
do2m do2m do2m do2m do2m do2m...
do2m do2m do2m do2m do2m do2m];
piano6 = [mi2m mi2m mi2m mi2m mi2m mi2m...
mi2m mi2m mi2m mi2m mi2m mi2m...
blkblock...
mi2m mi2m mi2m mi2m mi2m mi2m...
mi2m mi2m mi2m mi2m mi2m mi2m];
%% 主旋律钢琴1
v1 = [piano4 piano1];
%% 低音伴奏钢琴2
v2 = [blkblock1 piano2];
v3 = [piano6 piano1];
v4 = [piano5 piano1];
%% 高音伴奏钢琴2
v5 = [blkblock1 piano7];
v6 = [blkblock1 piano8];
v7 = [blkblock1 piano9];
%% 和音
s = v1+0.2*v2+v3+v4+0.4*v5+0.5*v6+0.5*v7;
s12m = max(max(s),abs(min(s)));
s = s./s12m;
sound(s,fs);
wavwrite(s,fs,32,'laputa.wav');