tract.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 50 行

% tract.m
% modifed by D. G. Childers 8/18/98
% update vocal tract parameter, e.g., area function and display waveforms
% tract(timeproc)
% This is used by syn_fire in synthesis.

espconst = 0.0000001;

wh=1.34651;
hk1=1.3431;
g1k=0.82081;

if (istartvt) 
	vfno = 1
%       arealeneqv(Ttemp);
%       arealeneqv(Tfsyn1);             /* OR arealeneqv2()--DING */
	nt1atemp=nt1a(vfno);
	for k=1:60,
		areatemp(k)=area(vfno,k);
		lengtemp(k)=leng(vfno,k);
	end
	syncaltf;
   if dpychoice==1
err_erg5_dpy=sagidpy(artset(vfno,:),wh,hk1,g1k,vfno,totfno,form_tar(vfno,:),syn_win2_h)/100;
	end     
else
	if vfno+1 <= totfno
		if (timeproc >= tartime(vfno+1)+espconst)
			vfno = vfno + 1
			nt1atemp=nt1a(vfno);
			for k=1:60,
				areatemp(k)=area(vfno,k);
				lengtemp(k)=leng(vfno,k);
			end
			syncaltf;
			if dpychoice==3
				figure(syn_win2_h);
				clf;
			end

			if dpychoice==1
err_erg5_dpy=sagidpy(artset(vfno,:),wh,hk1,g1k,vfno,totfno,form_tar(vfno,:),syn_win2_h)/100;
			end
		end
	end
	if vfno+1 <= totfno
		interpol;
	end
end