MATLAB中-27开3次方得不到-3的原因

原代码:

y=@(x) x.^(1/3);

运行结果:

  >> y(-27)

  ans =

  1.5000 + 2.5981i

原因:

用(-8)^(1/3)时,matlab其实是调用C = power(A,B) 命令的,即C = power(-8,1/3) 。
如果B的绝对值值小于1,函数返回的是复数根,如果想得到实数根,用nthroot 函数。

修改方案1:

y=@(x) nthroot(x,3);

修改方案2:

y=@(x) sign(x).*abs(x).^(1/3);
原文地址:https://www.cnblogs.com/qrlozte/p/7009149.html