| 查看: 3281 | 回复: 16 | |||
[求助]
1stopt微分方程组参数拟合 已有2人参与
|
|||
|
真心希望有人能帮我一下，这个问题已经困扰我很多天了，这10个金币是我所有的金币了。 具体的内容是这样 然后再得 c是光速为3*10^10，k=0.5，M、tsd和L是未知参数，需要拟合求解，现在只知道能得到t和LL的实验值。具体的计算应该是v和Eint都是不知道的，若知道M、tsd和L的值，我用的是ode45，参数v和Eint根据t来回算（初始值是0,0），然后带入到LL里面。现在反过来了不知道M、tsd和L的值我就不知道怎么写了，用MATLAB的ode45和lsqnonlin写了个程序，最后求出来不是很理想，而且跟初始值那些影响蛮大。用1stopt写了一个，运行不了（我觉得肯定是错的），不知道怎么写，能帮我改改，然后算出来吗？真的麻烦了，在线等啊。。。 Parameters M,tsd,L; Constant c=3*10^10,k=0.5; Variable t,v,E,LL; Function v'=E/(M*v*t); E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t; LL=E*t*(4*pi*c*v)/(3*k*M); Data; //t LL 502272 5.65432*10^42 762325 1.30331*10^43 1114354 9.00958*10^42 1364122 5.3601*10^42 1706100 3.56514*10^42 1964980 2.5611*10^42 2140916 2.27497*10^42 2313059 1.94825*10^42 2489034 1.50971*10^42 2742752 1.22059*10^42 3001044 1.04568*10^42 |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有206人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！！1stopt进行常微分方程组参数模拟
已经有5人回复
- 1stOpt软件微分方程参数拟合，运行后没有任何反应。
已经有5人回复
- 1stopt 微分方程数据拟合求助
已经有2人回复
- 1stopt微分方程组的参数估计的请教
已经有3人回复
- 用1stOpt求解微分方程
已经有2人回复
- [求助] 1stopt 微分方程代码求助
已经有5人回复
- 求帮忙运行一下下面的1stopt的程序，微分方程-参数估计
已经有6人回复
- 用1stopt软件做微分方程组的参数估计
已经有13人回复
- 求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
已经有9人回复
- 求助1stopt或者matlab 微分方程拟合
已经有0人回复
- 求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
已经有1人回复
- 求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 帮忙看看能否用1stopt计算出微分方程组
已经有0人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
10楼2016-03-04 17:26:28
2楼2016-03-03 19:41:21
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2016-03-03 21:12:00
4楼2016-03-03 21:50:59
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2016-03-04 11:37:12
6楼2016-03-04 13:04:23
|
Parameters M,tsd,L; Constant c=3*10^10,k=0.5; InitialODEValue t=0.1,E=0.1,v=0.1,LL=0; Variable t,LL； ODEFunction v'=E/(M*v*t); E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t; LL'=E*(4*pi*c*v)/(3*k*M) Data; //t LL 0.1 0.1 502272 1.12575E+37 762325 1.70965E+37 1114354 8.08503E+36 1364122 3.92934E+36 1706100 2.08964E+36 1964980 1.30337E+36 2140916 1.06261E+36 2313059 8.42281E+35 2489034 6.06545E+35 2742752 4.45025E+35 3001044 3.4844E+35 改成这样了，还是不能运行 |
7楼2016-03-04 16:15:27
8楼2016-03-04 16:17:08
|
Parameters M,tsd,L; Constant c=3*10^10,k=0.5; InitialODEValue E=0.1,v=0.1,LL=0.1; Variable t,LL； ODEFunction v'=E/(M*v*t); E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t; LL'=E*(4*pi*c*v)/(3*k*M)+E'*t*(4*pi*c*v)/(3*k*M)+E*t*(4*pi*c*v')/(3*k*M); Data; 0.1 0.1 502272 5.65432*10^42 762325 1.30331*10^43 1114354 9.00958*10^42 1364122 5.3601*10^42 1706100 3.56514*10^42 1964980 2.5611*10^42 2140916 2.27497*10^42 2313059 1.94825*10^42 2489034 1.50971*10^42 2742752 1.22059*10^42 3001044 1.04568*10^42 应该是这样，运行不了。。。。 |
9楼2016-03-04 16:48:12
28