版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

琚裕波

金蟲 (小有名氣)

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

[求助] 求助！這個程序錯在哪里？1stopt運行不了

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];                                     //對比數(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             //實驗平均溫度
T=TT(iter)+273.15
P_Ca=PP(iter)/0.101325
//反應平衡常數(shù)計算
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ù)計算
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ù)計算
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))
//反應速率常數(shù)計算
KT1=k1_a*Exp(k1_b/8.314/T)
KT2=k2_a*Exp(k2_b/8.314/T)
//反應速率計算
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)
//通過反應速率計算出口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
//計算值與輸入值進行對比
y1(iter)=yCO
y2(iter)=yCO2
Next
EndProgram;
Data;

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

回復此樓

» 猜你喜歡

一志愿西電085401數(shù)一英一299求調(diào)劑六級521 已經(jīng)有3人回復
生物考研337分求調(diào)劑已經(jīng)有5人回復
求調(diào)劑：085600材料與化工，考材科基，總分319 已經(jīng)有14人回復
285求調(diào)劑已經(jīng)有10人回復
化學工程085602 305分求調(diào)劑已經(jīng)有29人回復
理學07化學 303求調(diào)劑已經(jīng)有5人回復
081200-11408-276學碩求調(diào)劑已經(jīng)有3人回復
土木304求調(diào)劑已經(jīng)有5人回復
一志愿中國科學院大學265求調(diào)劑已經(jīng)有7人回復
286求調(diào)劑已經(jīng)有8人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

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

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

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

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

dingd

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

應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jjdg: 金幣+1, 辛苦了 2013-09-13 00:34:04
琚裕波: 回帖置頂 2013-09-13 14:28:13
琚裕波: 金幣+20, ★★★很有幫助, 繼續(xù)回答啊 2013-09-13 19:17:15

1.5破解版不行，5.0正版運行沒問題，不用任何修改。

不知道你的參數(shù)有無范圍限制，有無過擬合現(xiàn)象。

均方差(RMSE):0.00320856901844522
殘差平方和(SSE):0.00041179660584506
相關系數(shù)(R): 0.999147398619087
相關系數(shù)之平方(R^2): 0.998295524167289
確定系數(shù)(DC): 0.998293391304997
F統(tǒng)計(F-Statistic): 16.4900938090552

參數(shù) 最佳估算
-------------------- -------------
k1_a 16199.5274910722
k1_b -54640.7646512102
k2_a 194674.147413486
k2_b -70718.1386986521
kco_a -1.81845739744211
kco_b -23518.1117376316
kco2_a -1.95039686557898
kco2_b 15779.9670841516
kh2_a -0.851175550051731
kh2_b -780.551417096644

c1.jpg

c2.jpg

贊一下(1人)

回復此樓

7樓2013-09-11 11:52:01

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

查看全部 11 個回答

dingd

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

應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

jjdg

版主 (知名作家)

程序強帖: 9
應助: 223 (大學生)
貴賓: 3.926
金幣: 51005.2
紅花: 238
帖子: 7227
在線: 474.2小時
蟲號: 27000
注冊: 2003-10-21
專業(yè): 藥物代謝與藥物動力學
管轄: 程序語言

【答案】應助回帖

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

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

贊一下

回復此樓

努力學習！以正當途徑！獲得需要的知識！

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

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

琚裕波

金蟲 (小有名氣)

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

引用回帖:

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

用未注冊的1.0版可以運行，但用注冊版1.0和1.5都不行

贊一下

回復此樓

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

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

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

查看全部 11 個回答

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

[考研] 生物考研337分求調(diào)劑 +3	cgxin 2026-03-30	5/250	2026-03-31 11:47 by cgxin
[考研] 求調(diào)劑：085600材料與化工，考材科基，總分319 +13	678lucky 2026-03-31	14/700	2026-03-31 11:43 by JourneyLucky
[考研] 081200-11408-276學碩求調(diào)劑 +3	崔wj 2026-03-31	3/150	2026-03-31 11:32 by 求調(diào)劑zz
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +10	dongfang59 2026-03-27	10/500	2026-03-30 23:42 by 果果媽咪
[考研] 11408總分309，一志愿東南大學求調(diào)劑，不挑專業(yè) +5	天賦帶到THU 2026-03-29	6/300	2026-03-30 20:49 by dick_runner
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +7	Zuhui0306 2026-03-25	7/350	2026-03-30 19:01 by 源_2020
[考研] 一志愿南開大學0710生物學359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	7/350	2026-03-30 18:29 by 兔兔兔111223314
[考研] 328求調(diào)劑 +8	嗯滴的基本都 2026-03-27	8/400	2026-03-30 17:20 by Wang200018
[考研] 321求調(diào)劑 +7	璞玉~~ 2026-03-25	8/400	2026-03-29 06:41 by 544594351
[考研] 一志愿南京航空航天大學材料學碩求調(diào)劑 +3	@taotao 2026-03-28	3/150	2026-03-28 10:26 by JourneyLucky
[考研] 070300化學求調(diào)劑 +4	起個名咋這么難 2026-03-27	4/200	2026-03-27 21:39 by 83503孫老師
[考研] 化學調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考博] 26申博 +3	加油沖��！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 279 分求調(diào)劑 +4	睡個好覺_16 2026-03-24	4/200	2026-03-27 15:05 by 醉在風里
[考研] 085600，材料與化工321分，求調(diào)劑 +9	大饞小子 2026-03-27	9/450	2026-03-27 14:30 by mmm just
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 321求調(diào)劑 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 材料調(diào)劑 5+4	想要一壺桃花水 2026-03-25	10/500	2026-03-26 19:56 by 不吃魚的貓
[考研] 085602化學工程求調(diào)劑。 +4	平樂樂樂 2026-03-26	4/200	2026-03-26 17:57 by fmesaito
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

琚裕波

[求助] 求助！這個程序錯在哪里？1stopt運行不了

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

dingd

【答案】應助回帖

dingd

【答案】應助回帖

jjdg

【答案】應助回帖

琚裕波

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 求助！這個程序錯在哪里？1stopt運行不了

» 本主題相關價值貼推薦，對您同樣有幫助: