melp_init.asv

语音编码中的MELP的经典估计

%CONST
global FMCQ_CODEBOOK; 
global Wf;
global ham_win;
global stage1 stage2;
ham_win=hamming(200)';    			            %汉明窗系数
coeff;
stage;                        					%MSVQ四级本码(stage1(128),stage2(64,64,64))               
melp_wf;                      					%傅立叶幅度加权距离系数
codebook_fmcq1;	              					%傅立叶谱幅度码本(256)
codebook_fmcq2;

%source
%changed by jiang
%path=input('Please input drive[C:,D:...]\path..\path\filename(the type is .wav):\n    ','s');
%s=wavread(path)';        %%%%%%%%%%%%%%%
s=wavread('male1.wav')';
s=s*32767;               %%%%%%%%%%%%%%%
FRL=180;             %length of frame
Nframe=fix(length(s)/FRL);  %compute the frame number of the input file
%global sig_in
sig_in(1:FRL*2)=0;          					   %前一帧信号
%global cheb_s;
cheb_s(1:4)=0;  					               %60Hz四阶车比雪夫高通滤波器的初始状态

%global sig_1000
sig_1000(1:FRL*2)=0;
butter_s(1:6)=0;

%global melp_bands state_b state_e statet
melp_bands(1:5,1:FRL*2)=0;            		   %前一帧的五个子带信号
state_b(1:5,1:6)=0;                          %带通滤波器的状态
state_e(1:4,1:2)=0;                          %全波整形滤波器的初始状态
state_t(1:4,1:6)=0;        				      %包络检波中间过程的状态
%global melp_envelopes
melp_envelopes(1:4,1:FRL*2)=0;              	%前一幀的四个包络信号
%global pre_intp
pre_intp=40;                  					%前一帧的整数基音
%pre_intr=0.1;                 					%前一帧的整数基音对应的相关系数
frame_num=320;                 					%取的幀数
%global buffer
buffer=[50,50,50];                           %最近三个最强的基音值
%global pavg G2p
pavg=50;                                     %平均基音
G2p=20;                                      %前一帧的第2增益