| 查看: 789 | 回复: 3 | |||
BobbyYoung银虫 (初入文坛)
|
[求助]
matlab解非线性方程组 已有2人参与
|
|
本人matlab菜鸟一个，解非线性方程组时出现以下提示，请问是什么意思？该怎么解决？ 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 |
回复此楼» 猜你喜欢
德国亥姆霍兹Hereon中心招收两位医用镁合金腐蚀与LPSO相变方向2026公派博士生
已经有0人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，推荐试试
已经有11人回复
-
物理学I论文润色/翻译怎么收费?
已经有248人回复
-
2026-CJ开始申报了
已经有1人回复
-
西安电子科学大学杭州研究院刘丽香教授招收智能多模态传感器和微型储能器件方向博士
已经有13人回复
-
重庆交通大学光子学微结构与器件课题组招收2026年硕士研究生信息
已经有0人回复
-
一志愿郑大材料学硕298分，求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB解非线性方程组
已经有9人回复
- MATLAB求解五元非线性方程组，在线等
已经有7人回复
- matlab 非线性微分方程求解
已经有3人回复
- 求助解复杂非线性方程组的好的方法
已经有24人回复
- matlab非线性方程组该如何编程。求助
已经有10人回复
- 用matlab解非线性方程组
已经有4人回复
- 用matlab求解一个非线性方程组的解
已经有4人回复
- matlab 求解非线性方程组未果
已经有10人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- MATLAB求解非线性方程组
已经有5人回复
- 求Matlab解非线性方程代码
已经有12人回复
- matlab解非线性方程组解不开
已经有3人回复
- matlab 求解非线性方程 解出虚数解
已经有7人回复
- 10次代数方程怎么找解析解 最好详细给出matlab程序 谢谢
已经有5人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- Matlab：Newton-Raphson 解线非线性性方程组的问题
已经有12人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- matlab求解非线性方程组
已经有16人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-08-06 11:37:11
BobbyYoung
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 291
- 红花: 1
- 帖子: 31
- 在线: 24.7小时
- 虫号: 3155993
- 注册: 2014-04-22
- 性别: GG
- 专业: 无机化工
|
谢谢您，我更改初值试了一下，还是不行。因为我是想计算模拟水盐体系的溶解度数据，输入实验值的近似值，为 >>m0=[4.24,4.65,8.13]； >>[m,fval]=fsolve('FphaseEq',m0) m = 6.5588 7.2533 7.8359 fval = 1.0e+016 * 3.8464 0.0000 -0.0000 （以下每次函数值均相同，未列出） 然后我换了几个初值试了一下，结果如下（请见附件）。 由于是计算溶解度数据，计算值与实验值希望不会偏差太大，应该怎么解决matlab中提到的问题？小弟对于matlab所知实在十分有限，请前辈指点，非常感谢！  1.jpg |
3楼2014-08-06 15:51:49
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
【答案】应助回帖
感谢参与，应助指数 +1
|
fsolve归根到底还是迭代方法; 如果方程组自己没有0点,或者迭代法的Jacobian矩阵在某些步骤接近singular矩阵(数值不满秩)的情况下,会出现上面的情况. 如果一定要用fsolve,说明默认的trust region + dogleg 算法行不通; 可以改用: algorithm_name = 'trust-region-dogleg'; %这是默认算法 %'levenberg-marquardt' % 针对Jacobian矩阵非奇异效果较好的 options = optimset('Algorithm',algorithm_name); %或者通过设置下面两个选项用Gauss-Newton+Line search的方法: and NonlEqnAlgorithm options = optimset('LargeScale ','off',...%对规模较小问题适用 'NonlEqnAlgorithm','gn', ... % 选择G-N算法 'LineSearchType','cubicpoly'); % 用三次步长因子 然后把fsolve的options参数加上(参考help fsolve调用格式) 如果fsolve仍解决不了问题(非线性很强的,)可以考虑用最小二乘的方式把方程组改成残差平方和,用最优化的方式,粒子群,差分演化等(PSO或Differential evolution)求近似解即可. |
4楼2014-08-06 21:06:04
10