[交流] 用1stopt擬合動(dòng)力學(xué)方程參數(shù)無(wú)法運(yùn)行，新手不懂如何調(diào)試，請(qǐng)求各位高手幫忙，謝謝！已有2人參與

程序是利用龍格庫(kù)塔法模擬動(dòng)力學(xué)方程的相關(guān)參數(shù)，我用低版本的軟件運(yùn)行的時(shí)候一直報(bào)錯(cuò)，一直無(wú)法運(yùn)行，提示我的編寫有問(wèn)題，菜鳥不太懂，希望各位大神幫忙指正，謝謝啦。
下面是程序代碼：
Title "cyclohexanone";
Parameters k1,E1,k2,E2,kA,EA,kH2,EH2,kCNE,ECNE,kCOL,ECOL,DEN;
Variable V0,PP,TT,yAin,yH2in,yArin,yCNEin,yCOLin;
Variable yCNE[Output],yCOL[Output];
StartProgram [VB];
Dim iter As Integer
Dim N0 As Double
Dim RK_i As Integer
Dim h As Double
Dim i As Integer
dim j As Integer
dim c As Double
dim DM As Double
dim RCNE As Double
dim RCOL As Double
dim yAr As Double
dim yA As Double
dim yCNE As Double
dim yH2 As Double
dim yCOL As Double
Dim KT1 As Double
dim KT2 As Double
dim KA As Double
dim KH2 As Double
dim KCNE As Double
dim KCOL As Double
Dim T As Double
Dim BB As Double
Dim pH2 As Double
dim pCOL As Double
dim pCNE As Double
dim pA As Double
Dim Y_Y01 As Double
Dim Y_Y02 As Double
Dim Y_Y03 As Double
dim y_y1 As Double
dim y_y2 As Double
dim y_y3 As Double
Dim x_x1 As Double
Dim x_x2 As Double
Dim x_x3 As Double
dim d_d1 As Double
dim d_d2 As Double
dim d_d3 As Double
dim e_e1 As Double
dim e_e2 As Double
dim e_e3 As Double
dim e_e4 As Double
dim e_e5 As Double
for iter=0 to Datalength-1
N0=V0(iter)
h=0.35/20
Y_Y01=yCNEin(iter)
Y_Y02=yCOLin(iter)
Y_Y03=0
For RK_i=1 to 20
      e_e1=0.5
      e_e2=0.5
      e_e3=1
      e_e4=1
      e_e5=0.5
      x_x1=Y_Y01
      x_x2=Y_Y02
      x_x3=Y_Y03
      y_y1=Y_Y01
      y_y2=Y_Y02
      y_y3=Y_Y03
      For j=1 To 4
         yCNE=x_x1
         yCOL=x_x2

         BB=(1+2*yA)/(1+2*yAin(iter))
         yA=BB*(yAin(iter)+yCNEin(iter))-yCNE
         yH2=BB*(yH2in(iter)+yCOLin(iter)+2yCNEin(iter))-2yCNE-yCOL
         yAr=BB*yArin(iter)
         T=TT(iter)+273.15

         pA=PP(iter)*yA
         pH2=PP(iter)*yH2
         pCNE=PP(iter)*yCNE
         pCOL=PP(iter)*yCOL

         KT1=k1*Exp(-E1/(8.315*T))
         KT2=k2*Exp(-E2/(8.315*T))
         KA=kA*Exp(-EA/(8.315*T))
         KH2=kH2*Exp(-EH2/(8.315*T))
         KCNE=kCNE*Exp(-ECNE/(8.315*T))
         KCOL=kCOL*Exp(-ECOL/(8.315*T))

         DEN=1+KA*PA+KH2*PH2+KCNE*PCNE+KCOL*PCOL
         RCNE=(KT1*KA*PA*(KH2*PH2)^2)/(DEN^3)
         RCOL=(KT2*KCNE*PCNE*KH2*PH2)/(DEN^2)
         DM=(1+2*yCNE)/(N0*(1+2*yCNEin(iter)))
         d_d1=DM*(1+2*yCNE)*RCNE
         d_d2=DM*(RCOL+2*yCOL*RCNE)
         d_d3=1

         c=h*d_d1
         if j=1 then
            x_x1=Y_Y01+e_e1*c
            y_y1=y_y1+e_e2*c/3
         elseif j=2 then
            x_x1=Y_Y01+e_e2*c
            y_y1=y_y1+e_e3*c/3
         elseif j=3 then
            x_x1=Y_Y01+e_e3*c
            y_y1=y_y1+e_e4*c/3
         elseif j=4 then
            x_x1=Y_Y01+e_e4*c
            y_y1=y_y1+e_e5*c/3
         end if
         c=h*d_d2
         if j=1 then
            x_x2=Y_Y02+e_e1*c
            y_y2=y_y2+e_e2*c/3
         elseif j=2 then
            x_x2=Y_Y02+e_e2*c
            y_y2=y_y2+e_e3*c/3
         elseif j=3 then
            x_x2=Y_Y02+e_e3*c
            y_y2=y_y2+e_e4*c/3
         elseif j=4 then
            x_x2=Y_Y02+e_e4*c
            y_y2=y_y2+e_e5*c/3
         end if
         c=h*d_d3
         if j=1 then
            x_x3=Y_Y03+e_e1*c
            y_y3=y_y3+e_e2*c/3
         elseif j=2 then
            x_x3=Y_Y03+e_e2*c
            y_y3=y_y3+e_e3*c/3
         elseif j=3 then
            x_x3=Y_Y03+e_e3*c
            y_y3=y_y3+e_e4*c/3
         elseif j=4 then
            x_x3=Y_Y03+e_e4*c
            y_y3=y_y3+e_e5*c/3
         end if
      Next
      Y_Y01=y_y1
      Y_Y02=y_y2
      Y_Y03=y_y3
