subglimp.m

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

% subglimp.m
% Calculate the impedance of the subglottal system, a Foster-chain circuit.

Nfft = 64;


	rho=1.14e-3; c=3.53e4;

	for k=1:totfsub,
		Rsubg(k) = R(k)*0.1;
		Lsubg(k) = L(k)*1.0e-5;
		Csubg(k) = C(k)*1.0e-7;
	end

	df = 2*pi*5000.0/Nfft;
	Rg = 0.33*rho*c;
	Lg = 1.6e-3; % Ps=5 cm H2O; Ag=.25 cm^2 */
	w = df;
	j = 0;

	while j <= Nfft,
      Zsub(l) = 0.0+0.0*i;
      Zsubf(1)=Zsub(1);
		Glt = Rg+w*Lg*i;
		for k=1:totfsub,
			Rinv = 1/Rsubg(k);
			tmp = Rinv+w*Csubg(k)*i;
			Linv = 1/(0.0+w*Lsubg(k)*i);
			tmp = tmp+Linv;
			tmp = 1/tmp;
         Zsub(l) = Zsub(l)+tmp;
         Zsubf(l) = Zsub(l)+tmp;         
		end
      Zsub(l) = Zsub(l)+Glt;
      Zsubf(l) = Zsub(l)+Glt;      
		j = j + 1;
		w = w + df;
	end