24小时热门版块排行榜    

查看: 722  |  回复: 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算法试试,或者改变初值,或者检查下你的函数模 ...

主要就是解图片中的方程组x1,x2,x3.其它字母都给有值,求指点!!!
5楼2014-03-11 08:54:57
已阅   回复此楼   关注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算法试试,或者改变初值,或者检查下你的函数模 ...

嗯,解非线性方程组一般就是用fsolve函数,LM算法我不懂、、、
4楼2014-03-11 08:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +24 医学老男孩 2026-08-13 50/2500 2026-08-15 06:38 by 学员8dgXkO
[基金申请] filecode +7 cratir 2026-08-14 11/550 2026-08-15 06:19 by 学员8dgXkO
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见