當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

琚裕波

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 175.1
散金: 10
帖子: 51
在線: 20.7小時(shí)
蟲號(hào): 2094139
注冊(cè): 2012-10-29
性別: GG
專業(yè): 能源化工

[求助] 求助！這個(gè)程序錯(cuò)在哪里？1stopt運(yùn)行不了

Title "Kinetics";
Parameters k1_a,k1_b,k2_a,k2_b,KCO_a,KCO_b,KCO2_a,KCO2_b,KH2_a,KH2_b;  //需要擬合的參數(shù)
Variable TT,PP,Nin,yH2in,yCOin,yCO2in,yN2in,yN2out;                   //輸入數(shù)據(jù)
Variable y1[Output],y2[Output];                                     //對(duì)比數(shù)據(jù)
StartProgram [VB];
Dim iter As Integer
dim TBA As Double
dim RCO As Double
dim RCO2 As Double
Dim KT1 As Double
dim KF1 As Double
dim KT2 As Double
dim KF2 As Double
dim KH2 As Double
dim KCO As Double
dim KCO2 As Double
Dim T As Double
dim P_Ca As Double
Dim fH2 As Double
dim fCO2 As Double
dim fCO As Double
dim fM As Double
dim fH2O As Double
Dim Beta1 As Double
dim Beta2 As Double
dim Nt as double
dim yCO as double
dim yCO2 as double
for iter=0 to Datalength-1
TBA=503.2             //實(shí)驗(yàn)平均溫度
T=TT(iter)+273.15
P_Ca=PP(iter)/0.101325
//反應(yīng)平衡常數(shù)計(jì)算
KF1=Exp(13.1652+9203.26/T-5.92839*Ln(T)-0.00352404*T+0.0000102264*T*T-0.00000000769446*T^3+2.38583E-12*T^4)*(0.101325^(-2))
KF2=Exp(1.6654+4553.34/T-2.72613*Ln(T)-0.01422914*T+0.000017206*T*T-0.00000001106294*T^3+3.19698E-12*T^4)*(0.101325^(-2))
//逸度系數(shù)計(jì)算
fH2=PP(iter)*yH2in(iter)*Exp((0.110785+35.3324/T-5005.47/(T*T*T)-19.6109*yH2in(iter)/T-20.9799*yH2in(iter)*yH2in(iter)/T)*P_Ca/T)
fCO2=PP(iter)*yCO2in(iter)*Exp((-0.343605+428.452/T-69217700/(T*T*T)-327.402*yCO2in(iter)/T-374.954*yCO2in(iter)*yCO2in(iter)/T)*P_Ca/T)
fCO=PP(iter)*yCOin(iter)*Exp((-0.093261+189.156/T-399940/(T*T*T)-181.527*yCOin(iter)/T+140.001*yCOin(iter)*yCOin(iter)/T)*P_Ca/T)
Beta1=fM/(KF1*fCO*fH2^2)
Beta2=fM*fH2O/(KF2*fCO2*fH2^3)
//吸附平衡常數(shù)計(jì)算
KCO=Exp(KCO_a+KCO_b*(1/T-1/TBA))
KCO2=Exp(KCO2_a+KCO2_b*(1/T-1/TBA))
KH2=Exp(KH2_a+KH2_b*(1/T-1/TBA))
//反應(yīng)速率常數(shù)計(jì)算
KT1=k1_a*Exp(k1_b/8.314/T)
KT2=k2_a*Exp(k2_b/8.314/T)
//反應(yīng)速率計(jì)算
RCO=KT1*fCO*(fH2^2)*(1-Beta1)/((1+KCO*fCO+KCO2*fCO2+KH2*fH2)^3)
RCO2=KT2*fCO2*(fH2^3)*(1-Beta2)/((1+KCO*fCO+KCO2*fCO2+KH2*fH2)^4)
//通過反應(yīng)速率計(jì)算出口CO和CO2組成，其中1.9108為催化劑質(zhì)量
Nt=Nin(iter)*yN2in(iter)/yN2out(iter)
yCO=(Nin(iter)*yCOin(iter)-RCO*1.9108)/Nt
yCO2=(Nin(iter)*yCO2in(iter)-RCO2*1.9108)/Nt
//計(jì)算值與輸入值進(jìn)行對(duì)比
y1(iter)=yCO
y2(iter)=yCO2
Next
EndProgram;
Data;

[ Last edited by nono2009 on 2013-9-10 at 20:21 ]

回復(fù)此樓

» 猜你喜歡

