| 查看: 1950 | 回复: 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. |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
14楼2010-06-15 19:04:43
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