gen_aspi.m

这是一个用于语音信号处理的工具箱

% gen_aspi.m

function [ns]=gen_aspi(An_ah, An_afilt, pit_per)
% Generate aspiration noise.
% An_ah=60;An_afilt=10;pit_per=101;

scale = 1.0;
ptr_asp_delay=0;

% assign the aspiration gain parameter to frame parameter. 

ah = An_ah;
coeff = (An_afilt)*0.01;

for j=1:pit_per,
	rnd = randgau;
	per = rnd * ah;

	if (coeff > 0.0)
	  if (An_afilt > 0.0) 
		per = (1.0 - An_afilt/100.0) * per;
	  end        
	[per,ptr_asp_delay] = fir_iir(ptr_asp_delay,per,An_afilt/100.0,1);

	elseif (coeff < 0.0)
	  if (An_afilt > 0.0) 
		per = (1.0 - An_afilt/100.0) * per;
	  end        
	[per,ptr_asp_delay] = fir_iir(ptr_asp_delay,per, An_afilt/100.0,0);
	end

	if (coeff >= -1.0 & coeff < 0.0)
		scale = (1.0-coeff+(coeff*coeff))^0.5;
	elseif (coeff < 1.0 & coeff > 0.0)
		scale=(((1.0-coeff)*(1.0-coeff))/(1.0+coeff))^0.5;
	end

	ns(j) = per / scale;
end