| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 801 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
BobbyYoung銀蟲(chóng) (初入文壇)
|
[求助]
matlab解非線性方程組 已有2人參與
|
||
|
本人matlab菜鳥(niǎo)一個(gè),解非線性方程組時(shí)出現(xiàn)以下提示,請(qǐng)問(wèn)是什么意思?該怎么解決? Optimizer appears to be converging to a point which is not a root. Relative function value changing by less than max(options.TolFun^2,eps) but sum-of-squares of function values is greater than or equal to sqrt(options.TolFun) Try again with a new starting guess. m = 5.0482 5.9987 7.1051 fval = 1.0e+016 * 3.8464 0.0000 -0.0000 |

銀蟲(chóng) (初入文壇)

主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

鐵桿木蟲(chóng) (著名寫手)
ส็็็
|
fsolve歸根到底還是迭代方法; 如果方程組自己沒(méi)有0點(diǎn),或者迭代法的Jacobian矩陣在某些步驟接近singular矩陣(數(shù)值不滿秩)的情況下,會(huì)出現(xiàn)上面的情況. 如果一定要用fsolve,說(shuō)明默認(rèn)的trust region + dogleg 算法行不通; 可以改用: algorithm_name = 'trust-region-dogleg'; %這是默認(rèn)算法 %'levenberg-marquardt' % 針對(duì)Jacobian矩陣非奇異效果較好的 options = optimset('Algorithm',algorithm_name); %或者通過(guò)設(shè)置下面兩個(gè)選項(xiàng)用Gauss-Newton+Line search的方法: and NonlEqnAlgorithm options = optimset('LargeScale ','off',...%對(duì)規(guī)模較小問(wèn)題適用 'NonlEqnAlgorithm','gn', ... % 選擇G-N算法 'LineSearchType','cubicpoly'); % 用三次步長(zhǎng)因子 然后把fsolve的options參數(shù)加上(參考help fsolve調(diào)用格式) 如果fsolve仍解決不了問(wèn)題(非線性很強(qiáng)的,)可以考慮用最小二乘的方式把方程組改成殘差平方和,用最優(yōu)化的方式,粒子群,差分演化等(PSO或Differential evolution)求近似解即可. |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|