animate.m

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

	X9 = distance(P86_x,P86_y, R2_x,R2_y)*2.0;
	D9 = ar1(R2_x,R2_y, P86_x,P86_y, AS2_x,AS2_y, D9,G1_y,W_y);
	AX9 = distance(AS2_x,AS2_y, P86_x,P86_y)*2.0;
	AD9 = ar1(P86_x,P86_y, AS2_x,AS2_y, P16_x,P16_y, AD9,G1_y,W_y);
	X10 = distance(P16_x,P16_y, AS2_x,AS2_y)*2.0;
	D10 = ar1(AS2_x,AS2_y, P16_x,P16_y, AS4_x,AS4_y, D10,G1_y,W_y);
	AX10= distance(P88_x,P88_y, AS4_x,AS4_y)*2.0;
	AD10= ar1(P16_x,P16_y, AS4_x,AS4_y, P88_x,P88_y, AD10,G1_y,W_y);
	X11 = distance(P88_x,P88_y, AS4_x,AS4_y)*2.0;
	D11 = ar1(AS4_x,AS4_y, P88_x,P88_y, R4_x,R4_y, D11,G1_y,W_y);
	XR2 = distance(R4_x,R4_y, P88_x,P88_y)*2.0;
	RD2 = ar1(P88_x,P88_y, R4_x,R4_y, DLM_x,DLM_y, RD2,G1_y,W_y);
	X12 = distance(DLM_x,DLM_y, R4_x,R4_y)*2.0;
	D12 = ar1(R4_x,R4_y, DLM_x,DLM_y, R6_x,R6_y, D12,G1_y,W_y);
% AR9 */
	XR3 = distance(R6_x,R6_y, DLM_x,DLM_y)*2.0;
	RD3 = ar9(DLM_x,DLM_y, R6_x,R6_y, R8_x,R8_y, RD3, R5_x,R5_y,G1_y);
	XR4 = distance(R8_x,R8_y, R6_x,R6_y)*2.0;
	RD4 = ar9(R6_x,R6_y, R8_x,R8_y, R10_x,R10_y, RD4, R7_x,R7_y,G1_y);
	XR5 = distance(R10_x,R10_y, R8_x,R8_y)*2.0;
	RD5 = ar9(R8_x,R8_y, R10_x,R10_y, R12_x,R12_y, RD5, R9_x,R9_y,G1_y);
	XR6 = distance(R12_x,R12_y, R10_x,R10_y)*2.0;
	RD6 = ar9(R10_x,R10_y, R12_x,R12_y, P20_x,P20_y, RD6, R11_x,R11_y,G1_y);
	X13 = distance(P20_x,P20_y, R12_x,R12_y)*2.0;
	D13 = ar9(R12_x,R12_y, P20_x,P20_y, P22_x,P22_y, D13, P19_x,P19_y,G1_y);
	X14 = distance(P22_x,P22_y, P20_x,P20_y)*2.0;
	D14 = ar9(P20_x,P20_y, P22_x,P22_y, P24_x,P24_y, D14, P21_x,P21_y,G1_y);
	X15 = distance(P24_x,P24_y, P22_x,P22_y)*2.0;
	D15 = ar9(P22_x,P22_y, P24_x,P24_y, P26_x,P26_y, D15, P23_x,P23_y,G1_y);