291求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑已經(jīng)有8人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有12人回復(fù)
070300化學(xué)求調(diào)劑已經(jīng)有6人回復(fù)
考研化學(xué)308分求調(diào)劑已經(jīng)有7人回復(fù)
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑已經(jīng)有5人回復(fù)
工科0856求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) 已經(jīng)有4人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有16人回復(fù)
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）已經(jīng)有5人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助，用1stopt解方程組，兩次得到的解不唯一？什么原因，大神幫忙已經(jīng)有6人回復(fù)
sos!!1stopt解非線性方程組，一直出現(xiàn)常數(shù)定義錯(cuò)誤的問題！求大神幫助！已經(jīng)有12人回復(fù)
求大俠幫助跑一下1stopt 已經(jīng)有12人回復(fù)
求助用1stOpt的sharedmodel 運(yùn)行一組數(shù)據(jù)求參數(shù)，謝謝了啊，急需！已經(jīng)有4人回復(fù)
哪位朋友有注冊(cè)版1stOpt軟件？急需已經(jīng)有5人回復(fù)
1stOPT1.5版每次仿真結(jié)果不一樣，求助！急！已經(jīng)有3人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
求高版本1stopt，解決一個(gè)問題，微分方程組參數(shù)擬合已經(jīng)有3人回復(fù)
運(yùn)算過程有復(fù)數(shù)出現(xiàn)，求高版本1stOpt 已經(jīng)有6人回復(fù)
1stopt運(yùn)行程序，總是出現(xiàn)字符函數(shù)錯(cuò)誤的提示，請(qǐng)高手幫助已經(jīng)有5人回復(fù)
大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒有遇到過這種情況？已經(jīng)有8人回復(fù)
1stopt4階非線性常微分方程，幫小弟看一下，跪求。已經(jīng)有7人回復(fù)
1stOpt高版本求助！求高版本的幫忙走一下謝了！已經(jīng)有16人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
學(xué)1stopt要多久？已經(jīng)有12人回復(fù)
哪位大哥能幫我用1stopt3.0以上版本跑一個(gè)擬合程序已經(jīng)有11人回復(fù)
1stOpt的函數(shù)作圖不靠譜嗎？已經(jīng)有4人回復(fù)
【求助】此程序沒有錯(cuò)誤，怎么運(yùn)行不出來？已經(jīng)有6人回復(fù)
用1stopt進(jìn)行數(shù)據(jù)擬合時(shí)，無法運(yùn)行sharedmodel ,求助高手已經(jīng)有7人回復(fù)
【求助】安裝ansys時(shí)，運(yùn)行Keygen程序,但不會(huì)自動(dòng)生成一個(gè)文件License.dat,為什么? 已經(jīng)有12人回復(fù)

一切都會(huì)向著最好的方向發(fā)展！

1樓 2013-09-10 16:11:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

727786917

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 133.2
帖子: 149
在線: 38小時(shí)
蟲號(hào): 1895792
注冊(cè): 2012-07-16
性別: GG
專業(yè): 信息理論與信息系統(tǒng)

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝參與 2013-09-13 00:34:19
琚裕波: 金幣+1 2013-09-15 09:39:04

用未注冊(cè)的1.0版可以運(yùn)行，但用注冊(cè)版1.0和1.5都不行

贊一下(1人)

回復(fù)此樓

天道酬勤

9樓2013-09-12 22:03:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 11 個(gè)回答

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝參與 2013-09-11 00:50:34

1：沒見你的數(shù)據(jù)；
2：用的哪個(gè)版本��？網(wǎng)上1.5的處理這種編程模式估計(jì)有問題。

贊一下

回復(fù)此樓

2樓2013-09-10 20:55:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jjdg

版主 (知名作家)

程序強(qiáng)帖: 9
應(yīng)助: 223 (大學(xué)生)
貴賓: 3.926
金幣: 51005.2
紅花: 238
帖子: 7227
在線: 474.2小時(shí)
蟲號(hào): 27000
注冊(cè): 2003-10-21
專業(yè): 藥物代謝與藥物動(dòng)力學(xué)
管轄: 程序語言

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
ben_ladeng: 金幣+1, 謝謝指教 2013-09-21 22:14:41

不完整的代碼
你做kinetic可以去下個(gè)kinetica5.1來試試

贊一下

回復(fù)此樓

努力學(xué)習(xí)！以正當(dāng)途徑！獲得需要的知識(shí)！

3樓2013-09-11 00:51:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

琚裕波

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 175.1
散金: 10
帖子: 51
在線: 20.7小時(shí)
蟲號(hào): 2094139
注冊(cè): 2012-10-29
性別: GG
專業(yè): 能源化工

引用回帖:

3樓: Originally posted by jjdg at 2013-09-11 00:51:51
不完整的代碼
你做kinetic可以去下個(gè)kinetica5.1來試試

用未注冊(cè)的1.0版可以運(yùn)行，但用注冊(cè)版1.0和1.5都不行

贊一下

回復(fù)此樓

一切都會(huì)向著最好的方向發(fā)展！

4樓2013-09-11 10:31:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 070300化學(xué)求調(diào)劑 +6	苑豆豆 2026-03-20	6/300	2026-03-23 18:16 by YMU施老師
[考研] 招08考數(shù)學(xué) +6	laoshidan 2026-03-20	14/700	2026-03-23 14:37 by 15614359529
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 311求調(diào)劑 +6	冬十三 2026-03-18	6/300	2026-03-22 20:18 by edmund7
[考研] 求調(diào)劑院校信息 +6	CX 330 2026-03-21	6/300	2026-03-22 15:25 by 無懈可擊111
[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 求助 +5	夢(mèng)里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518