so.m

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

[so,fs,bits]=wavread('e:\so');
[soA,fs,bits]=wavread('e:\soA');
[sot,fs,bits]=wavread('e:\sot');
n=[0:(4.980-0.590)/(219630-26058):(4.980-0.590)];
e=exp(-n);
e=e';
soT=sot*ones(1,(219630-26058)/225);
soT=soT(:);
e=e*ones(1,2);
e(193501:193573,:)=[];
e(:,2)=[];
soS=e.*soT;
sof=soA;
sof(26059:(26058+193500))=soS;
sound(sof,fs);
wavwrite(sof,fs,'e:\sof');
subplot(2,1,1);
plot(so);
subplot(2,1,2);
plot(sof);