原代碼:
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);