matlab 没有精确解,Matlab求方程的解不是精确解

用solve求方程后的解代入到原式子，得到的结果不是0

syms a ZOe1 ZOo1 ZOe2 ZOo2 Z1 Z2 %自定义变量

H =

(Z1*sin(a)*(2*Z1 + ZOe1 + ZOo1)*(2*Z1 - ZOe1 - ZOo1 + 2*Z1*cos(a)^2 + ZOe1*cos(a)^2 + ZOo1*cos(a)^2)*(4*Z1^2*cos(a)^2 - 2*ZOe1*ZOo1 + ZOe1^2*cos(a)^2 + ZOo1^2*cos(a)^2 - ZOe1^2 - ZOo1^2 + 4*Z1*ZOe1*cos(a)^2 + 4*Z1*ZOo1*cos(a)^2 + 2*ZOe1*ZOo1*cos(a)^2)*i)/((ZOe1 + ZOo1)^2*(ZOe1*i + ZOo1*i)^2*(cos(2*a) + (4*ZOe2*ZOo2*(cos(2*a) + 1))/(ZOe2 - ZOo2)^2 - (cos(a)*(ZOe1 + ZOo1)*(ZOe2 + ZOo2)*(ZOo2 - ZOe2 + ZOe2*cos(a) + ZOo2*cos(a))*(ZOe2*i - ZOo2*i + ZOe2*cos(a)*i + ZOo2*cos(a)*i))/(2*ZOe1*ZOo1*sin(a)*(ZOe2 - ZOo2)^2)))

求得：

>> V=solve(H)

V =

pi + acos((ZOe1 + ZOo1)/(2*Z1 + ZOe1 + ZOo1))

pi - acos((ZOe1 + ZOo1)/(2*Z1 + ZOe1 + ZOo1))

pi + acos(((ZOe1 - 2*Z1 + ZOo1)/(2*Z1 + ZOe1 + ZOo1))^(1/2))

pi - acos(((ZOe1 - 2*Z1 + ZOo1)/(2*Z1 + ZOe1 + ZOo1))^(1/2))

将其中一个解代入原式子

syms a ZOe1 ZOo1 ZOe2 ZOo2 Z1 Z2 %自定义变量

a= pi + acos((ZOe1 + ZOo1)/(2*Z1 + ZOe1 + ZOo1));

H=(Z1*sin(a)*(2*Z1 + ZOe1 + ZOo1)*(2*Z1 - ZOe1 - ZOo1 + 2*Z1*cos(a)^2 + ZOe1*cos(a)^2 + ZOo1*cos(a)^2)*(4*Z1^2*cos(a)^2 - 2*ZOe1*ZOo1 + ZOe1^2*cos(a)^2 + ZOo1^2*cos(a)^2 - ZOe1^2 - ZOo1^2 + 4*Z1*ZOe1*cos(a)^2 + 4*Z1*ZOo1*cos(a)^2 + 2*ZOe1*ZOo1*cos(a)^2)*i)/((ZOe1 + ZOo1)^2*(ZOe1*i + ZOo1*i)^2*(cos(2*a) + (4*ZOe2*ZOo2*(cos(2*a) + 1))/(ZOe2 - ZOo2)^2 - (cos(a)*(ZOe1 + ZOo1)*(ZOe2 + ZOo2)*(ZOo2 - ZOe2 + ZOe2*cos(a) + ZOo2*cos(a))*(ZOe2*i - ZOo2*i + ZOe2*cos(a)*i + ZOo2*cos(a)*i))/(2*ZOe1*ZOo1*sin(a)*(ZOe2 - ZOo2)^2)));

>> H

H =

-(Z1*(1 - (ZOe1 + ZOo1)^2/(2*Z1 + ZOe1 + ZOo1)^2)^(1/2)*(2*Z1 + ZOe1 + ZOo1)*(2*Z1 - ZOe1 - ZOo1 + (2*Z1*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2 + (ZOe1*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2 + (ZOo1*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2)*((4*Z1^2*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2 - ZOe1^2 - ZOo1^2 - 2*ZOe1*ZOo1 + (ZOe1^2*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2 + (ZOo1^2*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2 + (4*Z1*ZOe1*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2 + (4*Z1*ZOo1*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2 + (2*ZOe1*ZOo1*(ZOe1 + ZOo1)^2)/(2*Z1 + ZOe1 + ZOo1)^2)*i)/((ZOe1 + ZOo1)^2*(ZOe1*i + ZOo1*i)^2*(cos(2*acos((ZOe1 + ZOo1)/(2*Z1 + ZOe1 + ZOo1))) + (4*ZOe2*ZOo2*(cos(2*acos((ZOe1 + ZOo1)/(2*Z1 + ZOe1 + ZOo1))) + 1))/(ZOe2 - ZOo2)^2 - ((ZOe1 + ZOo1)^2*(ZOe2 + ZOo2)*(ZOe2 - ZOo2 + (ZOe2*(ZOe1 + ZOo1))/(2*Z1 + ZOe1 + ZOo1) + (ZOo2*(ZOe1 + ZOo1))/(2*Z1 + ZOe1 + ZOo1))*(- ZOe2*i + ZOo2*i + (ZOe2*(ZOe1 + ZOo1)*i)/(2*Z1 + ZOe1 + ZOo1) + (ZOo2*(ZOe1 + ZOo1)*i)/(2*Z1 + ZOe1 + ZOo1)))/(2*ZOe1*ZOo1*(1 - (ZOe1 + ZOo1)^2/(2*Z1 + ZOe1 + ZOo1)^2)^(1/2)*(ZOe2 - ZOo2)^2*(2*Z1 + ZOe1 + ZOo1))))

结果不应该是0吗？