prevsubg.m

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

% prevsubg.m
% previous subglottal source frame

currsub=eval(get(curres_e_h,'String'));
%currsub=currsub-1;

if (currsub - 1) >= 1

	figure(subglot_h);	
	currsub=currsub-1;
	set(subresi_e_h,'String',round(R(currsub)));
	set(subresi_h,'Value',R(currsub));
	set(subindu_e_h,'String',round(L(currsub)));
	set(subindu_h,'Value',L(currsub));
	set(subcapa_e_h,'String',round(C(currsub)));
	set(subcapa_h,'Value',C(currsub));
	set(curres_e_h,'String',currsub);
	set(totres_e_h,'String',totfsub);
	
	Re=get(subresi_h,'Value')*0.1;
	In=get(subindu_h,'Value')*1e-5;
	Ca=get(subcapa_h,'Value')*1e-7;
	b=1.0/(Re*Ca);
	c=1.0/(In*Ca);
	d=1.0/Ca;
	pole=round(sqrt(c)/(2*pi));
	set(pole_e_h,'String',pole);
	A = 2*c + 2*d^2 / Re^2 - b^2;
	x1 =(A + sqrt(A^2-4*c^2)) / 2.0;
	x2 =(A - sqrt(A^2-4*c^2)) / 2.0;
	BW =round(abs(sqrt(x1)-sqrt(x2))/(2.0*pi));
	set(bw_e_h,'String',BW);

end