| 查看: 1952 | 回复: 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 求助!!谢谢了 ! |
» 猜你喜欢
聘U V热熔胶研究人员
已经有10人回复
职称评审没过,求安慰
已经有7人回复
求助文献
已经有3人回复
垃圾破二本职称评审标准
已经有10人回复
投稿返修后收到这样的回复,还有希望吗
已经有8人回复
三无产品还有机会吗
已经有6人回复
谈谈两天一夜的“延安行”
已经有13人回复
氨基封端PDMS和HDI反应快速固化
已经有11人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
» 本主题相关价值贴推荐,对您同样有帮助:
急!!!用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
- 专业: 可再生与替代能源利用中的
【答案】应助回帖
感谢参与,应助指数 +1
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与,应助指数 +1
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
w5113324: 金币+30 2016-03-17 11:19:17
w5113324: 金币+30 2016-03-17 11:19:17
送红花一朵













回复此楼