[討論] 數值解如何限定解為時數？

[seal46825 (Vanish)]

最近要跑的程式需要解一個函數f(x)很複雜
沒辦法用符號解去求解
所以我用數值的去解

他解有很多個 但我想要的是實數的解
我網路查如何得到實數解都只找到符號解的方法
想請問用數值解如何限定解是實數呢

我是用fsolve求解
CODE大概如下:

eq=@(x) f(x)
[ans,fval]=fsolve(eq,1)

不曉得如何讓ans只取實的解

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 180.217.158.18
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1525167137.A.55D.html
→ seal46825: 或是能不能讓他全部解都顯示出來
推 LiamIssac: real(x)
→ LiamIssac: 看看有沒有isreal這類函數 然後挑output是1的
推 sunev: 轉成optimized的問題，求f(x)^2+imag(x)*100000的最小值
→ seal46825: 感謝幫忙！