nt2rlc.m

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

% nt2rlc.m
% used by syn_fire.m in synthesis

clear i;
	glotoc = 1;
	exloc = 0;

	NING = 60;
	ISEC = 60/NING;
	NTS = round(29/ISEC);
	Nfft = 64;
	
	NTlen = 11;
	NTsno = 11;
	NTslen = NTlen/NTsno;
	
	NTArea(4) = 4.0;
	NTArea(5) = 6.0;
	NTArea(6) = 8.0;
	NTArea(7) = 8.0;
	NTArea(8) = 6.5;
	NTArea(9) = 4.0;
	NTArea(10) = 2.0;
	NTArea(11)= 2.0;
	NTArea(1) = nt1atemp + 0.0001;
	del_n = (NTArea(4) - NTArea(1))/3.0;
	NTArea(2) = NTArea(1) + del_n;
	NTArea(3) = NTArea(2) + del_n;

%	NTArea(2) = 2.0;
%	NTArea(3) = 3.0;

	NSIN = 0;
	Nsin(1)=0;
	Nsin(2)=0;



% function NTarea2RLC()

	w = 2*pi * 4000.0;
	rho = 1.14e-3;
	mu = 1.86e-4;
	c = 3.53e4;
	eta = 1.4;
	lambda = 5.5e-5;
	epsilon = 0.24;
	mass = 1.5;
	viscous = 1500;
	compilance = 0;
	rdem = sqrt(rho*mu*w/2.0);
	rc2 = rho*c*c;
	Gpc = (eta-1.0)*sqrt(lambda*w/(2.0*epsilon*rho))/rc2;

	for j=1:NTsno,
 
%   Perimeter = 2.0 * sqrt(pi * A) <----- circle
%            = 4.0 * sqrt(pi * A) <----- ellipse

		peri = 4.0 * sqrt(pi*NTArea(j));

		RNs(j) = peri * rdem/(NTArea(j)*NTArea(j));
		LNs(j) = rho / NTArea(j);
		CNp(j) = NTArea(j) / rc2;
		GNp(j) = peri * Gpc;

		wall = sqrt(pi)*NTArea(j)*NTslen; % Maeda */
%               wall = peri*NTslen; Flanagan */
		RNw(j) = viscous / wall;
		LNw(j) = mass / wall;
%               CNw(j) = wall / compilance;     */
	end
	RNr = (128.0*rho*c)/(9.0*pi*pi*NTArea(NTsno));
	LNr = 1.5 * ( 8.0*rho / (3.0*pi*sqrt(pi*NTArea(NTsno))) );