版塊導(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í)熱門版塊排行榜

返回列表

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ù)此樓

» 猜你喜歡

哈爾濱理工大學(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)有202人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(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的介紹，感覺略顯浮夸已經(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ā)表

[考研] 263求調(diào)劑 +4	yqdszhdap－ 2026-03-22	4/200	2026-03-22 21:20 by 1144970272
[考研] 284求調(diào)劑 +5	Zhao anqi 2026-03-22	5/250	2026-03-22 17:38 by barlinike
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 070300化學(xué)求調(diào)劑 +3	苑豆豆 2026-03-20	3/150	2026-03-22 14:48 by ColorlessPI
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +11	吃吃吃才有意義 2026-03-19	11/550	2026-03-21 18:23 by 學(xué)員8dgXkO
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 材料學(xué)碩297已過四六級(jí)求調(diào)劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao

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

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

» 猜你喜歡

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

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

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