| 查看: 1830 | 回复: 5 | |||
[求助]
[求助] 1stopt 微分方程代码求助 已有2人参与
|
|
这两天在用1stopt研究一个动力学微分方程的拟合求动力学参数，看起来应该是比较简单的方程，但是研究了很长时间，自己也根据网上能搜到的例子编了两个程序跑了下，但是结果相差比较大，我也不知道程序编的对不对，请懂的大神指点一二，跪谢了！我使用的是试用版的1stopt 5.0 方程是这样的： dx/dt=-p1*x; dy/dt=p1*x-p2*y; 初始值：t=0，此时x=3.5,y=0 测量的数据是t 和 y。具体数据请见程序中。p1和p2应为正。 程序1：（dx/dt=-p1*x 这个方程，我先积分成了x=x0*exp(-p1*t)，也不知道对不对？） //y'=dy/dt=3.5*p1*exp(-p1*t)-p2*y Constant V =0; //Initial Values Constant dt=1; Parameter p(2); Variable t, y; StartProgram [Pascal]; Procedure MainModel; var i,t: integer; y0 : double; Begin y0 := V; for i := 0 to DataLength - 1 do begin y0 := (3.5*p1*exp(-p1*t)-p2*y0)*dt+y0; y := y0; end; End; EndProgram; Data; //0 0 10 0.010515 14 0.023065 17 0.037836 20 0.062683 22 0.063288 24 0.064385 26 0.063359 28 0.062942 30 0.061384 32 0.060358 35 0.065102 45 0.059101 60 0.043515 这个程序的输出报告和图： ![[求助] 1stopt 微分方程代码求助](https://muchongimg.xmcimg.com/data/bcs/2015/0408/w76h3287083_1428458526_740.jpg#opennewwindow) ![[求助] 1stopt 微分方程代码求助-1](https://muchongimg.xmcimg.com/data/bcs/2015/0408/bw248h3287083_1428458553_282.jpg#opennewwindow) 关于程序1，我的主要疑惑是：1.不知道这个程序编的对不对？ 2.关于 dt 这个数值，在程序中的作用，我比较担心的是，我的t是不连续取值的，dt会不会使得t在程序中变成连续取值了？从图上看我的这一担心就更明显，因为横坐标只取到了12 程序2： InitialODEValue t=0,x=3.5,y=0; Parameters p1,p2; Variable t[Output],y[Output]; ODEFunction x'=-p1*x; y'=p1*x-p2*y; Data; //0 0 10 0.010514645 14 0.023065468 17 0.037836003 20 0.062682935 22 0.06328842 24 0.064385481 26 0.063358784 28 0.062942392 30 0.061383892 32 0.060358109 35 0.065102196 45 0.059100765 60 0.043515354 这个程序的输出报告和图： ![[求助] 1stopt 微分方程代码求助-2](https://muchongimg.xmcimg.com/data/bcs/2015/0408/w75h3287083_1428459505_278.jpg#opennewwindow) ![[求助] 1stopt 微分方程代码求助-3](https://muchongimg.xmcimg.com/data/bcs/2015/0408/bw201h3287083_1428459527_652.jpg#opennewwindow) 关于程序2，我的主要疑惑是：1.不知道这个程序编的对不对？ 2.关于ODEFunction 语句的使用是否正确？ODEFunction语句对微分函数的定义默认是怎么样的？在函数中我都没定义t，他是怎么计算的？也涉及到t的取值问题，ODEFunction是不是自动把t按照1,2,3,4,5……取值了？ 总结一下问题：1.请高手指导我这个动力学微分方程应该如何编程，能给代码就最好了；（授人以鱼，但是是最主要的。。。） 2.请指出我自己编写程序的问题，以及关于两个程序中t的取值，是否连续？还是间歇，如何确定（授人以渔） 3.关于ODEFunction的定义，关于程序1中dt的定义。 谢谢！！！ |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有138人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt自变量超多时是否有简单语法？AND 自变量多 数据量大的话1stopt都会算的特慢么
已经有7人回复
- 1stopt数据拟合求助
已经有8人回复
- 请教一下各位大神怎么求解 Runge–Kutta method的微分方程的各个系数
已经有12人回复
- 1stopt解方程组
已经有12人回复
- 求用高版本的1stOpt拟合方程
已经有8人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 一个二阶常微分方程的求解方法？
已经有5人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 求组！求解一个微分方程组的符号解，并画出图像！！！
已经有7人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-04-08 10:40:32
3楼2015-04-08 10:53:12
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
结果稳定唯一，但效果不是很好： Root of Mean Square Error (RMSE):0.0119983163868176 Sum of Square Error:0.00187147474953628 Correlation Coef. (R): 0.77730265287911 R-Square: 0.604199414172902 Determination Coef. (DC): 0.514455569529589 F-Statistic: 22.6549813151864 Parameter Best Estimate -------------------- ------------- p1 0.00131507520277041 p2 0.0675323910740055 ![[求助] 1stopt 微分方程代码求助-4](https://muchongimg.xmcimg.com/data/bcs/2015/0408/w214h291104_1428463918_135.jpg#opennewwindow) c2.jpg |
4楼2015-04-08 11:32:00
5楼2015-04-08 12:05:38
6楼2015-04-08 12:07:39
