| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2109 | 回復(fù): 16 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[交流]
【求助】幫忙看一個程序 已有3人參與
|
||||
|
%IAST(理想吸附溶液理論通過單組份的吸附數(shù)據(jù)獲得多組分的吸附數(shù)據(jù)) 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) %被積函數(shù):f %積分左端點(diǎn):a %積分右端點(diǎn):b %eps:精度 %積分結(jié)果:q %step 積分的子區(qū)間數(shù) if(nargin==3) eps=1.0e-4; %默認(rèn)精度為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次復(fù)合梯形積分 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. |
混合氣體吸附理論研究 |


鐵桿木蟲 (著名寫手)
方丈大師

木蟲 (小有名氣)
Matlab
|
大致了解一下 試試 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 ] |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +6 | 雙馬尾痞老板2 2026-04-01 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085602化學(xué)工程268分蹲調(diào)劑 +8 | 月照花林。 2026-04-01 | 8/400 |
|
|
[考研] 292求調(diào)劑 +17 | 木蟲er12138 2026-04-01 | 17/850 |
|
|
[考研] 286求調(diào)劑 +16 | PolarBear11 2026-03-26 | 16/800 |
|
|
[考研] 379求調(diào)劑 +3 | ?苦瓜不苦 2026-04-01 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
[考研] 0710生物學(xué)考研調(diào)劑 +3 | 李多米lee. 2026-03-27 | 4/200 |
|
|
[考研] 086000生物與醫(yī)藥 初試274求調(diào)劑 +5 | 小叮當(dāng)來了 2026-03-30 | 6/300 |
|
|
[考研] 348求調(diào)劑 +9 | zzzzyk123 2026-04-01 | 9/450 |
|
|
[考研] 262求調(diào)劑 +9 | 勵志一定發(fā)文章 2026-03-31 | 10/500 |
|
|
[考研] 材料調(diào)劑 +10 | Eujd1 2026-03-31 | 11/550 |
|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|
[考研] 一志愿華東理工大學(xué),080500學(xué)碩,317分,求調(diào)劑 +9 | s1145 2026-03-31 | 9/450 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 一志愿浙江大學(xué)工科動力工程370,數(shù)一121,專業(yè)課135,現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|