| 查看: 2247 | 回复: 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固定段从头到尾一模一样
已经有9人回复
是这周出结果还是下周出结果?
已经有3人回复
应该是93bebmhtak前后十一个字符比较关键
已经有27人回复
小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊
已经有9人回复
有时候,自然基金真的不能太认真 (我的申报经验)
已经有4人回复
帮忙看看fileCode
已经有13人回复
为什么网上很多人说本周 12号出结果
已经有7人回复
读博的好处
已经有3人回复
什么时候出结果,有咨询渠道???
已经有3人回复
分享一下我之前已中青C的计划书的filecode
已经有4人回复
» 本主题相关价值贴推荐,对您同样有帮助:
急!!!用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人回复
3楼2016-03-16 09:37:25
Flowers_Dai
银虫 (小有名气)
- 应助: 23 (小学生)
- 金币: 519.3
- 散金: 99
- 红花: 9
- 帖子: 212
- 在线: 50.6小时
- 虫号: 3575765
- 注册: 2014-12-04
- 专业: 可再生与替代能源利用中的
2楼2016-03-14 19:30:35
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










回复此楼