版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門版塊排行榜

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

返回列表

zhaoshazhu

新蟲 (小有名氣)

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

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

下面是我利用低版本的1stopt編寫四階龍格庫塔法計(jì)算動(dòng)力學(xué)方程，但是程序或許有問題我沒有算出結(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ù)此樓

» 猜你喜歡

歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有27人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有144人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有12人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)
山西大同大學(xué)物理學(xué)碩士研究生招收調(diào)劑生已經(jīng)有0人回復(fù)
[調(diào)劑信息]211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)
［調(diào)劑信息］211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)
[調(diào)劑信息]211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)

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

Table Curve 2D 數(shù)據(jù)擬合求助！已經(jīng)有5人回復(fù)
看了下1stOpt的介紹，感覺略顯浮夸已經(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擬合的過程中，程序編寫這一塊實(shí)在不會(huì)，求指教！已經(jīng)有10人回復(fù)
求高手幫忙擬合出方程已經(jīng)有7人回復(fù)
1stopt能進(jìn)行隱函數(shù)非線性擬合嗎？誰能幫擬合一下新型液體的粘度混合計(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的回帖

NicoleLas

木蟲 (正式寫手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時(shí)
蟲號(hào): 2872321
注冊(cè): 2013-12-15
性別: GG
專業(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

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時(shí)
蟲號(hào): 1812564
注冊(cè): 2012-05-12
專業(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

木蟲 (正式寫手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時(shí)
蟲號(hào): 2872321
注冊(cè): 2013-12-15
性別: GG
專業(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的回帖

zhaoshazhu

新蟲 (小有名氣)

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

引用回帖:

4樓: Originally posted by NicoleLas at 2015-01-06 19:19:39
估計(jì)1stOpt 4.0之后吧，直接支持微分方程擬合。...

老師說低版本的只要把四階龍格庫塔法編進(jìn)去也可以計(jì)算出來呢，可是我的就是沒有結(jié)果。

贊一下

回復(fù)此樓

5樓2015-01-06 19:37:53

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

NicoleLas

木蟲 (正式寫手)

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

引用回帖:

5樓: Originally posted by zhaoshazhu at 2015-01-06 19:37:53
老師說低版本的只要把四階龍格庫塔法編進(jìn)去也可以計(jì)算出來呢，可是我的就是沒有結(jié)果。...

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

贊一下

回復(fù)此樓

Be water, my friend.

6樓2015-01-06 19:47:53

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

zhaoshazhu

新蟲 (小有名氣)

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

引用回帖:

6樓: Originally posted by NicoleLas at 2015-01-06 19:47:53
手冊(cè)里確實(shí)有這樣的例子。
開頭幾句需要么？如果需要，去掉"//"，否則將被注釋掉，不起作用。
如果還是不可以用，哈哈，這個(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的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 zhaoshazhu 的主題更新

返回列表

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

[考研] 302求調(diào)劑 +4	zyx上岸！ 2026-04-02	4/200	2026-04-02 19:41 by 喬噠噠噠
[考研] 一志愿上海海洋大學(xué)083200食品學(xué)碩，求調(diào)劑，接受其他專業(yè) +6	what張 2026-04-01	7/350	2026-04-02 16:48 by zzsw+
[考研] 材料專碩322分 +9	哈哈哈吼吼吼哈 2026-04-02	9/450	2026-04-02 15:45 by zikuzi1853
[考研] 材料工程322分 +8	哈哈哈吼吼吼哈 2026-04-01	8/400	2026-04-02 11:53 by 3041
[考研] 求調(diào)劑！生物與醫(yī)藥專碩 +4	逆轉(zhuǎn)陸先生 2026-04-01	4/200	2026-04-02 11:51 by xiaoranmu
[考研] 265求調(diào)劑11408 +4	劉小鹿lu 2026-03-27	4/200	2026-04-02 11:27 by guyan1000
[考研] 292求調(diào)劑 +14	是妍子也是研子 2026-03-30	14/700	2026-04-02 10:17 by sanrepian
[基金申請(qǐng)] esi高被引論文是不是能對(duì)中標(biāo)有所加分和幫助呢 +4	redcom 2026-04-01	4/200	2026-04-02 10:16 by wqy8563
[考研] 【求調(diào)劑】新能源材料本科，一志愿211，初試321 +6	求調(diào)劑學(xué)校， 2026-04-02	6/300	2026-04-02 09:41 by 晴空210210
[考研] 343求調(diào)劑085601 +4	要努力學(xué)習(xí)x 2026-03-29	5/250	2026-04-02 09:01 by xuhui0306
[考研] 一志愿346上海大學(xué)生物學(xué) +3	上海大學(xué)346調(diào)劑 2026-04-01	3/150	2026-04-02 08:36 by w蟲蟲123
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14	秋二十二 2026-03-30	14/700	2026-04-01 11:45 by chemdavid
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +15	Doleres 2026-03-31	16/800	2026-04-01 09:49 by lfj11
[考研] 289求調(diào)劑 +7	BrightLL 2026-03-29	7/350	2026-03-31 22:05 by 544594351
[考研] 080500-315分復(fù)試調(diào)劑 +9	上岸3821 2026-03-31	9/450	2026-03-31 17:29 by 唐沐兒
[考研] 吉大生物學(xué)326分求調(diào)劑 +3	sunnyupup 2026-03-31	3/150	2026-03-31 09:28 by longlotian
[考研] 一志愿大連理工大學(xué)材料求調(diào)劑 +6	Gymno 2026-03-30	6/300	2026-03-31 07:26 by 無際的草原
[考研] 317分一志愿南理工材料工程本科湖工大求調(diào)劑 +12	芋泥小鈴鐺 2026-03-28	12/600	2026-03-30 17:06 by wangjy2002
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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

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

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