版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

[求助] 求1stopt指点已有1人参与

下面是我利用低版本的1stopt编写四阶龙格库塔法计算动力学方程，但是程序或许有问题我没有算出结果，请指点；
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

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

Table Curve 2D 数据拟合求助！已经有5人回复
看了下1stOpt的介绍，感觉略显浮夸已经有14人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有6人回复
Matlab用fsolve求解非线性方程，卡了三四天了跪求帮助！已经有35人回复
高手指点：如何使用matlab使数据行与行间空一行已经有6人回复
又会用proe5.0的吗求大神求教程已经有2人回复
1stopt解方程组已经有12人回复
在用origin拟合的过程中，程序编写这一块实在不会，求指教！已经有10人回复
求高手帮忙拟合出方程已经有7人回复
1stopt能进行隐函数非线性拟合吗？谁能帮拟合一下新型液体的粘度混合计算公式已经有5人回复
求助如果用SPSS对已知模型进行非线形回归？已经有7人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有9人回复
用1stopt或者origin非线性拟合磁导率已经有5人回复
SPSS求非线性方程参数已经有7人回复
1stopt4阶非线性常微分方程，帮小弟看一下，跪求。已经有7人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复
麻烦高人指点一下1stopt和spss的非线性回归中weibull的两参数的公式应该怎么表达啊？已经有9人回复
求SWMM软件高手指导已经有15人回复

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

4楼: Originally posted by NicoleLas at 2015-01-06 19:19:39
估计1stOpt 4.0之后吧，直接支持微分方程拟合。...

老师说低版本的只要把四阶龙格库塔法编进去也可以计算出来呢，可是我的就是没有结果。

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 7 个回答

NicoleLas

木虫 (正式写手)

应助: 66 (初中生)
金币: 4288.7
红花: 26
帖子: 514
在线: 177.2小时
虫号: 2872321
注册: 2013-12-15
性别: GG
专业: 计算机软件

【答案】应助回帖

感谢参与，应助指数 +1

1stOpt 1.5**版似乎不支持你的要求，而且你未去掉注释符号"//"

赞一下

回复此楼

Be water, my friend.

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

2楼: Originally posted by NicoleLas at 2015-01-06 18:57:32
1stOpt 1.5**版似乎不支持你的要求，而且你未去掉注释符号"//"

哪个版本支持啊？

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

应助: 66 (初中生)
金币: 4288.7
红花: 26
帖子: 514
在线: 177.2小时
虫号: 2872321
注册: 2013-12-15
性别: GG
专业: 计算机软件

引用回帖:

3楼: Originally posted by zhaoshazhu at 2015-01-06 19:11:33
哪个版本支持啊？...

估计1stOpt 4.0之后吧，直接支持微分方程拟合。

赞一下

回复此楼

Be water, my friend.

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 7 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 山东省面上项目限额评审 +3	石瑞0426 2026-03-19	3/150	2026-03-21 18:59 by RRRKKK
[考研] 286求调剂 +8	Faune 2026-03-21	8/400	2026-03-21 17:16 by 学员8dgXkO
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5	脱颖而出 2026-03-16	15/750	2026-03-21 10:16 by 脱颖而出
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 301求调剂 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 材料工程（专）一志愿985 初试335求调剂 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 初始318分求调剂（有工作经验） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 化学求调剂 +4	临泽境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 328求调剂，英语六级551，有科研经历 +4	生物工程调剂 2026-03-17	8/400	2026-03-21 02:12 by JourneyLucky
[考研] 一志愿华中科技大学，080502，354分求调剂 +5	守候夕阳CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 330求调剂 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 考研调剂求学校推荐 +3	伯乐29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 290求调剂 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 一志愿华中农业071010，总分320求调剂 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 学员8dgXkO
[考研] 260求调剂 +3	朱芷琳 2026-03-20	3/150	2026-03-20 20:35 by 学员8dgXkO
[考研] 材料专硕306英一数二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 301求调剂 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 一志愿，福州大学材料专硕339分求调剂 +3	木子momo青争 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 070303 总分349求调剂 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy

24小时热门版块排行榜

zhaoshazhu

[求助] 求1stopt指点 已有1人参与

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

zhaoshazhu

NicoleLas

【答案】应助回帖

zhaoshazhu

NicoleLas

[求助] 求1stopt指点已有1人参与