原代码：
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);
-
客服电话
-
APP下载
迪恩网络APP
随时随地掌握行业动态
-
官方微信
扫描二维码
关注迪恩网络微信公众号
