| 查看: 2110 | 回复: 1 | |||
[求助]
四个方程，5个未知数，求解其中一个未知数关于另外一个未知数的函数关系 已有1人参与
|
|
syms e syms par cL=286.4; R=8.314; Tm=865; k=0.814; rS=1.338; rL=1.278; p=(0.1e+6)+(par*10^6); cLb=(cL-p/R/Tm)/(1-(1-k)*(1-e)); cS=k*cLb*rS/rL; e=(cL*rS-cS*rL)*R*Tm/((cL*rS-cS*rL)*R*Tm+rL*p); 需要得到e与par之间的函数关系，并作图，请各位大神解答，谢谢！  |
回复此楼» 猜你喜欢
遇见不省心的家人很难过
已经有9人回复
-
博士延得我，科研能力直往上蹿
已经有3人回复
-
退学或坚持读
已经有24人回复
-
免疫学博士有名额，速联系
已经有14人回复
-
面上基金申报没有其他的参与者成吗
已经有4人回复
-
多组分精馏求助
已经有6人回复
【答案】应助回帖
感谢参与，应助指数 +1
|
syms e par cL=286.4; R=8.314; Tm=865; k=0.814; rS=1.338; rL=1.278; p=(0.1e+6)+(par*10^6); cLb=(cL-p/R/Tm)/(1-(1-k)*(1-e)); cS=k*cLb*rS/rL; Equ=e-(cL*rS-cS*rL)*R*Tm/((cL*rS-cS*rL)*R*Tm+rL*p); fsolv=solve(Equ==0,e) f1=eval(['@(par)' vectorize(fsolv(1))]) ; f2=eval(['@(par)' vectorize(fsolv(2))]) ; x=1:0.1:20; plot(x,f1(x),x,f2(x),'r') 算的有两个解，不知道对不对，楼主自己判断吧  |
2楼2015-09-16 16:50:52