IMA-ADPCM(16bit)语音压缩的matlab仿真

function [ yout ] = fun_adpcm( xin ) % adpcm压缩仿真-压缩过程 IncludePara; index = 0; prev_sample = 0; Len = length(xin); yout = zeros(size(xin)); for i = 1:Len cur_sample=xin(i); % 得到当前的采样数据 delta=cur_sample-prev_sample; % 计算出和上一个的增量 if (delta<0) delta = -delta; sb=8; else sb=0; % sb 保存的是符号位 end code =fix( 4*delta / step_table(index+1));%根据 steptable[] 得到一个 0~7 的值 if (code>7) code=7; % 它描述了声音强度的变化量 end yout(i) = bitor(code,sb); % 加上符号位保存起来 Edelta= fix((step_table(index+1)*code)/4) + fix(step_table(index+1)/8); if (sb==8) Edelta = -Edelta; end prev_sample = prev_sample + Edelta; index = index + index_adjust(code+1); % 根据声音强度调整下次取 steptable 的序号 if (index<0) index=0; % 便于下次得到更精确的变化量的描述 else if (index>88) index=88; end end end end