sbc2.asv

来自「语音信号子带编码源程序」· ASV 代码 · 共 31 行

function sbc2;
% function of subbands coding;
% accomplish 2 subbands coding of voice-signals;
% DPCM use 8,4 bits within eath subbands ;
clear;close all;
a=input('PLease input the bits vecter of the 2 subband-coding:');
a=input('请选择长度为4的字长向量，任一字长不能小于4。\n字长向量为：');
[m,n]=size(a);
if  n~=2||mean(a)<6||mean(a)>16
    display('You have input a error vection,try again!');
    [w,f,n]=wavread('E:/MATLAB6p5/work/w1.wav');
    wavplay(w/max(w),2*f);
    return; 
end
tic;
[x1,fs,nb]=wavread('E:\MATLAB6p5\work\2006\manspeech.wav');
wc=3400;wc=wc/fs;[m,n]=size(x1);
y=lowpass(x1,wc);
wavwrite(x1,fs/3,nb,'E:\manspeech1.wav');
[x,fs,nb]=wavread('E:\manspeech1.wav');
[y1,y2]=bandpass(x);
y1=code(y1,a(1));
y2=code(y2,a(2));
dec=decode(y1,a(1));y1=dec;
dec=decode(y2,a(2));y2=dec;
[yy1,yy2]=restructor(y1,y2);
f=result(yy1,yy2);
f=f/max(f);
max(f);[m1,n1]=size(f),
wavwrite(f,fs*3,max(a),'E:\manspeech2.wav');
toc,error=abs(sum((f(1:m)-x1))/sum(x1))