| 查看: 972 | 回复: 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 |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有118人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
- 金币: 68182.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.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
