| 查看: 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人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求助matlab一程序运行的问题，帮忙看看哪里不对
已经有4人回复- 拜托大家帮忙看个简单程序
已经有11人回复
- MC一个简单的程序，运行不了，哪位大侠帮看看怎么回事。多谢！
已经有7人回复
- 帮忙看一段程序
已经有3人回复
- 帮忙看个文件，非常感谢
已经有3人回复
- 帮忙看一下，觉得可以的话可以私下交流一下
已经有3人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】f90子程序怎么编译错误？请大家帮忙，谢谢
已经有3人回复
- 【求助】帮忙改写展宽程序
已经有11人回复
4楼2010-06-15 13:15:13
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
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
★
robert2020(金币+1):多谢专家意见，辛苦了！ 2010-06-15 23:08:20
robert2020(金币+1):多谢专家意见，辛苦了！ 2010-06-15 23:08:20
|
大致了解一下 试试 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); while (q1~=q2 & p*y1~=p1*x1 & p*(1-y1)~=p2*(1-x1)) f1=b1/(1+b1*p); f2=b2/(1+b2*p); q1=integral(f1,0,p1,eps); q2=integral(f2,0,p2,eps); end end [ Last edited by lijinfeng042 on 2010-6-15 at 15:40 ] |
5楼2010-06-15 15:38:54