应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 使用MATLAB求解3元3次方程组的问题
91/1返回列表
查看: 2006  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

w5113324

铁虫 (初入文坛)

[求助] 使用MATLAB求解3元3次方程组的问题 已有2人参与

方程如下：
T=1135;
Ptot=346;
P0=101;
Uf=0.704;
syms x y z;
Kr=exp(-2.63121*10^(-11)*T^4+1.24065*10^(-7)*T^3-2.25232*10^(-4)*T^2+1.95028*10^(-1)*T-66.1395);
Ks=exp(5.47301*10^(-12)*T^4-2.57479*10^(-8)*T^3+4.63742*10^(-5)*T^2-3.91500*10^(-2)*T+13.2097);
H2_in=16;CO_in=1.124;CO2_in=0.3;CH4_in=0.1;H2O_in=1.656;
n_tot=H2_in+CO_in+CO2_in+CH4_in+H2O_in;
eq1=(H2_in+3*x+y-z)^3*(CO_in-y+x)/((CH4_in-x)*(H2O_in-x-y+z)*n_tot^2)*(Ptot/P0)^2-Kr;
eq2=(H2_in+3*x+y-z)*(CO2_in-y)/((CO_in-y+x)*(H2O_in-x-y+z))-Ks;
eq3=z/(3*x+y)-Uf;
d=solve(eq1,eq2,eq3,'x','y','z');

求解后
警告: The solutions are parametrized by the symbols:
z = RootOf(z1^7 + (132884004*z1^6)/578125 +
(2119387706984217863273427132148046141508842279290362609727646198416856373*z1^5)/82123996273443879090756655219912020554359217704965320709505024000000 +
(6207335911728807795095662535116719152621565796678720774245761139522437154907*z1^4)/3798234827646779407947495303920930950639113818854646082814607360000000 +
(34919727432059595159713427078806550155318489300274916294518425509603094873863880435946079201671*z1^3)/580286932331463977597529434899502961023911489788027077179903793533437429678080000000000
+
(18286981928998639134663466563712586815724753527692813345702385641483117729280080289010590372106349*z1^2)/18133966635358249299922794840609467531997234055875846161871993547919919677440000000000000
+
(405661360233660275394850404499713117718113403204186863108828571787910707243798644736316320000947*z1)/70835807169368161327823417346130732546864195530765024069812474796562186240000000000000
+
33486843534555442984426309387901023909974465889584168019701785461544480007449703963860162330241/5534047435106887603736204480166463480223765275841017505454099593481420800000000000000,
z1) minus ({0} union RootOf(z1^6 +
(6774161150507946194092693822747561712566944460123370379580513419157650686440138290544*z1^5)/36570226538202626850896525893163958194276642217230977356667133758946469931248546875
+
(2178852024748883567751702897023743255495979447146828016997047008162749890257676724488951209399436414963086436281412407779986952156217*z1^4)/121199731338256297050878765449711016063934318578465466375010298420404532341066683284990608003352086435147857148336768496435200000
+
(13485725654430355568124649579869722309326655361776505617123219622688606479082343461173831946742834380971926372196792300676852340515308707*z1^3)/14013718935985884346507857255122836232392405585635069549610565754859274051935835254827039050387584994063970982776438857400320000000
+
(23201914853009842666708148721477926045546589776386973981940053329275189775752606484530375802591679479819136970010280811856492267885078387903926245623379121*z1^2)/856395493267256992263561212733654854689243716373301622458157945769443026862579406907434228013215001864658484357968938497480838051407069184000000000
+
(1490548085966046470167013641628008847628797008533819532435356182227702725505809000496370558016265616633485126307281732880338518064603684206021797103946104889*z1)/8363237238938056565073839968102098190324645667708023656817948689154717059204877020580412382941552752584555511308290415014461309095772160000000000000
+
12737781808759804950818611591386062116616220827244642646240037334836159306443888698155112948410337694678200008701493851386973122211063362030084087463521717/65337790929203566914639374750797642111911294278968934818890224134021227025038101723284471741730880879566839932096018867300478977310720000000000000,
z1) union RootOf(z1^6 +
(14813587029224799491896936188376888887216201516485969798208050019814634723659472290798776*z1^5)/46630401610458426065030877582261451931459709739395221764539286719747307128491766453125
+
(3418205991406512121115601666307370853450057566285240086969664593665591438671118189529152334887384258743430115980915458892169050219141373463*z1^4)/96588001401199456667774227513859821896222722624270116936266562893885990076832913159588751610248603461498876227638767426469888000000000
+
(114905440806174197433406381230300285280659040093211613577691129579621369445030477834892495145860375070832836755077097585253283296566565678269*z1^3)/55839938310068435886056975281450209533753761517156161353779106673027838013169027920387247024674973876179037819103662418427904000000000
+
(2195877091058608703927703773400965907325864796954260517983312358992652579896772630992787741838231486585883016915578379139586143876595790109946437*z1^2)/41321554349450642555682161708273155054977783522695559401796538938040600129745080661086562798259480668372487986136710189636648960000000000
+
(9357436690026630802021295981875470313479072178351230662537614951916965403393341239800666317912966337814780244936903113040550322048520727127536447*z1)/20176540209692696560391680521617751491688370860691191114158466278340136782102090166546173241337637045103753899480815522283520000000000000
+
182149648085692379586806621600313508984528014389778721435706293856207323374547175482285430739116551242027440808860142187582457833961489721547857/315258440776448383756120008150277367057630794698299861158726035599064637220345158852283956895900578829746154679387742535680000000000000,
z1))

