?*X2+a2*Y2+a3*Z2+a4-x2*X2*a9-x2*Y2*a10-x2*Z2*a11=x2,a5*X2+a6*Y2+a7*Z2+a8-y2*X2*a9-y2*Y2*a10-y2*Z2*a11=y2;
?*X3+a2*Y3+a3*Z3+a4-x3*X3*a9-x3*Y3*a10-x3*Z3*a11=x3,a5*X3+a6*Y3+a7*Z3+a8-y3*X3*a9-y3*Y3*a10-y3*Z3*a11=y3;
?*X4+a2*Y4+a3*Z4+a4-x4*X4*a9-x4*Y4*a10-x4*Z4*a11=x4,a5*X4+a6*Y4+a7*Z4+a8-y4*X4*a9-y4*Y4*a10-y4*Z4*a11=y4;
?*X5+a2*Y5+a3*Z5+a4-x5*X5*a9-x5*Y5*a10-x5*Z5*a11=x5,a5*X5+a6*Y5+a7*Z5+a8-y5*X5*a9-y5*Y5*a10-y5*Z5*a11=y5; %把已知量代入方程中进行求解
[a1,a2,a4,a5,a6,a8,a9,a10]=solve('a1*(-50)+a2*50+a4-51.6*50*a9+51.6*50*a10=-51.6','a5*(-50)+a6*50+a8-51.1*50*a9+51.5*50*a10=-51.1','a1*(-20)+a2*50+a4-20*24.9*a9+24.9*50*a10=-24.9','a5*(-20)+a6*50+a8-50.8*20*a9+50*50.8*a10=-50.8','a1*50+a2*50+a4-33.3*50*a9-33.3*50*a10=33.3','a5*50+a6*50+a8+42.9*50*a9-42.9*50*a10=-42.9','a1*50+a2*(-50)+a4-21.4*50*a9+21.4*50*a10=21.4','a5*50+a6*(-50)+a8-31.2*50*a9+31.2*50*a10=31.2','a1,a2,a4,a5,a6,a8,a9,a10')
21