| 9 | 1/1 | 返回列表 |
| 查看: 2044 | 回復: 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 求助!!謝謝了 ! |
銀蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
鐵桿木蟲 (職業(yè)作家)
鐵桿木蟲 (職業(yè)作家)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
[考研] 化工專碩求調(diào)劑 +5 | question挽風 2026-03-24 | 5/250 |
|
|
[考研] 本科211生物醫(yī)學工程085409求調(diào)劑339分 +3 | 里子木yy 2026-03-29 | 3/150 |
|
|
[碩博家園] 求調(diào)劑 有機化學考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
[考研] 一志愿北京工業(yè)大學,324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 070305高分子化學與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 312,生物學求調(diào)劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學材料學碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|