> In solve at 190

求助！！谢谢了 ！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-03-14 09:30:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Flowers_Dai

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你这个方程缺少条件吧，六个未知数，三个方程怎么解啊？是不是少给条件了，约束条件之类的，光这样应该解不出来的。
赞一下(1人) 回复此楼
2楼2016-03-14 19:30:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w5113324

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by Flowers_Dai at 2016-03-14 19:30:35
你这个方程缺少条件吧，六个未知数，三个方程怎么解啊？是不是少给条件了，约束条件之类的，光这样应该解不出来的。

我只有X,Y,Z三个未知数啊。其他的都是已知量~
赞一下 回复此楼
3楼2016-03-16 09:37:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
1stOpt数值解：

x: -1.87138571848282
y: -25.9960493233376
z: -22.2535793469822
赞一下(1人) 回复此楼
4楼2016-03-16 11:20:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

还有一组：

x: -3.76783079837869
y: -0.165068355947246
z: -8.07386676876265
赞一下(1人) 回复此楼
5楼2016-03-16 11:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w5113324

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2016-03-16 11:20:31
1stOpt数值解：

x: -1.87138571848282
y: -25.9960493233376
z: -22.2535793469822

大神！！怎么求出来的，求看看程序
赞一下 回复此楼
6楼2016-03-16 15:23:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w5113324

铁虫 (初入文坛)
引用回帖:
6楼: Originally posted by w5113324 at 2016-03-16 15:23:10
大神！！怎么求出来的，求看看程序...

为什么我用1stopt跑出来还是没结果啊
赞一下 回复此楼
7楼2016-03-17 09:57:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
w5113324: 金币+30 2016-03-17 11:19:17
盗版的也应该可以运行吧：
CODE:
Constant T=1135,Ptot=346,P0=101,Uf=0.704, H2_in=16,
         CO_in=1.124,CO2_in=0.3,CH4_in=0.1,H2O_in=1.656;
ConstStr Kr=exp(-2.63121*10^(-11)*T^4+1.24065*10^(-7)*T^3-2.25232*10^(-4)*T^2+1.95028*10^(-1)*T-66.1395),
         Ks=exp(5.47301*10^(-12)*T^4-2.57479*10^(-8)*T^3+4.63742*10^(-5)*T^2-3.91500*10^(-2)*T+13.2097),
         n_tot=H2_in+CO_in+CO2_in+CH4_in+H2O_in;
Function (H2_in+3*x+y-z)^3*(CO_in-y+x)/((CH4_in-x)*(H2O_in-x-y+z)*n_tot^2)*(Ptot/P0)^2-Kr;
         (H2_in+3*x+y-z)*(CO2_in-y)/((CO_in-y+x)*(H2O_in-x-y+z))-Ks;
         z/(3*x+y)-Uf;
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

w5113324
8楼2016-03-17 11:04:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w5113324

铁虫 (初入文坛)
送红花一朵
引用回帖:
8楼: Originally posted by dingd at 2016-03-17 11:04:12
盗版的也应该可以运行吧：

Constant T=1135,Ptot=346,P0=101,Uf=0.704, H2_in=16,
         CO_in=1.124,CO2_in=0.3,CH4_in=0.1,H2O_in=1.656;
ConstStr Kr=exp(-2.63121*10^(-11)*T^4+1.24065*10^(-7)*T^3-2. ...

太感谢！！！！！
回复此楼
9楼2016-03-17 11:19:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 w5113324 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