intersec.m

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

function [X_x,X_y]=intersec(p1_x,p1_y, p2_x,p2_y, p3_x,p3_y, p4_x,p4_y)

	dx1 = p2_x - p1_x;
	dx2 = p4_x - p3_x;
	if (dx1 == 0.0 & dx2 ~= 0.0)
		m2 = (p4_y - p3_y)/dx2;
		b2 = p3_y - m2*p3_x;
		X_x = p1_x;
		X_y = m2*X_x + b2;
		return;
	end
	if (dx1 ~= 0.0 & dx2 == 0.0)
		m1 = (p2_y - p1_y)/dx1;
		b1 = p1_y - m1*p1_x;
		X_x = p3_x;
		X_y = m1*X_x + b1;
		return;
	end
	if (dx1 == 0.0 & dx2 == 0.0)
		X_x = 0.000;
		X_y = 0.000;
		return;
	end
	m1 = (p2_y - p1_y)/dx1;
	b1 = p1_y - m1*p1_x;
	m2 = (p4_y - p3_y)/dx2;
	b2 = p3_y - m2*p3_x;
	if (m1 == m2)
		X_x = 0.000;
		X_y = 0.000;
		return;
	end
	X_x = -(b2 - b1)/(m2 - m1);
	X_y = m1*X_x + b1;
	return;