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

zhaoshazhu

應(yīng)助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時(shí)
蟲(chóng)號(hào): 1812564
注冊(cè): 2012-05-12
專(zhuān)業(yè): 環(huán)境化工

[求助] 求1stopt指點(diǎn) 已有1人參與

下面是我利用低版本的1stopt編寫(xiě)四階龍格庫(kù)塔法計(jì)算動(dòng)力學(xué)方程，但是程序或許有問(wèn)題我沒(méi)有算出結(jié)果，請(qǐng)指點(diǎn)；
Title "Kinetics";
//Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;
//Variable t,F0,F,yH2,yDMM,yMe,yH2in,yDMMin,yMein;
//Variable  y1[output],y2[output],y3[output];
StartProgram  [VB];
Dim iter As  Integer
dim yHPMin As  Double
dim yPDOin As  Double
dim yNPAin As  Double
dim F0 As  Double
dim F As  Double
Dim RK_i As  Integer
Dim h As  Double
Dim i As  Integer
dim j As  Integer
dim c As  Double
dim r1 As  Double
dim r2 As  Double
dim r3 As  Double
dim n_H2 As  Double
dim n_DMM As  Double
dim n_Me As  Double
dim y_H2 As  Double
dim y_DMM As  Double
dim y_Me As  Double
dim y_HPM As  Double
dim y_PDO As  Double
dim y_NPA As  Double
dim k1 As  Double
dim k2 As  Double
dim k3 As  Double
dim k4 As  Double
dim k5 As  Double
dim k6 As  Double
dim k7 As  Double
dim k8 As  Double
dim k9 As  Double
dim Y_Y01 As  Double
dim Y_Y02 As  Double
dim Y_Y03 As  Double
dim Y_Y04 As  Double
dim Y_Y05 As  Double
dim Y_Y06 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 to Datalenth-1
   yHPMin=0
   yPDOin=0
   yNPAin=0
   h=1.0008/20
   Y_Y01=yHPMin
   Y_Y02=yPDOin
   Y_Y03=yNPAin
   Y_Y04=yH2in(iter)
   Y_Y05=yDMMin(inter)
   Y_Y06=yMein(inter)

   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
      For j=1 To 4
            y_HPM=x_x1
            y_PDO=x_x2
            y_NPA=x_x3
            n_H2=F0*yH2- 2*F*(y_HPM+y_PDO+y_NPA )
            n_DMM=F0*yDMM- F*(y_HPM+y_PDO+y_NPA )
            n_Me=F0*yMe+F*(y_HPM+y_PDO)
            y_H2= n_H2/F
            y_DMM= n_DMM/F
            y_ME= n_Me/F
            r1 = k7*k1*k2*y_H2*y_DMM /(1+k1*y_H2+k2*y_DMM+k3*y_Me+k4*y_HPM+k5*y_PDO+k6*y_NPA)^2
            r2 = k8*k1*k4*y_H2*y_HPM/((1+k1*y_H2+k2*y_DMM+k3*y_Me+k4*y_HPM+k5*y_PDO+k6*y_NPA)^2
            r3 = k9*k1*k5*y_H2*y_PDO/((1+k1*y_H2+k2*y_DMM+k3*y_Me+k4*y_HPM+k5*y_PDO+k6*y_NPA)^2
            d_d1 = r1-r2
            d_d2= r2-r3
            d_d3 = r3

            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
   y1(inter)=y_y1
   y2(inter)=y_y2
   y3(inter)=y_y3

Next
EndProgram;
Data;
2.05 0.48911 0.48786 0.92932 0.00258 0.06810 0.000031 0.00116 0.001137
1.37 0.73291 0.73179 0.93029 0.00155 0.06817 0.000078 0.00135 0.000573
1.62 0.61922 0.61815 0.91758 0.00173 0.08068 0.000104 0.00173 0.000475
1.98 0.50545 0.50453 0.89928 0.00187 0.09884 0.000289 0.00236 0.000444
1.21 0.82549 0.82425 0.91772 0.00158 0.08070 0.000414 0.00190 0.000280
1.49 0.67418 0.67281 0.89895 0.00224 0.09881 0.000641 0.00202 0.000292
2.70 0.37162 0.36991 0.81544 0.00531 0.17925 0.001304 0.00232 0.000368
1.02 0.97819 0.97574 0.92936 0.00254 0.06810 0.000206 0.00178 0.000361
0.82 1.22297 1.21986 0.92918 0.00273 0.06809 0.000465 0.00161 0.000189
0.97 1.03243 1.03056 0.91722 0.00213 0.08065 0.000790 0.00138 0.000149
1.19 0.84404 0.84146 0.89756 0.00379 0.09865 0.000972 0.00173 0.000210
2.16 0.46514 0.46287 0.81436 0.00663 0.17902 0.001744 0.00212 0.000349

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有180人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

Table Curve 2D 數(shù)據(jù)擬合求助！已經(jīng)有5人回復(fù)
看了下1stOpt的介紹，感覺(jué)略顯浮夸已經(jīng)有14人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有6人回復(fù)
Matlab用fsolve求解非線性方程，卡了三四天了跪求幫助！已經(jīng)有35人回復(fù)
高手指點(diǎn)：如何使用matlab使數(shù)據(jù)行與行間空一行已經(jīng)有6人回復(fù)
又會(huì)用proe5.0的嗎求大神求教程已經(jīng)有2人回復(fù)
1stopt解方程組已經(jīng)有12人回復(fù)
在用origin擬合的過(guò)程中，程序編寫(xiě)這一塊實(shí)在不會(huì)，求指教！已經(jīng)有10人回復(fù)
求高手幫忙擬合出方程已經(jīng)有7人回復(fù)
1stopt能進(jìn)行隱函數(shù)非線性擬合嗎？誰(shuí)能幫擬合一下新型液體的粘度混合計(jì)算公式已經(jīng)有5人回復(fù)
求助如果用SPSS對(duì)已知模型進(jìn)行非線形回歸？已經(jīng)有7人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有9人回復(fù)
用1stopt或者origin非線性擬合磁導(dǎo)率已經(jīng)有5人回復(fù)
SPSS求非線性方程參數(shù) 已經(jīng)有7人回復(fù)
1stopt4階非線性常微分方程，幫小弟看一下，跪求。已經(jīng)有7人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
麻煩高人指點(diǎn)一下1stopt和spss的非線性回歸中weibull的兩參數(shù)的公式應(yīng)該怎么表達(dá)�。� 已經(jīng)有9人回復(fù)
求SWMM軟件高手指導(dǎo) 已經(jīng)有15人回復(fù)

1樓 2015-01-06 18:42:21

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

zhaoshazhu

新蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時(shí)
蟲(chóng)號(hào): 1812564
注冊(cè): 2012-05-12
專(zhuān)業(yè): 環(huán)境化工

引用回帖:

6樓: Originally posted by NicoleLas at 2015-01-06 19:47:53
手冊(cè)里確實(shí)有這樣的例子。
開(kāi)頭幾句需要么？如果需要，去掉"//"，否則將被注釋掉，不起作用。
如果還是不可以用，哈哈，這個(gè)我就不知道了，也許是因?yàn)椴皇钦娴木壒拾伞?img src="https://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/smile.gif" align="absmiddle" border="0">...

謝謝

回復(fù)此樓

7樓2015-01-06 20:59:59

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

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

NicoleLas

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時(shí)
蟲(chóng)號(hào): 2872321
注冊(cè): 2013-12-15
性別: GG
專(zhuān)業(yè): 計(jì)算機(jī)軟件

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

感謝參與，應(yīng)助指數(shù) +1

1stOpt 1.5**版似乎不支持你的要求，而且你未去掉注釋符號(hào)"//"

贊一下

回復(fù)此樓

Be water, my friend.

2樓2015-01-06 18:57:32

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

zhaoshazhu

新蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時(shí)
蟲(chóng)號(hào): 1812564
注冊(cè): 2012-05-12
專(zhuān)業(yè): 環(huán)境化工

引用回帖:

2樓: Originally posted by NicoleLas at 2015-01-06 18:57:32
1stOpt 1.5**版似乎不支持你的要求，而且你未去掉注釋符號(hào)"//"

哪個(gè)版本支持��？

回復(fù)此樓

3樓2015-01-06 19:11:33

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

NicoleLas

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時(shí)
蟲(chóng)號(hào): 2872321
注冊(cè): 2013-12-15
性別: GG
專(zhuān)業(yè): 計(jì)算機(jī)軟件

引用回帖:

3樓: Originally posted by zhaoshazhu at 2015-01-06 19:11:33
哪個(gè)版本支持啊？...

估計(jì)1stOpt 4.0之后吧，直接支持微分方程擬合。

贊一下

回復(fù)此樓

Be water, my friend.

4樓2015-01-06 19:19:39

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

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

[考研] 352求調(diào)劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +6	幸運(yùn)哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 293求調(diào)劑 +12	zjl的號(hào) 2026-03-16	17/850	2026-03-22 16:51 by i_cooler
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[考研] 306求0703調(diào)劑一志愿華中師范 +5	紙魚(yú)ly 2026-03-21	5/250	2026-03-21 17:11 by 學(xué)員8dgXkO
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專(zhuān)碩，總分308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢(mèng)人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