24小时热门版块排行榜    

查看: 721  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +8 a089 2026-08-07 8/400 2026-08-14 18:46 by beefly
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +6 Tide man 2026-08-14 7/350 2026-08-14 18:36 by 正能量1斤
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +22 医学老男孩 2026-08-13 48/2400 2026-08-14 18:33 by wuguocheng
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[硕博家园] 读博的好处 +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 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +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邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见