原代碼:
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);
原代碼:
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);
本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。