24小时热门版块排行榜    

查看: 716  |  回复: 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的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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的回帖
查看全部 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的回帖

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的回帖

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后的代码验证其规律 +10 医学老男孩 2026-08-13 27/1350 2026-08-13 21:43 by chemsynth
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 10/500 2026-08-13 20:42 by hmhminy
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 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 云上清扬
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见