% AR2 */
	X16 = distance(P26_x,P26_y, P24_x,P24_y)*2.0;
	D16 = ar2(P24_x,P24_y, P26_x,P26_y, P28_x,P28_y, D16, V_x,V_y);
	X17 = distance(P28_x,P28_y, P26_x,P26_y)*2.0;
	D17 = ar2(P26_x,P26_y, P28_x,P28_y, P30_x,P30_y, D17, P27_x,P27_y);
	X18 = distance(P30_x,P30_y, P28_x,P28_y)*2.0;
	D18 = ar2(P28_x,P28_y, P30_x,P30_y, P32_x,P32_y, D18, P29_x,P29_y);
	X19 = distance(P32_x,P32_y, P30_x,P30_y)*2.0;
	D19 = ar2(P30_x,P30_y, P32_x,P32_y, P34_x,P34_y, D19, P31_x,P31_y);
	X20 = distance(P34_x,P34_y, P32_x,P32_y)*2.0;
	D20 = ar2(P32_x,P32_y, P34_x,P34_y, P36_x,P36_y, D20, P33_x,P33_y);
	X21 = distance(P36_x,P36_y, P34_x,P34_y)*2.0;
	D21 = ar2(P34_x,P34_y, P36_x,P36_y, P38_x,P38_y, D21, P35_x,P35_y);
	X22 = distance(P38_x,P38_y, P36_x,P36_y)*2.0;
	D22 = ar2(P36_x,P36_y, P38_x,P38_y, P40_x,P40_y, D22, P37_x,P37_y);
	X23 = distance(P40_x,P40_y, P38_x,P38_y)*2.0;
	D23 = ar2(P38_x,P38_y, P40_x,P40_y, J2_x,J2_y, D23, P39_x,P39_y);
	XJ1 = distance(J2_x,J2_y, P40_x,P40_y)*2.0;
	JD1 = ar2(P40_x,P40_y, J2_x,J2_y, J4_x,J4_y, JD1, J1_x,J1_y);
	XJ2 = distance(J4_x,J4_y, J2_x,J2_y)*2.0;
	JD2 = ar2(J2_x,J2_y, J4_x,J4_y, J6_x,J6_y, JD2, J3_x,J3_y);
	XJ3 = distance(J6_x,J6_y, J4_x,J4_y)*2.0;
	JD3 = ar2(J4_x,J4_y, J6_x,J6_y, J8_x,J8_y, JD3, J5_x,J5_y);
	XJ4 = distance(J8_x,J8_y, J6_x,J6_y)*2.0;
	JD4 = ar2(J6_x,J6_y, J8_x,J8_y, P42_x,P42_y, JD4, J7_x,J7_y);
	X24 = distance(P42_x,P42_y, J8_x,J8_y)*2.0;
	D24 = ar2(J8_x,J8_y, P42_x,P42_y, P92_x,P92_y, D24, M_x,M_y);
% AR23 */
	X92 = distance(P92_x,P92_y, P42_x,P42_y)*2.0;
	D92 = ar23(P42_x,P42_y, P92_x,P92_y, P44_x,P44_y, D92, P91_x,P91_y);
	X25 = distance(P44_x,P44_y, P92_x,P92_y)*2.0;
	D25 = ar23(P92_x,P92_y, P44_x,P44_y, P46_x,P46_y, D25, P43_x,P43_y);
	X26 = distance(P46_x,P46_y, P44_x,P44_y)*2.0;
	D26 = ar23(P44_x,P44_y, P46_x,P46_y, P48_x,P48_y, D26, P45_x,P45_y);
	X27 = distance(P48_x,P48_y, P46_x,P46_y)*2.0;
	D27 = ar23(P46_x,P46_y, P48_x,P48_y, P50_x,P50_y, D27, P47_x,P47_y);
	X28 = distance(P50_x,P50_y, P48_x,P48_y)*2.0;
	D28 = ar23(P48_x,P48_y, P50_x,P50_y, P52_x,P52_y, D28, P49_x,P49_y);
	X29 = distance(P52_x,P52_y, P50_x,P50_y)*2.0;
	D29 = ar23(P50_x,P50_y, P52_x,P52_y, P54_x,P54_y, D29, P51_x,P51_y);
	X30 = distance(P54_x,P54_y, P52_x,P52_y)*2.0;
	D30 = ar23(P52_x,P52_y, P54_x,P54_y, Z56_x,Z56_y, D30, P53_x,P53_y);
	X31 = distance(Z56_x,Z56_y, P54_x,P54_y)*2.0;
	D31 = ar23(P54_x,P54_y, Z56_x,Z56_y, Z58_x,Z58_y, D31, Z55_x,Z55_y);
	X32 = distance(Z58_x,Z58_y, Z56_x,Z56_y)*2.0;
	D32 = ar23(Z56_x,Z56_y, Z58_x,Z58_y, Z60_x,Z60_y, D32, Z57_x,Z57_y);
	X33 = distance(Z60_x,Z60_y, Z58_x,Z58_y)*2.0;
	D33 = ar23(Z58_x,Z58_y, Z60_x,Z60_y, N0_x,N0_y, D33, Z59_x,Z59_y);
