main16.asv

语音信号子带编码源程序

% main function for test;
% 实现对四个子带的子带编码;
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);

[y11,y12]=bandpass(y1);[y13,y14]=bandpass(y2);[y15,y16]=bandpass(y3);[y17,y18]=bandpass(y4);
[y21,y22]=bandpass(y5);[y23,y24]=bandpass(y6);[y25,y26]=bandpass(y7);[y27,y28]=bandpass(y8);

y11=code(y11,8);y12=code(y12,8);y13=code(y13,8);y14=code(y14,4);
y15=code(y15,4);y16=code(y16,4);y17=code(y17,4);y18=code(y18,4);
y21=code(y21,4);y22=code(y22,4);y23=code(y23,4);y24=code(y24,4);
y25=code(y25,4);y26=code(y26,4);y27=code(y17,4);y28=code(y28,8);

dec=decode(y11,8);y11=dec;dec=decode(y12,8);y12=dec;dec=decode(y13,8);y13=dec;dec=decode(y14,4);y14=dec;
dec=decode(y15,4);y15=dec;dec=decode(y16,4);y16=dec;dec=decode(y17,4);y17=dec;dec=decode(y18,4);y18=dec;
dec=decode(y21,4);y21=dec;dec=decode(y22,4);y22=dec;dec=decode(y23,4);y23=dec;dec=decode(y24,4);y24=dec;
dec=decode(y25,4);y25=dec;dec=decode(y26,4);y26=dec;dec=decode(y27,4);y27=dec;dec=decode(y28,8);y28=dec;

[yy11,yy12]=restructor(y11,y12);[yy13,yy14]=restructor(y13,y14);[yy15,yy16]=restructor(y15,y16);[yy17,yy18]=restructor(y17,y18);
f11=result(yy11,yy12);f12=result(yy13,yy14);f13=result(yy15,yy16);f14=result(yy17,yy18);
[yy21,yy22]=restructor(y21,y22);[yy23,yy24]=restructor(y23,y24);[yy25,yy26]=restructor(y25,y26);[yy27,yy28]=restructor(y27,y28);
f21=result(yy21,yy22);f22=result(yy23,yy24);f23=result(yy25,yy26);f24=result(yy27,yy28);

[g1,g2]=restructor(f11,f12);[g3,g4]=restructor(f13,f14);[g5,g6]=restructor(f21,f22);[g7,g8]=restructor(f23,f24);
g11=result(g1,g2);g12=result(g3,g4);g13=result(g5,g6);g14=result(g7,g8);

[yy1,yy2]=restructor(g11,g12);[yy3,yy4]=restructor(g13,g14);
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:\manspeech2.wav');