main8.asv

语音信号子带编码源程序

 function main8
% 实现对四个子带的子带编码;
clear;close all;
[x1,fs,nb]=wavread('E:\MATLAB6p5\work\2006\manspeech.wav');
wc=3400;wc=wc/fs;
y=lowpass(x1,wc);
wavwrite(x1,fs/3,nb,'E:\manspeech1.wav');
[x,fs,nb]=wavread('E:\manspeech1.wav');
[y1,y2]=bandpass(x);m1=y1;m2=y2;

[a1,a2]=bandpass(m1);[a3,a4]=bandpass(m2);
[y1,y2]=bandpass(a1);[y3,y4]=bandpass(a2);[y5,y6]=bandpass(a3);[y7,y8]=bandpass(a4);
y1=code(y1,8);y2=code(y2,8);y3=code(y3,4);y4=code(y4,4);
y5=code(y5,4);y6=code(y6,4);y7=code(y7,16);y8=code(y8,16);
dec=decode(y1,8);y1=dec;dec=decode(y2,8);y2=dec;
dec=decode(y3,4);y3=dec;dec=decode(y4,4);y4=dec;
dec=decode(y5,4);y5=dec;dec=decode(y6,4);y6=dec;
dec=decode(y7,16);y7=dec;dec=decode(y8,16);y8=dec;

[yy1,yy2]=restructor(y1,y2);[yy3,yy4]=restructor(y3,y4);
f1=result(yy1,yy2);f2=result(yy3,yy4);
[yy1,yy2]=restructor(y5,y6);[yy3,yy4]=restructor(y7,y8);
f3=result(yy1,yy2);f4=result(yy3,yy4);

[yy1,yy2]=restructor(f1,f2);[yy3,yy4]=restructor(f3,f4);
ff1=result(yy1,yy2);ff2=result(yy3,yy4);

[yy1,yy2]=restructor(ff1,ff2);
re=result(yy1,yy2);

re=re/max(re);
wavwrite(re,3*fs,nb,'E:\manspeech_8.wav');