Next
yCNE(iter)=y_y1
yCOL(iter)=y_y2
Next
EndProgram;
Data;
4.06 0.2 120 0.113 0.788 0.099 0.0001 0.0001 0.00796539 7.23E-05
2.03 0.2 120 0.114 0.787 0.099 0.0001 0.0001 0.018276944 0.000165986
1.02 0.2 120 0.117 0.784 0.099 0.0001 0.0001 0.027001782 0.000272745
1.02 0.2 120 0.115 0.786 0.099 0.0001 0.0001 0.028928857 0.000292211
0.51 0.2 120 0.114 0.787 0.099 0.0001 0.0001 0.05419921 0.000602822
1.66 0.2 141 0.133 0.771 0.096 0.0001 0.0001 0.024173136 0.000244173
4.07 0.2 140 0.115 0.787 0.098 0.0001 0.0001 0.018688673 0.000169726
2.04 0.2 140 0.117 0.785 0.098 0.0001 0.0001 0.041576064 0.00041996
1.03 0.2 140 0.127 0.775 0.098 0.0001 0.0001 0.077775178 0.001668313
1.03 0.2 140 0.12 0.781 0.099 0.0001 0.0001 0.067879227 0.000824444
0.52 0.2 141 0.125 0.777 0.098 0.0001 0.0001 0.124489331 0.002152918
4.46 0.2 162 0.122 0.776 0.102 0.0001 0.0001 0.011771229 9.49E-05
3.32 0.2 162 0.121 0.774 0.105 0.0001 0.0001 0.016481008 0.000149676
1.66 0.2 162 0.133 0.771 0.096 0.0001 0.0001 0.039389391 0.000397873
4.08 0.2 161 0.117 0.784 0.099 0.0001 0.0001 0.035950835 0.000326496
2.05 0.2 161 0.12 0.781 0.099 0.0001 0.0001 0.071190475 0.000791805
1.04 0.2 161 0.137 0.766 0.097 0.0001 0.0001 0.137103964 0.002229333
1.04 0.2 161 0.132 0.771 0.097 0.0001 0.0001 0.134541091 0.00218766
4.45 0.2 182 0.119 0.778 0.103 0.0001 0.0001 0.018492489 0.000149133
3.32 0.2 182 0.122 0.772 0.106 0.0001 0.0001 0.029740384 0.000270094
1.66 0.2 182 0.13 0.773 0.097 0.0001 0.0001 0.060335598 0.000671073
4.09 0.2 181 0.12 0.782 0.098 0.0001 0.0001 0.058478171 0.000590689
2.06 0.2 181 0.126 0.777 0.097 0.0001 0.0001 0.107592821 0.001527687
1.05 0.2 181 0.141 0.763 0.096 0.0001 0.0001 0.186241228 0.006156735
1.05 0.2 182 0.142 0.762 0.096 0.0001 0.0001 0.183768277 0.005488445
1.04 0.2 181 0.138 0.766 0.096 0.0001 0.0001 0.179658491 0.005747588
1.04 0.2 182 0.136 0.767 0.097 0.0001 0.0001 0.176490707 0.005458475

文件中是我需要模擬的動(dòng)力學(xué)方程

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 動(dòng)力學(xué)方程.docx

2017-01-10 14:23:50, 13.85 K

» 猜你喜歡

重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有186人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有21人回復(fù)

1樓 2017-01-10 14:24:52

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

也許有用：http://www.gaoyang168.com/bbs/viewthread.php?tid=10962975

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

2樓2017-01-10 15:28:10

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

lipenggg

專家顧問(wèn) (正式寫手)

專家經(jīng)驗(yàn): +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時(shí)
蟲號(hào): 3929201
注冊(cè): 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計(jì)算模擬

1.5版本不支持編程

回復(fù)此樓

一團(tuán)和氣

3樓2017-01-16 10:43:17

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

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

[考研] 0805求調(diào)劑 +5	是水分 2026-03-31	5/250	2026-03-31 20:48 by yuq
[考研] 359求調(diào)劑 +6	鄧鄧鄧書書 2026-03-25	6/300	2026-03-31 20:16 by 無(wú)際的草原
[考研] 0856求調(diào)劑 +9	楒桉 2026-03-28	9/450	2026-03-31 19:06 by 暮澤12
[考研] 324分 085600材料與化工 +6	呆鵝oor 2026-03-27	6/300	2026-03-31 19:03 by xixiangyue6220
[考研] 285求調(diào)劑 +6	AZMK 2026-03-30	11/550	2026-03-31 17:33 by AZMK
[考研] 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。 +12	nnnnnnn5 2026-03-25	12/600	2026-03-31 16:21 by 記事本2026
[考研] 322求調(diào)劑 +10	舊吢 2026-03-24	10/500	2026-03-31 13:31 by 無(wú)際的草原
[考研] 085600 295分求調(diào)劑 +13	W55j 2026-03-30	15/750	2026-03-31 13:29 by 王亮_大連醫(yī)科大
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-31	4/200	2026-03-31 11:56 by jp9609
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑六級(jí)521 +4	愛(ài)吃大鴨梨 2026-03-31	4/200	2026-03-31 11:51 by 搏擊518
[考研] 一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 +3	西瓜田的守望者 2026-03-30	3/150	2026-03-31 11:08 by asdfzly
[考研] 086000生物與醫(yī)藥調(diào)劑 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv