大侠看下什么错误
程序如下：<BR>syms theta Ls Lt Ms;<BR>D=;<BR>L=;<BR>H=inv(D)<BR>A=D*L*H<BR>其中D是派克变换矩阵，L是一个电机的电感矩阵，H是派克变换矩阵的逆矩阵，这个程序只要是讲L矩阵进行派克变换，其中电机的电感矩阵派克变换后应该是与theta无关的一个常数矩阵才对，但是这个程序的计算结果却不是常数矩阵，计算结果如下：<BR>H =<BR> <BR>[ 3/2*(sin(theta+1/3*pi)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3/2*(cos(theta+1/3*pi)-sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3*(cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))]<BR>[ -3/2*(-sin(theta)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3/2*(cos(theta)+sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3*(cos(theta)*cos(theta+1/6*pi)+sin(theta)*sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))]<BR>[ -3/2*(sin(theta)+sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), -3/2*(cos(theta)+cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3*(cos(theta)*sin(theta+1/3*pi)-sin(theta)*cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))]<BR> <BR> <BR> <BR>A =<BR> <BR>[ 3/2*(2/3*cos(theta)*(Ls+Lt*cos(2*theta))-2/3*cos(theta+1/3*pi)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*sin(theta+1/6*pi)*(-Ms-Lt*sin(2*theta+1/6*pi)))*(sin(theta+1/3*pi)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(2/3*cos(theta)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*cos(theta+1/3*pi)*(Ls-Lt*sin(2*theta+1/6*pi))-2/3*sin(theta+1/6*pi)*(-Ms+Lt*cos(2*theta)))*(-sin(theta)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(2/3*cos(theta)*(-Ms-Lt*sin(2*theta+1/6*pi))-2/3*cos(theta+1/3*pi)*(-Ms+Lt*cos(2*theta))-2/3*sin(theta+1/6*pi)*(Ls-Lt*cos(2*theta+1/3*pi)))*(sin(theta)+sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3/2*(2/3*cos(theta)*(Ls+Lt*cos(2*theta))-2/3*cos(theta+1/3*pi)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*sin(theta+1/6*pi)*(-Ms-Lt*sin(2*theta+1/6*pi)))*(cos(theta+1/3*pi)-sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3/2*(2/3*cos(theta)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*cos(theta+1/3*pi)*(Ls-Lt*sin(2*theta+1/6*pi))-2/3*sin(theta+1/6*pi)*(-Ms+Lt*cos(2*theta)))*(cos(theta)+sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(2/3*cos(theta)*(-Ms-Lt*sin(2*theta+1/6*pi))-2/3*cos(theta+1/3*pi)*(-Ms+Lt*cos(2*theta))-2/3*sin(theta+1/6*pi)*(Ls-Lt*cos(2*theta+1/3*pi)))*(cos(theta)+cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3*(2/3*cos(theta)*(Ls+Lt*cos(2*theta))-2/3*cos(theta+1/3*pi)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*sin(theta+1/6*pi)*(-Ms-Lt*sin(2*theta+1/6*pi)))*(cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3*(2/3*cos(theta)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*cos(theta+1/3*pi)*(Ls-Lt*sin(2*theta+1/6*pi))-2/3*sin(theta+1/6*pi)*(-Ms+Lt*cos(2*theta)))*(cos(theta)*cos(theta+1/6*pi)+sin(theta)*sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3*(2/3*cos(theta)*(-Ms-Lt*sin(2*theta+1/6*pi))-2/3*cos(theta+1/3*pi)*(-Ms+Lt*cos(2*theta))-2/3*sin(theta+1/6*pi)*(Ls-Lt*cos(2*theta+1/3*pi)))*(cos(theta)*sin(theta+1/3*pi)-sin(theta)*cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))]<BR>[ 3/2*(-2/3*sin(theta)*(Ls+Lt*cos(2*theta))+2/3*sin(theta+1/3*pi)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*cos(theta+1/6*pi)*(-Ms-Lt*sin(2*theta+1/6*pi)))*(sin(theta+1/3*pi)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(-2/3*sin(theta)*(-Ms-Lt*cos(2*theta+1/3*pi))+2/3*sin(theta+1/3*pi)*(Ls-Lt*sin(2*theta+1/6*pi))-2/3*cos(theta+1/6*pi)*(-Ms+Lt*cos(2*theta)))*(-sin(theta)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(-2/3*sin(theta)*(-Ms-Lt*sin(2*theta+1/6*pi))+2/3*sin(theta+1/3*pi)*(-Ms+Lt*cos(2*theta))-2/3*cos(theta+1/6*pi)*(Ls-Lt*cos(2*theta+1/3*pi)))*(sin(theta)+sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3/2*(-2/3*sin(theta)*(Ls+Lt*cos(2*theta))+2/3*sin(theta+1/3*pi)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*cos(theta+1/6*pi)*(-Ms-Lt*sin(2*theta+1/6*pi)))*(cos(theta+1/3*pi)-sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3/2*(-2/3*sin(theta)*(-Ms-Lt*cos(2*theta+1/3*pi))+2/3*sin(theta+1/3*pi)*(Ls-Lt*sin(2*theta+1/6*pi))-2/3*cos(theta+1/6*pi)*(-Ms+Lt*cos(2*theta)))*(cos(theta)+sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(-2/3*sin(theta)*(-Ms-Lt*sin(2*theta+1/6*pi))+2/3*sin(theta+1/3*pi)*(-Ms+Lt*cos(2*theta))-2/3*cos(theta+1/6*pi)*(Ls-Lt*cos(2*theta+1/3*pi)))*(cos(theta)+cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3*(-2/3*sin(theta)*(Ls+Lt*cos(2*theta))+2/3*sin(theta+1/3*pi)*(-Ms-Lt*cos(2*theta+1/3*pi))-2/3*cos(theta+1/6*pi)*(-Ms-Lt*sin(2*theta+1/6*pi)))*(cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3*(-2/3*sin(theta)*(-Ms-Lt*cos(2*theta+1/3*pi))+2/3*sin(theta+1/3*pi)*(Ls-Lt*sin(2*theta+1/6*pi))-2/3*cos(theta+1/6*pi)*(-Ms+Lt*cos(2*theta)))*(cos(theta)*cos(theta+1/6*pi)+sin(theta)*sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3*(-2/3*sin(theta)*(-Ms-Lt*sin(2*theta+1/6*pi))+2/3*sin(theta+1/3*pi)*(-Ms+Lt*cos(2*theta))-2/3*cos(theta+1/6*pi)*(Ls-Lt*cos(2*theta+1/3*pi)))*(cos(theta)*sin(theta+1/3*pi)-sin(theta)*cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))]<BR>[ 3/2*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(sin(theta+1/3*pi)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(-sin(theta)+cos(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(sin(theta)+sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3/2*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(cos(theta+1/3*pi)-sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3/2*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(cos(theta)+sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))-3/2*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(cos(theta)+cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi)), 3*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(cos(theta)*cos(theta+1/6*pi)+sin(theta)*sin(theta+1/6*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))+3*(1/3*Ls+1/3*Lt*cos(2*theta)-2/3*Ms-1/3*Lt*cos(2*theta+1/3*pi)-1/3*Lt*sin(2*theta+1/6*pi))*(cos(theta)*sin(theta+1/3*pi)-sin(theta)*cos(theta+1/3*pi))/(cos(theta)*sin(theta+1/3*pi)+cos(theta)*cos(theta+1/6*pi)-sin(theta)*cos(theta+1/3*pi)+sin(theta)*sin(theta+1/6*pi)+cos(theta+1/3*pi)*cos(theta+1/6*pi)+sin(theta+1/6*pi)*sin(theta+1/3*pi))]<BR><BR>请大侠指点下到底是什么错误啊，多谢了！
页:
[1]