关于Matlab的roots函数

很显然，如果是无参数的式子可以直接用 ‘double()’ 显示出结果。但如果像上图式子中含有参数，就很烦。
然后找到一个人介绍了 ‘roots’ 函数。 原文
但是按照对方的方法尝试后发现结果不太对劲，就去了Matlab官网上去查了查，果然不对劲。

前面都是废话。下面介绍roots。

roots

函数名称：多项式根，以列向量的形式多项式的根

p=[ ];
r=roots(p);

p 是一个包含 n+1 多项式系数的向量。
上图中的问题就可以写为

p=[1 -13*b 40*b^2 -32*b^3];
roots(p);

得到结果

ans =(13*b)/3 + ((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3) + (49*b^2)/(9*((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3))(13*b)/3 - ((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3)/2 - (49*b^2)/(18*((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3)) - (3^(1/2)*(((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3) - (49*b^2)/(9*((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3)))*1i)/2(13*b)/3 - ((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3)/2 - (49*b^2)/(18*((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3)) + (3^(1/2)*(((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3) - (49*b^2)/(9*((-(1264*b^6)/27)^(1/2) + (289*b^3)/27)^(1/3)))*1i)/2

遗憾的是，由于我只需要第三个根，而matlab里输出的根好像没什么顺序，我的问题还是没解决orz