24小时热门版块排行榜    

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

xwl196

银虫 (小有名气)

引用回帖:
6楼: Originally posted by dingd at 2015-04-09 16:33:16
版本问题。...

高版本的1stopt找不到,只能回到matlab中继续奋斗了
感谢!
7楼2015-04-09 16:40:33
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 【长期招聘】复旦大学 AI 制药实验室 博士后招聘 +3 第三种态度 2026-10-28 9/450 2026-11-02 13:06 by Grace_PepLib
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 12:20 by n0P4fnsXKbBK
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 BCdvBAEpt3zm 2026-11-01 4/200 2026-11-02 10:23 by dcc2yOIRGGTT
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 10:02 by dcc2yOIRGGTT
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:13 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 05:54 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:37 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 X5bHxOULxFqb 2026-10-30 27/1350 2026-11-02 05:28 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
信息提示
请填处理意见