matlab怎么算多元多项式,如何提取多元多项式的系数？

Q=-1+1/b^2*y^2+1/a^2*z0^2-2/a^2*z*z0-1/c^2*x0^2*cos(theta_z)^2-1/c^2*x0^2*cos(theta_x)^2-1/b^2*y0^2*cos(theta_y)^2-1/c^2*x^2*cos(theta_x)^2-1/c^2*x^2*cos(theta_z)^2-1/a^2*z0^2*cos(theta_y)^2-1/b^2*y^2*cos(theta_y)^2-1/b^2*y^2*cos(theta_x)^2-1/b^2*y0^2*cos(theta_x)^2+1/b^2*cos(theta_x)^2*x0^2+1/c^2*cos(theta_x)^2*y0^2+1/b^2*cos(theta_x)^2*x^2+1/c^2*cos(theta_x)^2*y^2-1/a^2*z^2*cos(theta_y)^2+1/c^2*y^2*cos(theta_z)^2+1/c^2*y0^2*cos(theta_z)^2+1/b^2*x0^2*cos(theta_z)^2+1/b^2*cos(theta_y)^2*z0^2+1/b^2*x^2*cos(theta_z)^2+1/a^2*cos(theta_y)^2*y^2+1/b^2*cos(theta_y)^2*z^2+1/a^2*cos(theta_y)^2*y0^2-1/b^2*y0^2*cos(theta_z)^2-1/b^2*y^2*cos(theta_z)^2-2/c^2*x*x0-2/b^2*y*y0+4/b^2*cos(theta_x)^2*x0*y*cos(theta_z)*sin(theta_z)-4/b^2*cos(theta_x)^2*y0*x0*cos(theta_z)*sin(theta_z)+4/b^2*cos(theta_x)^2*y0*x*cos(theta_z)*sin(theta_z)-4/b^2*cos(theta_x)^2*y*x*cos(theta_z)*sin(theta_z)+2/b^2*cos(theta_x)^2*y0^2*cos(theta_z)^2+2/b^2*cos(theta_x)^2*y^2*cos(theta_z)^2+1/c^2*cos(theta_y)^2*cos(theta_x)^2*z0^2+1/c^2*cos(theta_y)^2*cos(theta_x)^2*z^2-2/c^2*cos(theta_y)^2*cos(theta_x)^2*z0*z+2/b^2*y0*y*cos(theta_z)^2-2/b^2*x0*y*cos(theta_z)*sin(theta_z)+1/a^2*cos(theta_y)^2*x^2*cos(theta_z)^2+1/c^2*x^2+2/c^2*x0^2*cos(theta_x)^2*cos(theta_z)^2-2/c^2*y*x*cos(theta_z)*sin(theta_z)+4/c^2*y*x*cos(theta_z)*sin(theta_z)*cos(theta_x)^2-2/c^2*y*x*cos(theta_z)*sin(theta_z)*cos(theta_y)^2*cos(theta_x)^2+1/a^2*cos(theta_y)^2*x0^2*cos(theta_z)^2+1/b^2*y0^2*cos(theta_z)^2*cos(theta_y)^2+1/b^2*y0^2*cos(theta_y)^2*cos(theta_x)^2-1/b^2*y0^2*cos(theta_z)^2*cos(theta_y)^2*cos(theta_x)^2+2/c^2*x^2*cos(theta_x)^2*cos(theta_z)^2-2/a^2*cos(theta_y)*y0*z*sin(theta_y)*sin(theta_z)-2/a^2*cos(theta_y)*z0*y*sin(theta_y)*sin(theta_z)-2/b^2*cos(theta_y)*cos(theta_x)*y0*sin(theta_x)*z*cos(theta_z)+2/a^2*z*z0*cos(theta_y)^2-4/b^2*cos(theta_x)*sin(theta_x)*x0*x*cos(theta_z)*sin(theta_y)*sin(theta_z)-1/c^2*cos(theta_y)^2*x^2*cos(theta_x)^2*cos(theta_z)^2+2/b^2*cos(theta_x)*sin(theta_x)*y0*x0*sin(theta_y)-1/b^2*cos(theta_y)^2*cos(theta_x)^2*y^2*cos(theta_z)^2+2/a^2*cos(theta_y)*z*y*sin(theta_y)*sin(theta_z)+4/c^2*cos(theta_x)*sin(theta_x)*x0*x*cos(theta_z)*sin(theta_y)*sin(theta_z)-2/b^2*cos(theta_y)^2*cos(theta_x)^2*x0*x*cos(theta_z)^2+4/b^2*cos(theta_x)*sin(theta_x)*y0*y*cos(theta_z)*sin(theta_y)*sin(theta_z)-2/b^2*cos(theta_x)^2*x0*y*cos(theta_z)*sin(theta_z)*cos(theta_y)^2+2/b^2*cos(theta_x)^2*y0*y*cos(theta_z)^2*cos(theta_y)^2-4/c^2*cos(theta_x)*sin(theta_x)*y0*y*cos(theta_z)*sin(theta_y)*sin(theta_z)+1/b^2*cos(theta_y)^2*cos(theta_x)^2*x^2*cos(theta_z)^2+2/b^2*cos(theta_y)*cos(theta_x)*z*sin(theta_x)*y*cos(theta_z)-2/a^2*cos(theta_y)*z0*x*cos(theta_z)*sin(theta_y)+2/a^2*cos(theta_y)*z0*x0*cos(theta_z)*sin(theta_y)+2/c^2*cos(theta_y)*cos(theta_x)*sin(theta_x)*y0*z*cos(theta_z)-2/c^2*cos(theta_y)*cos(theta_x)^2*z*x*cos(theta_z)*sin(theta_y)-4/c^2*cos(theta_x)^2*x0*y*cos(theta_z)*sin(theta_z)-2/c^2*cos(theta_y)*cos(theta_x)^2*z0*x0*cos(theta_z)*sin(theta_y)-2/c^2*cos(theta_x)^2*y^2*cos(theta_z)^2+2/c^2*x0*y*cos(theta_z)*sin(theta_z)-2/b^2*cos(theta_y)^2*z0*z+2/b^2*cos(theta_y)^2*z0*z*cos(theta_x)^2+2/c^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*y*cos(theta_z)+2/c^2*x*x0*cos(theta_z)^2-2/c^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*y0*cos(theta_z)-2/b^2*cos(theta_y)*cos(theta_x)*z*sin(theta_x)*x*sin(theta_z)-2/b^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*y*cos(theta_z)+2/c^2*cos(theta_y)*cos(theta_x)^2*x0*z*cos(theta_z)*sin(theta_y)+1/b^2*cos(theta_y)^2*cos(theta_x)^2*x0^2*cos(theta_z)^2+2/c^2*x*x0*cos(theta_x)^2+2/c^2*cos(theta_x)*sin(theta_x)*y0*x*sin(theta_y)-2/b^2*cos(theta_x)*sin(theta_x)*x0*y*sin(theta_y)+2/b^2*y*y0*cos(theta_y)^2-2/b^2*cos(theta_y)*z0*y0*sin(theta_y)*sin(theta_z)+2/b^2*cos(theta_y)*z0*y0*sin(theta_y)*sin(theta_z)*cos(theta_x)^2+2/c^2*cos(theta_y)*cos(theta_x)^2*z0*x*cos(theta_z)*sin(theta_y)-2/b^2*cos(theta_y)*z*x*cos(theta_z)*sin(theta_y)-2/b^2*y*y0*cos(theta_z)^2*cos(theta_y)^2+2/b^2*y*y0*cos(theta_x)^2-2/b^2*y*y0*cos(theta_y)^2*cos(theta_x)^2-2/b^2*x0*x*cos(theta_z)^2+2/b^2*x0*x*cos(theta_z)^2*cos(theta_y)^2+2/b^2*cos(theta_y)*z*x*cos(theta_z)*sin(theta_y)*cos(theta_x)^2-2/a^2*cos(theta_y)*x0*z*cos(theta_z)*sin(theta_y)+2/b^2*cos(theta_y)*z0*x*cos(theta_z)*sin(theta_y)-2/b^2*cos(theta_y)*z0*x*cos(theta_z)*sin(theta_y)*cos(theta_x)^2-2/c^2*y0*x0*cos(theta_z)*sin(theta_z)+4/c^2*cos(theta_x)*sin(theta_x)*y0*x0*cos(theta_z)^2*sin(theta_y)-2/c^2*y0*y*cos(theta_z)^2+2/b^2*x0*y*cos(theta_z)*sin(theta_z)*cos(theta_y)^2+2/b^2*cos(theta_y)*cos(theta_x)*x0*sin(theta_x)*z*sin(theta_z)-2/c^2*cos(theta_y)*cos(theta_x)*z*sin(theta_x)*y*cos(theta_z)-2/c^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*x*sin(theta_z)-2/b^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*x0*sin(theta_z)+2/b^2*cos(theta_x)*sin(theta_x)*y*x*sin(theta_y)-2/b^2*cos(theta_x)*sin(theta_x)*y0*x*sin(theta_y)+2/a^2*cos(theta_y)*z0*y0*sin(theta_y)*sin(theta_z)-2/c^2*cos(theta_y)*cos(theta_x)*sin(theta_x)*x0*z*sin(theta_z)+1/c^2*x0^2+2/c^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*x0*sin(theta_z)-2/c^2*cos(theta_x)*sin(theta_x)*y0*x0*sin(theta_y)+2/c^2*cos(theta_y)*cos(theta_x)*z*sin(theta_x)*x*sin(theta_z)+1/b^2*y^2*cos(theta_z)^2*cos(theta_y)^2-2/a^2*cos(theta_y)^2*x0*x*cos(theta_z)^2-1/c^2*cos(theta_x)^2*y^2*cos(theta_y)^2+2/a^2*cos(theta_y)^2*y*x*cos(theta_z)*sin(theta_z)-4/c^2*cos(theta_x)*sin(theta_x)*y0*x*cos(theta_z)^2*sin(theta_y)+1/b^2*y^2*cos(theta_y)^2*cos(theta_x)^2-2/b^2*cos(theta_x)^2*x^2*cos(theta_z)^2-1/c^2*cos(theta_x)^2*y0^2*cos(theta_y)^2-2/c^2*cos(theta_x)^2*y0^2*cos(theta_z)^2-2/b^2*cos(theta_x)^2*x0^2*cos(theta_z)^2-4/c^2*cos(theta_x)*sin(theta_x)*x0*y*cos(theta_z)^2*sin(theta_y)-2/b^2*cos(theta_x)*sin(theta_x)*y^2*cos(theta_z)*sin(theta_y)*sin(theta_z)+2/c^2*cos(theta_x)*sin(theta_x)*y^2*cos(theta_z)*sin(theta_y)*sin(theta_z)-2/c^2*cos(theta_x)*sin(theta_x)*x^2*cos(theta_z)*sin(theta_y)*sin(theta_z)+4/b^2*cos(theta_x)*sin(theta_x)*x0*y*cos(theta_z)^2*sin(theta_y)-4/b^2*cos(theta_x)*sin(theta_x)*y*x*cos(theta_z)^2*sin(theta_y)+4/c^2*cos(theta_x)*sin(theta_x)*y*x*cos(theta_z)^2*sin(theta_y)-4/b^2*cos(theta_x)*sin(theta_x)*y0*x0*cos(theta_z)^2*sin(theta_y)-2/b^2*cos(theta_x)*sin(theta_x)*y0^2*cos(theta_z)*sin(theta_y)*sin(theta_z)+2/b^2*cos(theta_x)*sin(theta_x)*x0^2*cos(theta_z)*sin(theta_y)*sin(theta_z)+2/c^2*cos(theta_x)*sin(theta_x)*y0^2*cos(theta_z)*sin(theta_y)*sin(theta_z)+2/b^2*cos(theta_x)*sin(theta_x)*x^2*cos(theta_z)*sin(theta_y)*sin(theta_z)+2/a^2*cos(theta_y)^2*y0*x0*cos(theta_z)*sin(theta_z)-2/a^2*cos(theta_y)^2*x0*y*cos(theta_z)*sin(theta_z)+2/a^2*cos(theta_y)^2*y0*y*cos(theta_z)^2+2/c^2*cos(theta_y)^2*x0*y*cos(theta_z)*sin(theta_z)*cos(theta_x)^2-2/b^2*cos(theta_y)*z*y*sin(theta_y)*sin(theta_z)+2/b^2*cos(theta_y)*z*y*sin(theta_y)*sin(theta_z)*cos(theta_x)^2+4/c^2*cos(theta_x)^2*y0*x0*cos(theta_z)*sin(theta_z)-2/c^2*cos(theta_x)^2*y0*x0*cos(theta_z)*sin(theta_z)*cos(theta_y)^2+2/b^2*cos(theta_x)^2*y0*x0*cos(theta_z)*sin(theta_z)*cos(theta_y)^2-2/b^2*y0*x0*cos(theta_z)*sin(theta_z)*cos(theta_y)^2+2/b^2*y0*x*cos(theta_z)*sin(theta_z)*cos(theta_y)^2-2/b^2*y0*x*cos(theta_z)*sin(theta_z)*cos(theta_y)^2*cos(theta_x)^2+2/c^2*cos(theta_x)*sin(theta_x)*x0*y*sin(theta_y)+2/c^2*y0*x*cos(theta_z)*sin(theta_z)-4/c^2*y0*x*cos(theta_z)*sin(theta_z)*cos(theta_x)^2+2/c^2*y0*x*cos(theta_z)*sin(theta_z)*cos(theta_y)^2*cos(theta_x)^2+2/b^2*y*x*cos(theta_z)*sin(theta_z)+2/c^2*cos(theta_y)*cos(theta_x)^2*z0*y*sin(theta_y)*sin(theta_z)-2/c^2*cos(theta_y)*cos(theta_x)^2*z*y*sin(theta_y)*sin(theta_z)+2/c^2*cos(theta_y)*cos(theta_x)^2*y0*z*sin(theta_y)*sin(theta_z)+2/a^2*cos(theta_y)*z*x*cos(theta_z)*sin(theta_y)-2/a^2*cos(theta_y)^2*y0*x*cos(theta_z)*sin(theta_z)-2/c^2*cos(theta_x)*sin(theta_x)*x0^2*cos(theta_z)*sin(theta_y)*sin(theta_z)-4/b^2*cos(theta_x)^2*y0*y*cos(theta_z)^2-2/b^2*y0*x*cos(theta_z)*sin(theta_z)+2/b^2*y0*x0*cos(theta_z)*sin(theta_z)-2/c^2*cos(theta_y)*cos(theta_x)^2*z0*y0*sin(theta_y)*sin(theta_z)-1/b^2*x0^2*cos(theta_z)^2*cos(theta_y)^2-1/b^2*x^2*cos(theta_z)^2*cos(theta_y)^2-2/b^2*cos(theta_y)*x0*z*cos(theta_z)*sin(theta_y)*cos(theta_x)^2+4/b^2*cos(theta_x)*sin(theta_x)*y0*x*cos(theta_z)^2*sin(theta_y)-4/c^2*cos(theta_x)^2*x0*x*cos(theta_z)^2+2/c^2*cos(theta_x)^2*x0*x*cos(theta_z)^2*cos(theta_y)^2+1/b^2*y0^2-1/b^2*cos(theta_y)^2*z0^2*cos(theta_x)^2-1/b^2*cos(theta_y)^2*z^2*cos(theta_x)^2-1/a^2*cos(theta_y)^2*y^2*cos(theta_z)^2-1/a^2*cos(theta_y)^2*y0^2*cos(theta_z)^2+2/b^2*cos(theta_y)*x0*z*cos(theta_z)*sin(theta_y)+1/a^2*z^2-2/c^2*cos(theta_x)^2*y0*y+2/c^2*cos(theta_x)^2*y0*y*cos(theta_y)^2+4/c^2*cos(theta_x)^2*y0*y*cos(theta_z)^2-2/c^2*cos(theta_x)^2*y0*y*cos(theta_z)^2*cos(theta_y)^2-2/b^2*cos(theta_x)^2*x0*x+4/b^2*cos(theta_x)^2*x0*x*cos(theta_z)^2+1/c^2*cos(theta_y)^2*y^2*cos(theta_z)^2*cos(theta_x)^2+1/c^2*cos(theta_y)^2*y0^2*cos(theta_z)^2*cos(theta_x)^2-1/c^2*cos(theta_y)^2*x0^2*cos(theta_x)^2*cos(theta_z)^2-2/a^2*cos(theta_y)^2*y0*y+2/b^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*x*sin(theta_z)-2/b^2*y*x*cos(theta_z)*sin(theta_z)*cos(theta_y)^2+2/b^2*y*x*cos(theta_z)*sin(theta_z)*cos(theta_y)^2*cos(theta_x)^2+2/b^2*cos(theta_y)*z0*y*sin(theta_y)*sin(theta_z)-2/b^2*cos(theta_y)*y0*z*sin(theta_y)*sin(theta_z)*cos(theta_x)^2-2/c^2*cos(theta_x)*sin(theta_x)*y*x*sin(theta_y)+2/b^2*cos(theta_y)*cos(theta_x)*z0*sin(theta_x)*y0*cos(theta_z)-2/b^2*cos(theta_y)*z0*y*sin(theta_y)*sin(theta_z)*cos(theta_x)^2-2/b^2*cos(theta_y)*z0*x0*cos(theta_z)*sin(theta_y)+2/b^2*cos(theta_y)*z0*x0*cos(theta_z)*sin(theta_y)*cos(theta_x)^2+2/b^2*cos(theta_y)*y0*z*sin(theta_y)*sin(theta_z)，

如何提取x^2,y^2,z^2,x*y,y*z,z*x,x,y,z的系数？a,b,c,a0,b0,c0,theta_x,theta_y,theta_z都是常量，之后会赋值。

顺便再问个小问题，定义了多项式x^2*A+2*x*y*D+2*x*z*F+2*x*G+y^2*B+2*y*z*E+2*y*H+z^2*C+2*z*I+J，各个系数也是稍后会赋值，但是怎么使得赋值后系数就直接显示数值，而不是原来的A、B、C……呢？