% AR 4 */
	X34 = distance(N0_x,N0_y, Z60_x,Z60_y)*2.0;
	D34 = ar4(Z60_x,Z60_y, N0_x,N0_y, N2_x,N2_y, D34, N_x,N_y);
	X35 = distance(N2_x,N2_y, N0_x,N0_y)*2.0;
	D35 = ar4(N0_x,N0_y, N2_x,N2_y, N4_x,N4_y, D35, N1_x,N1_y);
	X36 = distance(N4_x,N4_y, N2_x,N2_y)*2.0;
	D36 = ar4(N2_x,N2_y, N4_x,N4_y, N6_x,N6_y, D36, N3_x,N3_y);
	X37 = distance(N6_x,N6_y, N4_x,N4_y)*2.0;
	D37 = ar4(N4_x,N4_y, N6_x,N6_y, N8_x,N8_y, D37, N5_x,N5_y);
	X38 = distance(N8_x,N8_y, N6_x,N6_y)*2.0;
	D38 = ar4(N6_x,N6_y, N8_x,N8_y, N10_x,N10_y, D38, N7_x,N7_y);
	X39 = distance(N10_x,N10_y, N8_x,N8_y)*2.0;
	D39 = ar4(N8_x,N8_y, N10_x,N10_y, N12_x,N12_y, D39, N9_x,N9_y);
	X40 = distance(N12_x,N12_y, N10_x,N10_y)*2.0;
	D40 = ar4(N10_x,N10_y, N12_x,N12_y, N14_x,N14_y, D40, N11_x,N11_y);
	X41 = distance(N14_x,N14_y, N12_x,N12_y)*2.0;
	D41 = ar4(N12_x,N12_y, N14_x,N14_y, P62_x,P62_y, D41, N13_x,N13_y);
% AR 5 */
	X42 = distance(P62_x,P62_y, N14_x,N14_y)*2.0;
	D42 = ar5(N14_x,N14_y, P62_x,P62_y, P64_x,P64_y, D42, L5_x,L5_y, L6_x,L6_y, L7_x,L7_y);
% area */
	area=zeros(1,60);
	leng=zeros(1,60);
	area(1)  = DD1;
	area(2)  = DD2;
	area(3)  = DD3;
	area(4)  = DD4;
	area(5)  = DD5;
	area(6)  = DD6;
	area(7)  = DD7;
	area(8)  = DD8;
	area(9)  = ND9;
	area(10)  = D5;
	area(11) = D6;
	area(12) = D7;
	area(13) = D8;
	area(14) = RD1;
	area(15) = D9;
	area(16) = AD9;
	area(17) = D10;
	area(18) = AD10;
	area(19) = D11;
	area(20) = RD2;
	area(21) = D12;
	area(22) = RD3;
	area(23) = RD4;
	area(24) = RD5;
	area(25) = RD6;
	area(26) = D13;
	area(27) = D14;
	area(28) = D15;
	area(29) = D16;
	area(30) = D17;
	area(31) = D18;
	area(32) = D19;
	area(33) = D20;
	area(34) = D21;
	area(35) = D22;
	area(36) = D23;
	area(37) = JD1;
	area(38) = JD2;
	area(39) = JD3;
	area(40) = JD4;
	area(41) = D24;
	area(42) = D92;
	area(43) = D25;
	area(44) = D26;
	area(45) = D27;
	area(46) = D28;
	area(47) = D29;
	area(48) = D30;
	area(49) = D31;
	area(50) = D32;
	area(51) = D33;
	area(52) = D34;
	area(53) = D35;
	area(54) = D36;
	area(55) = D37;
	area(56) = D38;
	area(57) = D39;
	area(58) = D40;
	area(59) = D41;
	area(60) = D42;
	NT1A = 2.0*distance(Vs_x,Vs_y, V_x,V_y);        % refer to Ding diss. P.82 */
