24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【求助】帮忙看一个程序已有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.

回复此楼

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有9人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

好好学习，天天向上。

1楼2010-06-14 21:05:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

引用回帖:

Originally posted by change0618 at 2010-06-15 18:14:03:
对于单组份吸附方程f=b/(1+b*p); 当b*p~=-1时, 积分式子为F=ln(1+b*p) ;
你的以下三行命令q1=integral(f1,0,p1,eps); q2=integral(f2,0,p2,eps);q1=q2;
可以变成了ln(1+b1*p1)= ln(1+b2*p2);
即b1*p1= b2*p2;
...

有时候，是非线性LANGMUIR公式了，
f=(m1*p/(1+m1*p))+(m2*p/(1+m2*p))
甚至三个这样的公式连加。
我就是想看看IAST理论该怎么编，手算毕竟很慢。