| 查看: 3058 | 回复: 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 |
回复此楼» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有272人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！！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
5