cint.m

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

function [cin_x,cin_y]=cint(c1_x,c1_y, rt, p2_x,p2_y, p1_x,p1_y)

	alf = angle(c1_x,c1_y, p2_x,p2_y) - angle(p1_x,p1_y, p2_x,p2_y);
	s = abs(sin(alf));
	so = (distance(c1_x,c1_y, p2_x,p2_y) * s)/rt;
	if (so == 1.0)
		om = pi/2.0;
	end
	if (so < 1.0)
		om = atan(so/sqrt(abs(1.0 - so*so)));
		tet = om - abs(alf);
	end
	if (so < 1.0 & alf >= 0.0)
		ang = abs(angle(p2_x,p2_y, c1_x,c1_y) + tet);
		[cin_x,cin_y] = polar2re(c1_x,c1_y, rt, ang);
	end
	if (so < 1.0 & alf < 0.0)
	[cin_x,cin_y] = polar2re(c1_x,c1_y, rt, angle(p2_x,p2_y,c1_x,c1_y)-tet);
	end
	if (so > 1.0)
		cin_x = 0.000;
		cin_y = 0.000;
	end