| 查看: 2389 | 回复: 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. |
回复此楼» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有13人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组-博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
13楼2010-06-15 18:51:27
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