我求解符号方程:
>> A=sym('[34618-150*w^2,-16502,0,0,0;-16502,31144-150*w^2,-14642,0,0;0,-14642,27018-150*w^2,-12376,0;0,0,-12376,22447-150*w^2,-10071;0,0,0,-10071,10071-150*w^2]')
A =
[ 34618-150*w^2, -16502, 0, 0, 0]
[ -16502, 31144-150*w^2, -14642, 0, 0]
[ 0, -14642, 27018-150*w^2, -12376, 0]
[ 0, 0, -12376, 22447-150*w^2, -10071]
[ 0, 0, 0, -10071, 10071-150*w^2]
>> solve(det(A),'w')
ans =
solvelib::Union({[-z1^(1/2), 1], [z1^(1/2), 1]}, z1, RootOf(X19^5 - (20883*X19^4)/25 + (17855633*X19^3)/75 - (11413518430603*X19^2)/421875 + (33558897596166236*X19)/31640625 - 631449249611223316/87890625, X19))
如何将此符号答案转换为数字?
我使用'double'但答案不正确?
我使用vpa,但是当行列式替代时不能为零!!!
请帮助我!
谢谢