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

登录注册

应《网络安全法》要求，自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
专业: 环境化工

引用回帖:

6楼: Originally posted by NicoleLas at 2015-01-06 19:47:53
手册里确实有这样的例子。
开头几句需要么？如果需要，去掉"//"，否则将被注释掉，不起作用。
如果还是不可以用，哈哈，这个我就不知道了，也许是因为不是正版的缘故吧。

...

谢谢

回复此楼

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

已阅回复此楼关注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 个回答

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

[考研] 328求调剂，英语六级551，有科研经历 +5	生物工程调剂 2026-03-17	9/450	2026-03-21 23:32 by zhujy1982
[考研] 材料学硕301分求调剂 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 一志愿211，0703化学310分求调剂 +3	努力奋斗112 2026-03-15	3/150	2026-03-21 22:21 by peike
[考研] 0703化学调剂 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 学员8dgXkO
[考研] 298求调剂 +4	上岸6666@ 2026-03-20	4/200	2026-03-21 17:14 by 学员8dgXkO
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4	顺顺顺mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 考研调剂求学校推荐 +3	伯乐29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 350求调剂 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] 一志愿西南交通专硕材料355 本科双非求调剂 +5	西南交通专材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 295复试调剂 +8	简木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 0817 化学工程 299分求调剂有科研经历有二区文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿南京航空航天大学大学，080500材料科学与工程学硕 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 320求调剂0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 学员8dgXkO
[考研] 0703化学调剂 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 0854，计算机类招收调剂 +3	胡辣汤放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 考研求调剂 +3	橘颂. 2026-03-17	4/200	2026-03-17 21:43 by 有只狸奴
[考研] 326求调剂 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 无懈可击111
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷

24小时热门版块排行榜

zhaoshazhu

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

» 猜你喜欢

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

zhaoshazhu

NicoleLas

【答案】应助回帖

zhaoshazhu

NicoleLas

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