是关于解含有贝塞尔函数的一个方程，求教各位高手
   clear
syms k a b
a=0.04;
b=0.1841;
z=((k^3+(k-1)/a^2).*(k^3-(k-1)/a^2)*besseli(1,k.*a).*besselj(1,k.*a))/(k^2*(k^3-(k-1)/a^2)*besseli(1,k.*a).*besselj(0,k.*a)+k^2*(k^3+(k-1)/a^2)*besseli(0,k.*a).*besselj(1,k.*a)+2*k^4*(b-1)/a*besseli(1,k.*a).*besselj(1,k.*a))
 >> solve(z)

        
 应该不影响解值的，a=0.05.  应该为多个解，是个关于频率的值，我求出来居然有复数。感觉不大对，是不是我解错了。