| 查看: 2236 | 回复: 8 | ||
[求助]
使用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 求助！！谢谢了 ！ |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急！！！用matlab求解微分方程的问题
已经有10人回复
- Matlab 求解非线性方程组
已经有6人回复
- 用1stOpt解出非线性方程组的解作为初值用Matlab求解，解不出来
已经有4人回复
- 求助使用Matlab求解定积分的问题，望高手帮忙解答。
已经有1人回复
- MATLAB求解取车的最高车速
已经有0人回复
- 使用MatLab求解给定函数的最小值（变量使用我给定的值）
已经有1人回复
- 使用MATLAB对Fick第二定律编程求解
已经有4人回复
- matlab求解固定床二维均相模型求解轴向和径向温度梯度的实例加代码
已经有1人回复
- 使用matlab中LMI工具箱求解线性矩阵不等式程序
已经有3人回复
- matlab 隐函数求解
已经有6人回复
- 关于matlab求解非线性指数方程组出现问题
已经有3人回复
- matlab 求解riccati 方程
已经有4人回复
- 重金求助高手解决二元高次方程组MatlAB求解问题
已经有3人回复
- 使用matlab求解riccati 方程
已经有8人回复
- 求助matlab求解三角函数方程
已经有9人回复
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- matlab求解方程组的问题
已经有2人回复
- 使用Matlab求解数学方程
已经有1人回复
Flowers_Dai
银虫 (小有名气)
- 应助: 23 (小学生)
- 金币: 519.3
- 散金: 99
- 红花: 9
- 帖子: 212
- 在线: 50.6小时
- 虫号: 3575765
- 注册: 2014-12-04
- 专业: 可再生与替代能源利用中的
2楼2016-03-14 19:30:35
3楼2016-03-16 09:37:25
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2016-03-16 11:20:31
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2016-03-16 11:21:07
6楼2016-03-16 15:23:10
7楼2016-03-17 09:57:16
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
w51133248楼2016-03-17 11:04:12
9楼2016-03-17 11:19:34
20