| 查看: 2431 | 回复: 16 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】帮忙看一个程序 已有3人参与
|
|||
|
%IAST(理想吸附溶液理论通过单组份的吸附数据获得多组分的吸附数据） function [p1,p2,x1]=IAST(b1,b2,p,y1) f1=b1/(1+b1*p); f2=b2/(1+b2*p); q1=integral(f1,0,p1,eps); q2=integral(f2,0,p2,eps); q1=q2; p*y1=p1*x1; p*(1-y1)=p2*(1-x1); end function q=integral(f,a,b,eps) %被积函数：f %积分左端点：a %积分右端点：b %eps:精度 %积分结果:q %step 积分的子区间数 if(nargin==3) eps=1.0e-4; %默认精度为0.0001 end n=1; h=(b-a)/2; q1=0; q2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h; while abs(q2-q1)>eps n=n+1; h=(b-a)/n; q1=q2; q2=0; for i=0:n-1 %第n次复合梯形积分 x=a+h*i; x1=x+h; q2=q2+(h/2)*(subs(sym(f),findsym(sym(f)),x)+... subs(sym(f),findsym(sym(f)),x1)); end end q=q2; end report 报告： 2: Function return value 'p1' appears to never be set. 2: Function return value 'p2' appears to never be set. 2: Function return value 'x1' appears to never be set. 2: Input argument 'y1' appears never to be used. 5: The value assigned here to variable 'q1' might never be used. 7: The value assigned here to variable 'q1' might never be used. 8: Parse error at '=': usage appears to be invalid MATLAB syntax. 9: Parse error at '=': usage appears to be invalid MATLAB syntax. |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
16楼2010-06-15 20:01:37
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2010-06-14 21:57:30
3楼2010-06-15 00:41:23
4楼2010-06-15 13:15:13
200