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);

 

posted @ 2017-06-14 15:34  rldts  阅读(1981)  评论(0)    收藏  举报