大侠看下什么错误

linux853

程序如下：<BR>syms theta Ls Lt Ms;<BR>D=[2/3*cos(theta),2/3*cos(theta-2*pi/3),2/3*cos(theta+2*pi/3);<BR>     (-2/3)*sin(theta),(-2/3)*sin(theta-2*pi/3),(-2/3)*sin(theta+2*pi/3);<BR>      1/3,1/3,1/3];<BR>L=[Ls+Lt*cos(2*theta),-Ms-Lt*cos(2*(theta+pi/6)),-Ms-Lt*cos(2*(theta+5*pi/6));<BR>   -Ms-Lt*cos(2*(theta+pi/6)),Ls+Lt*cos(2*(theta-2*pi/3)),-Ms-Lt*cos(2*(theta-pi/2));<BR>   -Ms-Lt*cos(2*(theta+5*pi/6)),-Ms-Lt*cos(2*(theta-pi/2)),Ls+Lt*cos(2*(theta+2*pi/3))];<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>请大侠指点下到底是什么错误啊，多谢了！