griddpy.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 502 行 · 第 1/2 页

	ry = distance(pme_x,pme_y, H_x,H_y);
	dltv = 0.25*(tdl - angle(H_x,H_y, pme_x,pme_y));
	if (hy > 0.0 & pme_y <= PP_y & pme_y > H_y)
		dltv = dltv + pi/2.0;
	end
	[P84_x,P84_y] = polar2re(pme_x,pme_y, ry, tdl-(dltv*1.0));
	[P12_x,P12_y] = polar2re(pme_x,pme_y, ry, tdl-(dltv*2.0));
	[P82_x,P82_y] = polar2re(pme_x,pme_y, ry, tdl-(dltv*3.0));
	[P16_x,P16_y] = midpoint(PP_x,PP_y, DL_x,DL_y);
	[P86_x,P86_y] = midpoint(PP_x,PP_y, P16_x,P16_y);
	[R2_x,R2_y] = midpoint(PP_x,PP_y, P86_x,P86_y);
	[AS2_x,AS2_y] = midpoint(P16_x,P16_y, P86_x,P86_y);
	[P88_x,P88_y] = midpoint(DL_x,DL_y, P16_x,P16_y);
	[AS4_x,AS4_y] = midpoint(P16_x,P16_y, P88_x,P88_y);
	[R4_x,R4_y] = midpoint(DL_x,DL_y, P88_x,P88_y);
	P81_x = W_x;
	P81_y = P82_y;
	P11_x = W_x;
	P11_y = P12_y;
	P83_x = W_x;
	P83_y = P84_y;
	R1_x  = W_x;
	R1_y  = R2_y;
	P85_x = W_x;
	P85_y = P86_y;
	AS1_x = W_x;
	AS1_y = AS2_y;
	P15_x = W_x;
	P15_y = P16_y;
	AS3_x = W_x;
	AS3_y = AS4_y;
	P87_x = W_x;
	P87_y = P88_y;
	R3_x  = W_x;
	R3_y  = R4_y;
	Z = angle(V_x,V_y, cv_x,cv_y);
	Q = (Z - angle(M_x,M_y, cv_x,cv_y))/12.0;
	E = distance(cv_x,cv_y, M_x,M_y);
	[P27_x,P27_y] = polar2re(cv_x,cv_y, E, Z-(Q*1.0));
	[P29_x,P29_y] = polar2re(cv_x,cv_y, E, Z-(Q*2.0));
	[P31_x,P31_y] = polar2re(cv_x,cv_y, E, Z-(Q*3.0));
	[P33_x,P33_y] = polar2re(cv_x,cv_y, E, Z-(Q*4.0));
	[P35_x,P35_y] = polar2re(cv_x,cv_y, E, Z-(Q*5.0));
	[P37_x,P37_y] = polar2re(cv_x,cv_y, E, Z-(Q*6.0));
	[P39_x,P39_y] = polar2re(cv_x,cv_y, E, Z-(Q*7.0));
	[J1_x,J1_y] = polar2re(cv_x,cv_y, E, Z-(Q*8.0));
	[J3_x,J3_y] = polar2re(cv_x,cv_y, E, Z-(Q*9.0));
	[J5_x,J5_y] = polar2re(cv_x,cv_y, E, Z-(Q*10.0));
	[J7_x,J7_y] = polar2re(cv_x,cv_y, E, Z-(Q*11.0));

	[P26_x,P26_y] = cint(tongc_x,tongc_y, rt, V_x,V_y, cv_x,cv_y);
	ccv_x = cv_x;
	ccv_y = cv_y;
	if (P26_y < DL_y)
	ccv_x = tongc_x;
	ccv_y = tongc_y;
	[P26_x,P26_y] = polar2re(tongc_x,tongc_y, rt, angle(V_x,V_y, tongc_x,tongc_y));
	end
	tdl = angle(DL_x,DL_y, tongc_x,tongc_y);
	dltv = (tdl - angle(P26_x,P26_y, tongc_x,tongc_y)) * 0.125;
	[R6_x,R6_y] = polar2re(tongc_x,tongc_y, rt, tdl-(dltv*1.0));
	[R8_x,R8_y] = polar2re(tongc_x,tongc_y, rt, tdl-(dltv*2.0));
	[R10_x,R10_y] = polar2re(tongc_x,tongc_y, rt, tdl-(dltv*3.0));
	[R12_x,R12_y] = polar2re(tongc_x,tongc_y, rt, tdl-(dltv*4.0));
	[P20_x,P20_y] = polar2re(tongc_x,tongc_y, rt, tdl-(dltv*5.0));
	[P22_x,P22_y] = polar2re(tongc_x,tongc_y, rt, tdl-(dltv*6.0));
	[Z24_x,Z24_y] = polar2re(tongc_x,tongc_y, rt, tdl-(dltv*7.0));

	wv_x = W_x;
	wv_y = V_y;
	[R11_x,R11_y] = midpoint(wv_x,wv_y, P17_x,P17_y);
	[R7_x,R7_y] = midpoint(R11_x,R11_y, P17_x,P17_y);
	[R5_x,R5_y] = midpoint(R7_x,R7_y, P17_x,P17_y);
	[R9_x,R9_y] = midpoint(R11_x,R11_y, R7_x,R7_y);
	[P21_x,P21_y] = midpoint(wv_x,wv_y, R11_x,R11_y);
	[P19_x,P19_y] = midpoint(P21_x,P21_y, R11_x,R11_y);
	[P23_x,P23_y] = midpoint(wv_x,wv_y, P21_x,P21_y);
	if (P21_y > W_y)
		[P21_x,P21_y] = intersec(P21_x,P21_y, P22_x,P22_y, W_x,W_y, V_x,V_y);
	end
	if (P23_y > W_y)
		[P23_x,P23_y] = intersec(P23_x,P23_y, Z24_x,Z24_y, W_x,W_y, V_x,V_y);
	end
	[P28_x,P28_y] = cint(tongc_x,tongc_y, rt, P27_x,P27_y, ccv_x,ccv_y);
	[P30_x,P30_y] = cint(tongc_x,tongc_y, rt, P29_x,P29_y, ccv_x,ccv_y);
	[P32_x,P32_y] = cint(tongc_x,tongc_y, rt, P31_x,P31_y, ccv_x,ccv_y);
	[P34_x,P34_y] = cint(tongc_x,tongc_y, rt, P33_x,P33_y, ccv_x,ccv_y);
	[P36_x,P36_y] = cint(tongc_x,tongc_y, rt, P35_x,P35_y, ccv_x,ccv_y);
	[P38_x,P38_y] = cint(tongc_x,tongc_y, rt, P37_x,P37_y, ccv_x,ccv_y);
	[P40_x,P40_y] = cint(tongc_x,tongc_y, rt, P39_x,P39_y, ccv_x,ccv_y);
	[J2_x,J2_y] = cint(tongc_x,tongc_y, rt, J1_x,J1_y, ccv_x,ccv_y);
	[J4_x,J4_y] = cint(tongc_x,tongc_y, rt, J3_x,J3_y, ccv_x,ccv_y);
	[J6_x,J6_y] = cint(tongc_x,tongc_y, rt, J5_x,J5_y, ccv_x,ccv_y);
	[J8_x,J8_y] = cint(tongc_x,tongc_y, rt, J7_x,J7_y, ccv_x,ccv_y);
	[Z42_x,Z42_y] = cint(tongc_x,tongc_y, rt, M_x,M_y, ccv_x,ccv_y);
	c1b = distance(c1_x,c1_y, B_x,B_y);
	c2t = distance(c2_x,c2_y, T_x,T_y);
	if (B_x < P34_x)
		[P34_x,P34_y] = cint(c1_x,c1_y, c1b, P33_x,P33_y, ccv_x,ccv_y);
	end
	if (B_x < P36_x)
		[P36_x,P36_y] = cint(c1_x,c1_y, c1b, P35_x,P35_y, ccv_x,ccv_y);
	end
	if (B_x < P38_x)
		[P38_x,P38_y] = cint(c1_x,c1_y, c1b, P37_x,P37_y, ccv_x,ccv_y);
	end
	if (B_x < P40_x)
		[P40_x,P40_y] = cint(c1_x,c1_y, c1b, P39_x,P39_y, ccv_x,ccv_y);
	end
	if (B_x < Z42_x)
		[Z42_x,Z42_y] = cint(c1_x,c1_y, c1b, M_x,M_y, ccv_x,ccv_y);
	end
	if (B_x < J8_x)
		[J8_x,J8_y] = cint(c1_x,c1_y, c1b, J7_x,J7_y, ccv_x,ccv_y);
	end
	if (B_x < J6_x)
		[J6_x,J6_y] = cint(c1_x,c1_y, c1b, J5_x,J5_y, ccv_x,ccv_y);
	end
	if (B_x < J4_x)
		[J4_x,J4_y] = cint(c1_x,c1_y, c1b, J3_x,J3_y, ccv_x,ccv_y);
	end
	if (B_x < J2_x)
		[J2_x,J2_y] = cint(c1_x,c1_y, c1b, J1_x,J1_y, ccv_x,ccv_y);
	end
	if (T_x < Z42_x)
	[Z42_x,Z42_y] = pgc1(M_x,M_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	end
	[P92_x,P92_y] = pgc1(P91_x,P91_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[P44_x,P44_y] = pgc1(P43_x,P43_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[P46_x,P46_y] = pgc1(P45_x,P45_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[P48_x,P48_y] = pgc1(P47_x,P47_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[P50_x,P50_y] = pgc1(P49_x,P49_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[P52_x,P52_y] = pgc1(P51_x,P51_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[P54_x,P54_y] = pgc1(P53_x,P53_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[Z56_x,Z56_y] = pgc1(Z55_x,Z55_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[Z58_x,Z58_y] = pgc1(Z57_x,Z57_y, cmn_x,cmn_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[Z60_x,Z60_y] = pgc1(Z59_x,Z59_y, L60_x,L60_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	if (Z60_x <= T_x)
	[Z60_x,Z60_y] = pgc1(Z59_x,Z59_y, cv_x,cv_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	end
	[NN_x,NN_y] = pgc1(N_x,N_y,  L0_x,L0_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[N2_x,N2_y] = pgc1(N1_x,N1_y, L2_x,L2_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[N4_x,N4_y] = pgc1(N3_x,N3_y, LN4_x,LN4_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[N6_x,N6_y] = pgc1(N5_x,N5_y, LN6_x,LN6_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[N8_x,N8_y] = pgc1(N7_x,N7_y, L8_x,L8_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[N10_x,N10_y] = pgc1(N9_x,N9_y, L10_x,L10_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[N12_x,N12_y] = pgc1(N11_x,N11_y, L12_x,L12_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	[N14_x,N14_y] = pgc1(N13_x,N13_y, L14_x,L14_y, c1_x,c1_y, c1b, c2_x,c2_y, c2t,PF_x,PF_y,PS_x,PS_y,JAW_x,JAW_y,L6_x,L6_y,L7_x,L7_y,T_x);
	if (L7_x < N13_x)
		N14_x = N13_x;
		N14_y = L7_y;
	end
	P61_x = L7_x;
	P61_y = L3_y;
	P62_x = L7_x;
	P62_y = L7_y;
	if (T_x >= P61_x)
		[P62_x,P62_y] = cint(c1_x,c1_y, c1b, P61_x,P61_y, P62_x,P62_y);
	end
	if (P92_x == 0.000 | B_x > P92_x)
		[P92_x,P92_y] = cint(tongc_x,tongc_y, rt, P91_x,P91_y, cmn_x,cmn_y);
	end
	if (P44_x == 0.000 | B_x > P44_x)
		[P44_x,P44_y] = cint(tongc_x,tongc_y, rt, P43_x,P43_y, cmn_x,cmn_y);
	end
	if (P46_x == 0.000 | B_x > P46_x)
		[P46_x,P46_y] = cint(tongc_x,tongc_y, rt, P45_x,P45_y, cmn_x,cmn_y);
	end
	if (P48_x == 0.000 | B_x > P48_x | P48_x < P46_x)
		[P48_x,P48_y] = cint(tongc_x,tongc_y, rt, P47_x,P47_y, cmn_x,cmn_y);
	end
	if (P50_x == 0.000 | B_x > P50_x | P50_x < P48_x)
		[P50_x,P50_y] = cint(tongc_x,tongc_y, rt, P49_x,P49_y, cmn_x,cmn_y);
	end
	if (P52_x == 0.000 | B_x > P52_x | P52_x < P50_x)
		[P52_x,P52_y] = cint(tongc_x,tongc_y, rt, P51_x,P51_y, cmn_x,cmn_y);
	end
	if (P54_x == 0.000 | B_x > P54_x | P54_x < P52_x)
		[P54_x,P54_y] = cint(tongc_x,tongc_y, rt, P53_x,P53_y, cmn_x,cmn_y);
	end
	if (Z56_x == 0.000 | B_x > Z56_x | Z56_x < P54_x)
		[Z56_x,Z56_y] = cint(tongc_x,tongc_y, rt, Z55_x,Z55_y, cmn_x,cmn_y);
	end
	if (Z58_x == 0.000 | B_x > Z58_x | Z58_x < Z56_x)
		[Z58_x,Z58_y] = cint(tongc_x,tongc_y, rt, Z57_x,Z57_y, cmn_x,cmn_y);
	end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	figure(syn_win2_h);
	subplot(1,3,2);
	hold;
	two2seg(Z4_x,Z4_y, Z3_x,Z3_y, ox, oy, hei, xrat, yrat);
	two2seg(Z6_x,Z6_y, Z5_x,Z5_y, ox, oy, hei, xrat, yrat);
	two2seg(Z8_x,Z8_y, Z7_x,Z7_y, ox, oy, hei, xrat, yrat);
	two2seg(Z10_x,Z10_y, H3_x,H3_y, ox, oy, hei, xrat, yrat);
	two2seg(Z12_x,Z12_y, Z11_x,Z11_y, ox, oy, hei, xrat, yrat);
	two2seg(Z14_x,Z14_y, H2_x,H2_y, ox, oy, hei, xrat, yrat);
	two2seg(Z16_x,Z16_y, Z15_x,Z15_y, ox, oy, hei, xrat, yrat);
	two2seg(Z18_x,Z18_y, H1_x,H1_y, ox, oy, hei, xrat, yrat);
	two2seg(P81_x,P81_y, P82_x,P82_y, ox, oy, hei, xrat, yrat);
	two2seg(P11_x,P11_y, P12_x,P12_y, ox, oy, hei, xrat, yrat);
	two2seg(P83_x,P83_y, P84_x,P84_y, ox, oy, hei, xrat, yrat);
	two2seg(P13_x,P13_y, PP_x,PP_y, ox, oy, hei, xrat, yrat);
	two2seg(R1_x,R1_y, R2_x,R2_y, ox, oy, hei, xrat, yrat);
	two2seg(P85_x,P85_y, P86_x,P86_y, ox, oy, hei, xrat, yrat);
	two2seg(AS1_x,AS1_y, AS2_x,AS2_y, ox, oy, hei, xrat, yrat);
	two2seg(P15_x,P15_y, P16_x,P16_y, ox, oy, hei, xrat, yrat);
	two2seg(AS3_x,AS3_y, AS4_x,AS4_y, ox, oy, hei, xrat, yrat);
	two2seg(P87_x,P87_y, P88_x,P88_y, ox, oy, hei, xrat, yrat);
	two2seg(R3_x,R3_y, R4_x,R4_y, ox, oy, hei, xrat, yrat);
	two2seg(P17_x,P17_y, DL_x,DL_y, ox, oy, hei, xrat, yrat);
	two2seg(R5_x,R5_y, R6_x,R6_y, ox, oy, hei, xrat, yrat);
	two2seg(R7_x,R7_y, R8_x,R8_y, ox, oy, hei, xrat, yrat);
	two2seg(R9_x,R9_y, R10_x,R10_y, ox, oy, hei, xrat, yrat);
	two2seg(R11_x,R11_y, R12_x,R12_y, ox, oy, hei, xrat, yrat);
	two2seg(P19_x,P19_y, P20_x,P20_y, ox, oy, hei, xrat, yrat);
	two2seg(P21_x,P21_y, P22_x,P22_y, ox, oy, hei, xrat, yrat);
	two2seg(P23_x,P23_y, Z24_x,Z24_y, ox, oy, hei, xrat, yrat);
	two2seg(V_x,V_y, P26_x,P26_y, ox, oy, hei, xrat, yrat);
	two2seg(P27_x,P27_y, P28_x,P28_y, ox, oy, hei, xrat, yrat);
	two2seg(P29_x,P29_y, P30_x,P30_y, ox, oy, hei, xrat, yrat);
	two2seg(P31_x,P31_y, P32_x,P32_y, ox, oy, hei, xrat, yrat);
	two2seg(P33_x,P33_y, P34_x,P34_y, ox, oy, hei, xrat, yrat);
	two2seg(P35_x,P35_y, P36_x,P36_y, ox, oy, hei, xrat, yrat);
	two2seg(P37_x,P37_y, P38_x,P38_y, ox, oy, hei, xrat, yrat);
	two2seg(P39_x,P39_y, P40_x,P40_y, ox, oy, hei, xrat, yrat);
	two2seg(J1_x,J1_y, J2_x,J2_y, ox, oy, hei, xrat, yrat);
	two2seg(J3_x,J3_y, J4_x,J4_y, ox, oy, hei, xrat, yrat);
	two2seg(J5_x,J5_y, J6_x,J6_y, ox, oy, hei, xrat, yrat);
	two2seg(J7_x,J7_y, J8_x,J8_y, ox, oy, hei, xrat, yrat);
	two2seg(M_x,M_y, Z42_x,Z42_y, ox, oy, hei, xrat, yrat);
	two2seg(P91_x,P91_y, P92_x,P92_y, ox, oy, hei, xrat, yrat);
	two2seg(P43_x,P43_y, P44_x,P44_y, ox, oy, hei, xrat, yrat);
	two2seg(P45_x,P45_y, P46_x,P46_y, ox, oy, hei, xrat, yrat);
	two2seg(P47_x,P47_y, P48_x,P48_y, ox, oy, hei, xrat, yrat);
	two2seg(P49_x,P49_y, P50_x,P50_y, ox, oy, hei, xrat, yrat);
	two2seg(P51_x,P51_y, P52_x,P52_y, ox, oy, hei, xrat, yrat);
	two2seg(P53_x,P53_y, P54_x,P54_y, ox, oy, hei, xrat, yrat);
	two2seg(Z55_x,Z55_y, Z56_x,Z56_y, ox, oy, hei, xrat, yrat);
	two2seg(Z57_x,Z57_y, Z58_x,Z58_y, ox, oy, hei, xrat, yrat);
	two2seg(Z59_x,Z59_y, Z60_x,Z60_y, ox, oy, hei, xrat, yrat);
	two2seg(N_x,N_y, NN_x,NN_y, ox, oy, hei, xrat, yrat);
	two2seg(N1_x,N1_y, N2_x,N2_y, ox, oy, hei, xrat, yrat);
	two2seg(N3_x,N3_y, N4_x,N4_y, ox, oy, hei, xrat, yrat);
	two2seg(N5_x,N5_y, N6_x,N6_y, ox, oy, hei, xrat, yrat);
	two2seg(N7_x,N7_y, N8_x,N8_y, ox, oy, hei, xrat, yrat);
	two2seg(N9_x,N9_y, N10_x,N10_y, ox, oy, hei, xrat, yrat);
	two2seg(N11_x,N11_y, N12_x,N12_y, ox, oy, hei, xrat, yrat);
	two2seg(N13_x,N13_y, N14_x,N14_y, ox, oy, hei, xrat, yrat);
	two2seg(P61_x,P61_y, P62_x,P62_y, ox, oy, hei, xrat, yrat);
	hold;