si.m

基于Matlab的八度音阶发音演奏程序

[si,fs,bits]=wavread('e:\si');
[siA,fs,bits]=wavread('e:\siA');
[sit,fs,bits]=wavread('e:\sit');
n=[0:(4.984-0.722)/(219799-31881):(4.984-0.722)];
e=exp(-n);
e=e';
siT=sit*ones(1,(219799-31881)/179);
siT=siT(:);
e=e*ones(1,2);
e(187772:187919,:)=[];
e(:,2)=[];
siS=e.*siT;
sif=siA;
sif(31882:(31881+187771))=siS;
sound(sif,fs);
wavwrite(sif,fs,'e:\sif');
subplot(2,1,1);
plot(si);
subplot(2,1,2);
plot(sif);