% length */
	leng(1)  = XX1;
	leng(2)  = XX2;
	leng(3)  = XX3;
	leng(4)  = XX4;
	leng(5)  = XX5;
	leng(6)  = XX6;
	leng(7)  = XX7;
	leng(8)  = XX8;
	leng(9)  = XX9;
	leng(10)  = X5;
	leng(11) = X6;
	leng(12) = X7;
	leng(13) = X8;
	leng(14) = XR1;
	leng(15) = X9;
	leng(16) = AX9;
	leng(17) = X10;
	leng(18) = AX10;
	leng(19) = X11;
	leng(20) = XR2;
	leng(21) = X12;
	leng(22) = XR3;
	leng(23) = XR4;
	leng(24) = XR5;
	leng(25) = XR6;
	leng(26) = X13;
	leng(27) = X14;
	leng(28) = X15;
	leng(29) = X16;
	leng(30) = X17;
	leng(31) = X18;
	leng(32) = X19;
	leng(33) = X20;
	leng(34) = X21;
	leng(35) = X22;
	leng(36) = X23;
	leng(37) = XJ1;
	leng(38) = XJ2;
	leng(39) = XJ3;
	leng(40) = XJ4;
	leng(41) = X24;
	leng(42) = X92;
	leng(43) = X25;
	leng(44) = X26;
	leng(45) = X27;
	leng(46) = X28;
	leng(47) = X29;
	leng(48) = X30;
	leng(49) = X31;
	leng(50) = X32;
	leng(51) = X33;
	leng(52) = X34;
	leng(53) = X35;
	leng(54) = X36;
	leng(55) = X37;
	leng(56) = X38;
	leng(57) = X39;
	leng(58) = X40;
	leng(59) = X41;
	leng(60) = X42;

	p1_x=5;
	p1_y=6;
	p2_x=3;
	p2_y=9;
	center_x=0;
	center_y=0;
	ox=4;
	oy=3;
	wid=3;
	hei=4;
	xrat=0.5;
	yrat=0.5;

% First draw fix structure
 	figure(animate_h);
	clf;
	Wh_x = W_x;
	Wh_y = H_y;
	draw_lin(Wh_x, Wh_y, G2_x, G2_y, ox, oy, wid, hei, xrat, yrat);
	hold;
	draw_lin(G2_x, G2_y, G_x, G_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(G_x, G_y, G1_x, G1_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(H1_x,H1_y, H2_x, H2_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(H2_x, H2_y, K_x, K_y, ox, oy, wid, hei, xrat, yrat);

% Draw lips

	draw_lin(N_x,N_y, U_x,U_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(JAW_x,JAW_y, PS_x,PS_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(PS_x,PS_y, PF_x,PF_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(U_x,U_y, L3_x,L3_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(L3_x,L3_y, L5_x,L5_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(JAW_x,JAW_y, L6_x,L6_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(L6_x,L6_y, L7_x,L7_y, ox, oy, wid, hei, xrat, yrat);

% Draw hyoid 

	if (hy < 0.0)
		draw_arc(PP_x,PP_y, H_x,H_y, pme_x,pme_y, ox, oy, wid, hei, xrat, yrat);
	else
		draw_arc(H_x,H_y, PP_x,PP_y, pme_x,pme_y, ox, oy, wid, hei, xrat, yrat);
	end
	draw_lin(DL_x,DL_y, PP_x,PP_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(H_x,H_y, H1_x,H1_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(W_x,W_y, Wh_x,Wh_y, ox, oy, wid, hei, xrat, yrat);

%  Draw tonge body & tonge blade

	draw_arc(N_x,N_y, M_x,M_y, cmn_x,cmn_y, ox, oy, wid, hei, xrat, yrat);
	draw_arc(T_x,T_y, B_x,B_y, c1_x,c1_y, ox, oy, wid, hei, xrat, yrat);
	draw_arc(PF_x,PF_y, T_x,T_y, c2_x,c2_y, ox, oy, wid, hei, xrat, yrat);
% Draw velum
	
	draw_arc(M_x,M_y, V_x,V_y, cv_x,cv_y, ox, oy, wid, hei, xrat, yrat);
	draw_lin(V_x,V_y, W_x,W_y, ox, oy, wid, hei, xrat, yrat);
	draw_arc(B_x,B_y, DL_x,DL_y, tongc_x,tongc_y, ox, oy, wid, hei, xrat, yrat);
%	title('Midsagittal profile','color',TEXT_COLOR);
	
	ss=sprintf('Midsagittal profile');
	text(5,4.65,ss,'color',TEXT_COLOR,'FontSize',10);
	ss=sprintf('frame #: %4.0f%', currf);
	text(5.2,4.4,ss,'color',TEXT_COLOR,'FontSize',10);

	axis([floor(ox+xrat*G2_x-1.0) ceil(ox+xrat*L3_x+1.0) 1 4.5]);
%       axis([4.5 7.5 1 4.5]);

	axis('off');

	drawnow;