24小时热门版块排行榜    

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

xwl196

银虫 (小有名气)

[求助] matlab求解含中间变量的非线性方程组 已有1人参与

方程组如附件所示,其中:
T=328.2;
ρmc=1.0254;
Eii=0.432*195.2;
C=1.23;
Eis=-1250;
ρads为中间变量
求大神指点:1.如何求解想x和y;2.或者如何从方程(1)中把ρads表示成关于x的表达式。

matlab求解含中间变量的非线性方程组
1.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
5楼: Originally posted by xwl196 at 2015-04-09 16:29:27
你用的1stopt是哪个版本的,我用1stopt1.5运行的结果如下:
"model nex"

====== 结果 ======

迭代数: 16
计算用时(时:分:秒:毫秒): 00:00:00:50
计算中止原因: 达到收敛判定标准
优化算法: 麦 ...

版本问题。
6楼2015-04-09 16:33:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
nex有数值吗,多少?
2楼2015-04-09 14:36:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwl196

银虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2015-04-09 14:36:36
nex有数值吗,多少?

ρb=x=[[0.0116
0.0240
0.0502
0.0804
0.1170
0.1608
0.2152
0.2972
0.4123
0.5257
0.6100];
所以问题转化为在方程(1)中用matlab的solve求解ρads,然后将ρads带入方程(2)中求解nex就好了!
matlab程序:
>> T=328.2;
rhomc=1.0254;
Eii=0.432*195.2;
C=1.23;
Eis=-1250;
x=[0.0116
0.0240
0.0502
0.0804
0.1170
0.1608
0.2152
0.2972
0.4123
0.5257
0.6100];
>> syms rhoads
>> [rhoads]=solve(log(rhoads*(rhomc-x)/(x*(rhomc-rhoads)))+(7*rhoads-8*x)/rhomc*Eii/T+Eis/T==0)
运行后出现的提示是:
??? Error using ==> mupadmex
Error in MuPAD command: Array sizes must match.

Error in ==> sym.plus at 14
    X = mupadmex('symobj::zip',A.s,B.s,'_plus');
没找到问题出现在哪?
3楼2015-04-09 15:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xwl196: 金币+10, ★★★很有帮助 2015-04-09 16:40:44
Solve命令是求符号解,不一定有的。可以是数值解。
下面是用1stOpt求解,很简单:
CODE:
Constant T=328.2,rhomc=1.0254,Eii=0.432*195.2,C=1.23,Eis=-1250;
LoopConstant x=[0.0116,0.0240,0.0502,0.0804,0.1170,0.1608,0.2152,0.2972,0.4123,0.5257,0.6100];
PassParameter nex=2*C*(rhoads-x)/rhomc;
Function
ln(rhoads*(rhomc-x)/(x*(rhomc-rhoads)))+(7*rhoads-8*x)/rhomc*Eii/T+Eis/T;

Ñ-»·3£êy x        rhoads        nex
0.0116        0.257851015773962        0.590771892728639
0.024        0.378121213925542        0.849559378054255
0.0502        0.520370658449439        1.1279693970993
0.0804        0.618855300376744        1.29178860827656
0.117        0.699976313181814        1.3985973575456
0.1608        0.769244989810605        1.45969833716997
0.2152        0.831497991791605        1.47853818978677
0.2972        0.896153280504329        1.43692712116311
0.4123        0.952359132344292        1.2956363034591
0.5257        0.984696573619486        1.10116205490924
0.61        0.999729281534294        0.934985403329788
4楼2015-04-09 16:05:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 12/600 2026-08-12 01:27 by 南京中西医
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 29/1450 2026-08-12 00:24 by zhanghaozhu
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 9/450 2026-08-11 22:12 by 医学老男孩
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见