24小时热门版块排行榜    

查看: 719  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ThomsonW

铁虫 (初入文坛)

[求助] matlab 已有2人参与

我需要解一个非线性三元方程组,见图片。其中P&/P=0.1/3,m=50,k1=6.7
,k2=5.56E+13,k3=1.24E+11。编辑的matlab程序如下,可出现了些问题,我不知道怎么回事?求高人指点
function dd
x0 = [0.01  0.4  0.2]';
x = fsolve(@NonlinEqs,x0);
%               
function f = NonlinEqs(x);
tmp= 50 - 2*x(1) - x(2) - x(3);
f(1) =(51-x(1))*(x(1)-x(2))*(x(1)-x(3))*0.1/(3*(1-x(1))*tmp^2) -6.7;
f(2) = x(2)*(x(2)+x(3))/((x(1)-x(2))*tmp) -5.56E+13;
f(3) = x(3)*(x(2)+x(3))/((x(1)-x(3))*tmp) -1.24E+11;
计算结果如下:
>> dd

No solution found.

fsolve stopped because the relative size of the current step is less than the
default value of the step size tolerance squared, but the vector of function values
is not near zero as measured by the default value of the function tolerance.

<stopping criteria details>


fsolve stopped because the relative norm of the current step, 6.778107e-013, is less than
max(options.TolX^2,eps) = 1.000000e-012. However, the sum of squared function values,
r = 3.091375e+027, exceeds sqrt(options.TolFun) = 1.000000e-003.

Optimization Metric                                               Options
relative norm(step) =  6.78e-013                max(TolX^2,eps) =  1e-012 (default)
r =  3.09e+027                                    sqrt(TolFun) = 1.0e-003 (default)

我不知道上面说的是怎么回事,求高人帮助!!!!!!
matlab
CIYO$ZU85C~TZ2%NP599D0Q.jpg



[ Last edited by ThomsonW on 2014-3-10 at 21:00 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ThomsonW

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by 戴钢盔的猪头 at 2014-03-11 00:52:18
看你跑完程序的提示,求解采用的算法应该是信赖域折线(狗腿)法。程序在你给的初值的情况下开始迭代,发现迭代步已经缩小到阈值以下,目标函数仍未收敛。你可以改用LM算法试试,或者改变初值,或者检查下你的函数模 ...

嗯,解非线性方程组一般就是用fsolve函数,LM算法我不懂、、、
4楼2014-03-11 08:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

AV兄

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
ThomsonW: 金币+1, 有帮助 2014-03-11 08:52:58
x0 = [0.01  0.4  0.2]';  
这句话后面有个撇号  去掉试试  我也是新手  呵呵
欢迎光临我的 AI 社区 http://keepdl.org
2楼2014-03-10 22:47:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ThomsonW: 金币+10, ★★★很有帮助 2014-03-11 08:52:44
看你跑完程序的提示,求解采用的算法应该是信赖域折线(狗腿)法。程序在你给的初值的情况下开始迭代,发现迭代步已经缩小到阈值以下,目标函数仍未收敛。你可以改用LM算法试试,或者改变初值,或者检查下你的函数模型。

[ 发自手机版 http://muchong.com/3g ]
3楼2014-03-11 00:52:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ThomsonW

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by 戴钢盔的猪头 at 2014-03-11 00:52:18
看你跑完程序的提示,求解采用的算法应该是信赖域折线(狗腿)法。程序在你给的初值的情况下开始迭代,发现迭代步已经缩小到阈值以下,目标函数仍未收敛。你可以改用LM算法试试,或者改变初值,或者检查下你的函数模 ...

主要就是解图片中的方程组x1,x2,x3.其它字母都给有值,求指点!!!
5楼2014-03-11 08:54:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +12 医学老男孩 2026-08-13 32/1600 2026-08-14 08:19 by 医学老男孩
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +71 AI8RaGaPaSCR 2026-08-07 74/3700 2026-08-14 08:15 by am8QyKH2Q7LX
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 30/1500 2026-08-14 07:55 by 医学老男孩
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见