前面已经有syms Dx Jx Jy Jz real 为什么还有这一项 conj((4*Dx^2 + Jy^2 + Jz^2 + 2*Jy*Jz)^(1/2)) 求大神啊!!
syms Jx Jy Jz Dx B real;H=[Jz i*Dx -i*Dx Jx-Jy;-i*Dx -Jz Jx+Jy i*Dx;i*Dx Jx+Jy -Jz -i*Dx;Jx-Jy -i*Dx i*Dx Jz];[V,D]=eig(H);D1=diag(exp(diag(-B*D)));Pt=V*D1*inv(V),S=[0 0 0 -1;0 0 1 0;0 1 0 0;-1 0 0 0];L=Pt*S*conj(Pt)*S,Q=eig(L)
这是整个程序,哪里出错了啊!谢谢啊