load?matlab1.mat;
%sound（y11116000）;

y=zeros（2700001）;

[y111y222y333]=gcc（y111y222y333）;%?通過(guò)GCC來(lái)講三個(gè)信號(hào)對(duì)齊?

sigma=zeros（4801）;??%sigma是指數(shù)窗?。20ms一幀，16000*20ms=320個(gè)點(diǎn)為一幀。1/2重疊，即160點(diǎn)重疊，160個(gè)點(diǎn)不重疊
for?i=1:1:480;
sigma（i）=0.999999^i;
end

rn=zeros（4801）;?%rn?函數(shù)
rn（1）=1;
rn（241）=1;
rn（2:240）=2;

wlow=zeros（4801）;%低通函數(shù)?
wlow（1:200）=1;

for?i=1:840
????
????low=320*（i-1）+1;
????high=low+479;
????
????y1=y111（low:high）;
????y2=y222（low:high）;
????y3=y333（low:high）;
????
????y1=y1.*sigma;
????y2=y2.*sigma;
????y3=y3.*sigma;
????
????dcty1=dct（y1）;
????dcty2=dct（y2）;
????dcty3=dct（y3）;
????
????dcty1=log（abs（dcty1））;
????dcty2=log（abs（dcty2））;
????dcty3=log（abs（dcty3））;
????
????idcty1=idct（dcty1）;
????idcty2=idct（dcty2）;
????idcty3=idct（dcty3）;
????
????y0=（idcty1+idcty2+idcty3）/3;
????y0=y0.*wlow;
????y0=idct（exp（dct（y0）））;
????y0=y0./sigma;
????y（low:high）=y（low:high）+y0;
end

f=ones（201）;
f=f/20;



ylst=conv（yf）;
?sound（ylst16000）;
?
?
?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1058??2014-07-31?18:58??quhunxiangDCT.m